The Amsterdam Ultra-high field adult lifespan database (AHEAD): A freely available multimodal 7 Tesla submillimeter magnetic resonance imaging database

ATOMMIC: An Advanced Toolbox for Multitask Medical Imaging Consistency to facilitate Artificial Intelligence applications from acquisition to analysis in Magnetic Resonance Imaging

BACKGROUND AND OBJECTIVES

Artificial intelligence (AI) is revolutionizing Magnetic Resonance Imaging (MRI) along the acquisition and processing chain. Advanced AI frameworks have been applied in various successive tasks, such as image reconstruction, quantitative parameter map estimation, and image segmentation. However, existing frameworks are often designed to perform tasks independently of each other or are focused on specific models or single datasets, limiting generalization. This work introduces the Advanced Toolbox for Multitask Medical Imaging Consistency (ATOMMIC), a novel open-source toolbox that streamlines AI applications for accelerated MRI reconstruction and analysis. ATOMMIC implements several tasks using deep learning (DL) models and enables MultiTask Learning (MTL) to perform related tasks in an integrated manner, targeting generalization in the MRI domain.

METHODS

We conducted a comprehensive literature review and analyzed 12,479 GitHub repositories to assess the current landscape of AI frameworks for MRI. Subsequently, we demonstrate how ATOMMIC standardizes workflows and improves data interoperability, enabling effective benchmarking of various DL models across MRI tasks and datasets. To showcase ATOMMIC's capabilities, we evaluated twenty-five DL models on eight publicly available datasets, focusing on accelerated MRI reconstruction, segmentation, quantitative parameter map estimation, and joint accelerated MRI reconstruction and segmentation using MTL.

RESULTS

ATOMMIC's high-performance training and testing capabilities, utilizing multiple GPUs and mixed precision support, enable efficient benchmarking of multiple models across various tasks. The framework's modular architecture implements each task through a collection of data loaders, models, loss functions, evaluation metrics, and pre-processing transformations, facilitating seamless integration of new tasks, datasets, and models. Our findings demonstrate that ATOMMIC supports MTL for multiple MRI tasks with harmonized complex-valued and real-valued data support while maintaining active development and documentation. Task-specific evaluations demonstrate that physics-based models outperform other approaches in reconstructing highly accelerated acquisitions. These high-quality reconstruction models also show superior accuracy in estimating quantitative parameter maps. Furthermore, when combining high-performing reconstruction models with robust segmentation networks through MTL, performance is improved in both tasks.

CONCLUSIONS

ATOMMIC advances MRI reconstruction and analysis by leveraging MTL and ensuring consistency across tasks, models, and datasets. This comprehensive framework serves as a versatile platform for researchers to use existing AI methods and develop new approaches in medical imaging.

Probabilistic Tractography-Based Tremor Network Connectivity in Tremor Dominant Parkinson's Disease and Essential Tremor plus

Background: The basal ganglia-thalamocortical (BGTC) and cerebello-thalamocortical (CTC) networks are implicated in tremor genesis; however, exact contributions across disorders have not been studied. Objective: Evaluate the structural connectivity of BGTC and CTC in tremor-dominant Parkinson's disease (TDPD) and essential tremor plus (ETP) with the aid of probabilistic tractography and graph theory analysis. Methods: Structural connectomes of the BGTC and CTC were generated by probabilistic tractography for TDPD (n = 25), ETP (ET with rest tremor, n = 25), and healthy control (HC, n = 22). The Brain Connectivity Toolbox was used for computing standard topological graph measures of segregation, integration, and centrality. Tremor severity was ascertained using the Fahn-Tolosa-Marin tremor rating scale (FTMRS). Results: There was no difference in total FTMRS scores. Compared with HC, TDPD had a lower global efficiency and characteristic path length. Abnormality in segregation, integration, and centrality of bilateral putamen, globus pallidus externa (GPe), and GP interna (GPi), with reduction of centrality of right caudate and cerebellar lobule 8, was observed. ETP showed reduction in segregation and integration of right GPe and GPi, ventrolateral posterior nucleus, and centrality of right putamen, compared with HC. Differences between TDPD and ETP were a reduction of strength of the right putamen, and lower clustering coefficient, local efficiency, and strength of the left GPi in TDPD. Conclusions: Contrary to expectations, TDPD and ETP may not be significantly different with regard to tremor pathogenesis, with definite overlaps. There may be fundamental similarities in network disruption across different tremor disorders with the same tremor activation patterns, along with disease-specific changes.

Quantitative MRI at 7-Tesla reveals novel frontocortical myeloarchitecture anomalies in major depressive disorder

Whereas meta-analytical data highlight abnormal frontocortical macrostructure (thickness/surface area/volume) in Major Depressive Disorder (MDD), the underlying microstructural processes remain uncharted, due to the use of conventional MRI scanners and acquisition techniques. We uniquely combined Ultra-High Field MRI at 7.0 Tesla with Quantitative Imaging to map intracortical myelin (proxied by longitudinal relaxation time T1) and iron concentration (proxied by transverse relaxation time T2*), microstructural processes deemed particularly germane to cortical macrostructure. Informed by meta-analytical evidence, we focused specifically on orbitofrontal and rostral anterior cingulate cortices among adult MDD patients (N = 48) and matched healthy controls (HC; N = 10). Analyses probed the association of MDD diagnosis and clinical profile (severity, medication use, comorbid anxiety disorders, childhood trauma) with aforementioned microstructural properties. MDD diagnosis (p's < 0.05, Cohen's D = 0.55-0.66) and symptom severity (p's < 0.01, r = 0.271-0.267) both related to decreased intracortical myelination (higher T1 values) within the lateral orbitofrontal cortex, a region tightly coupled to processing negative affect and feelings of sadness in MDD. No relations were found with local iron concentrations. These findings allow uniquely fine-grained insights on frontocortical microstructure in MDD, and cautiously point to intracortical demyelination as a possible driver of macroscale cortical disintegrity in MDD.

Cortical folding correlates to aging and Alzheimer's Disease's cognitive and CSF biomarkers

This manuscript presents the quantification and correlation of three aspects of Alzheimer's Disease evolution, including structural, biochemical, and cognitive assessments. We aimed to test a novel structural biomarker for neurodegeneration based on a cortical folding model for mammals. Our central hypothesis is that the cortical folding variable, representative of axonal tension in white matter, is an optimal discriminator of pathological aging and correlates with altered loadings in Cerebrospinal Fluid samples and a decline in cognition and memory. We extracted morphological features from T1w 3T MRI acquisitions using FreeSurfer from 77 Healthy Controls (age = 66 ± 8.4, 69% females), 31 Mild Cognitive Impairment (age = 72 ± 4.8, 61% females), and 13 Alzheimer's Disease patients (age = 77 ± 6.1, 62% females) of recruited volunteers in Brazil to test its discriminative power using optimal cut-point analysis. Cortical folding distinguishes the groups with reasonable accuracy (Healthy Control-Alzheimer's Disease, accuracy = 0.82; Healthy Control-Mild Cognitive Impairment, accuracy = 0.56). Moreover, Cerebrospinal Fluid biomarkers (total Tau, A[Formula: see text]1-40, A[Formula: see text]1-42, and Lipoxin) and cognitive scores (Cognitive Index, Rey's Auditory Verbal Learning Test, Trail Making Test, Digit Span Backward) were correlated with the global neurodegeneration in MRI aiming to describe health, disease, and the transition between the two states using morphology.

Architecture of the subthalamic nucleus

The subthalamic nucleus (STN) is a major neuromodulation target for the alleviation of neurological and neuropsychiatric symptoms using deep brain stimulation (DBS). STN-DBS is today applied as treatment in Parkinson´s disease, dystonia, essential tremor, and obsessive-compulsive disorder (OCD). STN-DBS also shows promise as a treatment for refractory Tourette syndrome. However, the internal organization of the STN has remained elusive and challenges researchers and clinicians: How can this small brain structure engage in the multitude of functions that renders it a key hub for therapeutic intervention of a variety of brain disorders ranging from motor to affective to cognitive? Based on recent gene expression studies of the STN, a comprehensive view of the anatomical and cellular organization, including revelations of spatio-molecular heterogeneity, is now possible to outline. In this review, we focus attention to the neurobiological architecture of the STN with specific emphasis on molecular patterns discovered within this complex brain area. Studies from human, non-human primate, and rodent brains now reveal anatomically defined distribution of specific molecular markers. Together their spatial patterns indicate a heterogeneous molecular architecture within the STN. Considering the translational capacity of targeting the STN in severe brain disorders, the addition of molecular profiling of the STN will allow for advancement in precision of clinical STN-based interventions.

WarpDrive: Improving spatial normalization using manual refinements

Spatial normalization-the process of mapping subject brain images to an average template brain-has evolved over the last 20+ years into a reliable method that facilitates the comparison of brain imaging results across patients, centers & modalities. While overall successful, sometimes, this automatic process yields suboptimal results, especially when dealing with brains with extensive neurodegeneration and atrophy patterns, or when high accuracy in specific regions is needed. Here we introduce WarpDrive, a novel tool for manual refinements of image alignment after automated registration. We show that the tool applied in a cohort of patients with Alzheimer's disease who underwent deep brain stimulation surgery helps create more accurate representations of the data as well as meaningful models to explain patient outcomes. The tool is built to handle any type of 3D imaging data, also allowing refinements in high-resolution imaging, including histology and multiple modalities to precisely aggregate multiple data sources together.

Echoes from Intrinsic Connectivity Networks in the Subcortex

Decades of research have greatly improved our understanding of intrinsic human brain organization in terms of functional networks and the transmodal hubs within the cortex at which they converge. However, substrates of multinetwork integration in the human subcortex are relatively uncharted. Here, we leveraged recent advances in subcortical atlasing and ultra-high field (7 T) imaging optimized for the subcortex to investigate the functional architecture of 14 individual structures in healthy adult males and females with a fully data-driven approach. We revealed that spontaneous neural activity in subcortical regions can be decomposed into multiple independent subsignals that correlate with, or "echo," the activity in functional networks across the cortex. Distinct subregions of the thalamus, striatum, claustrum, and hippocampus showed a varied pattern of echoes from attention, control, visual, somatomotor, and default mode networks, demonstrating evidence for a heterogeneous organization supportive of functional integration. Multiple network activity furthermore converged within the globus pallidus externa, substantia nigra, and ventral tegmental area but was specific to one subregion, while the amygdala and pedunculopontine nucleus preferentially affiliated with a single network, showing a more homogeneous topography. Subregional connectivity of the globus pallidus interna, subthalamic nucleus, red nucleus, periaqueductal gray, and locus coeruleus did not resemble patterns of cortical network activity. Together, these finding describe potential mechanisms through which the subcortex participates in integrated and segregated information processing and shapes the spontaneous cognitive dynamics during rest.SIGNIFICANCE STATEMENT Despite the impact of subcortical dysfunction on brain health and cognition, large-scale functional mapping of subcortical structures severely lags behind that of the cortex. Recent developments in subcortical atlasing and imaging at ultra-high field provide new avenues for studying the intricate functional architecture of the human subcortex. With a fully data-driven analysis, we reveal subregional connectivity profiles of a large set of noncortical structures, including those rarely studied in fMRI research. The results have implications for understanding how the functional organization of the subcortex facilitates integrative processing through cross-network information convergence, paving the way for future work aimed at improving our knowledge of subcortical contributions to intrinsic brain dynamics and spontaneous cognition.

Fine scale hippocampus morphology variation cross 552 healthy subjects from age 20 to 80

The cerebral cortex varies over the course of a person's life span: at birth, the surface is smooth, before becoming more bumpy (deeper sulci and thicker gyri) in middle age, and thinner in senior years. In this work, a similar phenomenon was observed on the hippocampus. It was previously believed the fine-scale morphology of the hippocampus could only be extracted only with high field scanners (7T, 9.4T); however, recent studies show that regular 3T MR scanners can be sufficient for this purpose. This finding opens the door for the study of fine hippocampal morphometry for a large amount of clinical data. In particular, a characteristic bumpy and subtle feature on the inferior aspect of the hippocampus, which we refer to as hippocampal dentation, presents a dramatic degree of variability between individuals from very smooth to highly dentated. In this report, we propose a combined method joining deep learning and sub-pixel level set evolution to efficiently obtain fine-scale hippocampal segmentation on 552 healthy subjects. Through non-linear dentation extraction and fitting, we reveal that the bumpiness of the inferior surface of the human hippocampus has a clear temporal trend. It is bumpiest between 40 and 50 years old. This observation should be aligned with neurodevelopmental and aging stages.

Pupil size reflects activation of subcortical ascending arousal system nuclei during rest

Neuromodulatory nuclei that are part of the ascending arousal system (AAS) play a crucial role in regulating cortical state and optimizing task performance. Pupil diameter, under constant luminance conditions, is increasingly used as an index of activity of these AAS nuclei. Indeed, task-based functional imaging studies in humans have begun to provide evidence of stimulus-driven pupil-AAS coupling. However, whether there is such a tight pupil-AAS coupling during rest is not clear. To address this question, we examined simultaneously acquired resting-state fMRI and pupil-size data from 74 participants, focusing on six AAS nuclei: the locus coeruleus, ventral tegmental area, substantia nigra, dorsal and median raphe nuclei, and cholinergic basal forebrain. Activation in all six AAS nuclei was optimally correlated with pupil size at 0-2 s lags, suggesting that spontaneous pupil changes were almost immediately followed by corresponding BOLD-signal changes in the AAS. These results suggest that spontaneous changes in pupil size that occur during states of rest can be used as a noninvasive general index of activity in AAS nuclei. Importantly, the nature of pupil-AAS coupling during rest appears to be vastly different from the relatively slow canonical hemodynamic response function that has been used to characterize task-related pupil-AAS coupling.

Affine transformation edited and refined deep neural network for quantitative susceptibility mapping

Deep neural networks have demonstrated great potential in solving dipole inversion for Quantitative Susceptibility Mapping (QSM). However, the performances of most existing deep learning methods drastically degrade with mismatched sequence parameters such as acquisition orientation and spatial resolution. We propose an end-to-end AFfine Transformation Edited and Refined (AFTER) deep neural network for QSM, which is robust against arbitrary acquisition orientation and spatial resolution up to 0.6 mm isotropic at the finest. The AFTER-QSM neural network starts with a forward affine transformation layer, followed by a Unet for dipole inversion, then an inverse affine transformation layer, followed by a Residual Dense Network (RDN) for QSM refinement. Simulation and in-vivo experiments demonstrated that the proposed AFTER-QSM network architecture had excellent generalizability. It can successfully reconstruct susceptibility maps from highly oblique and anisotropic scans, leading to the best image quality assessments in simulation tests and suppressed streaking artifacts and noise levels for in-vivo experiments compared with other methods. Furthermore, ablation studies showed that the RDN refinement network significantly reduced image blurring and susceptibility underestimation due to affine transformations. In addition, the AFTER-QSM network substantially shortened the reconstruction time from minutes using conventional methods to only a few seconds.

Spatial profiles provide sensitive MRI measures of the midbrain micro- and macrostructure

The midbrain is the rostral-most part of the brainstem. It contains numerous nuclei and white matter tracts, which are involved in motor, auditory and visual processing, and changes in their structure and function have been associated with aging, as well as neurodegenerative disorders. Current tools for estimating midbrain subregions and their structure with MRI require high resolution and multi-parametric quantitative MRI measures. We propose an approach that relies on morphology to calculate profiles along the midbrain and show these profiles are sensitive to the underlying macrostructure of the midbrain. First, we show that the midbrain structure can be sampled, within subject space, along three main axes of the left and right midbrain, producing profiles that are similar across subjects. We use two data sets with different field strengths, that contain R1, R2* and QSM maps and show that the profiles are highly correlated both across subjects and between datasets. Next, we compare profiles of the midbrain that sample ROIs, and show that the profiles along the first two axes sample the midbrain in a way that reliably separates the main structures, i.e., the substantia nigra, the red nucleus, and periaqueductal gray. We further show that age differences which are localized to specific nuclei, are reflected in the profiles. Finally, we generalize the same approach to calculate midbrain profiles on a third clinically relevant dataset using HCP subjects, with metrics such as the diffusion tensor and semi-quantitative data such as T1w/T2w maps. Our results suggest that midbrain profiles, both of quantitative and semi-quantitative estimates are sensitive to the underlying macrostructure of the midbrain. The midbrain profiles are calculated in native space, and rely on simple measurements. We show that it is robust and can be easily expanded to different datasets, and as such we hope that it will be of great use to the community and to the study of the midbrain in particular.

A unified model for reconstruction and R2* mapping of accelerated 7T data using the quantitative recurrent inference machine

Quantitative MRI (qMRI) acquired at the ultra-high field of 7 Tesla has been used in visualizing and analyzing subcortical structures. qMRI relies on the acquisition of multiple images with different scan settings, leading to extended scanning times. Data redundancy and prior information from the relaxometry model can be exploited by deep learning to accelerate the imaging process. We propose the quantitative Recurrent Inference Machine (qRIM), with a unified forward model for joint reconstruction and R₂^*-mapping from sparse data, embedded in a Recurrent Inference Machine (RIM), an iterative inverse problem-solving network. To study the dependency of the proposed extension of the unified forward model to network architecture, we implemented and compared a quantitative End-to-End Variational Network (qE2EVN). Experiments were performed with high-resolution multi-echo gradient echo data of the brain at 7T of a cohort study covering the entire adult life span. The error in reconstructed R₂^* from undersampled data relative to reference data significantly decreased for the unified model compared to sequential image reconstruction and parameter fitting using the RIM. With increasing acceleration factor, an increasing reduction in the reconstruction error was observed, pointing to a larger benefit for sparser data. Qualitatively, this was following an observed reduction of image blurriness in R₂^*-maps. In contrast, when using the U-Net as network architecture, a negative bias in R₂^* in selected regions of interest was observed. Compressed Sensing rendered accurate, but less precise estimates of R₂^*. The qE2EVN showed slightly inferior reconstruction quality compared to the qRIM but better quality than the U-Net and Compressed Sensing. Subcortical maturation over age measured by a linearly increasing interquartile range of R₂^* in the striatum was preserved up to an acceleration factor of 9. With the integrated prior of the unified forward model, the proposed qRIM can exploit the redundancy among repeated measurements and shared information between tasks, facilitating relaxometry in accelerated MRI.

Susceptibility networks reveal independent patterns of brain iron abnormalities in multiple sclerosis

Brain iron homeostasis is necessary for healthy brain function. MRI and histological studies have shown altered brain iron levels in the brains of patients with multiple sclerosis (MS), particularly in the deep gray matter (DGM). Previous studies were able to only partially separate iron-modifying effects because of incomplete knowledge of iron-modifying processes and influencing factors. It is therefore unclear to what extent and at which stages of the disease different processes contribute to brain iron changes. We postulate that spatially covarying magnetic susceptibility networks determined with Independent Component Analysis (ICA) reflect, and allow for the study of, independent processes regulating iron levels. We applied ICA to quantitative susceptibility maps for 170 individuals aged 9-81 years without neurological disease ("Healthy Aging" (HA) cohort), and for a cohort of 120 patients with MS and 120 age- and sex-matched healthy controls (HC; together the "MS/HC" cohort). Two DGM-associated "susceptibility networks" identified in the HA cohort (the Dorsal Striatum and Globus Pallidus Interna Networks) were highly internally reproducible (i.e. "robust") across multiple ICA repetitions on cohort subsets. DGM areas overlapping both robust networks had higher susceptibility levels than DGM areas overlapping only a single robust network, suggesting that these networks were caused by independent processes of increasing iron concentration. Because MS is thought to accelerate brain aging, we hypothesized that associations between age and the two robust DGM-associated networks would be enhanced in patients with MS. However, only one of these networks was altered in patients with MS, and it had a null age association in patients with MS rather than a stronger association. Further analysis of the MS/HC cohort revealed three additional disease-related networks (the Pulvinar, Mesencephalon, and Caudate Networks) that were differentially altered between patients with MS and HCs and between MS subtypes. Exploratory regression analyses of the disease-related networks revealed differential associations with disease duration and T2 lesion volume. Finally, analysis of ROI-based disease effects in the MS/HC cohort revealed an effect of disease status only in the putamen ROI and exploratory regression analysis did not show associations between the caudate and pulvinar ROIs and disease duration or T2 lesion volume, showing the ICA-based approach was more sensitive to disease effects. These results suggest that the ICA network framework increases sensitivity for studying patterns of brain iron change, opening a new avenue for understanding brain iron physiology under normal and disease conditions.

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.

Freezing of gait in Parkinson’s disease is associated with the microstructural and functional changes of globus pallidus internus

Background

Freezing of gait (FOG) is a common motor symptom in advanced Parkinson’s disease (PD). However, the pathophysiology mechanism of FOG is not fully understood. The purpose of this study was to investigate microstructural abnormalities in subcortical gray matter and alterations in functional connectivity of the nuclei with microstructural changes. In addition, the correlations between these microstructural and functional changes and the severity of FOG were measured.

Materials and methods

Twenty-four patients with FOG (PD-FOG), 22 PD patients without FOG (PD-nFOG), and 27 healthy controls (HC) were recruited. FOG Questionnaire (FOGQ) and Gait and Falling Questionnaire (GFQ) were assessed, and Timed Up and Go (TUG) tests were performed in PD-FOG patients. All subjects underwent diffusion tensor imaging (DTI) and resting-state functional MRI scanning. The DTI measures, including fractional anisotropy (FA), mean diffusivity (MD), radial diffusivity (RD), and axial diffusivity (AD), were extracted and measured from basal ganglia, thalamus, and substantia nigra. The nuclei with microstructural alterations were selected as seed regions to perform the seed-based resting-state functional connectivity.

Results

The MD and RD values of the right globus pallidus internus (GPi) were significantly higher in patients with PD-FOG compared with PD-nFOG patients and HC. In PD-FOG patients, the MD and RD values of the right GPi were significantly correlated with the time of the TUG test in both ON and OFF states. The MD values were also correlated with the GFQ scores in PD-FOG patients. Resting-state functional connectivity between the right GPi and left middle occipital gyri decreased significantly in PD-FOG patients compared to PD-nFOG patients, and was negatively correlated with GFQ scores as well as the time of ON state TUG in PD-FOG patients.

Conclusion

Microstructural alterations in the right GPi and functional connectivity between the right GPi and visual cortex may be associated with the pathophysiological mechanisms of FOG in PD patients.

Quantitative susceptibility mapping as an imaging biomarker for Alzheimer’s disease: The expectations and limitations

Alzheimer’s disease (AD) is the most common type of dementia and a distressing diagnosis for individuals and caregivers. Researchers and clinical trials have mainly focused on β-amyloid plaques, which are hypothesized to be one of the most important factors for neurodegeneration in AD. Meanwhile, recent clinicopathological and radiological studies have shown closer associations of tau pathology rather than β-amyloid pathology with the onset and progression of Alzheimer’s symptoms. Toward a biological definition of biomarker-based research framework for AD, the 2018 National Institute on Aging–Alzheimer’s Association working group has updated the ATN classification system for stratifying disease status in accordance with relevant pathological biomarker profiles, such as cerebral β-amyloid deposition, hyperphosphorylated tau, and neurodegeneration. In addition, altered iron metabolism has been considered to interact with abnormal proteins related to AD pathology thorough generating oxidative stress, as some prior histochemical and histopathological studies supported this iron-mediated pathomechanism. Quantitative susceptibility mapping (QSM) has recently become more popular as a non-invasive magnetic resonance technique to quantify local tissue susceptibility with high spatial resolution, which is sensitive to the presence of iron. The association of cerebral susceptibility values with other pathological biomarkers for AD has been investigated using various QSM techniques; however, direct evidence of these associations remains elusive. In this review, we first briefly describe the principles of QSM. Second, we focus on a large variety of QSM applications, ranging from common applications, such as cerebral iron deposition, to more recent applications, such as the assessment of impaired myelination, quantification of venous oxygen saturation, and measurement of blood– brain barrier function in clinical settings for AD. Third, we mention the relationships among QSM, established biomarkers, and cognitive performance in AD. Finally, we discuss the role of QSM as an imaging biomarker as well as the expectations and limitations of clinically useful diagnostic and therapeutic implications for AD.

Nigral neuropathology of Parkinson's motor subtypes coincide with circuitopathies: a scoping review

The neuropathological substrates of Parkinson’s disease (PD) patients with motor subtypes tremor-dominance (TD), non-tremor dominance (nTD), postural instability and gait difficulty (PIGD), and akinetic-rigid (AR) are not completely differentiated. While extensive pathological research has been conducted on neuronal tissue of PD patients, data have not been discussed in the context of mechanistic circuitry theories differentiating motor subtypes. It is, therefore, expected that a more specific and tailored management of PD symptoms can be accomplished by understanding symptom-specific neuropathological mechanisms with the detail histology can provide. This scoping review gives an overview of the literature comparing TD and nTD PD motor subtypes by clarify observed pathology with underlying physiological circuitry theories. Studies using an array of pathological examination techniques have shown significant differences between TD and nTD PD subtypes. nTD PD patients show higher neuronal loss, gliosis, extraneuronal melanin deposits, and neuroaxonal dystrophy in multiple subregions of the substantia nigra (SN) related to the overactivity of the indirect motor loop. TD patients show more severe cell loss specifically in medial SN subdivisions, and have damage in the retrorubral field A-8 that projects to the dorsolateral striatum and ventromedial thalamus in the direct motor loop. Pathological studies are consistent with neuroimaging data and support contemporary mechanistic circuitry theories of PD motor symptom genesis. Further multimodal neuroimaging and histological studies are required to validate and expand upon these findings.

Establishing a Baseline for Human Cortical Folding Morphological Variables: A Multisite Study

Differences in the way human cerebral cortices fold have been correlated to health, disease, development, and aging. However, to obtain a deeper understanding of the mechanisms that generate such differences, it is useful to derive one's morphometric variables from the first principles. This study explores one such set of variables that arise naturally from a model for universal self-similar cortical folding that was validated on comparative neuroanatomical data. We aim to establish a baseline for these variables across the human lifespan using a heterogeneous compilation of cross-sectional datasets as the first step to extending the model to incorporate the time evolution of brain morphology. We extracted the morphological features from structural MRI of 3,650 subjects: 3,095 healthy controls (CTL) and 555 patients with Alzheimer's Disease (AD) from 9 datasets, which were harmonized with a straightforward procedure to reduce the uncertainty due to heterogeneous acquisition and processing. The unprecedented possibility of analyzing such a large number of subjects in this framework allowed us to compare CTL and AD subjects' lifespan trajectories, testing if AD is a form of accelerated aging at the brain structural level. After validating this baseline from development to aging, we estimate the variables' uncertainties and show that Alzheimer's Disease is similar to premature aging when measuring global and local degeneration. This new methodology may allow future studies to explore the structural transition between healthy and pathological aging and may be essential to generate data for the cortical folding process simulations.

Mapping microstructural gradients of the human striatum in normal aging and Parkinson's disease

Mapping structural spatial change (i.e., gradients) in the striatum is essential for understanding the function of the basal ganglia in both health and disease. We developed a method to identify and quantify gradients of microstructure in the single human brain in vivo. We found spatial gradients in the putamen and caudate nucleus of the striatum that were robust across individuals, clinical conditions, and datasets. By exploiting multiparametric quantitative MRI, we found distinct, spatially dependent, aging-related alterations in water content and iron concentration. Furthermore, we found cortico-striatal microstructural covariation, showing relations between striatal structural gradients and cortical hierarchy. In Parkinson's disease (PD) patients, we found abnormal gradients in the putamen, revealing changes in the posterior putamen that explain patients' dopaminergic loss and motor dysfunction. Our work provides a noninvasive approach for studying the spatially varying, structure-function relationship in the striatum in vivo, in normal aging and PD.

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.

A high-resolution multi-shell 3T diffusion magnetic resonance imaging dataset as part of the Amsterdam Ultra-high field adult lifespan database (AHEAD)

In order to further our understanding of brain function and the underlying networks, more advanced diffusion weighted magnetic resonance imaging (DWI MRI) data are essential. Here we present freely available high-resolution multi-shell multi-directional 3 Tesla (T) DWI MRI data as part of the ‘Amsterdam Ultra-high field adult lifespan database’ (AHEAD). The 3T DWI AHEAD dataset include 1.28mm isotropic whole brain DWI data of 49 healthy adult participants between 18 and 90 years old. The acquired data include DWIs at three non-zero b-values (48 directions, b-value 700 s/mm²; 56 directions, b-value 1000 s/mm²; 64 directions, b-value 1600 s/mm²) including a total of twelve volumes with a b-value of 0 s/mm² (b0 volumes). In addition, eight b0 volumes with a reversed phase encoding direction were acquired to correct for distortions. To facilitate future use, the DWI data have been denoised, corrected for eddy currents, susceptibility-induced off-resonance field distortions, bias fields, and are skull stripped.

Catching Wandering Minds with Tapping Fingers: Neural and Behavioral Insights into Task-unrelated Cognition

When the human mind wanders, it engages in episodes during which attention is focused on self-generated thoughts rather than on external task demands. Although the sustained attention to response task is commonly used to examine relationships between mind wandering and executive functions, limited executive resources are required for optimal task performance. In the current study, we aimed to investigate the relationship between mind wandering and executive functions more closely by employing a recently developed finger-tapping task to monitor fluctuations in attention and executive control through task performance and periodical experience sampling during concurrent functional magnetic resonance imaging (fMRI) and pupillometry. Our results show that mind wandering was preceded by increases in finger-tapping variability, which was correlated with activity in dorsal and ventral attention networks. The entropy of random finger-tapping sequences was related to activity in frontoparietal regions associated with executive control, demonstrating the suitability of this paradigm for studying executive functioning. The neural correlates of behavioral performance, pupillary dynamics, and self-reported attentional state diverged, thus indicating a dissociation between direct and indirect markers of mind wandering. Together, the investigation of these relationships at both the behavioral and neural level provided novel insights into the identification of underlying mechanisms of mind wandering.

Charting human subcortical maturation across the adult lifespan with in vivo 7 T MRI

The human subcortex comprises hundreds of unique structures. Subcortical functioning is crucial for behavior, and disrupted function is observed in common neurodegenerative diseases. Despite their importance, human subcortical structures continue to be difficult to study in vivo. Here we provide a detailed account of 17 prominent subcortical structures and ventricles, describing their approximate iron and myelin contents, morphometry, and their age-related changes across the normal adult lifespan. The results provide compelling insights into the heterogeneity and intricate age-related alterations of these structures. They also show that the locations of many structures shift across the lifespan, which is of direct relevance for the use of standard magnetic resonance imaging atlases. The results further our understanding of subcortical morphometry and neuroimaging properties, and of normal aging processes which ultimately can improve our understanding of neurodegeneration.

CEREBRUM-7T: Fast and Fully Volumetric Brain Segmentation of 7 Tesla MR Volumes

Ultra-high-field magnetic resonance imaging (MRI) enables sub-millimetre resolution imaging of the human brain, allowing the study of functional circuits of cortical layers at the meso-scale. An essential step in many functional and structural neuroimaging studies is segmentation, the operation of partitioning the MR images in anatomical structures. Despite recent efforts in brain imaging analysis, the literature lacks in accurate and fast methods for segmenting 7-tesla (7T) brain MRI. We here present CEREBRUM-7T, an optimised end-to-end convolutional neural network, which allows fully automatic segmentation of a whole 7T T1_w MRI brain volume at once, without partitioning the volume, pre-processing, nor aligning it to an atlas. The trained model is able to produce accurate multi-structure segmentation masks on six different classes plus background in only a few seconds. The experimental part, a combination of objective numerical evaluations and subjective analysis, confirms that the proposed solution outperforms the training labels it was trained on and is suitable for neuroimaging studies, such as layer functional MRI studies. Taking advantage of a fine-tuning operation on a reduced set of volumes, we also show how it is possible to effectively apply CEREBRUM-7T to different sites data. Furthermore, we release the code, 7T data, and other materials, including the training labels and the Turing test.

Manual delineation approaches for direct imaging of the subcortex

The growing interest in the human subcortex is accompanied by an increasing number of parcellation procedures to identify deep brain structures in magnetic resonance imaging (MRI) contrasts. Manual procedures continue to form the gold standard for parcellating brain structures and is used for the validation of automated approaches. Performing manual parcellations is a tedious process which requires a systematic and reproducible approach. For this purpose, we created a series of protocols for the anatomical delineation of 21 individual subcortical structures. The intelligibility of the protocols was assessed by calculating Dice similarity coefficients for ten healthy volunteers. In addition, dilated Dice coefficients showed that manual parcellations created using these protocols can provide high-quality training data for automated algorithms. Here, we share the protocols, together with three example MRI datasets and the created manual delineations. The protocols can be applied to create high-quality training data for automated parcellation procedures, as well as for further validation of existing procedures and are shared without restrictions with the research community.

Quantitative Susceptibility Mapping of the Basal Ganglia and Thalamus at 9.4 Tesla

The thalamus (Th) and basal ganglia (BG) are central subcortical connectivity hubs of the human brain, whose functional anatomy is still under intense investigation. Nevertheless, both substructures contain a robust and reproducible functional anatomy. The quantitative susceptibility mapping (QSM) at ultra-high field may facilitate an improved characterization of the underlying functional anatomy in vivo. We acquired high-resolution QSM data at 9.4 Tesla in 21 subjects, and analyzed the thalamic and BG by using a prior defined functional parcellation. We found a more substantial contribution of paramagnetic susceptibility sources such as iron in the pallidum in contrast to the caudate, putamen, and Th in descending order. The diamagnetic susceptibility sources such as myelin and calcium revealed significant contributions in the Th parcels compared with the BG. This study presents a detailed nuclei-specific delineation of QSM-provided diamagnetic and paramagnetic susceptibility sources pronounced in the BG and the Th. We also found a reasonable interindividual variability as well as slight hemispheric differences. The results presented here contribute to the microstructural knowledge of the Th and the BG. In specific, the study illustrates QSM values (myelin, calcium, and iron) in functionally similar subregions of the Th and the BG.

Cortical and subcortical contributions to interference resolution and inhibition - An fMRI ALE meta-analysis

Interacting with our environment requires the selection of appropriate responses and the inhibition of others. Such effortful inhibition is achieved by a number of interference resolution and global inhibition processes. This meta-analysis including 57 studies and 73 contrasts revisits the overlap and differences in brain areas supporting interference resolution and global inhibition in cortical and subcortical brain areas. Activation likelihood estimation was used to discern the brain regions subserving each type of cognitive control. Individual contrast analysis revealed a common activation of the bilateral insula and supplementary motor areas. Subtraction analyses demonstrated the voxel-wise differences in recruitment in a number of areas including the precuneus in the interference tasks and the frontal pole and dorsal striatum in the inhibition tasks. Our results display a surprising lack of subcortical involvement within these types of cognitive control, a finding that is likely to reflect a systematic gap in the field of functional neuroimaging.

Quantity and quality: Normative open-access neuroimaging databases

The focus of this article is to compare twenty normative and open-access neuroimaging databases based on quantitative measures of image quality, namely, signal-to-noise (SNR) and contrast-to-noise ratios (CNR). We further the analysis through discussing to what extent these databases can be used for the visualization of deeper regions of the brain, such as the subcortex, as well as provide an overview of the types of inferences that can be drawn. A quantitative comparison of contrasts including T1-weighted (T1w) and T2-weighted (T2w) images are summarized, providing evidence for the benefit of ultra-high field MRI. Our analysis suggests a decline in SNR in the caudate nuclei with increasing age, in T1w, T2w, qT1 and qT2* contrasts, potentially indicative of complex structural age-dependent changes. A similar decline was found in the corpus callosum of the T1w, qT1 and qT2* contrasts, though this relationship is not as extensive as within the caudate nuclei. These declines were accompanied by a declining CNR over age in all image contrasts. A positive correlation was found between scan time and the estimated SNR as well as a negative correlation between scan time and spatial resolution. Image quality as well as the number and types of contrasts acquired by these databases are important factors to take into account when selecting structural data for reuse. This article highlights the opportunities and pitfalls associated with sampling existing databases, and provides a quantitative backing for their usage.

A probabilistic atlas of the human ventral tegmental area (VTA) based on 7 Tesla MRI data

Functional magnetic resonance imaging (fMRI) BOLD signal is commonly localized by using neuroanatomical atlases, which can also serve for region of interest analyses. Yet, the available MRI atlases have serious limitations when it comes to imaging subcortical structures: only 7% of the 455 subcortical nuclei are captured by current atlases. This highlights the general difficulty in mapping smaller nuclei deep in the brain, which can be addressed using ultra-high field 7 Tesla (T) MRI. The ventral tegmental area (VTA) is a subcortical structure that plays a pivotal role in reward processing, learning and memory. Despite the significant interest in this nucleus in cognitive neuroscience, there are currently no available, anatomically precise VTA atlases derived from 7 T MRI data that cover the full region of the VTA. Here, we first provide a protocol for multimodal VTA imaging and delineation. We then provide a data description of a probabilistic VTA atlas based on in vivo 7 T MRI data.

Multi-contrast anatomical subcortical structures parcellation

The human subcortex is comprised of more than 450 individual nuclei which lie deep in the brain. Due to their small size and close proximity, up until now only 7% have been depicted in standard MRI atlases. Thus, the human subcortex can largely be considered as terra incognita. Here, we present a new open-source parcellation algorithm to automatically map the subcortex. The new algorithm has been tested on 17 prominent subcortical structures based on a large quantitative MRI dataset at 7 Tesla. It has been carefully validated against expert human raters and previous methods, and can easily be extended to other subcortical structures and applied to any quantitative MRI dataset. In sum, we hope this novel parcellation algorithm will facilitate functional and structural neuroimaging research into small subcortical nuclei and help to chart terra incognita.

3 versus 7 Tesla magnetic resonance imaging for parcellations of subcortical brain structures in clinical settings

7 Tesla (7T) magnetic resonance imaging holds great promise for improved visualization of the human brain for clinical purposes. To assess whether 7T is superior regarding localization procedures of small brain structures, we compared manual parcellations of the red nucleus, subthalamic nucleus, substantia nigra, globus pallidus interna and externa. These parcellations were created on a commonly used clinical anisotropic clinical 3T with an optimized isotropic (o)3T and standard 7T scan. The clinical 3T MRI scans did not allow delineation of an anatomically plausible structure due to its limited spatial resolution. o3T and 7T parcellations were directly compared. We found that 7T outperformed the o3T MRI as reflected by higher Dice scores, which were used as a measurement of interrater agreement for manual parcellations on quantitative susceptibility maps. This increase in agreement was associated with higher contrast to noise ratios for smaller structures, but not for the larger globus pallidus segments. Additionally, control-analyses were performed to account for potential biases in manual parcellations by assessing semi-automatic parcellations. These results showed a higher consistency for structure volumes for 7T compared to optimized 3T which illustrates the importance of the use of isotropic voxels for 3D visualization of the surgical target area. Together these results indicate that 7T outperforms c3T as well as o3T given the constraints of a clinical setting.

Methodological Considerations for Neuroimaging in Deep Brain Stimulation of the Subthalamic Nucleus in Parkinson's Disease Patients

Deep brain stimulation (DBS) of the subthalamic nucleus is a neurosurgical intervention for Parkinson's disease patients who no longer appropriately respond to drug treatments. A small fraction of patients will fail to respond to DBS, develop psychiatric and cognitive side-effects, or incur surgery-related complications such as infections and hemorrhagic events. In these cases, DBS may require recalibration, reimplantation, or removal. These negative responses to treatment can partly be attributed to suboptimal pre-operative planning procedures via direct targeting through low-field and low-resolution magnetic resonance imaging (MRI). One solution for increasing the success and efficacy of DBS is to optimize preoperative planning procedures via sophisticated neuroimaging techniques such as high-resolution MRI and higher field strengths to improve visualization of DBS targets and vasculature. We discuss targeting approaches, MRI acquisition, parameters, and post-acquisition analyses. Additionally, we highlight a number of approaches including the use of ultra-high field (UHF) MRI to overcome limitations of standard settings. There is a trade-off between spatial resolution, motion artifacts, and acquisition time, which could potentially be dissolved through the use of UHF-MRI. Image registration, correction, and post-processing techniques may require combined expertise of traditional radiologists, clinicians, and fundamental researchers. The optimization of pre-operative planning with MRI can therefore be best achieved through direct collaboration between researchers and clinicians.

Sharpness in motion corrected quantitative imaging at 7T

Sub-millimeter imaging at 7T has opened new possibilities for qualitatively and quantitatively studying brain structure as it evolves throughout the life span. However, subject motion introduces image blurring on the order of magnitude of the spatial resolution and is thus detrimental to image quality. Such motion can be corrected for, but widespread application has not yet been achieved and quantitative evaluation is lacking. This raises a need to quantitatively measure image sharpness throughout the brain. We propose a method to quantify sharpness of brain structures at sub-voxel resolution, and use it to assess to what extent limited motion is related to image sharpness. The method was evaluated in a cohort of 24 healthy volunteers with a wide and uniform age range, aiming to arrive at results that largely generalize to larger populations. Using 3D fat-excited motion navigators, quantitative R₁, R₂^* and Quantitative Susceptibility Maps and T₁-weighted images were retrospectively corrected for motion. Sharpness was quantified in all modalities for selected regions of interest (ROI) by fitting the sigmoidally shaped error function to data within locally homogeneous clusters. A strong, almost linear correlation between motion and sharpness improvement was observed, and motion correction significantly improved sharpness. Overall, the Full Width at Half Maximum reduced from 0.88 mm to 0.70 mm after motion correction, equivalent to a 2.0 times smaller voxel volume. Motion and sharpness were not found to correlate with the age of study participants. We conclude that in our data, motion correction using fat navigators is overall able to restore the measured sharpness to the imaging resolution, irrespective of the amount of motion observed during scanning.