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

Prior-guided individualized thalamic parcellation based on local diffusion characteristics

Comprising numerous subnuclei, the thalamus intricately interconnects the cortex and subcortex, orchestrating various facets of brain functions. Extracting personalized parcellation patterns for these subnuclei is crucial, as different thalamic nuclei play varying roles in cognition and serve as therapeutic targets for neuromodulation. However, accurately delineating the thalamic nuclei boundary at the individual level is challenging due to intersubject variability. In this study, we proposed a prior-guided parcellation (PG-par) method to achieve robust individualized thalamic parcellation based on a central-boundary prior. We first constructed probabilistic atlas of thalamic nuclei using high-quality diffusion MRI datasets based on the local diffusion characteristics. Subsequently, high-probability voxels in the probabilistic atlas were utilized as prior guidance to train unique multiple classification models for each subject based on a multilayer perceptron. Finally, we employed the trained model to predict the parcellation labels for thalamic voxels and construct individualized thalamic parcellation. Through a test-retest assessment, the proposed prior-guided individualized thalamic parcellation exhibited excellent reproducibility and the capacity to detect individual variability. Compared with group atlas registration and individual clustering parcellation, the proposed PG-par demonstrated superior parcellation performance under different scanning protocols and clinic settings. Furthermore, the prior-guided individualized parcellation exhibited better correspondence with the histological staining atlas. The proposed prior-guided individualized thalamic parcellation method contributes to the personalized modeling of brain parcellation.

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.

Altered dynamic functional connectivity of the thalamus subregions in patients with schizophrenia

BACKGROUND

Previous neuroimaging studies indicated that patients with schizophrenia showed impaired thalamus and thalamo-cortical circuits. However, the dynamic functional connectivity (dFC) patterns of the thalamus remain unclear. In this study, we explored the dFC of the thalamus in SZ patients and whether clinical features are correlated with altered dFC.

METHODS

Forty-three patients with schizophrenia and 31 healthy controls underwent 3.0 T rs-fMRI. Based on the human Brainnetome atlas, the thalamus is divided into 8 subregions. Subsequently, we performed flexible least squares method to calculate the dFC of each thalamus subregions.

RESULTS

Compared with healthy controls, patients with schizophrenia exhibited increased dFC between the thalamus and cerebellar, visual-related cortex, sensorimotor-related cortex, and frontal lobe. In addition, we found that the dFC of the thalamus and the right fusiform gyrus was negatively associated with age of onset.

CONCLUSIONS

Our findings demonstrated that the dFC of specific thalamus sub-regions is altered in patients with schizophrenia. Our results further suggested the dysconnectivity of thalamus plays an important role in the pathophysiology of schizophrenia.

Characterizing functional modules in the human thalamus: coactivation-based parcellation and systems-level functional decoding

The human thalamus relays sensory signals to the cortex and facilitates brain-wide communication. The thalamus is also more directly involved in sensorimotor and various cognitive functions but a full characterization of its functional repertoire, particularly in regard to its internal anatomical structure, is still outstanding. As a putative hub in the human connectome, the thalamus might reveal its functional profile only in conjunction with interconnected brain areas. We therefore developed a novel systems-level Bayesian reverse inference decoding that complements the traditional neuroinformatics approach towards a network account of thalamic function. The systems-level decoding considers the functional repertoire (i.e., the terms associated with a brain region) of all regions showing co-activations with a predefined seed region in a brain-wide fashion. Here, we used task-constrained meta-analytic connectivity-based parcellation (MACM-CBP) to identify thalamic subregions as seed regions and applied the systems-level decoding to these subregions in conjunction with functionally connected cortical regions. Our results confirm thalamic structure-function relationships known from animal and clinical studies and revealed further associations with language, memory, and locomotion that have not been detailed in the cognitive neuroscience literature before. The systems-level decoding further uncovered large systems engaged in autobiographical memory and nociception. We propose this novel decoding approach as a useful tool to detect previously unknown structure-function relationships at the brain network level, and to build viable starting points for future studies.

In humans, striato-pallido-thalamic projections are largely segregated by their origin in either the striosome-like or matrix-like compartments

Cortico-striato-thalamo-cortical (CSTC) loops are fundamental organizing units in mammalian brains. CSTCs process limbic, associative, and sensorimotor information in largely separated but interacting networks. CTSC loops pass through paired striatal compartments, striosome (aka patch) and matrix, segregated pools of medium spiny projection neurons with distinct embryologic origins, cortical/subcortical structural connectivity, susceptibility to injury, and roles in behaviors and diseases. Similarly, striatal dopamine modulates activity in striosome and matrix in opposite directions. Routing CSTCs through one compartment may be an anatomical basis for regulating discrete functions. We used differential structural connectivity, identified through probabilistic diffusion tractography, to distinguish the striatal compartments (striosome-like and matrix-like voxels) in living humans. We then mapped compartment-specific projections and quantified structural connectivity between each striatal compartment, the globus pallidus interna (GPi), and 20 thalamic nuclei in 221 healthy adults. We found that striosome-originating and matrix-originating streamlines were segregated within the GPi: striosome-like connectivity was significantly more rostral, ventral, and medial. Striato-pallido-thalamic streamline bundles that were seeded from striosome-like and matrix-like voxels transited spatially distinct portions of the white matter. Matrix-like streamlines were 5.7-fold more likely to reach the GPi, replicating animal tract-tracing studies. Striosome-like connectivity dominated in six thalamic nuclei (anteroventral, central lateral, laterodorsal, lateral posterior, mediodorsal-medial, and medial geniculate). Matrix-like connectivity dominated in seven thalamic nuclei (centromedian, parafascicular, pulvinar-anterior, pulvinar-lateral, ventral lateral-anterior, ventral lateral-posterior, ventral posterolateral). Though we mapped all thalamic nuclei independently, functionally-related nuclei were matched for compartment-level bias. We validated these results with prior thalamostriate tract tracing studies in non-human primates and other species; where reliable data was available, all agreed with our measures of structural connectivity. Matrix-like connectivity was lateralized (left > right hemisphere) in 18 thalamic nuclei, independent of handedness, diffusion protocol, sex, or whether the nucleus was striosome-dominated or matrix-dominated. Compartment-specific biases in striato-pallido-thalamic structural connectivity suggest that routing CSTC loops through striosome-like or matrix-like voxels is a fundamental mechanism for organizing and regulating brain networks. Our MRI-based assessments of striato-thalamic connectivity in humans match and extend the results of prior tract tracing studies in animals. Compartment-level characterization may improve localization of human neuropathologies and improve neurosurgical targeting in the GPi and thalamus.

A phylogenetically-conserved axis of thalamocortical connectivity in the human brain

The thalamus enables key sensory, motor, emotive, and cognitive processes via connections to the cortex. These projection patterns are traditionally considered to originate from discrete thalamic nuclei, however recent work showing gradients of molecular and connectivity features in the thalamus suggests the organisation of thalamocortical connections occurs along a continuous dimension. By performing a joint decomposition of densely sampled gene expression and non-invasive diffusion tractography in the adult human thalamus, we define a principal axis of genetic and connectomic variation along a medial-lateral thalamic gradient. Projections along this axis correspond to an anterior-posterior cortical pattern and are aligned with electrophysiological properties of the cortex. The medial-lateral axis demonstrates phylogenetic conservation, reflects transitions in neuronal subtypes, and shows associations with neurodevelopment and common brain disorders. This study provides evidence for a supra-nuclear axis of thalamocortical organisation characterised by a graded transition in molecular properties and anatomical connectivity.

Psilocybin induces spatially constrained alterations in thalamic functional organizaton and connectivity

BACKGROUND

Classic psychedelics, such as psilocybin and LSD, and other serotonin 2A receptor (5-HT2AR) agonists evoke acute alterations in perception and cognition. Altered thalamocortical connectivity has been hypothesized to underlie these effects, which is supported by some functional MRI (fMRI) studies. These studies have treated the thalamus as a unitary structure, despite known differential 5-HT2AR expression and functional specificity of different intrathalamic nuclei. Independent Component Analysis (ICA) has been previously used to identify reliable group-level functional subdivisions of the thalamus from resting-state fMRI (rsfMRI) data. We build on these efforts with a novel data-maximizing ICA-based approach to examine psilocybin-induced changes in intrathalamic functional organization and thalamocortical connectivity in individual participants.

METHODS

Baseline rsfMRI data (n=38) from healthy individuals with a long-term meditation practice was utilized to generate a statistical template of thalamic functional subdivisions. This template was then applied in a novel ICA-based analysis of the acute effects of psilocybin on intra- and extra-thalamic functional organization and connectivity in follow-up scans from a subset of the same individuals (n=18). We examined correlations with subjective reports of drug effect and compared with a previously reported analytic approach (treating the thalamus as a single functional unit).

RESULTS

Several intrathalamic components showed significant psilocybin-induced alterations in spatial organization, with effects of psilocybin largely localized to the mediodorsal and pulvinar nuclei. The magnitude of changes in individual participants correlated with reported subjective effects. These components demonstrated predominant decreases in thalamocortical connectivity, largely with visual and default mode networks. Analysis in which the thalamus is treated as a singular unitary structure showed an overall numerical increase in thalamocortical connectivity, consistent with previous literature using this approach, but this increase did not reach statistical significance.

CONCLUSIONS

We utilized a novel analytic approach to discover psilocybin-induced changes in intra- and extra-thalamic functional organization and connectivity of intrathalamic nuclei and cortical networks known to express the 5-HT2AR. These changes were not observed using whole-thalamus analyses, suggesting that psilocybin may cause widespread but modest increases in thalamocortical connectivity that are offset by strong focal decreases in functionally relevant intrathalamic nuclei.

Multi-atlas thalamic nuclei segmentation on standard T1-weighed MRI with application to normal aging

Specific thalamic nuclei are implicated in healthy aging and age-related neurodegenerative diseases. However, few methods are available for robust automated segmentation of thalamic nuclei. The threefold aims of this study were to validate the use of a modified thalamic nuclei segmentation method on standard T1 MRI data, to apply this method to quantify age-related volume declines, and to test functional meaningfulness by predicting performance on motor testing. A modified version of THalamus Optimized Multi-Atlas Segmentation (THOMAS) generated 22 unilateral thalamic nuclei. For validation, we compared nuclear volumes obtained from THOMAS parcellation of white-matter-nulled (WMn) MRI data to T1 MRI data in 45 participants. To examine the effects of age/sex on thalamic nuclear volumes, T1 MRI available from a second data set of 121 men and 117 women, ages 20-86 years, were segmented using THOMAS. To test for functional ramifications, composite regions and constituent nuclei were correlated with Grooved Pegboard test scores. THOMAS on standard T1 data showed significant quantitative agreement with THOMAS from WMn data, especially for larger nuclei. Sex differences revealing larger volumes in men than women were accounted for by adjustment with supratentorial intracranial volume (sICV). Significant sICV-adjusted correlations between age and thalamic nuclear volumes were detected in 20 of the 22 unilateral nuclei and whole thalamus. Composite Posterior and Ventral regions and Ventral Anterior/Pulvinar nuclei correlated selectively with higher scores from the eye-hand coordination task. These results support the use of THOMAS for standard T1-weighted data as adequately robust for thalamic nuclear parcellation.

In-vivo data-driven parcellation of Heschl's gyrus using structural connectivity

The human auditory cortex around Heschl’s gyrus (HG) exhibits diverging patterns across individuals owing to the heterogeneity of its substructures. In this study, we investigated the subregions of the human auditory cortex using data-driven machine-learning techniques at the individual level and assessed their structural and functional profiles. We studied an openly accessible large dataset of the Human Connectome Project and identified the subregions of the HG in humans using data-driven clustering techniques with individually calculated imaging features of cortical folding and structural connectivity information obtained via diffusion magnetic resonance imaging tractography. We characterized the structural and functional profiles of each HG subregion according to the cortical morphology, microstructure, and functional connectivity at rest. We found three subregions. The first subregion (HG₁) occupied the central portion of HG, the second subregion (HG₂) occupied the medial-posterior-superior part of HG, and the third subregion (HG₃) occupied the lateral-anterior-inferior part of HG. The HG₃ exhibited strong structural and functional connectivity to the association and paralimbic areas, and the HG₁ exhibited a higher myelin density and larger cortical thickness than other subregions. A functional gradient analysis revealed a gradual axis expanding from the HG₂ to the HG₃. Our findings clarify the individually varying structural and functional organization of human HG subregions and provide insights into the substructures of the human auditory cortex.

Convolutional Neural Network Based Frameworks for Fast Automatic Segmentation of Thalamic Nuclei from Native and Synthesized Contrast Structural MRI

Thalamic nuclei have been implicated in several neurological diseases. Thalamic nuclei parcellation from structural MRI is challenging due to poor intra-thalamic nuclear contrast while methods based on diffusion and functional MRI are affected by limited spatial resolution and image distortion. Existing multi-atlas based techniques are often computationally intensive and time-consuming. In this work, we propose a 3D convolutional neural network (CNN) based framework for thalamic nuclei parcellation using T1-weighted Magnetization Prepared Rapid Gradient Echo (MPRAGE) images. Transformation of images to an efficient representation has been proposed to improve the performance of subsequent classification tasks especially when working with limited labeled data. We investigate this by transforming the MPRAGE images to White-Matter-nulled MPRAGE (WMn-MPRAGE) contrast, previously shown to exhibit good intra-thalamic nuclear contrast, prior to the segmentation step. We trained two 3D segmentation frameworks using MPRAGE images (n = 35 subjects): (a) a native contrast segmentation (NCS) on MPRAGE images and (b) a synthesized contrast segmentation (SCS) where synthesized WMn-MPRAGE representation generated by a contrast synthesis CNN were used. Thalamic nuclei labels were generated using THOMAS, a multi-atlas segmentation technique proposed for WMn-MPRAGE images. The segmentation accuracy and clinical utility were evaluated on a healthy cohort (n = 12) and a cohort (n = 45) comprising of healthy subjects and patients with alcohol use disorder (AUD), respectively. Both the segmentation CNNs yielded comparable performances on most thalamic nuclei with Dice scores greater than 0.84 for larger nuclei and at least 0.7 for smaller nuclei. However, for some nuclei, the SCS CNN yielded significant improvements in Dice scores (medial geniculate nucleus, P = 0.003, centromedian nucleus, P = 0.01) and percent volume difference (ventral anterior, P = 0.001, ventral posterior lateral, P = 0.01) over NCS. In the AUD cohort, the SCS CNN demonstrated a significant atrophy in ventral lateral posterior nucleus in AUD patients compared to healthy age-matched controls (P = 0.01), agreeing with previous studies on thalamic atrophy in alcoholism, whereas the NCS CNN showed spurious atrophy of the ventral posterior lateral nucleus. CNN-based segmentation of thalamic nuclei provides a fast and automated technique for thalamic nuclei prediction in MPRAGE images. The transformation of images to an efficient representation, such as WMn-MPRAGE, can provide further improvements in segmentation performance.

A unified 3D map of microscopic architecture and MRI of the human brain

We present the first three-dimensional (3D) concordance maps of cyto- and fiber architecture of the human brain, combining histology, immunohistochemistry, and 7-T quantitative magnetic resonance imaging (MRI), in two individual specimens. These 3D maps each integrate data from approximately 800 microscopy sections per brain, showing neuronal and glial cell bodies, nerve fibers, and interneuronal populations, as well as ultrahigh-field quantitative MRI, all coaligned at the 200-μm scale to the stacked blockface images obtained during sectioning. These unprecedented 3D multimodal datasets are shared without any restrictions and provide a unique resource for the joint study of cell and fiber architecture of the brain, detailed anatomical atlasing, or modeling of the microscopic underpinnings of MRI contrasts.

The anterior thalamic nuclei: core components of a tripartite episodic memory system

Standard models of episodic memory focus on hippocampal-parahippocampal interactions, with the neocortex supplying sensory information and providing a final repository of mnemonic representations. However, recent advances have shown that other regions make distinct and equally critical contributions to memory. In particular, there is growing evidence that the anterior thalamic nuclei have a number of key cognitive functions that support episodic memory. In this article, we describe these findings and argue for a core, tripartite memory system, comprising a 'temporal lobe' stream (centred on the hippocampus) and a 'medial diencephalic' stream (centred on the anterior thalamic nuclei) that together act on shared cortical areas. We demonstrate how these distributed brain regions form complementary and necessary partnerships in episodic memory formation.

Ultrasonic Deep Brain Neuromodulation in Acute Disorders of Consciousness: A Proof-of-Concept

The promotion of recovery in patients who have entered a disorder of consciousness (DOC; e.g., coma or vegetative states) following severe brain injury remains an enduring medical challenge despite an ever-growing scientific understanding of these conditions. Indeed, recent work has consistently implicated altered cortical modulation by deep brain structures (e.g., the thalamus and the basal ganglia) following brain damage in the arising of, and recovery from, DOCs. The (re)emergence of low-intensity focused ultrasound (LIFU) neuromodulation may provide a means to selectively modulate the activity of deep brain structures noninvasively for the study and treatment of DOCs. This technique is unique in its combination of relatively high spatial precision and noninvasive implementation. Given the consistent implication of the thalamus in DOCs and prior results inducing behavioral recovery through invasive thalamic stimulation, here we applied ultrasound to the central thalamus in 11 acute DOC patients, measured behavioral responsiveness before and after sonication, and applied functional MRI during sonication. With respect to behavioral responsiveness, we observed significant recovery in the week following thalamic LIFU compared with baseline. With respect to functional imaging, we found decreased BOLD signals in the frontal cortex and basal ganglia during LIFU compared with baseline. In addition, we also found a relationship between altered connectivity of the sonicated thalamus and the degree of recovery observed post-LIFU.

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.

In Vivo Super-Resolution Track-Density Imaging for Thalamic Nuclei Identification

The development of novel techniques for the in vivo, non-invasive visualization and identification of thalamic nuclei has represented a major challenge for human neuroimaging research in the last decades. Thalamic nuclei have important implications in various key aspects of brain physiology and many of them show selective alterations in various neurologic and psychiatric disorders. In addition, both surgical stimulation and ablation of specific thalamic nuclei have been proven to be useful for the treatment of different neuropsychiatric diseases. The present work aimed at describing a novel protocol for histologically guided delineation of thalamic nuclei based on short-tracks track-density imaging (stTDI), which is an advanced imaging technique exploiting high angular resolution diffusion tractography to obtain super-resolved white matter maps. We demonstrated that this approach can identify up to 13 distinct thalamic nuclei bilaterally with very high inter-subject (ICC: 0.996, 95% CI: 0.993-0.998) and inter-rater (ICC:0.981; 95% CI:0.963-0.989) reliability, and that both subject-based and group-level thalamic parcellation show a fair share of similarity to a recent standard-space histological thalamic atlas. Finally, we showed that stTDI-derived thalamic maps can be successfully employed to study structural and functional connectivity of the thalamus and may have potential implications both for basic and translational research, as well as for presurgical planning purposes.