A joint and integrated analysis of multi-site diffusion MRI (dMRI) datasets can dramatically increase the statistical power of neuroimaging studies and enable comparative studies pertaining to ...several brain disorders. However, dMRI data sets acquired on multiple scanners cannot be naively pooled for joint analysis due to scanner specific nonlinear effects as well as differences in acquisition parameters. Consequently, for joint analysis, the dMRI data has to be harmonized, which involves removing scanner-specific differences from the raw dMRI signal. In this work, we propose a dMRI harmonization method that is capable of removing scanner-specific effects, while accounting for minor differences in acquisition parameters such as b-value, spatial resolution and number of gradient directions. We validate our algorithm on dMRI data acquired from two sites: Philadelphia Neurodevelopmental Cohort (PNC) with 800 healthy adolescents (ages 8–22 years) and Brigham and Women's Hospital (BWH) with 70 healthy subjects (ages 14–54 years). In particular, we show that gender and age-related maturation differences in different age groups are preserved after harmonization, as measured using effect sizes (small, medium and large), irrespective of the test sample size. Since we use matched control subjects from different scanners to estimate scanner-specific effects, our goal in this work is also to determine the minimum number of well-matched subjects needed from each site to achieve best harmonization results. Our results indicate that at-least 16 to 18 well-matched healthy controls from each site are needed to reliably capture scanner related differences. The proposed method can thus be used for retrospective harmonization of raw dMRI data across sites despite differences in acquisition parameters, while preserving inter-subject anatomical variability.
•A multi-site diffusion MRI harmonization method that can remove scanner-specific effects across sites.•The proposed method accounts for minor differences in acquisition parameters such as b-value, spatial resolution and number of gradient directions.•Biological differences due to gender and age in different groups are preserved.•At least 16 to 18 well-matched healthy controls from each site are needed to reliably capture scanner related differences.
Several prominent theories of schizophrenia suggest that structural white matter pathologies may follow a developmental, maturational, and/or degenerative process. However, a lack of lifespan studies ...has precluded verification of these theories. Here, we analyze the largest sample of carefully harmonized diffusion MRI data to comprehensively characterize age-related white matter trajectories, as measured by fractional anisotropy (FA), across the course of schizophrenia. Our analysis comprises diffusion scans of 600 schizophrenia patients and 492 healthy controls at different illness stages and ages (14-65 years), which were gathered from 13 sites. We determined the pattern of age-related FA changes by cross-sectionally assessing the timing of the structural neuropathology associated with schizophrenia. Quadratic curves were used to model between-group FA differences across whole-brain white matter and fiber tracts at each age; fiber tracts were then clustered according to both the effect-sizes and pattern of lifespan white matter FA differences. In whole-brain white matter, FA was significantly lower across the lifespan (up to 7%; p < 0.0033) and reached peak maturation younger in patients (27 years) compared to controls (33 years). Additionally, three distinct patterns of neuropathology emerged when investigating white matter fiber tracts in patients: (1) developmental abnormalities in limbic fibers, (2) accelerated aging and abnormal maturation in long-range association fibers, (3) severe developmental abnormalities and accelerated aging in callosal fibers. Our findings strongly suggest that white matter in schizophrenia is affected across entire stages of the disease. Perhaps most strikingly, we show that white matter changes in schizophrenia involve dynamic interactions between neuropathological processes in a tract-specific manner.
•A novel deep learning tract segmentation method with a new 2D multi-channel fiber feature descriptor.•Consistent tractography segmentation of subjects across the full lifespan.•Good generalization ...to tractography data from multiple different fiber tracking methods.
Display omitted
White matter tract segmentation, i.e. identifying tractography fibers (streamline trajectories) belonging to anatomically meaningful fiber tracts, is an essential step to enable tract quantification and visualization. In this study, we present a deep learning tractography segmentation method (DeepWMA) that allows fast and consistent identification of 54 major deep white matter fiber tracts from the whole brain. We create a large-scale training tractography dataset of 1 million labeled fiber samples, and we propose a novel 2D multi-channel feature descriptor (FiberMap) that encodes spatial coordinates of points along each fiber. We learn a convolutional neural network (CNN) fiber classification model based on FiberMap and obtain a high fiber classification accuracy of 90.99% on the training tractography data with ground truth fiber labels. Then, the method is evaluated on a test dataset of 597 diffusion MRI scans from six independently acquired populations across genders, the lifespan (1 day - 82 years), and different health conditions (healthy control, neuropsychiatric disorders, and brain tumor patients). We perform comparisons with two state-of-the-art tract segmentation methods. Experimental results show that our method obtains a highly consistent tract segmentation result, where on average over 99% of the fiber tracts are successfully identified across all subjects under study, most importantly, including neonates and patients with space-occupying brain tumors. We also demonstrate good generalization of the method to tractography data from multiple different fiber tracking methods. The proposed method leverages deep learning techniques and provides a fast and efficient tool for brain white matter segmentation in large diffusion MRI tractography datasets.
The heterogeneity of schizophrenia has defied efforts to derive reproducible and definitive anatomical maps of structural brain changes associated with the disorder. We aimed to map deviations from ...normative ranges of brain structure for individual patients and evaluate whether the loci of individual deviations recapitulated group-average brain maps of schizophrenia pathology. For each of 48 white matter tracts and 68 cortical regions, normative percentiles of variation in fractional anisotropy (FA) and cortical thickness (CT) were established using diffusion-weighted and structural MRI from healthy adults (n = 195). Individuals with schizophrenia (n = 322) were classified as either within the normative range for healthy individuals of the same age and sex (5-95% percentiles), infra-normal (<5% percentile) or supra-normal (>95% percentile). Repeating this classification for each tract and region yielded a deviation map for each individual. Compared to the healthy comparison group, the schizophrenia group showed widespread reductions in FA and CT, involving virtually all white matter tracts and cortical regions. Paradoxically, however, no more than 15-20% of patients deviated from the normative range for any single tract or region. Furthermore, 79% of patients showed infra-normal deviations for at least one locus (healthy individuals: 59 ± 2%, p < 0.001). Thus, while infra-normal deviations were common among patients, their anatomical loci were highly inconsistent between individuals. Higher polygenic risk for schizophrenia associated with a greater number of regions with infra-normal deviations in CT (r = -0.17, p = 0.006). We conclude that anatomical loci of schizophrenia-related changes are highly heterogeneous across individuals to the extent that group-consensus pathological maps are not representative of most individual patients. Normative modeling can aid in parsing schizophrenia heterogeneity and guiding personalized interventions.
Diffusion MRI studies consistently report group differences in white matter between individuals diagnosed with schizophrenia and healthy controls. Nevertheless, the abnormalities found at the ...group‐level are often not observed at the individual level. Among the different approaches aiming to study white matter abnormalities at the subject level, normative modeling analysis takes a step towards subject‐level predictions by identifying affected brain locations in individual subjects based on extreme deviations from a normative range. Here, we leveraged a large harmonized diffusion MRI dataset from 512 healthy controls and 601 individuals diagnosed with schizophrenia, to study whether normative modeling can improve subject‐level predictions from a binary classifier. To this aim, individual deviations from a normative model of standard (fractional anisotropy) and advanced (free‐water) dMRI measures, were calculated by means of age and sex‐adjusted z‐scores relative to control data, in 18 white matter regions. Even though larger effect sizes are found when testing for group differences in z‐scores than are found with raw values (p < .001), predictions based on summary z‐score measures achieved low predictive power (AUC < 0.63). Instead, we find that combining information from the different white matter tracts, while using multiple imaging measures simultaneously, improves prediction performance (the best predictor achieved AUC = 0.726). Our findings suggest that extreme deviations from a normative model are not optimal features for prediction. However, including the complete distribution of deviations across multiple imaging measures improves prediction, and could aid in subject‐level classification.
We use a large (n = 1113) harmonized diffusion MRI dataset of individuals diagnosed with schizophrenia and healthy controls, in order to study the predictive performance of a normative modeling approach. This approach compares each individual with distribution derived from healthy controls to derive measures of deviation from the normal distribution. Our results show that combining deviation information from different white matter tracts, while using multiple imaging measures simultaneously, improves prediction performance, and could therefore aid in subject‐level classification.
White matter (WM) abnormalities are repeatedly demonstrated across the schizophrenia time-course. However, our understanding of how demographic and clinical variables interact, influence, or are ...dependent on WM pathologies is limited. The most well-known barriers to progress are heterogeneous findings due to small sample sizes and the confounding influence of age on WM. The present study leverages access to the harmonized diffusion magnetic-resonance-imaging data and standardized clinical data from 13 international sites (597 schizophrenia patients (SCZ)). Fractional anisotropy (FA) values for all major WM structures in patients were predicted based on FA models estimated from a healthy population (n = 492). We utilized the deviations between predicted and real FA values to answer three essential questions. (1) "Which clinical variables explain WM abnormalities?". (2) "Does the degree of WM abnormalities predict symptom severity?". (3) "Does sex influence any of those relationships?". Regression and mediator analyses revealed that a longer duration-of-illness is associated with more severe WM abnormalities in several tracts. In addition, they demonstrated that a higher antipsychotic medication dose is related to more severe corpus callosum abnormalities. A structural equation model revealed that patients with more WM abnormalities display higher symptom severity. Last, the results exhibited sex-specificity. Males showed a stronger association between duration-of-illness and WM abnormalities. Females presented a stronger association between WM abnormalities and symptom severity, with IQ impacting this relationship. Our findings provide clear evidence for the interaction of demographic, clinical, and behavioral variables with WM pathology in SCZ. Our results also point to the need for longitudinal studies, directly investigating the casualty and sex-specificity of these relationships, as well as the impact of cognitive resiliency on structure-function relationships.
Characterization of anisotropy via diffusion MRI reveals fiber crossings in a substantial portion of voxels within the white-matter (WM) regions of the human brain. A considerable number of such ...voxels could exhibit asymmetric features such as bends and junctions. However, widely employed reconstruction methods yield symmetric Orientation Distribution Functions (ODFs) even when the underlying geometry is asymmetric. In this paper, we employ inter-voxel directional filtering approaches through a cone model to reveal more information regarding the cytoarchitectural organization within the voxel. The cone model facilitates a sharpening of the ODFs in some directions while suppressing peaks in other directions, thus yielding an Asymmetric ODF (AODF) field. We also show that a scalar measure of AODF asymmetry can be employed to obtain new contrast within the human brain. The feasibility of the technique is demonstrated on in vivo data obtained from the MGH-USC Human Connectome Project (HCP) and Parkinson's Progression Markers Initiative (PPMI) Project database. Characterizing asymmetry in neural tissue cytoarchitecture could be important for localizing and quantitatively assessing specific neuronal pathways.
Research suggests electroconvulsive therapy (ECT) induces an acute neuroinflammatory response and changes in white matter (WM) structural connectivity. However, whether these processes are related, ...either to each other or to eventual treatment outcomes, has yet to be determined. We examined the relationship between levels of peripheral pro-inflammatory cytokines and diffusion imaging-indexed changes in WM microstructure in individuals with treatment-resistant depression (TRD) who underwent ECT. Forty-two patients were assessed at baseline, after their second ECT (T2), and after completion of ECT (T3). A Montgomery Åsberg Depression Rating Scale improvement of >50% post-ECT defined ECT-responders (n = 19) from non-responders (n = 23). Thirty-four controls were also examined. Tissue-specific fractional anisotropy (FAt) was estimated using diffusion imaging data and the Free-Water method in 17 WM tracts. Inflammatory panels were evaluated from peripheral blood. Cytokines were examined to characterize the association between potential ECT-induced changes in an inflammatory state and WM microstructure. Longitudinal trajectories of both measures were also examined separately for ECT-responders and non-responders. Patients exhibited elevated Interleukin-8 (IL-8) levels at baseline compared to controls. In patients, correlations between IL-8 and FAt changes from baseline to T2 were significant in the positive direction in the right superior longitudinal fasciculus (R-SLF) and right cingulum (R-CB) (p
= 0.003). In these tracts, linear mixed-effects models revealed that trajectories of IL-8 and FAt were significantly positively correlated across all time points in responders, but not non-responders (R-CB-p = .001; R-SLF-p = 0.008). Our results suggest that response to ECT in TRD may be mediated by IL-8 and WM microstructure.
Ketamine is increasingly being used as a therapeutic for treatment-resistant depression (TRD), yet the effects of ketamine on the human brain remain largely unknown. This pilot study employed ...diffusion magnetic resonance imaging (dMRI) to examine relationships between ketamine treatment and white matter (WM) microstructure, with the aim of increasing the current understanding of ketamine's neural mechanisms of action in humans. Longitudinal dMRI data were acquired from 13 individuals with TRD two hours prior to (pre-infusion), and four hours following (post-infusion), an intravenous ketamine infusion. Free-water imaging was employed to quantify cerebrospinal fluid-corrected mean fractional anisotropy (FA) in 15 WM bundles pre- and post-infusion. Analyses revealed that higher pre-infusion FA in the left cingulum bundle and the left superior longitudinal fasciculus was associated with greater depression symptom improvement 24 h post-ketamine. Moreover, four hours after intravenous administration of ketamine, FA rapidly increased in numerous WM bundles in the brain; this increase was significantly associated with 24 h symptom improvement in select bundles. Overall, the results of this preliminary study suggest that WM properties, as measured by dMRI, may have a potential impact on clinical improvement following ketamine. Ketamine administration additionally appears to be associated with rapid WM diffusivity changes, suggestive of rapid changes in WM microstructure. This study thus points to pre-treatment WM structure as a potential factor associated with ketamine's clinical efficacy, and to post-treatment microstructural changes as a candidate neuroimaging marker of ketamine's cellular mechanisms.
Subtle alterations in white matter microstructure are observed in youth at clinical high risk (CHR) for psychosis. However, the timing of these changes and their relationships to the emergence of ...psychosis remain unclear. Here, we track the evolution of white matter abnormalities in a large, longitudinal cohort of CHR individuals comprising the North American Prodrome Longitudinal Study (NAPLS-3). Multi-shell diffusion magnetic resonance imaging data were collected across multiple timepoints (1-5 over 1 year) in 286 subjects (aged 12-32 years): 25 CHR individuals who transitioned to psychosis (CHR-P; 61 scans), 205 CHR subjects with unknown transition outcome after the 1-year follow-up period (CHR-U; 596 scans), and 56 healthy controls (195 scans). Linear mixed effects models were fitted to infer the impact of age and illness-onset on variation in the fractional anisotropy of cellular tissue (FA
) and the volume fraction of extracellular free water (FW). Baseline measures of white matter microstructure did not differentiate between HC, CHR-U and CHR-P individuals. However, age trajectories differed between the three groups in line with a developmental effect: CHR-P and CHR-U groups displayed higher FA
in adolescence, and 4% lower FA
by 30 years of age compared to controls. Furthermore, older CHR-P subjects (20+ years) displayed 4% higher FW in the forceps major (p < 0.05). Prospective analysis in CHR-P did not reveal a significant impact of illness onset on regional FA
or FW, suggesting that transition to psychosis is not marked by dramatic change in white matter microstructure. Instead, clinical high risk for psychosis-regardless of transition outcome-is characterized by subtle age-related white matter changes that occur in tandem with development.