Michael Dayan

The role of iron in gray matter degeneration in Huntington's disease: A magnetic resonance imaging study

In Huntingtons disease, iron accumulation in basal ganglia accompanies neuronal loss. However, if iron content changes with disease progression and how it relates to gray matter atrophy is not clear yet. We explored iron content in basal ganglia and cortex and its relationship with gray matter volume in 77 mutation carriers [19 presymptomatic, 8 with soft symptoms (SS), and 50 early‐stage patients) and 73 matched‐controls by T2*relaxometry and T1‐weighted imaging on a 3T scanner. The ANCOVA model showed that iron accumulates in the caudate in presymptomatic subjects (P = 0.004) and remains relatively stable along disease stages in this nucleus; while increases in putamen and globus pallidus (P < 0.05). Volume instead decreases in basal ganglia, starting from the caudate (P < 0.0001) and extending to the putamen and globus pallidus (P ≤ 0.001). The longer the disease duration and the higher the CAG repeats, the higher the iron accumulation and the smaller the volume. In the cortex, iron decreases in parieto‐occipital areas in SS (P < 0.027); extending to premotor and parieto‐temporo‐occipital areas in patients (P < 0.003); while volume declines in frontoparietal and temporal areas in presymptomatic (P < 0.023) and SS (P < 0.045), and extends throughout the cortex, with the exception of anterior frontal regions, in patients (P < 0.023). There is an inverse correlation between volume and iron levels in putamen, globus pallidus and the anterior cingulate; and a direct correlation in cortical structures (SMA‐sensoriomotor and temporo‐occipital). Iron homeostasis is affected in the disease; however, there appear to be differences in the role played by iron in basal ganglia and in cortex. Hum Brain Mapp, 36:–66, 2015.

Modeling the Relationship among Gray Matter Atrophy, Abnormalities in Connecting White Matter, and Cognitive Performance in Early Multiple Sclerosis

BACKGROUND AND PURPOSE: Quantitative assessment of clinical and pathologic consequences of white matter abnormalities in multiple sclerosis is critical in understanding the pathways of disease. This study aimed to test whether gray matter atrophy was related to abnormalities in connecting white matter and to identify patterns of imaging biomarker abnormalities that were related to patient processing speed. MATERIALS AND METHODS: Image data and Symbol Digit Modalities Test scores were collected from a cohort of patients with early multiple sclerosis. The Network Modification Tool was used to estimate connectivity irregularities by projecting white matter abnormalities onto connecting gray matter regions. Partial least-squares regression quantified the relationship between imaging biomarkers and processing speed as measured by the Symbol Digit Modalities Test. RESULTS: Atrophy in deep gray matter structures of the thalami and putamen had moderate and significant correlations with abnormalities in connecting white matter (r = 0.39–0.41, P < .05 corrected). The 2 models of processing speed, 1 for each of the WM imaging biomarkers, had goodness-of-fit (R2) values of 0.42 and 0.30. A measure of the impact of white matter lesions on the connectivity of occipital and parietal areas had significant nonzero regression coefficients. CONCLUSIONS: We concluded that deep gray matter regions may be susceptible to inflammation and/or demyelination in white matter, possibly having a higher sensitivity to remote degeneration, and that lesions affecting visual processing pathways were related to processing speed. The Network Modification Tool may be used to quantify the impact of early white matter abnormalities on both connecting gray matter structures and processing speed.

Consensus between pipelines in structural brain networks

Structural brain networks may be reconstructed from diffusion MRI tractography data and have great potential to further our understanding of the topological organisation of brain structure in health and disease. Network reconstruction is complex and involves a series of processesing methods including anatomical parcellation, registration, fiber orientation estimation and whole-brain fiber tractography. Methodological choices at each stage can affect the anatomical accuracy and graph theoretical properties of the reconstructed networks, meaning applying different combinations in a network reconstruction pipeline may produce substantially different networks. Furthermore, the choice of which connections are considered important is unclear. In this study, we assessed the similarity between structural networks obtained using two independent state-of-the-art reconstruction pipelines. We aimed to quantify network similarity and identify the core connections emerging most robustly in both pipelines. Similarity of network connections was compared between pipelines employing different atlases by merging parcels to a common and equivalent node scale. We found a high agreement between the networks across a range of fiber density thresholds. In addition, we identified a robust core of highly connected regions coinciding with a peak in similarity across network density thresholds, and replicated these results with atlases at different node scales. The binary network properties of these core connections were similar between pipelines but showed some differences in atlases across node scales. This study demonstrates the utility of applying multiple structural network reconstrution pipelines to diffusion data in order to identify the most important connections for further study.

MRI Analysis of White Matter Myelin Water Content in Multiple Sclerosis: A Novel Approach Applied to Finding Correlates of Cortical Thinning

A novel lesion-mask free method based on a gamma mixture model was applied to myelin water fraction (MWF) maps to estimate the association between cortical thickness and myelin content, and how it differs between relapsing-remitting (RRMS) and secondary-progressive multiple sclerosis (SPMS) groups (135 and 23 patients, respectively). It was compared to an approach based on lesion masks. The gamma mixture distribution of whole brain, white matter (WM) MWF was characterized with three variables: the mode (most frequent value) m1 of the gamma component shown to relate to lesion, the mode m2 of the component shown to be associated with normal appearing (NA) WM, and the mixing ratio (λ) between the two distributions. The lesion-mask approach relied on the mean MWF within lesion and within NAWM. A multivariate regression analysis was carried out to find the best predictors of cortical thickness for each group and for each approach. The gamma-mixture method was shown to outperform the lesion-mask approach in terms of adjusted R2, both for the RRMS and SPMS groups. The predictors of the final gamma-mixture models were found to be m1 (β = 1.56, p < 0.005), λ (β = −0.30, p < 0.0005) and age (β = −0.0031, p < 0.005) for the RRMS group (adjusted R2 = 0.16), and m2 (β = 4.72, p < 0.0005) for the SPMS group (adjusted R2 = 0.45). Further, a DICE coefficient analysis demonstrated that the lesion mask had more overlap to an ROI associated with m1, than to an ROI associated with m2 (p < 0.00001), and vice versa for the NAWM mask (p < 0.00001). These results suggest that during the relapsing phase, focal WM damage is associated with cortical thinning, yet in SPMS patients, global WM deterioration has a much stronger influence on secondary degeneration. Through these findings, we demonstrate the potential contribution of myelin loss on neuronal degeneration at different disease stages and the usefulness of our statistical reduction technique which is not affected by the typical bias associated with approaches based on lesion masks.

Profilometry: A new statistical framework for the characterization of white matter pathways, with application to multiple sclerosis

Aims: describe a new “profilometry” framework for the multimetric analysis of white matter tracts, and demonstrate its application to multiple sclerosis (MS) with radial diffusivity (RD) and myelin water fraction (MWF). Methods: A cohort of 15 normal controls (NC) and 141 MS patients were imaged with T1, T2 FLAIR, T2 relaxometry and diffusion MRI (dMRI) sequences. T1 and T2 FLAIR allowed for the identification of patients having lesion(s) on the tracts studied, with a special focus on the forceps minor. T2 relaxometry provided MWF maps, while dMRI data yielded RD maps and the tractography required to compute MWF and RD tract profiles. The statistical framework combined a multivariate analysis of covariance (MANCOVA) and a linear discriminant analysis (LDA) both accounting for age and gender, with multiple comparison corrections. Results: In the single‐case case study the profilometry visualization showed a clear departure of MWF and RD from the NC normative data at the lesion location(s). Group comparison from MANCOVA demonstrated significant differences at lesion locations, and a significant age effect in several tracts. The follow‐up LDA analysis suggested MWF better discriminates groups than RD. Discussion and conclusion: While progress has been made in both tract‐profiling and metrics for white matter characterization, no single framework for a joint analysis of multimodality tract profiles accounting for age and gender is known to exist. The profilometry analysis and visualization appears to be a promising method to compare groups using a single score from MANCOVA while assessing the contribution of each metric with LDA. Hum Brain Mapp 37:989–1004, 2016.

Functional brain connectivity is predictable from anatomic network's Laplacian eigen-structure

&NA; How structural connectivity (SC) gives rise to functional connectivity (FC) is not fully understood. Here we mathematically derive a simple relationship between SC measured from diffusion tensor imaging, and FC from resting state fMRI. We establish that SC and FC are related via (structural) Laplacian spectra, whereby FC and SC share eigenvectors and their eigenvalues are exponentially related. This gives, for the first time, a simple and analytical relationship between the graph spectra of structural and functional networks. Laplacian eigenvectors are shown to be good predictors of functional eigenvectors and networks based on independent component analysis of functional time series. A small number of Laplacian eigenmodes are shown to be sufficient to reconstruct FC matrices, serving as basis functions. This approach is fast, and requires no time‐consuming simulations. It was tested on two empirical SC/FC datasets, and was found to significantly outperform generative model simulations of coupled neural masses.

Estimating function from structure in epileptics using graph diffusion model

The relationship between anatomic and resting state functional connectivity (FC) of large-scale brain networks has been of interest and has been investigated in a number of articles. In a recent article we introduced a graph diffusion model which predicts the functional network from the structural network in healthy brains. In this work we apply the graph diffusion model to two types of epilepsy, medial temporal sclerosis epilepsy (TLE-MTS), and MRI-normal temporal lobe epilepsy (TLE-no). We show that it is possible to estimate function from structure in non-healthy brains. We conclude that TLE-MTS on average requires a higher graph diffusion depth to estimate FC than both the healthy or the TLE-no groups. This suggests that an overly strong FC/SC relationship might be a sign of poor brain condition.

Estimating brain's functional graph from the structural graph's Laplacian

The interplay between the brain’s function and structure has been of immense interest to the neuroscience and connectomics communities. In this work we develop a simple linear model relating the structural network and the functional network. We propose that the two networks are related by the structural network’s Laplacian up to a shift. The model is simple to implement and gives accurate prediction of function’s eigenvalues at the subject level and its eigenvectors at group level.