Bayesian Parameter Estimation and Segmentation in the Multi-Atlas Random Orbit Model

On the optimization of 3D inflow-based vascular-space-occupancy (iVASO) MRI for the quantification of arterial cerebral blood volume (CBVa)

PURPOSE

The inflow-based vascular-space-occupancy (iVASO) MRI was originally developed in a single-slice mode to measure arterial cerebral blood volume (CBVa). When vascular crushers are applied in iVASO, the signals can be sensitized predominantly to small pial arteries and arterioles. The purpose of this study is to perform a systematic optimization and evaluation of a 3D iVASO sequence on both 3 T and 7 T for the quantification of CBVa values in the human brain.

METHODS

Three sets of experiments were performed in three separate cohorts. (1) 3D iVASO MRI protocols were compared to single-slice iVASO, and the reproducibility of whole-brain 3D iVASO MRI was evaluated. (2) The effects from different vascular crushers in iVASO were assessed. (3) 3D iVASO MRI results were evaluated in arterial and venous blood vessels identified using ultrasmall-superparamagnetic-iron-oxides-enhanced MRI to validate its arterial origin.

RESULTS

3D iVASO scans showed signal-to-noise ratio (SNR) and CBVa measures consistent with single-slice iVASO with reasonable intrasubject reproducibility. Among the iVASO scans performed with different vascular crushers, the whole-brain 3D iVASO scan with a motion-sensitized-driven-equilibrium preparation with two binomial refocusing pulses and an effective TE of 50 ms showed the best suppression of macrovascular signals, with a relatively low specific absorption rate. When no vascular crusher was applied, the CBVa maps from 3D iVASO scans showed large CBVa values in arterial vessels but well-suppressed signals in venous vessels.

CONCLUSION

A whole-brain 3D iVASO MRI scan was optimized for CBVa measurement in the human brain. When only microvascular signals are desired, a motion-sensitized-driven-equilibrium-based vascular crusher with binomial refocusing pulses can be applied in 3D iVASO.

Baseline functional connectivity predicts who will benefit from neuromodulation: evidence from primary progressive aphasia

Background

Identifying the characteristics of individuals who demonstrate response to an intervention allows us to predict who is most likely to benefit from certain interventions. Prediction is challenging in rare and heterogeneous diseases, such as primary progressive aphasia (PPA), that have varying clinical manifestations. We aimed to determine the characteristics of those who will benefit most from transcranial direct current stimulation (tDCS) of the left inferior frontal gyrus (IFG) using a novel heterogeneity and group identification analysis.

Methods

We compared the predictive ability of demographic and clinical patient characteristics (e.g., PPA variant and disease progression, baseline language performance) vs. functional connectivity alone (from resting-state fMRI) in the same cohort.

Results

Functional connectivity alone had the highest predictive value for outcomes, explaining 62% and 75% of tDCS effect of variance in generalization (semantic fluency) and in the trained outcome of the clinical trial (written naming), contrasted with <15% predicted by clinical characteristics, including baseline language performance. Patients with higher baseline functional connectivity between the left IFG (opercularis and triangularis), and between the middle temporal pole and posterior superior temporal gyrus, were most likely to benefit from tDCS.

Conclusions

We show the importance of a baseline 7-minute functional connectivity scan in predicting tDCS outcomes, and point towards a precision medicine approach in neuromodulation studies. The study has important implications for clinical trials and practice, providing a statistical method that addresses heterogeneity in patient populations and allowing accurate prediction and enrollment of those who will most likely benefit from specific interventions.

Regional differences in the link between water exchange rate across the blood-brain barrier and cognitive performance in normal aging

The blood-brain barrier (BBB) undergoes functional changes with aging which may contribute to cognitive decline. A novel, diffusion prepared arterial spin labeling-based MRI technique can measure the rate of water exchange across the BBB (kw) and may thus be sensitive to age-related alterations in water exchange at the BBB. However, studies investigating relationships between kw and cognition have reported different directions of association. Here, we begin to investigate the direction of associations between kw and cognition in different brain regions, and their possible underpinnings, by evaluating links between kw, cognitive performance, and MRI markers of cerebrovascular dysfunction and/or damage. Forty-seven healthy older adults (age range 61-84) underwent neuroimaging to obtain whole-brain measures of kw, cerebrovascular reactivity (CVR), and white matter hyperintensity (WMH) volumes. Additionally, participants completed uniform data set (Version 3) neuropsychological tests of executive function (EF) and episodic memory (MEM). Voxel-wise linear regressions were conducted to test associations between kw and cognitive performance, CVR, and WMH volumes. We found that kw in the frontoparietal brain regions was positively associated with cognitive performance but not with CVR or WMH volumes. Conversely, kw in the basal ganglia was negatively associated with cognitive performance and CVR and positively associated with regional, periventricular WMH volume. These regionally dependent associations may relate to different physiological underpinnings in the relationships between kw and cognition in neocortical versus subcortical brain regions in older adults.

A Majorization-Minimization Algorithm for Neuroimage Registration

Intensity-based image registration is critical for neuroimaging tasks, such as 3D reconstruction, times-series alignment, and common coordinate mapping. The gradient-based optimization methods commonly used to solve this problem require a careful selection of step-length. This limitation imposes substantial time and computational costs. Here we propose a gradient-independent rigid-motion registration algorithm based on the majorization-minimization (MM) principle. Each iteration of our intensity-based MM algorithm reduces to a simple point-set rigid registration problem with a closed form solution that avoids the step-length issue altogether. The details of the algorithm are presented, and an error bound for its more practical truncated form is derived. The performance of the MM algorithm is shown to be more effective than gradient descent on simulated images and Nissl stained coronal slices of mouse brain. We also compare and contrast the similarities and differences between the MM algorithm and another gradient-free registration algorithm called the block-matching method. Finally, extensions of this algorithm to more complex problems are discussed.

Accurate measurement of magnetic resonance parkinsonism index by a fully automatic and deep learning quantification pipeline

OBJECTIVES

This study aims at a fully automatic pipeline for measuring the magnetic resonance parkinsonism index (MRPI) using deep learning methods.

METHODS

MRPI is defined as the product of the pons area to the midbrain area ratio and the middle cerebellar peduncle (MCP) width to the superior cerebellar peduncle (SCP) width ratio. In our proposed pipeline, we first used nnUNet to segment the brainstem and then employed HRNet to identify two key boundary points so as to sub-divide the whole brainstem into midbrain and pons. HRNet was also employed to predict the MCP endpoints for measuring the MCP width. Finally, we segmented the SCP on an oblique coronal plane and calculated its width. A total of 400 T1-weighted magnetic resonance images (MRIs) were used to train the nnUNet and HRNet models. Five-fold cross-validation was conducted to evaluate our proposed pipeline's performance on the training dataset. We also evaluated the performance of our proposed pipeline on three external datasets. Two of them had two raters manually measuring the MRPI values, providing insights into automatic accuracy versus inter-rater variability.

RESULTS

We obtained average absolute percentage errors (APEs) of 17.21%, 18.17%, 20.83%, and 22.83% on the training dataset and the three external validation datasets, while the inter-rater average APE measured on the first two external validation datasets was 11.31%. Our proposed pipeline significantly improved the MRPI quantification accuracy over a representative state-of-the-art traditional approach (p < 0.001).

CONCLUSION

The proposed automatic pipeline can accurately predict MRPI that is comparable with manual measurement.

CLINICAL RELEVANCE STATEMENT

This study presents an automated magnetic resonance parkinsonism index measurement tool that can analyze large amounts of magnetic resonance images, enhance the efficiency of Parkinsonism-Plus syndrome diagnosis, reduce the workload of clinicians, and minimize the impact of human factors on diagnosis.

KEY POINTS

• We propose an automatic pipeline for measuring the magnetic resonance parkinsonism index from magnetic resonance images. • The effectiveness of the proposed pipeline is successfully established on multiple datasets and comparisons with inter-rater measurements. • The proposed pipeline significantly outperforms a state-of-the-art quantification approach, being much closer to ground truth.

Identifying brain hierarchical structures associated with Alzheimer's disease using a regularized regression method with tree predictors

Brain segmentation at different levels is generally represented as hierarchical trees. Brain regional atrophy at specific levels was found to be marginally associated with Alzheimer's disease outcomes. In this study, we propose an ℓ1 -type regularization for predictors that follow a hierarchical tree structure. Considering a tree as a directed acyclic graph, we interpret the model parameters from a path analysis perspective. Under this concept, the proposed penalty regulates the total effect of each predictor on the outcome. With regularity conditions, it is shown that under the proposed regularization, the estimator of the model coefficient is consistent in ℓ2 -norm and the model selection is also consistent. When applied to a brain sMRI dataset acquired from the Alzheimer's Disease Neuroimaging Initiative (ADNI), the proposed approach identifies brain regions where atrophy in these regions demonstrates the declination in memory. With regularization on the total effects, the findings suggest that the impact of atrophy on memory deficits is localized from small brain regions, but at various levels of brain segmentation. Data used in preparation of this paper were obtained from the ADNI database.

Longitudinal imaging highlights preferential basal ganglia circuit atrophy in Huntington's disease

Huntington's disease is caused by a CAG repeat expansion in the Huntingtin gene (HTT), coding for polyglutamine in the Huntingtin protein, with longer CAG repeats causing earlier age of onset. The variable 'Age' × ('CAG'-L), where 'Age' is the current age of the individual, 'CAG' is the repeat length and L is a constant (reflecting an approximation of the threshold), termed the 'CAG Age Product' (CAP) enables the consideration of many individuals with different CAG repeat expansions at the same time for analysis of any variable and graphing using the CAG Age Product score as the X axis. Structural MRI studies have showed that progressive striatal atrophy begins many years prior to the onset of diagnosable motor Huntington's disease, confirmed by longitudinal multicentre studies on three continents, including PREDICT-HD, TRACK-HD and IMAGE-HD. However, previous studies have not clarified the relationship between striatal atrophy, atrophy of other basal ganglia structures, and atrophy of other brain regions. The present study has analysed all three longitudinal datasets together using a single image segmentation algorithm and combining data from a large number of subjects across a range of CAG Age Product score. In addition, we have used a strategy of normalizing regional atrophy to atrophy of the whole brain, in order to determine which regions may undergo preferential degeneration. This made possible the detailed characterization of regional brain atrophy in relation to CAG Age Product score. There is dramatic selective atrophy of regions involved in the basal ganglia circuit-caudate, putamen, nucleus accumbens, globus pallidus and substantia nigra. Most other regions of the brain appear to have slower but steady degeneration. These results support (but certainly do not prove) the hypothesis of circuit-based spread of pathology in Huntington's disease, possibly due to spread of mutant Htt protein, though other connection-based mechanisms are possible. Therapeutic targets related to prion-like spread of pathology or other mechanisms may be suggested. In addition, they have implications for current neurosurgical therapeutic approaches, since delivery of therapeutic agents solely to the caudate and putamen may miss other structures affected early, such as nucleus accumbens and output nuclei of the striatum, the substantia nigra and the globus pallidus.

Mixed longitudinal and cross-sectional analyses of deep gray matter and white matter using diffusion weighted images in premanifest and manifest Huntington's disease

Changes in the brain of patients with Huntington's disease (HD) begin years before clinical onset, so it remains critical to identify biomarkers to track these early changes. Metrics derived from tensor modeling of diffusion-weighted MRIs (DTI), that indicate the microscopic brain structure, can add important information to regional volumetric measurements. This study uses two large-scale longitudinal, multicenter datasets, PREDICT-HD and IMAGE-HD, to trace changes in DTI of HD participants with a broad range of CAP scores (a product of CAG repeat expansion and age), including those with pre-manifest disease (i.e., prior to clinical onset). Utilizing a fully automated data-driven approach to study the whole brain divided in regions of interest, we traced changes in DTI metrics (diffusivity and fractional anisotropy) versus CAP scores, using sigmoidal and linear regression models. We identified points of inflection in the sigmoidal regression using change-point analysis. The deep gray matter showed more evident and earlier changes in DTI metrics over CAP scores, compared to the deep white matter. In the deep white matter, these changes were more evident and occurred earlier in superior and posterior areas, compared to anterior and inferior areas. The curves of mean diffusivity vs. age of HD participants within a fixed CAP score were different from those of controls, indicating that the disease has an additional effect to age on the microscopic brain structure. These results show the regional and temporal vulnerability of the white matter and deep gray matter in HD, with potential implications for experimental therapeutics.

Brain regions associated with olfactory dysfunction in first episode psychosis patients

OBJECTIVES

Olfactory dysfunction is reproducibly reported in psychotic disorders, particularly in association with negative symptoms. The superior frontal gyrus (SFG) has been frequently studied in patients with psychotic disorders, in particular with their associations with negative symptoms. The relationship between olfactory functions and brain structure has been studied in healthy controls (HCs). Nevertheless, the studies with patients with psychotic disorders are limited. Here we report the olfactory-brain relationship in a first episode psychosis (FEP) cohort through both hypothesis-driven (centred on the SFG) and data-driven approaches.

METHODS

Using data from 88 HCs and 76 FEP patients, we evaluated the correlation between olfactory functions and structural/resting-state functional magnetic resonance imaging (MRI) data.

RESULTS

We found a significant correlation between the left SFG volume and odour discrimination in FEP patients, but not in HCs. We also observed a significant correlation between rs-fMRI connectivity involving the left SFG and odour discrimination in FEP patients, but not in HCs. The data-driven approach didn't observe any significant correlations, possibly due to insufficient statistical power.

CONCLUSION

The left SFG may be a promising brain region in the context of olfactory dysfunction and negative symptoms in FEP.

Multivariate associations between behavioural dimensions and white matter across children and adolescents with and without attention-deficit/hyperactivity disorder

BACKGROUND

Attention deficit/hyperactivity disorder (ADHD) is a heterogeneous neurodevelopmental disorder. Integrity of white matter microstructure plays a key role in the neural mechanism of ADHD presentations. However, the relationships between specific behavioural dimensions and white matter microstructure are less well known. This study aimed to identify associations between white matter and a broad set of clinical features across children and adolescent with and without ADHD using a data-driven multivariate approach.

METHOD

We recruited a total of 130 children (62 controls and 68 ADHD) and employed regularized generalized canonical correlation analysis to characterize the associations between white matter and a comprehensive set of clinical measures covering three domains, including symptom, cognition and behaviour. We further applied linear discriminant analysis to integrate these associations to explore potential developmental effects.

RESULTS

We delineated two brain-behaviour dimensional associations in each domain resulting a total of six multivariate patterns of white matter microstructural alterations linked to hyperactivity-impulsivity and mild affected; executive functions and working memory; externalizing behaviour and social withdrawal, respectively. Apart from executive function and externalizing behaviour sharing similar white matter patterns, all other dimensions linked to a specific pattern of white matter microstructural alterations. The multivariate dimensional association scores showed an overall increase and normalization with age in ADHD group while remained stable in controls.

CONCLUSIONS

We found multivariate neurobehavioral associations exist across ADHD and controls, which suggested that multiple white matter patterns underlie ADHD heterogeneity and provided neural bases for more precise diagnosis and individualized treatment.

Default mode and fronto-parietal network associations with IQ development across childhood in autism

Background

Intellectual disability affects approximately one third of individuals with autism spectrum disorder (autism). Yet, a major unresolved neurobiological question is what differentiates autistic individuals with and without intellectual disability. Intelligence quotients (IQs) are highly variable during childhood. We previously identified three subgroups of autistic children with different trajectories of intellectual development from early (2–3½ years) to middle childhood (9–12 years): (a) persistently high: individuals whose IQs remained in the normal range; (b) persistently low: individuals whose IQs remained in the range of intellectual disability (IQ < 70); and (c) changers: individuals whose IQs began in the range of intellectual disability but increased to the normal IQ range. The frontoparietal (FPN) and default mode (DMN) networks have established links to intellectual functioning. Here, we tested whether brain regions within the FPN and DMN differed volumetrically between these IQ trajectory groups in early childhood.

Methods

We conducted multivariate distance matrix regression to examine the brain regions within the FPN (11 regions x 2 hemispheres) and the DMN (12 regions x 2 hemispheres) in 48 persistently high (18 female), 108 persistently low (32 female), and 109 changers (39 female) using structural MRI acquired at baseline. FPN and DMN regions were defined using networks identified in Smith et al. (Proc Natl Acad Sci U S A 106:13040–5, 2009). IQ trajectory groups were defined by IQ measurements from up to three time points spanning early to middle childhood (mean age time 1: 3.2 years; time 2: 5.4 years; time 3: 11.3 years).

Results

The changers group exhibited volumetric differences in the DMN compared to both the persistently low and persistently high groups at time 1. However, the persistently high group did not differ from the persistently low group, suggesting that DMN structure may be an early predictor for change in IQ trajectory. In contrast, the persistently high group exhibited differences in the FPN compared to both the persistently low and changers groups, suggesting differences related more to concurrent IQ and the absence of intellectual disability.

Conclusions

Within autism, volumetric differences of brain regions within the DMN in early childhood may differentiate individuals with persistently low IQ from those with low IQ that improves through childhood. Structural differences in brain networks between these three IQ-based subgroups highlight distinct neural underpinnings of these autism sub-phenotypes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s11689-022-09460-y.

Neurometabolic Dysfunction in SPG11 Hereditary Spastic Paraplegia

BACKGROUND

Pathogenic variants in SPG11 cause the most common autosomal recessive complicated hereditary spastic paraplegia. Besides the prototypical combination of spastic paraplegia with a thin corpus callosum, obesity has increasingly been reported in this multisystem neurodegenerative disease. However, a detailed analysis of the metabolic state is lacking.

METHODS

In order to characterize metabolic alterations, a cross-sectional analysis was performed comparing SPG11 patients (n = 16) and matched healthy controls (n = 16). We quantified anthropometric parameters, body composition as determined by bioimpedance spectroscopy, and serum metabolic biomarkers, and we measured hypothalamic volume by high-field MRI.

RESULTS

Compared to healthy controls, SPG11 patients exhibited profound changes in body composition, characterized by increased fat tissue index, decreased lean tissue index, and decreased muscle mass. The presence of lymphedema correlated with increased extracellular fluid. The serum levels of the adipokines leptin, resistin, and progranulin were significantly altered in SPG11 while adiponectin and C1q/TNF-related protein 3 (CTRP-3) were unchanged. MRI volumetry revealed a decreased hypothalamic volume in SPG11 patients.

CONCLUSIONS

Body composition, adipokine levels, and hypothalamic volume are altered in SPG11. Our data indicate a link between obesity and hypothalamic neurodegeneration in SPG11 and imply that specific metabolic interventions may prevent obesity despite severely impaired mobility in SPG11.

Differences in functional connectivity distribution after transcranial direct-current stimulation: A connectivity density point of view

In this manuscript, we consider the problem of relating functional connectivity measurements viewed as statistical distributions to outcomes. We demonstrate the utility of using the distribution of connectivity on a study of resting-state functional magnetic resonance imaging association with an intervention. The method uses the estimated density of connectivity between nodes of interest as a functional covariate. Moreover, we demonstrate the utility of the procedure in an instance where connectivity is naturally considered an outcome by reversing the predictor/response relationship using case/control methodology. The method utilizes the density quantile, the density evaluated at empirical quantiles, instead of the empirical density directly. This improved the performance of the method by highlighting tail behavior, though we emphasize that by being flexible and non-parametric, the technique can detect effects related to the central portion of the density. To demonstrate the method in an application, we consider 47 primary progressive aphasia patients with various levels of language abilities. These patients were randomly assigned to two treatment arms, transcranial direct-current stimulation and language therapy versus sham (language therapy only), in a clinical trial. We use the method to analyze the effect of direct stimulation on functional connectivity. As such, we estimate the density of correlations among the regions of interest and study the difference in the density post-intervention between treatment arms. We discover that it is the tail of the density, rather than the mean or lower order moments of the distribution, that demonstrates a significant impact in the classification. The new approach has several benefits. Among them, it drastically reduces the number of multiple comparisons compared with edge-wise analysis. In addition, it allows for the investigation of the impact of functional connectivity on the outcomes where the connectivity is not geometrically localized.

A Machine Learning-Based Approach to Discrimination of Tauopathies Using [18 F]PM-PBB3 PET Images

BACKGROUND

We recently developed a positron emission tomography (PET) probe, [¹⁸ F]PM-PBB3, to detect tau lesions in diverse tauopathies, including mixed three-repeat and four-repeat (3R + 4R) tau fibrils in Alzheimer's disease (AD) and 4R tau aggregates in progressive supranuclear palsy (PSP). For wider availability of this technology for clinical settings, bias-free quantitative evaluation of tau images without a priori disease information is needed.

OBJECTIVE

We aimed to establish tau PET pathology indices to characterize PSP and AD using a machine learning approach and test their validity and tracer capabilities.

METHODS

Data were obtained from 50 healthy control subjects, 46 patients with PSP Richardson syndrome, and 37 patients on the AD continuum. Tau PET data from 114 regions of interest were subjected to Elastic Net cross-validation linear classification analysis with a one-versus-the-rest multiclass strategy to obtain a linear function that discriminates diseases by maximizing the area under the receiver operating characteristic curve. We defined PSP- and AD-tau scores for each participant as values of the functions optimized for differentiating PSP (4R) and AD (3R + 4R), respectively, from others.

RESULTS

The discriminatory ability of PSP- and AD-tau scores assessed as the area under the receiver operating characteristic curve was 0.98 and 1.00, respectively. PSP-tau scores correlated with the PSP rating scale in patients with PSP, and AD-tau scores correlated with Mini-Mental State Examination scores in healthy control-AD continuum patients. The globus pallidus and amygdala were highlighted as regions with high weight coefficients for determining PSP- and AD-tau scores, respectively.

CONCLUSIONS

These findings highlight our technology's unbiased capability to identify topologies of 3R + 4R versus 4R tau deposits. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Altered Development of Amygdala-Connected Brain Regions in Males and Females with Autism

Altered amygdala development is implicated in the neurobiology of autism, but little is known about the coordinated development of the brain regions directly connected with the amygdala. Here we investigated the volumetric development of an amygdala-connected network, defined as the set of brain regions with monosynaptic connections with the amygdala, in autism from early to middle childhood. A total of 950 longitudinal structural MRI scans were acquired from 282 children (93 female) with autism and 128 children with typical development (61 female) at up to four time points (mean ages: 39, 52, 64, and 137 months, respectively). Volumes from 32 amygdala-connected brain regions were examined using mixed effects multivariate distance matrix regression. The Social Responsiveness Scale-2 was administered to assess degree of autistic traits and social impairments. The amygdala-connected network exhibited persistent diagnostic differences (p values ≤ 0.03) that increased over time (p values ≤ 0.02). These differences were most prominent in autistics with more impacted social functioning at baseline. This pattern was not observed across regions without monosynaptic amygdala connection. We observed qualitative sex differences. In males, the bilateral subgenual anterior cingulate cortices were most affected, while in females the left fusiform and superior temporal gyri were most affected. In conclusion, (1) autism is associated with widespread alterations to the development of brain regions connected with the amygdala, which were associated with autistic social behaviors; and (2) autistic males and females exhibited different patterns of alterations, adding to a growing body of evidence of sex differences in the neurobiology of autism.SIGNIFICANCE STATEMENT Global patterns of development across brain regions with monosynaptic connection to the amygdala differentiate autism from typical development, and are modulated by social functioning in early childhood. Alterations to brain regions within the amygdala-connected network differed in males and females with autism. Results also indicate larger volumetric differences in regions having monosynaptic connection with the amygdala than in regions without monosynaptic connection.

Cerebral hemodynamic response to caffeine: effect of dietary caffeine consumption

Caffeine has a significant effect on cerebrovascular systems, and the dual action of caffeine on both neural and vascular responses leads to concerns for the interpretation of blood oxygenation level-dependent (BOLD) functional MRI. However, potential differences in the brain response to caffeine with regard to consumption habits have not been fully elucidated, as BOLD responses may vary with the dietary caffeine consumption history. The main aim of this study was to characterize the acute effect of caffeine on cerebral hemodynamic responses in participants with different patterns of caffeine consumption habits. Fifteen non-habitual and 11 habitual volunteers were included in this study. The cerebral blood flow (CBF) and cerebrovascular reactivity (CVR) to the breath-hold challenge were measured before and after 200 mg caffeine administration. The non-habitual individuals exhibited a pattern of progressive reduction in CBF with time. The CVR was diminished in the caffeinated condition (P < 0.05). In the habitual group, the pattern of CBF decrease was smaller and homogeneous across the brain, and reached steady state rapidly. The CVR was not affected in the presence of caffeine (P > 0.05). Our results demonstrated that the cerebral hemodynamic response to caffeine was subject to the habitual consumption patterns of the participants. The compromised CVR following caffeine administration in the non-habitual group may partially explain the suppressed BOLD response to a visual stimulation in low-caffeine-level users.

Superficial and deep white matter diffusion abnormalities in focal epilepsies

OBJECTIVE

This study was undertaken to evaluate superficial-white matter (WM) and deep-WM magnetic resonance imaging diffusion tensor imaging (DTI) metrics and identify distinctive patterns of microstructural abnormalities in focal epilepsies of diverse etiology, localization, and response to antiseizure medication (ASM).

METHODS

We examined DTI data for 113 healthy controls and 113 patients with focal epilepsies: 51 patients with temporal lobe epilepsy (TLE) and hippocampal sclerosis (HS) refractory to ASM, 27 with pharmacoresponsive TLE-HS, 15 with temporal lobe focal cortical dysplasia (FCD), and 20 with frontal lobe FCD. To assess WM microstructure, we used a multicontrast multiatlas parcellation of DTI. We evaluated fractional anisotropy (FA), mean diffusivity (MD), radial diffusivity (RD), and axial diffusivity (AD), and assessed within-group differences ipsilateral and contralateral to the epileptogenic lesion, as well as between-group differences, in regions of interest (ROIs).

RESULTS

The TLE-HS groups presented more widespread superficial- and deep-WM diffusion abnormalities than both FCD groups. Concerning superficial WM, TLE-HS groups showed multilobar ipsilateral and contralateral abnormalities, with less extensive distribution in pharmacoresponsive patients. Both the refractory TLE-HS and pharmacoresponsive TLE-HS groups also presented pronounced changes in ipsilateral frontotemporal ROIs (decreased FA and increased MD, RD, and AD). Conversely, FCD patients showed diffusion changes almost exclusively adjacent to epileptogenic areas.

SIGNIFICANCE

Our findings add further evidence of widespread abnormalities in WM diffusion metrics in patients with TLE-HS compared to other focal epilepsies. Notably, superficial-WM microstructural damage in patients with FCD is more restricted around the epileptogenic lesion, whereas TLE-HS groups showed diffuse WM damage with ipsilateral frontotemporal predominance. These findings suggest the potential of superficial-WM analysis for better understanding the biological mechanisms of focal epilepsies, and identifying dysfunctional networks and their relationship with the clinical-pathological phenotype. In addition, lobar superficial-WM abnormalities may aid in the diagnosis of subtle FCDs.

A Multi-Atlas-Based [18F]9-Fluoropropyl-(+)-Dihydrotetrabenazine Positron Emission Tomography Image Segmentation Method for Parkinson's Disease Quantification

Objectives

[18F]9-fluoropropyl-(+)-dihydrotetrabenazine ([18F]-FP-DTBZ) positron emission tomography (PET) provides reliable information for the diagnosis of Parkinson's disease (PD). In this study, we proposed a multi-atlas-based [18F]-FP-DTBZ PET image segmentation method for PD quantification assessment.

Methods

A total of 99 subjects from Xuanwu Hospital of Capital Medical University were included in this study, and both brain PET and magnetic resonance (MR) scans were conducted. Data from 20 subjects were used to generate atlases, based on which a multi-atlas-based [18F]-FP-DTBZ PET segmentation method was developed especially for striatum and its subregions. The proposed method was compared with the template-based method through striatal subregion parcellation performance and the standard uptake value ratio (SUVR) quantification accuracy. Discriminant analysis between healthy controls (HCs) and PD patients was further performed.

Results

Segmentation results of the multi-atlas-based method showed better consistency than the template-based method with the ground truth, yielding a dice coefficient of 0.81 over 0.73 on the full striatum. The SUVRs calculated by the multi-atlas-based method had an average interclass correlation coefficient (ICC) of 0.953 with the standardized result, whereas the template-based method only reached 0.815. The SUVRs of HCs were generally higher than that of patients with PD and showed significant differences in all of the striatal subregions (all p < 0.001). The median and posterior putamen performed best in discriminating patients with PD from HCs.

Conclusion

The proposed multi-atlas-based [18F]-FP-DTBZ PET image segmentation method achieved better performance than the template-based method, indicating great potential in improving accuracy and efficiency for PD diagnosis in clinical routine.

Dataset of relationship between longitudinal change in cognitive performance and functional connectivity in cognitively normal older individuals

The data show an association between measured and predicted changes in cognitive performance in older adults who are cognitively normal. Changes in cognitive performance over two years were assessed using the Cognitive Composite Score. The prediction of change in cognitive function was based on changes in pairwise functional connectivity between 80 gray matter regions examined by resting-state functional magnetic resonance imaging. A feature extraction process based on the Variable Importance Testing Approach (VITA) identified changes in 11 pairs of functional connections associated with the default mode network as features related to changes in cognitive performance. Linear and elastic net regression models were applied to these 11 features to predict changes in cognitive performance over two years. A relationship between the 11 features and the geriatric depression score was also shown. The dataset supplements the research findings in the "Changes in pairwise functional connectivity associated with changes in cognitive performance in cognitively normal older individuals: a two-year observational study" published in Oishi et al. (2022). The raw rs-fMRI correlation matrix and associated clinical data can be accessed upon request from the BIOCARD website (www.biocard-se.org) and can be reused for predictive model building.

The Contribution of Working Memory Areas to Verbal Learning and Recall in Primary Progressive Aphasia

Recent evidence of domain-specific working memory (WM) systems has identified the areas and networks which are involved in phonological, orthographic, and semantic WM, as well as in higher level domain-general WM functions. The contribution of these areas throughout the process of verbal learning and recall is still unclear. In the present study, we asked, what is the contribution of domain-specific specialized WM systems in the course of verbal learning and recall? To answer this question, we regressed the perfusion data from pseudo-continuous arterial spin labeling (pCASL) MRI with all the immediate, consecutive, and delayed recall stages of the Rey Auditory Verbal Learning Test (RAVLT) from a group of patients with Primary Progressive Aphasia (PPA), a neurodegenerative syndrome in which language is the primary deficit. We found that the early stages of verbal learning involve the areas with subserving phonological processing (left superior temporal gyrus), as well as semantic WM memory (left angular gyrus, AG_L). As learning unfolds, areas with subserving semantic WM (AG_L), as well as lexical/semantic (inferior temporal and fusiform gyri, temporal pole), and episodic memory (hippocampal complex) become more involved. Finally, a delayed recall depends entirely on semantic and episodic memory areas (hippocampal complex, temporal pole, and gyri). Our results suggest that AG_L subserving domain-specific (semantic) WM is involved only during verbal learning, but a delayed recall depends only on medial and cortical temporal areas.

Clinical and Molecular Correlates of Abnormal Changes in the Cerebellum and Globus Pallidus in Fragile X Premutation

BACKGROUND

Fragile X premutation carriers (55-200 CGG triplets) may develop a progressive neurodegenerative disorder, fragile X-associated tremor/ataxia syndrome (FXTAS), after the age of 50. The neuroradiologic markers of FXTAS are hyperintense T2-signals in the middle cerebellar peduncle-the MCP sign. We recently noticed abnormal T2-signals in the globus pallidus in male premutation carriers and controls but the prevalence and clinical significance were unknown.

METHODS

We estimated the prevalence of the MCP sign and pallidal T2-abnormalities in 230 male premutation carriers and 144 controls (aged 8-86), and examined the associations with FXTAS symptoms, CGG repeat length, and iron content in the cerebellar dentate nucleus and globus pallidus.

RESULTS

Among participants aged ≥45 years (175 premutation carriers and 82 controls), MCP sign was observed only in premutation carriers (52 vs. 0%) whereas the prevalence of pallidal T2-abnormalities approached significance in premutation carriers compared with controls after age-adjustment (25.1 vs. 13.4%, p = 0.069). MCP sign was associated with impaired motor and executive functioning, and the additional presence of pallidal T2-abnormalities was associated with greater impaired executive functioning. Among premutation carriers, significant iron accumulation was observed in the dentate nucleus, and neither pallidal or MCP T2-abnormalities affected measures of the dentate nucleus. While the MCP sign was associated with CGG repeat length >75 and dentate nucleus volume correlated negatively with CGG repeat length, pallidal T2-abnormalities did not correlate with CGG repeat length. However, pallidal signal changes were associated with age-related accelerated iron depletion and variability and having both MCP and pallidal signs further increased iron variability in the globus pallidus.

CONCLUSIONS

Only the MCP sign, not pallidal abnormalities, revealed independent associations with motor and cognitive impairment; however, the occurrence of combined MCP and pallidal T2-abnormalities may present a risk for greater cognitive impairment and increased iron variability in the globus pallidus.

A multimodal study of a first episode psychosis cohort: potential markers of antipsychotic treatment resistance

Treatment resistant (TR) psychosis is considered to be a significant cause of disability and functional impairment. Numerous efforts have been made to identify the clinical predictors of TR. However, the exploration of molecular and biological markers is still at an early stage. To understand the TR condition and identify potential molecular and biological markers, we analyzed demographic information, clinical data, structural brain imaging data, and molecular brain imaging data in 7 Tesla magnetic resonance spectroscopy from a first episode psychosis cohort that includes 136 patients. Age, gender, race, smoking status, duration of illness, and antipsychotic dosages were controlled in the analyses. We found that TR patients had a younger age at onset, more hospitalizations, more severe negative symptoms, a reduction in the volumes of the hippocampus (HP) and superior frontal gyrus (SFG), and a reduction in glutathione (GSH) levels in the anterior cingulate cortex (ACC), when compared to non-TR patients. The combination of multiple markers provided a better classification between TR and non-TR patients compared to any individual marker. Our study shows that ACC-GSH, HP and SFG volumes, and age at onset, could potentially be biomarkers for TR diagnosis, while hospitalization and negative symptoms could be used to evaluate the progression of the disease. Multimodal cohorts are essential in obtaining a comprehensive understanding of brain disorders.

Factorisation-Based Image Labelling

Segmentation of brain magnetic resonance images (MRI) into anatomical regions is a useful task in neuroimaging. Manual annotation is time consuming and expensive, so having a fully automated and general purpose brain segmentation algorithm is highly desirable. To this end, we propose a patched-based labell propagation approach based on a generative model with latent variables. Once trained, our Factorisation-based Image Labelling (FIL) model is able to label target images with a variety of image contrasts. We compare the effectiveness of our proposed model against the state-of-the-art using data from the MICCAI 2012 Grand Challenge and Workshop on Multi-Atlas Labelling. As our approach is intended to be general purpose, we also assess how well it can handle domain shift by labelling images of the same subjects acquired with different MR contrasts.

Application of the dual stream model to neurodegenerative disease: Evidence from a multivariate classification tool in primary progressive aphasia

BACKGROUND

A clinical diagnosis of primary progressive aphasia relies on behavioral characteristics and patterns of atrophy to determine a variant: logopenic; nonfluent/agrammatic; or semantic. The dual stream model (Hickok & Poeppel, 2000; 2004; 2007; 2015) is a contemporary paradigm that has been applied widely to understand brain-behavior relationships; however, applications to neurodegenerative diseases like primary progressive aphasia are limited.

AIMS

The primary aim of this study is to determine if the dual stream model can be applied to a neurodegenerative disease, such as primary progressive aphasia, using both behavioral and neuroimaging data.

METHODS & PROCEDURES

We analyzed behavioral and neuroimaging data to apply a multivariate classification tool (support vector machines) to determine if the dual stream model extends to primary progressive aphasia. Sixty-four individuals with primary progressive aphasia were enrolled (26 logopenic variant, 20 nonfluent/agrammatic variant, and 18 semantic variant) and administered four behavioral tasks to assess three linguistic domains (naming, repetition, and semantic knowledge). We used regions of interest from the dual stream model and calculated the cortical volume for gray matter regions and white matter structural volumes and fractional anisotropy. We applied a multivariate classification tool (support vector machines) to distinguish variants based on behavioral performance and patterns of atrophy.

OUTCOMES & RESULTS

Behavioral performance discriminates logopenic from semantic variant and nonfluent/agrammatic from semantic variant. Cortical volume distinguishes all three variants. White matter structural volumes and fractional anisotropy primarily distinguish nonfluent/agrammatic from semantic variant. Regions of interest that contribute to each classification in cortical and white matter analyses demonstrate alignment of logopenic and nonfluent/agrammatic variants to the dorsal stream, while the semantic variant aligns with the ventral stream.

CONCLUSIONS

A novel implementation of an automated multivariate classification suggests that the dual stream model can be extended to primary progressive aphasia. Variants are distinguished by behavioral and neuroanatomical patterns and align to the dorsal and ventral streams of the dual stream model.

Machine Learning for Predicting Motor Improvement After Acute Subcortical Infarction Using Baseline Whole Brain Volumes

Background. Neuroimaging biomarkers are valuable predictors of motor improvement after stroke, but there is a gap between published evidence and clinical usage. Objective. In this work, we aimed to investigate whether machine learning techniques, when applied to a combination of baseline whole brain volumes and clinical data, can accurately predict individual motor outcome after stroke. Methods. Upper extremity Fugl-Meyer Assessments (FMA-UE) were conducted 1 week and 12 weeks, and structural MRI was performed 1 week, after onset in 56 patients with subcortical infarction. Proportional recovery model residuals were employed to assign patients to proportional and poor recovery groups (34 vs 22). A sophisticated machine learning scheme, consisting of conditional infomax feature extraction, synthetic minority over-sampling technique for nominal and continuous, and bagging classification, was employed to predict motor outcomes, with the input features being a combination of baseline whole brain volumes and clinical data (FMA-UE scores). Results. The proposed machine learning scheme yielded an overall balanced accuracy of 87.71% in predicting proportional vs poor recovery outcomes, a sensitivity of 93.77% in correctly identifying poor recovery outcomes, and a ROC AUC of 89.74%. Compared with only using clinical data, adding whole brain volumes can significantly improve the classification performance, especially in terms of the overall balanced accuracy (from 80.88% to 87.71%) and the sensitivity (from 92.23% to 93.77%). Conclusions. Experimental results suggest that a combination of baseline whole brain volumes and clinical data, when equipped with appropriate machine learning techniques, may provide valuable information for personalized rehabilitation planning after subcortical infarction.

A Hybrid Learning Pipeline for Automated Diagnosis of First-Episode Schizophrenia Utilizing T1-weighted Images

In this work, we proposed and validated a hybrid learning pipeline for automated diagnosis of first-episode schizophrenia (FES) utilizing T1-weighted images. Amygdalar and hippocampal shape abnormalities in FES have been observed in previous studies. In this work, we jointly made use of two types of features, together with advanced machine learning techniques, for an automated discrimination of FES and healthy control (96 versus 102). Specifically, we first employed a ResNet34 model to extract convolutional neural network (CNN) features. We then combined these CNN features with shape features of the bilateral hippocampi and the bilateral amygdalas, before being inputted to advanced classification algorithms such as the Gradient Boosting Decision Tree (GBDT) for classifying between FES and healthy control. Shape features were represented using log Jacobian determinants, through a well-established statistical shape analysis pipeline. When combining CNN with hippocampal shape, the best results came from utilizing GBDT as the classifier, with an overall accuracy of 75.15%, a sensitivity of 69.35%, a specificity of 80.19%, an F1 of 72.16%, and an AUC of 79.68%. When combing CNN and amygdalar shape, the best results came from utilizing Bagging as the classifier, with an overall accuracy of 74.39%, a sensitivity of 67.93%, a specificity of 80%, an F1 of 71.11%, and an AUC of 80.98%. Compared with using each single set of features, either CNN or shape, significant improvements have been observed, in terms of FES discrimination. To the best of our knowledge, this is the first work that has tried to combine CNN features and hippocampal/amygdalar shape features for automated FES identification.

Down-sampling template curve to accelerate LDDMM-curve with application to shape analysis of the corpus callosum

Large deformation diffeomorphic metric mapping for curve (LDDMM-curve) has been widely used in deformation based statistical shape analysis of the mid-sagittal corpus callosum. A main limitation of LDDMM-curve is that it is time-consuming and computationally complex. In this study, down-sampling strategies for accelerating LDDMM-curve are investigated and tested on two large datasets, one on Alzheimer's disease (155 Alzheimer's disease, 325 mild cognitive impairment and 185 healthy controls) and the other on first-episode schizophrenia (92 first-episode schizophrenia and 106 healthy controls). For both datasets a variety of down-sampling factors are tested in terms of registration accuracy, registration speed, and most importantly disease-related patterns. Experimental results revealed that down-sampling template curve by a factor of 2 can significantly reduce the running time of LDDMM-curve without sacrificing the registration accuracy. Also, the disease-induced patterns, or more specifically the group comparison results, were almost identical before and after down-sampling. It is also shown that there was no need to down-sample the target population curves but only the single template curve of the study of interest. Comprehensive analyses were conducted.

Neuroanatomic Volume Differences in Tinnitus and Hearing Loss

OBJECTIVES

To investigate neuroanatomic volume differences in tinnitus and hearing loss.

STUDY DESIGN

Cross-sectional.

METHODS

Sixteen regions of interest (ROIs) in adults (43 male, 29 female) were examined using 3Tesla structural magnetic resonance imaging in four cohorts: 1) tinnitus with moderate hearing loss (N = 31), 2) moderate hearing loss only (N = 15), 3) tinnitus with normal hearing (N = 17), and 4) normal hearing only (N = 13). ROI volumes were corrected for brain size, age, and sex variations. Analysis of covariance (ANCOVA) and post hoc Tukey's test were used to isolate the effects of tinnitus and hearing loss on volume differences. Effect sizes were calculated as the fraction of total variance (η² ) in ANCOVA models and percent of mean volume difference relative to mean total volume.

RESULTS

The four cohort ANCOVA revealed tinnitus and hearing loss cohorts to have increased volume in the corona radiata (η²  = 0.192; P = .0018) and decreased volume in the nucleus accumbens (η²  = 0.252; P < .0001), caudate nucleus (η²  = 0.188; P = .002), and inferior fronto-occipital fasciculus (η²  = 0.250; P = .0001). Tinnitus with normal hearing showed decreased volume in the nucleus accumbens (22.0%; P = .001) and inferior fronto-occipital fasciculus (18.1%; P = .002), and hearing loss only showed increased volume in the corona radiata (10.7%; P = .01) and decreased volume in the nucleus accumbens (22.1%; P = .001), caudate nucleus (16.1%; P = .004), and inferior fronto-occipital fasciculus (18.3%; P = .003).

CONCLUSION

Tinnitus and hearing loss have overlapping effects on neurovolumetric alterations, especially impacting the nucleus accumbens and inferior fronto-occipital fasciculus. Neurovolumetric studies on tinnitus or hearing loss can be more complete by accounting for those two clinical dimensions separately and jointly.

LEVEL OF EVIDENCE

3 Laryngoscope, 131:1863-1868, 2021.

Brain-derived neurotrophic factor Val66Met variant on brain volumes in infants

The brain-derived neurotrophic factor (BDNF) has many important roles in neurogenesis and neuronal health. BDNF is also involved in learning and memory. Individuals with BDNF-Val66Met variant (Met +) are at higher risk for neuropsychiatric disorders and have smaller hippocampi and amgydalae compared to those without this variant (Met -). Whether these smaller brain volumes are already present at birth is unknown and were evaluated. 66 newborn infants were genotyped for BDNF-rs6265 and had brain MRI scans. The T1-weighted images were automatically parcellated for hippocampus and amygdala, as well as the intracranial volume (ICV), total brain volume, total gray and white matter, using a multi-atlas label fusion method implemented in the MRICloud ( https://braingps.anatomyworks.org ). The segmented brain volumes were normalized to the ICV for group comparisons. The two infant groups were not different in their demographics and birth characteristics. However, compared to Met - infants, the Met + infants had smaller hippocampi (p = 0.013), smaller amygdalae (p = 0.041), and less steep age-related declines in total brain volume and % white matter volume. The smaller relative hippocampal and amygdala volumes in Met + infants suggest that the Met + genotype affected prenatal developmental processes. In addition, the slower age-dependent declines in the relative total brain and white matter volumes of the Met + group in this cross-sectional dataset suggest the BDNF-Val66Met variant might have an ongoing negative influence on the postnatal developmental processes.

Hippocampal and Amygdalar Morphological Abnormalities in Alzheimer's Disease Based on Three Chinese MRI Datasets

BACKGROUND

Early detection of Alzheimer's disease (AD) and its early stage, the mild cognitive impairment (MCI), has important scientific, clinical and social significance. Magnetic resonance imaging (MRI) based statistical shape analysis provides an opportunity to detect regional structural abnormalities of brain structures caused by AD and MCI.

OBJECTIVE

In this work, we aimed to employ a well-established statistical shape analysis pipeline, in the framework of large deformation diffeomorphic metric mapping, to identify and quantify the regional shape abnormalities of the bilateral hippocampus and amygdala at different prodromal stages of AD, using three Chinese MRI datasets collected from different domestic hospitals.

METHODS

We analyzed the region-specific shape abnormalities at different stages of the neuropathology of AD by comparing the localized shape characteristics of the bilateral hippocampi and amygdalas between healthy controls and two disease groups (MCI and AD). In addition to group comparison analyses, we also investigated the association between the shape characteristics and the Mini Mental State Examination (MMSE) of each structure of interest in the disease group (MCI and AD combined) as well as the discriminative power of different morphometric biomarkers.

RESULTS

We found the strongest disease pathology (regional atrophy) at the subiculum and CA1 subregions of the hippocampus and the basolateral, basomedial as well as centromedial subregions of the amygdala. Furthermore, the shape characteristics of the hippocampal and amygdalar subregions exhibiting the strongest AD related atrophy were found to have the most significant positive associations with the MMSE. Employing the shape deformation marker of the hippocampus or the amygdala for automated MCI or AD detection yielded a significant accuracy boost over the corresponding volume measurement.

CONCLUSION

Our results suggested that the amygdalar and hippocampal morphometrics, especially those of shape morphometrics, can be used as auxiliary indicators for monitoring the disease status of an AD patient.

A Deformation-Based Shape Study of the Corpus Callosum in First Episode Schizophrenia

Background: Previous first-episode schizophrenia (FES) studies have reported abnormalities in the volume and mid-sagittal size of the corpus callosum (CC), but findings have been inconsistent. Besides, the CC shape has rarely been analyzed in FES. Therefore, in this study, we investigated FES-related CC shape abnormalities using 198 participants [92 FES patients and 106 healthy controls (HCs)]. Methods: We conducted statistical shape analysis of the mid-sagittal CC curve in a large deformation diffeomorphic metric mapping framework. The CC was divided into the genu, body, and splenium (gCC, bCC, and sCC) to target the key CC sub-regions affected by the FES pathology. Gender effects have been investigated. Results: There were significant area differences between FES and HC in the entire CC and gCC but not in bCC nor sCC. In terms of the localized shape morphometrics, significant region-specific shape inward-deformations were detected in the superior portion of gCC and the anterosuperior portion of bCC in FES. These global area and local shape morphometric abnormalities were restricted to female FES but not male FES. Conclusions: gCC was significantly affected in the neuropathology of FES and this finding was specific to female FES. This study suggests that gCC may be a key sub-region that is vulnerable to the neuropathology of FES, specifically in female patients. The morphometrics of gCC may serve as novel and efficient biomarkers for screening female FES patients.

A 7 Tesla Amygdalar-Hippocampal Shape Analysis of Lithium Response in Bipolar Disorder

Research to discover clinically useful predictors of lithium response in patients with bipolar disorder has largely found them to be elusive. We demonstrate here that detailed neuroimaging may have the potential to fill this important gap in mood disorder therapeutics. Lithium treatment and bipolar disorder have both been shown to affect anatomy of the hippocampi and amygdalae but there is no consensus on the nature of their effects. We aimed to investigate structural surface anatomy changes in amygdala and hippocampus correlated with treatment response in bipolar disorder. Patients with bipolar disorder (N = 14) underwent lithium treatment, were classified by response status at acute and long-term time points, and scanned with 7 Tesla structural MRI. Large Deformation Diffeomorphic Metric Mapping was applied to detect local differences in hippocampal and amygdalar anatomy between lithium responders and non-responders. Anatomy was also compared to 21 healthy comparison participants. A patch of the ventral surface of the left hippocampus was found to be significantly atrophied in non-responders as compared to responders at the acute time point and was associated at a trend-level with long-term response status. We did not detect an association between response status and surface anatomy of the right hippocampus or amygdala. To the best of our knowledge, this is the first shape analysis of hippocampus and amygdala in bipolar disorder using 7 Tesla MRI. These results can inform future work investigating possible neuroimaging predictors of lithium response in bipolar disorder.

Multimodal MRI assessment for first episode psychosis: A major change in the thalamus and an efficient stratification of a subgroup

Multi‐institutional brain imaging studies have emerged to resolve conflicting results among individual studies. However, adjusting multiple variables at the technical and cohort levels is challenging. Therefore, it is important to explore approaches that provide meaningful results from relatively small samples at institutional levels. We studied 87 first episode psychosis (FEP) patients and 62 healthy subjects by combining supervised integrated factor analysis (SIFA) with a novel pipeline for automated structure‐based analysis, an efficient and comprehensive method for dimensional data reduction that our group recently established. We integrated multiple MRI features (volume, DTI indices, resting state fMRI—rsfMRI) in the whole brain of each participant in an unbiased manner. The automated structure‐based analysis showed widespread DTI abnormalities in FEP and rs‐fMRI differences between FEP and healthy subjects mostly centered in thalamus. The combination of multiple modalities with SIFA was more efficient than the use of single modalities to stratify a subgroup of FEP (individuals with schizophrenia or schizoaffective disorder) that had more robust deficits from the overall FEP group. The information from multiple MRI modalities and analytical methods highlighted the thalamus as significantly abnormal in FEP. This study serves as a proof‐of‐concept for the potential of this methodology to reveal disease underpins and to stratify populations into more homogeneous sub‐groups.

Gray matter alteration in heroin-dependent men: An atlas-based magnetic resonance imaging study

Previous imaging studies on heroin addiction have reported brain morphological alterations. However, the effects of heroin exposure on gray matter volume varied among different studies due to different factors such as substitution treatment or mandatory abstinence. Meanwhile, the relationship between gray matter and heroin use history remains unknown. Thirty-three male heroin-dependent (HD) individuals who are not under any substitution treatment or mandatory abstinence and 40 male healthy controls (HC) were included in this structural magnetic resonance imaging study. With an atlas-based approach, gray matter structures up to individual functional area were delineated, and the differences in their volumes between the HD and HC groups were analyzed. In addition, the relationship between gray matter volume and duration of heroin use was explored. The HD group demonstrated significantly lower cortical volume mainly in the prefrontal cortex and mesolimbic dopaminergic regions across different parcellation levels, whereas several visual and somatosensory cortical regions in the HD group had greater volume relative to the HC group at a more detailed parcellation level. The duration of heroin use was negatively correlated with the gray matter volume of prefrontal cortex. These findings suggest that heroin addiction be related to gray matter alteration and might be related to damage/maladaption of the inhibitory control, reward, visual, and somatosensory functions of the brain, although cognitive correlates are warranted in future study. In addition, the atlas-based morphology analysis is a potential tool to help researchers search biomarkers of heroin addiction.

A deep learning framework for pancreas segmentation with multi-atlas registration and 3D level-set

In this paper, we propose and validate a deep learning framework that incorporates both multi-atlas registration and level-set for segmenting pancreas from CT volume images. The proposed segmentation pipeline consists of three stages, namely coarse, fine, and refine stages. Firstly, a coarse segmentation is obtained through multi-atlas based 3D diffeomorphic registration and fusion. After that, to learn the connection feature, a 3D patch-based convolutional neural network (CNN) and three 2D slice-based CNNs are jointly used to predict a fine segmentation based on a bounding box determined from the coarse segmentation. Finally, a 3D level-set method is used, with the fine segmentation being one of its constraints, to integrate information of the original image and the CNN-derived probability map to achieve a refine segmentation. In other words, we jointly utilize global 3D location information (registration), contextual information (patch-based 3D CNN), shape information (slice-based 2.5D CNN) and edge information (3D level-set) in the proposed framework. These components form our cascaded coarse-fine-refine segmentation framework. We test the proposed framework on three different datasets with varying intensity ranges obtained from different resources, respectively containing 36, 82 and 281 CT volume images. In each dataset, we achieve an average Dice score over 82%, being superior or comparable to other existing state-of-the-art pancreas segmentation algorithms.

Whole brain volume and cortical thickness abnormalities in Wilson's disease: a clinical correlation study

Wilson's disease (WD) is an inherited autosomal recessive disorder of copper metabolism, and its neurological and neuropsychiatric manifestations are associated with copper accumulation in brain. A few neuroimaging studies have shown that gray matter atrophy in WD affects both subcortical structures and cortex. This study aims to quantitatively evaluate the morphometric brain abnormalities in patients with WD in terms of whole brain volume and cortical thickness and their associations with clinical severity of WD. Thirty patients clinically diagnosed as WD with neurological manifestations and 25 healthy controls (HC) were recruited. 3D T1-weighted images were segmented into 276 whole-brain regions of interest (ROIs) and 68 cortical ROIs. WD-vs-HC group comparisons were then conducted for each ROI. The associations between those morphometric measurements and the Global Assessment Scale (GAS) score for WD were analyzed. Compared with HC, significant WD-related volumetric decreases were found in the bilateral subcortical nuclei (putamen, globus pallidus, caudate nucleus, substantia nigra, red nucleus and thalamus), diffuse white matter and several gray matter regions. WD patients showed reduced cortical thickness in the left precentral gyrus and the left insula. Further, the volumes of the right globus pallidus, bilateral putamen, right external capsule and left superior longitudinal fasciculus were negatively correlated with GAS. Our results indicated that significant WD-related morphometric abnormalities were quantified in terms of whole-brain volumes and cortical thicknesses, some of which correlated significantly to the clinical severity of WD. Those morphometrics may provide a potentially effective biomarker of WD.

Effects of spinal cord stimulation on voxel-based brain morphometry in patients with failed back surgery syndrome

OBJECTIVE

Despite the clinical effectiveness of Spinal Cord Stimulation (SCS), potential structural brain modifications have not been explored. Our aim was to identify structural volumetric changes during subsensory SCS, in patients with Failed Back Surgery Syndrome (FBSS).

METHODS

In this cohort study, twenty-two FBSS patients underwent a magnetic resonance imaging protocol before SCS and 3 months after SCS. Clinical parameters were correlated with volumetric changes, calculated with voxel-based morphometry.

RESULTS

After 3 months, a significant volume decrease was found in the inferior frontal gyrus, precuneus, cerebellar posterior lobe and middle temporal gyrus. Significant increases were found in the inferior temporal gyrus, precentral gyrus and the middle frontal gyrus after SCS. Additionally, significant increases in volume of superior frontal and parietal white matter and a significant decrease in volume of white matter underlying the premotor/middle frontal gyrus were revealed after SCS. A significant correlation was highlighted between white matter volume underlying premotor/middle frontal gyrus and leg pain relief.

CONCLUSIONS

This study revealed for the first time that SCS is able to induce volumetric changes in gray and white matter, suggesting the reversibility of brain alterations after chronic pain treatment.

SIGNIFICANCE

Volumetric brain alterations are observable after 3 months of subsensory SCS in FBSS patients.

Left ventricle segmentation using a Bayesian approach with distance dependent shape priors

We propose a method for segmentation of the left ventricle in magnetic resonance cardiac images. The framework consists of an initial Bayesian segmentation of the central slice of the volume. This segmentation is used to locate a shape prior for the LV myocardial tissue. This shape prior is determined using the fact that the myocardium is approximately annular as seen in the short-axis. Then a second Bayesian segmentation is performed to obtain the final result. This procedure is repeated for the rest of the slices. An extrapolation of the area of the LV is used to determine a stopping criterion. The method was evaluated on the databases of the Cardiac Atlas project. Our results demonstrate a suitable accuracy for myocardial segmentation (≈0.8 Dice's coefficient). For the endocardium and the epicardium the Dice's coefficients are 0.94 and 0.9 respectively. The accuracy was also evaluated in terms of the Hausdorff distance and the average distance. For the myocardium we obtain 8 mm and 2 mm respectively. Our results demonstrate the capability and merits of the proposed method to estimate the structure of the LV. The method requires minimal user input and generates results with quality comparable to more complex approaches. This paper suggests a new efficient approach for automatic LV quantification based on a Bayesian technique with shape priors with errors comparable to state-of-the-art techniques.

Short- and long-term reproducibility of BOLD signal change induced by breath-holding at 1.5 and 3 T

Cerebrovascular reactivity (CVR) can give insight into the cerebrovascular function. CVR can be estimated by measuring a blood-oxygen-level-dependent (BOLD) response combined with breath-holding (BH). The reproducibility of this technique has been addressed and existing studies have focused on short-term reproducibility using a 3 T magnetic resonance imaging (MRI) system. However, little is known about the long-term reproducibility of this procedure and the corresponding reproducibility using a 1.5 T MRI system. Here, we systematically examined the short- and long-term reproducibility of BOLD responses to BH across field strengths. Nine subjects participated in three MRI sessions separated by 30 minutes (sessions 1 and 2: short term) and 68-92 days (sessions 1 and 3, long term) at both 1.5 and 3 T MRI. Our findings revealed that significant differences between field strengths were detected in the activated gray matter volume and BOLD signal change (both P < 0.001), with smaller magnitudes at 1.5 T. However, activation patterns were reproducible, independent of the time interval, brain region or field strength. All interscan coefficient of variation values were below the 33% fiducial limit, and the intraclass correlation coefficient values were above 0.4, which is usually considered the acceptability limit in functional studies. These findings suggest that the response of BOLD signal to BH for assessing CVR is reproducible over time at 1.5 and 3 T. This technique can be considered a tool for monitoring longitudinal changes in patients with cerebrovascular diseases, and its use should be encouraged for clinical 1.5 T MRI systems.

Diffeomorphic Registration With Intensity Transformation and Missing Data: Application to 3D Digital Pathology of Alzheimer's Disease

This paper examines the problem of diffeomorphic image registration in the presence of differing image intensity profiles and sparsely sampled, missing, or damaged tissue. Our motivation comes from the problem of aligning 3D brain MRI with 100-micron isotropic resolution to histology sections at 1 × 1 × 1,000-micron resolution with multiple varying stains. We pose registration as a penalized Bayesian estimation, exploiting statistical models of image formation where the target images are modeled as sparse and noisy observations of the atlas. In this injective setting, there is no assumption of symmetry between atlas and target. Cross-modality image matching is achieved by jointly estimating polynomial transformations of the atlas intensity. Missing data is accommodated via a multiple atlas selection procedure where several atlas images may be of homogeneous intensity and correspond to "background" or "artifact." The two concepts are combined within an Expectation-Maximization algorithm, where atlas selection posteriors and deformation parameters are updated iteratively and polynomial coefficients are computed in closed form. We validate our method with simulated images, examples from neuropathology, and a standard benchmarking dataset. Finally, we apply it to reconstructing digital pathology and MRI in standard atlas coordinates. By using a standard convolutional neural network to detect tau tangles in histology slices, this registration method enabled us to quantify the 3D density distribution of tauopathy throughout the medial temporal lobe of an Alzheimer's disease postmortem specimen.

Interaction between ventricular expansion and structural changes in the corpus callosum and putamen in males with FMR1 normal and premutation alleles

Ventricular enlargement (VE) is commonly observed in aging and fragile X-associated tremor/ataxia syndrome (FXTAS), a late-onset neurodegenerative disorder. VE may generate a mechanical force causing structural deformation. In this longitudinal study, we examined the relationships between VE and structural changes in the corpus callosum (CC) and putamen. MRI scans (2-7/person over 0.2-7.5 years) were acquired from 22 healthy controls, 26 unaffected premutation carriers (PFX-), and 39 carriers affected with FXTAS (PFX+). Compared with controls, PFX- demonstrated enlarged fourth ventricles, whereas PFX+ displayed enlargement in both third and fourth ventricles, CC thinning, putamen atrophy/deformation (thinning and increased distance), and accelerated expansions in lateral ventricles. Common for all groups, baseline VE predicted accelerated CC thinning and putamen atrophy/deformation and conversely, baseline CC and putamen atrophy/deformation and enlarged third and fourth ventricles predicted accelerated lateral ventricular expansion. The results suggest a progressive VE within the 4 ventricles as FXTAS develops and a deleterious cycle between VE and brain deformation that may commonly occur during aging and FXTAS progression but become accelerated in FXTAS.

Developmental aspects of FXAND in a man with the FMR1 premutation

BACKGROUND

Fragile X mental retardation 1 (FMR1) premutation can cause developmental problems including autism spectrum disorder (ASD), social anxiety, depression, and attention deficit hyperactivity disorder (ADHD). These problems fall under an umbrella term of Fragile X-associated Neuropsychiatric Disorders (FXAND) and is separate from Fragile X-associated Tremor/Ataxia syndrome (FXTAS), a neurodegenerative disorder.

METHODS/CLINICAL CASE

A 26-year-old Caucasian male with the Fragile X premutation who presented with multiple behavior and emotional problems including depression and anxiety at 10 years of age. He was evaluated at 13, 18, and 26 years old with age-appropriate cognitive assessments, psychiatric evaluations, and an MRI of the brain.

RESULTS

The Autism Diagnostic Observation Scale (ADOS) was done at 13 years old and showed the patient has autism spectrum disorder (ASD). An evaluation at 18 years old showed a full-scale IQ of 64. A Kiddie Schedule for Affective Disorders and Schizophrenia (K-SADS) performed at 26 years old confirmed the previous impression of social anxiety disorder, agoraphobia disorder, and selective mutism. His MRI acquired at 26 years old showed enlarged ventricles, increased frontal subarachnoid spaces, and hypergyrification.

CONCLUSION

This is an exemplary case of an FMR1 premutation carrier with significant psychiatric and cognitive issues that demonstrates Fragile X-associated Neuropsychiatric Disorders (FXAND) as separate from the other well-known premutation disorders.

Baseline MRI associates with later naming status in primary progressive aphasia

Advanced imaging studies in neurodegenerative disease have yielded new insights into subtypes of disease, progression of disease in various brain regions, and changes in structural and functional connectivity between brain regions related to symptom progression. However, few studies have revealed imaging markers at baseline that correlate with rate or degree of decline in function. Here we tested the hypothesis that imaging features at baseline correlate with outcome of naming in primary progressive aphasia. We obtained longitudinal multimodal imaging in 15 individuals with primary progressive aphasia at the same time points as assessment of naming. We found that functional connectivity between particular brain regions (measured with resting state functional connectivity magnetic resonance imaging) is strongly associated with accuracy of naming 21 months later, independently of baseline severity of naming impairment. These data indicate that functional connectivity may carry information about later performance in naming, and is potentially useful for refining prognosis.

Anatomical and Functional Characterization in Children With Unilateral Cerebral Palsy: An Atlas-Based Analysis

Background. Variability in hand function among children with unilateral cerebral palsy (UCP) might reflect the type of brain injury and resulting anatomical sequelae. Objective. We used atlas-based analysis of structural images to determine whether children with periventricular (PV) versus middle cerebral artery (MCA) injuries might exhibit unique anatomical characteristics that account for differences in hand function. Methods. Forty children with UCP underwent structural brain imaging using 3-T magnetic resonance imaging. Brain lesions were classified as PV or MCA. A group of 40 typically developing (TD) children served as comparison controls. Whole brains were parcellated into 198 structures (regions of interest) to obtain volume estimates. Dexterity and bimanual hand function were assessed. Unbiased, differential expression analysis was performed to determine volumetric differences between PV and MCA groups. Principal component analysis (PCA) was performed and the top 3 components were extracted to perform regression on hand function. Results. Children with PV had significantly better hand function than children with MCA. Multidimensional scaling analysis of volumetric data revealed separate clustering of children with MCA, PV, and TD children. PCA extracted anatomical components that comprised the 2 types of brain injury. In the MCA group, reductions of volume were concentrated in sensorimotor structures of the injured hemisphere. Models using PCA predicted hand function with greater accuracy than models based on qualitative brain injury type. Conclusions. Our results highlight unique quantitative differences in children with UCP that also predict differences in hand function. The systematic discrimination between groups found in our study reveals future questions about the potential prognostic utility of this approach.

Amygdalar and Hippocampal Morphometry Abnormalities in First-Episode Schizophrenia Using Deformation-Based Shape Analysis

In this study, we investigated and quantified the amygdalar and hippocampal morphometry abnormalities exerted by first-episode schizophrenia using a total of 92 patients and 106 healthy control participants. Magnetic resonance imaging (MRI) based automated segmentation was conducted to obtain the amygdalar and hippocampal segmentations. Disease-versus-control volume differences of the bilateral amygdalas and hippocampi were quantified. In addition, deformation-based statistical shape analysis was employed to quantify the region-specific shape abnormalities of each structure of interest. To better identify the key relevant areas in the pathology of first-episode schizophrenia, each structure was divided into four subregions; CA1, CA2, CA3 combined with dentate gyrus for the hippocampus in each hemisphere and basolateral, basomedial, centromedial, and lateral nucleus for the amygdala in each hemisphere. We observed significant global volume reduction and localized shape atrophy in each of the four structures of interest. The amygdalar shape abnormalities mainly occurred at the basolateral and centromedial subregions, whereas the hippocampal shape abnormalities mainly concentrated on the CA1 and CA2 subregions. For the same structure, the one on the right hemisphere was affected more by the disease pathology than that on the left hemisphere. To conclude, we have successfully quantified the global and local morphometric abnormalities of the bilateral amygdalas and hippocampi using a sophisticated statistical analysis pipeline and high-field subregion segmentations, with MRI data of a considerable sample size. This study is one of the very first of such kind in first-episode schizophrenia analyses.

Brain volumes as predictors of tDCS effects in primary progressive aphasia

The current study aims to determine the brain areas critical for response to anodal transcranial direct current stimulation (tDCS) in PPA. Anodal tDCS and sham were administered over the left inferior frontal gyrus (IFG), combined with written naming/spelling therapy. Thirty people with PPA were included in this study, and assessed immediately, 2 weeks, and 2 months post-therapy. We identified anatomical areas whose volumes significantly predicted the additional tDCS effects. For trained words, the volumes of the left Angular Gyrus and left Posterior Cingulate Cortex predicted the additional tDCS gain. For untrained words, the volumes of the left Middle Frontal Gyrus, left Supramarginal Gyrus, and right Posterior Cingulate Cortex predicted the additional tDCS gain. These findings show that areas involved in language, attention and working memory contribute to the maintenance and generalization of stimulation effects. The findings highlight that tDCS possibly affects areas anatomically or functionally connected to stimulation targets.

Using deep Siamese neural networks for detection of brain asymmetries associated with Alzheimer's Disease and Mild Cognitive Impairment

In recent studies, neuroanatomical volume and shape asymmetries have been seen during the course of Alzheimer's Disease (AD) and could potentially be used as preclinical imaging biomarkers for the prediction of Mild Cognitive Impairment (MCI) and AD dementia. In this study, a deep learning framework utilizing Siamese neural networks trained on paired lateral inter-hemispheric regions is used to harness the discriminative power of whole-brain volumetric asymmetry. The method uses the MRICloud pipeline to yield low-dimensional volumetric features of pre-defined atlas brain structures, and a novel non-linear kernel trick to normalize these features to reduce batch effects across datasets and populations. By working with the low-dimensional features, Siamese networks were shown to yield comparable performance to studies that utilize whole-brain MR images, with the advantage of reduced complexity and computational time, while preserving the biological information density. Experimental results also show that Siamese networks perform better in certain metrics by explicitly encoding the asymmetry in brain volumes, compared to traditional prediction methods that do not use the asymmetry, on the ADNI and BIOCARD datasets.

Significant Dimension Reduction of 3D Brain MRI using 3D Convolutional Autoencoders

Content-based image retrieval (CBIR) is a technology designed to retrieve images from a database based on visual features. While the CBIR is highly desired, it has not been applied to clinical neuroradiology, because clinically relevant neuroradiological features are swamped by a huge number of noisy and unrelated voxel information. Thus, effective dimension reduction is the key to successful CBIR. We propose a novel dimensional compression method based on 3D convolutional autoencoders (3D-CAE), which was applied to the ADNI2 3D brain MRI dataset. Our method succeeded in compressing 5 million voxel information to only 150 dimensions, while preserving clinically relevant neuroradiological features. The RMSE per voxel was as low as 8.4%, suggesting a promise of our method toward the application to the CBIR.

Test-retest reproducibility of a multi-atlas automated segmentation tool on multimodality brain MRI

INTRODUCTION

The increasing use of large sample sizes for population and personalized medicine requires high-throughput tools for imaging processing that can handle large amounts of data with diverse image modalities, perform a biologically meaningful information reduction, and result in comprehensive quantification. Exploring the reproducibility of these tools reveals the specific strengths and weaknesses that heavily influence the interpretation of results, contributing to transparence in science.

METHODS

We tested-retested the reproducibility of MRICloud, a free automated method for whole-brain, multimodal MRI segmentation and quantification, on two public, independent datasets of healthy adults.

RESULTS

The reproducibility was extremely high for T1-volumetric analysis, high for diffusion tensor images (DTI) (however, regionally variable), and low for resting-state fMRI.

CONCLUSION

In general, the reproducibility of the different modalities was slightly superior to that of widely used software. This analysis serves as a normative reference for planning samples and for the interpretation of structure-based MRI studies.