Retrospective harmonization of multi-site diffusion MRI data acquired with different acquisition parameters

MidRISH: Unbiased harmonization of rotationally invariant harmonics of the diffusion signal

Data harmonization is necessary for removing confounding effects in multi-site diffusion image analysis. One such harmonization method, LinearRISH, scales rotationally invariant spherical harmonic (RISH) features from one site ("target") to the second ("reference") to reduce confounding scanner effects. However, reference and target site designations are not arbitrary and resultant diffusion metrics (fractional anisotropy, mean diffusivity) are biased by this choice. In this work we propose MidRISH: rather than scaling reference RISH features to target RISH features, we project both sites to a mid-space. We validate MidRISH with the following experiments: harmonizing scanner differences from 37 matched patients free of cognitive impairment, and harmonizing acquisition and study differences on 117 matched patients free of cognitive impairment. We find that MidRISH reduces bias of reference selection while preserving harmonization efficacy of LinearRISH. Users should be cautious when performing LinearRISH harmonization. To select a reference site is to choose diffusion metric effect-size. Our proposed method eliminates the bias-inducing site selection step.

Voxel-Wise Fusion of 3T and 7T Diffusion MRI Data to Extract more Accurate Fiber Orientations

While 7T diffusion magnetic resonance imaging (dMRI) has high spatial resolution, its diffusion imaging quality is usually affected by signal loss due to B1 inhomogeneity, T2 decay, susceptibility, and chemical shift. In contrast, 3T dMRI has relative higher diffusion angular resolution, but lower spatial resolution. Combination of 3T and 7T dMRI, thus, may provide more detailed and accurate information about the voxel-wise fiber orientations to better understand the structural brain connectivity. However, this topic has not yet been thoroughly explored until now. In this study, we explored the feasibility of fusing 3T and 7T dMRI data to extract voxel-wise quantitative parameters at higher spatial resolution. After 3T and 7T dMRI data was preprocessed, respectively, 3T dMRI volumes were coregistered into 7T dMRI space. Then, 7T dMRI data was harmonized to the coregistered 3T dMRI B0 (b = 0) images. Last, harmonized 7T dMRI data was fused with 3T dMRI data according to four fusion rules proposed in this study. We employed high-quality 3T and 7T dMRI datasets (N = 24) from the Human Connectome Project to test our algorithms. The diffusion tensors (DTs) and orientation distribution functions (ODFs) estimated from the 3T-7T fused dMRI volumes were statistically analyzed. More voxels containing multiple fiber populations were found from the fused dMRI data than from 7T dMRI data set. Moreover, extra fiber directions were extracted in temporal brain regions from the fused dMRI data at Otsu's thresholds of quantitative anisotropy, but could not be extracted from 7T dMRI dataset. This study provides novel algorithms to fuse intra-subject 3T and 7T dMRI data for extracting more detailed information of voxel-wise quantitative parameters, and a new perspective to build more accurate structural brain networks.

Impact of harmonization on the reproducibility of MRI radiomic features when using different scanners, acquisition parameters, and image pre-processing techniques: a phantom study

This study investigated the impact of ComBat harmonization on the reproducibility of radiomic features extracted from magnetic resonance images (MRI) acquired on different scanners, using various data acquisition parameters and multiple image pre-processing techniques using a dedicated MRI phantom. Four scanners were used to acquire an MRI of a nonanatomic phantom as part of the TCIA RIDER database. In fast spin-echo inversion recovery (IR) sequences, several inversion durations were employed, including 50, 100, 250, 500, 750, 1000, 1500, 2000, 2500, and 3000 ms. In addition, a 3D fast spoiled gradient recalled echo (FSPGR) sequence was used to investigate several flip angles (FA): 2, 5, 10, 15, 20, 25, and 30 degrees. Nineteen phantom compartments were manually segmented. Different approaches were used to pre-process each image: Bin discretization, Wavelet filter, Laplacian of Gaussian, logarithm, square, square root, and gradient. Overall, 92 first-, second-, and higher-order statistical radiomic features were extracted. ComBat harmonization was also applied to the extracted radiomic features. Finally, the Intraclass Correlation Coefficient (ICC) and Kruskal-Wallis's (KW) tests were implemented to assess the robustness of radiomic features. The number of non-significant features in the KW test ranged between 0-5 and 29-74 for various scanners, 31-91 and 37-92 for three times tests, 0-33 to 34-90 for FAs, and 3-68 to 65-89 for IRs before and after ComBat harmonization, with different image pre-processing techniques, respectively. The number of features with ICC over 90% ranged between 0-8 and 6-60 for various scanners, 11-75 and 17-80 for three times tests, 3-83 to 9-84 for FAs, and 3-49 to 3-63 for IRs before and after ComBat harmonization, with different image pre-processing techniques, respectively. The use of various scanners, IRs, and FAs has a great impact on radiomic features. However, the majority of scanner-robust features is also robust to IR and FA. Among the effective parameters in MR images, several tests in one scanner have a negligible impact on radiomic features. Different scanners and acquisition parameters using various image pre-processing might affect radiomic features to a large extent. ComBat harmonization might significantly impact the reproducibility of MRI radiomic features.

Perfusion imaging by arterial spin labeling in migraine: A literature review

Arterial spin labeling (ASL) is a non-invasive magnetic resonance imaging (MRI) method for the assessment of cerebral blood flow (CBF). This review summarizes recent ASL-based investigations in adult and pediatric patients with migraine with aura, migraine without aura, and chronic migraine. A systematic search according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines was conducted within PubMed and reference sections of articles identified from April 2014 to November 2022. Out of 236 initial articles, 20 remained after filtering, encompassing data from 1155 subjects in total. Cross-sectional studies in adults showed inconsistent results, while longitudinal studies demonstrated that cerebral perfusion changes over the migraine cycle can be tracked using ASL. The most consistent findings were observed in ictal states among pediatric migraine patients, where studies showed hypoperfusion matching aura symptoms during early imaging followed by hyperperfusion. Overall, ASL is a useful but currently underutilized modality for evaluating cerebral perfusion in patients with migraine. The generalizability of results is currently limited by heterogeneities regarding study design and documentation of clinical variables (e.g., relation of attacks to scanning timepoint, migraine subtypes). Future MRI studies should consider augmenting imaging protocols with ASL to further elucidate perfusion dynamics in migraine.

An Introduction to the Human Connectome Project for Early Psychosis

BACKGROUND

The time following a recent onset of psychosis is a critical period during which intervention may be maximally effective. Studying individuals in this period also offers an opportunity to investigate putative brain biomarkers of illness prior to the long-term effects of chronicity and medication. The Human Connectome Project for Early Psychosis (HCP-EP) was funded by the National Institutes of Mental Health (NIMH) as an extension of the original Human Connectome Project's approach to understanding the human brain and its structural and functional connections.

DESIGN

The HCP-EP data were collected at 3 sites in Massachusetts (Beth Israel Deaconess Medical Center, McLean Hospital, and Massachusetts General Hospital), and one site in Indiana (Indiana University). Brigham and Women's Hospital served as the data coordination center and as an imaging site.

RESULTS

The HCP-EP dataset includes high-quality clinical, cognitive, functional, neuroimaging, and blood specimen data acquired from 303 individuals between the ages of 16-35 years old with affective psychosis (n = 75), non-affective psychosis (n = 148), and healthy controls (n = 80). Participants with early psychosis were within 5 years of illness onset (mean duration = 1.9 years, standard deviation = 1.4 years). All data and novel or modified analytic tools developed as part of the study are publicly available to the research community through the NIMH Data Archive (NDA) or GitHub (https://github.com/pnlbwh).

CONCLUSIONS

This paper provides an overview of the specific HCP-EP procedures, assessments, and protocols, as well as a brief characterization of the study participants to make it easier for researchers to use this rich dataset. Although we focus here on discussing and comparing affective and non-affective psychosis groups, the HCP-EP dataset also provides sufficient information for investigators to group participants differently.

A diffusion MRI tractography atlas for concurrent white matter mapping across Eastern and Western populations

The study of brain differences across Eastern and Western populations provides vital insights for understanding potential cultural and genetic influences on cognition and mental health. Diffusion MRI (dMRI) tractography is an important tool in assessing white matter (WM) connectivity and brain tissue microstructure across different populations. However, a comprehensive investigation into WM fiber tracts between Eastern and Western populations is challenged due to the lack of a cross-population WM atlas and the large site-specific variability of dMRI data. This study presents a dMRI tractography atlas, namely the East-West WM Atlas, for concurrent WM mapping between Eastern and Western populations and creates a large, harmonized dMRI dataset (n=306) based on the Human Connectome Project and the Chinese Human Connectome Project. The curated WM atlas, as well as subject-specific data including the harmonized dMRI data, the whole brain tractography data, and parcellated WM fiber tracts and their diffusion measures, are publicly released. This resource is a valuable addition to facilitating the exploration of brain commonalities and differences across diverse cultural backgrounds.

Tractometry of Human Visual White Matter Pathways in Health and Disease

Diffusion-weighted MRI (dMRI) provides a unique non-invasive view of human brain tissue properties. The present review article focuses on tractometry analysis methods that use dMRI to assess the properties of brain tissue within the long-range connections comprising brain networks. We focus specifically on the major white matter tracts that convey visual information. These connections are particularly important because vision provides rich information from the environment that supports a large range of daily life activities. Many of the diseases of the visual system are associated with advanced aging, and tractometry of the visual system is particularly important in the modern aging society. We provide an overview of the tractometry analysis pipeline, which includes a primer on dMRI data acquisition, voxelwise model fitting, tractography, recognition of white matter tracts, and calculation of tract tissue property profiles. We then review dMRI-based methods for analyzing visual white matter tracts: the optic nerve, optic tract, optic radiation, forceps major, and vertical occipital fasciculus. For each tract, we review background anatomical knowledge together with recent findings in tractometry studies on these tracts and their properties in relation to visual function and disease. Overall, we find that measurements of the brain's visual white matter are sensitive to a range of disorders and correlate with perceptual abilities. We highlight new and promising analysis methods, as well as some of the current barriers to progress toward integration of these methods into clinical practice. These barriers, such as variability in measurements between protocols and instruments, are targets for future development.

Reduced cross-scanner variability using vendor-agnostic sequences for single-shell diffusion MRI

PURPOSE

To reduce the inter-scanner variability of diffusion MRI (dMRI) measures between scanners from different vendors by developing a vendor-neutral dMRI pulse sequence using the open-source vendor-agnostic Pulseq platform.

METHODS

We implemented a standard EPI based dMRI sequence in Pulseq. We tested it on two clinical scanners from different vendors (Siemens Prisma and GE Premier), systematically evaluating and comparing the within- and inter-scanner variability across the vendors, using both the vendor-provided and Pulseq dMRI sequences. Assessments covered both a diffusion phantom and three human subjects, using standard error (SE) and Lin's concordance correlation to measure the repeatability and reproducibility of standard DTI metrics including fractional anisotropy (FA) and mean diffusivity (MD).

RESULTS

Identical dMRI sequences were executed on both scanners using Pulseq. On the phantom, the Pulseq sequence showed more than a 2.5× reduction in SE (variability) across Siemens and GE scanners. Furthermore, Pulseq sequences exhibited markedly reduced SE in-vivo, maintaining scan-rescan repeatability while delivering lower variability in FA and MD (more than 50% reduction in cortical/subcortical regions) compared to vendor-provided sequences.

CONCLUSION

The Pulseq diffusion sequence reduces the cross-scanner variability for both phantom and in-vivo data, which will benefit multi-center neuroimaging studies and improve the reproducibility of neuroimaging studies.

TractoSCR: a novel supervised contrastive regression framework for prediction of neurocognitive measures using multi-site harmonized diffusion MRI tractography

Neuroimaging-based prediction of neurocognitive measures is valuable for studying how the brain's structure relates to cognitive function. However, the accuracy of prediction using popular linear regression models is relatively low. We propose a novel deep regression method, namely TractoSCR, that allows full supervision for contrastive learning in regression tasks using diffusion MRI tractography. TractoSCR performs supervised contrastive learning by using the absolute difference between continuous regression labels (i.e., neurocognitive scores) to determine positive and negative pairs. We apply TractoSCR to analyze a large-scale dataset including multi-site harmonized diffusion MRI and neurocognitive data from 8,735 participants in the Adolescent Brain Cognitive Development (ABCD) Study. We extract white matter microstructural measures using a fine parcellation of white matter tractography into fiber clusters. Using these measures, we predict three scores related to domains of higher-order cognition (general cognitive ability, executive function, and learning/memory). To identify important fiber clusters for prediction of these neurocognitive scores, we propose a permutation feature importance method for high-dimensional data. We find that TractoSCR obtains significantly higher accuracy of neurocognitive score prediction compared to other state-of-the-art methods. We find that the most predictive fiber clusters are predominantly located within the superficial white matter and projection tracts, particularly the superficial frontal white matter and striato-frontal connections. Overall, our results demonstrate the utility of contrastive representation learning methods for regression, and in particular for improving neuroimaging-based prediction of higher-order cognitive abilities. Our code will be available at: https://github.com/SlicerDMRI/TractoSCR.

Topological abnormalities of the morphometric similarity network of the cerebral cortex in schizophrenia

A morphometric similarity (MS) network can be constructed using multiple magnetic resonance imaging parameters of each cortical region. An MS network can be used to assess the similarity between cortical regions. Although MS networks can detect microstructural alterations and capture connections between histologically similar cortical areas, the influence of schizophrenia on the topological characteristics of MS networks remains unclear. We obtained T1- and diffusion-weighted images of 239 healthy controls and 190 individuals with schizophrenia to construct the MS network. Group comparisons of the mean MS of the cortical regions and subnetworks were performed. The strengths of the connections between the cortical regions and the global and nodal network indices were compared between the groups. Clinical associations with the network indices were tested using Spearman's rho. Compared with healthy controls, individuals with schizophrenia had significant group differences in the mean MS of several cortical regions and subnetworks. Individuals with schizophrenia had both superior and inferior strengths of connections between cortical regions compared with those of healthy controls. We observed regional abnormalities of the MS network in individuals with schizophrenia regarding lower centrality values of the pars opercularis, superior frontal, and superior temporal areas. Specific nodal network measures of the right pars opercularis and left superior temporal areas were associated with illness duration in individuals with schizophrenia. We identified regional abnormalities of the MS network in schizophrenia with the left superior temporal area possibly being a key region in topological organization and cortical connections.

Subcortical-cortical white matter connectivity in adults with autism spectrum disorder and schizophrenia patients

Autism spectrum disorder (ASD) and schizophrenia (SZ) are neuropsychiatric disorders that overlap in symptoms associated with social-cognitive impairment. Alterations of the cingulate cortex, subcortical, medial-temporal, and orbitofrontal structures are frequently reported in both disorders. In this study, we examined white-matter connectivity between these structures in adults with ASD and SZ patients compared with their respective neurotypical controls and indirectly with each other, using probabilistic and local DTI tractography. This exploratory study utilized publicly available neuroimaging databases, of adults with ASD (ABIDE II; n = 28) and SZ (COBRE; n = 38), age-gender matched neurotypicals (NT) and associated phenotypic data. Tractography was performed using Freesurfer and MRtrix software, and diffusion metrics of white-matter tracts between cingulate-, orbitofrontal- cortices, subcortical structures, parahippocampal, entorhinal cortex were assessed. In ASD, atypical diffusivity parameters were found in the isthmus cingulate and parahippocampal connectivity to subcortical and rostral-anterior cingulate, which were also associated with IQ and social skills (SRS). In contrast, atypical diffusivity parameters were observed between the medial-orbitofrontal cortex and subcortical structures in SZ, and were associated with executive function (i.e., IQ, processing speed) and emotional regulation. Overall, the results suggest that defects in the isthmus cingulate, medial-orbitofrontal, and striato-limbic white matter connectivity may help unravel the neural underpinnings of executive and social-emotional dysfunction at the core of neuropsychiatric disorders.

Protocol for the Chinese Real-World Evidence for Acute Spinal Cord Injury (ChiRES) study: a prospective, observational, multicentre cohort study of acute spinal cord injury

INTRODUCTION

Spinal cord injury (SCI) is a catastrophic event with devastating physical, social and occupational consequences for patients and their families. The number of patients with acute SCI in China continues to grow rapidly, but there have been no large prospective cohort studies of patients with acute SCI. This proposed study aims to establish a multicentre, extensive sample cohort of clinical data and biological samples of patients in China, which would aid the systematisation and standardisation of clinical research and treatment of acute SCI, thus reducing the heavy burden of acute SCI on patients and society.

METHODS AND ANALYSIS

The Chinese Real-World Evidence for Acute Spinal Cord Injury (ChiRES) study is an observational, multicentre cohort study of patients with acute SCI admitted to the Qilu Hospital of Shandong University and other participating centres with prospective collection of their clinical data and biological samples. We aim to recruit 2097 patients in this study. Demographics, disease history, emergency intervention information, motor and sensory examinations, surgical information, medication information and rehabilitation evaluation will be recorded. This will facilitate the development of a prediction model for complications and prognosis of patients with acute SCI and an evaluation of the current management of acute SCI. Among these variables, detailed information on surgical treatment will also be used to assess procedures for acute SCI treatment. Outcome measurements, including the International Standard for Neurological Classification of Spinal Cord Injury examinations, the occurrence of complications and death, will be performed repeatedly during follow-up. We will analyse imaging data and blood samples to develop SCI imaging markers and biomarkers.

ETHICS AND DISSEMINATION

This study protocol has been approved by the Medical Ethics Committee of the Qilu Hospital of Shandong University and all other participating centres. The findings will be disseminated in peer-reviewed journals and academic conferences.

Accelerating Medicines Partnership® Schizophrenia (AMP® SCZ): Rationale and Study Design of the Largest Global Prospective Cohort Study of Clinical High Risk for Psychosis

This article describes the rationale, aims, and methodology of the Accelerating Medicines Partnership® Schizophrenia (AMP® SCZ). This is the largest international collaboration to date that will develop algorithms to predict trajectories and outcomes of individuals at clinical high risk (CHR) for psychosis and to advance the development and use of novel pharmacological interventions for CHR individuals. We present a description of the participating research networks and the data processing analysis and coordination center, their processes for data harmonization across 43 sites from 13 participating countries (recruitment across North America, Australia, Europe, Asia, and South America), data flow and quality assessment processes, data analyses, and the transfer of data to the National Institute of Mental Health (NIMH) Data Archive (NDA) for use by the research community. In an expected sample of approximately 2000 CHR individuals and 640 matched healthy controls, AMP SCZ will collect clinical, environmental, and cognitive data along with multimodal biomarkers, including neuroimaging, electrophysiology, fluid biospecimens, speech and facial expression samples, novel measures derived from digital health technologies including smartphone-based daily surveys, and passive sensing as well as actigraphy. The study will investigate a range of clinical outcomes over a 2-year period, including transition to psychosis, remission or persistence of CHR status, attenuated positive symptoms, persistent negative symptoms, mood and anxiety symptoms, and psychosocial functioning. The global reach of AMP SCZ and its harmonized innovative methods promise to catalyze the development of new treatments to address critical unmet clinical and public health needs in CHR individuals.

Diffusion MRI harmonization via personalized template mapping

One fundamental challenge in diffusion magnetic resonance imaging (dMRI) harmonization is to disentangle the contributions of scanner-related effects from the variable brain anatomy for the observed imaging signals. Conventional harmonization methods rely on establishing an atlas space to resolve anatomical variability and generate a unified inter-site mapping function. However, this approach is limited in accounting for the misalignment of neuroanatomy that still widely persists even after registration, especially in regions close to cortical boundaries. To overcome this challenge, we propose a personalized framework in this paper to more effectively address the confounding from the misalignment of neuroanatomy in dMRI harmonization. Instead of using a common template representing site-effects for all subjects, the main novelty of our method is the adaptive computation of personalized templates for both source and target scanning sites to estimate the inter-site mapping function. We integrate our method with the rotation invariant spherical harmonics (RISH) features to achieve the harmonization of dMRI signals between sites. In our experiments, the proposed approach is applied to harmonize the dMRI data acquired from two scanning platforms: Siemens Prisma and GE MR750 from the Adolescent Brain Cognitive Development dataset and compared with a state-of-the-art method based on RISH features. Our results indicate that the proposed harmonization framework achieves superior performance not only in reducing inter-site variations due to scanner differences but also in preserving sex-related biological variability in original cohorts. Moreover, we assess the impact of harmonization on the estimation of fiber orientation distributions and show the robustness of the personalized harmonization procedure in preserving the fiber orientation of original dMRI signals.

Executive functioning, behavior, and white matter microstructure in the chronic phase after pediatric mild traumatic brain injury: results from the adolescent brain cognitive development study

BACKGROUND

Mild traumatic brain injury (mTBI) is common in children. Long-term cognitive and behavioral outcomes as well as underlying structural brain alterations following pediatric mTBI have yet to be determined. In addition, the effect of age-at-injury on long-term outcomes is largely unknown.

METHODS

Children with a history of mTBI (n = 406; Mage = 10 years, SDage = 0.63 years) who participated in the Adolescent Brain Cognitive Development (ABCD) study were matched (1:2 ratio) with typically developing children (TDC; n = 812) and orthopedic injury (OI) controls (n = 812). Task-based executive functioning, parent-rated executive functioning and emotion-regulation, and self-reported impulsivity were assessed cross-sectionally. Regression models were used to examine the effect of mTBI on these domains. The effect of age-at-injury was assessed by comparing children with their first mTBI at either 0-3, 4-7, or 8-10 years to the respective matched TDC controls. Fractional anisotropy (FA) and mean diffusivity (MD), both MRI-based measures of white matter microstructure, were compared between children with mTBI and controls.

RESULTS

Children with a history of mTBI displayed higher parent-rated executive dysfunction, higher impulsivity, and poorer self-regulation compared to both control groups. At closer investigation, these differences to TDC were only present in one respective age-at-injury group. No alterations were found in task-based executive functioning or white matter microstructure.

CONCLUSIONS

Findings suggest that everyday executive function, impulsivity, and emotion-regulation are affected years after pediatric mTBI. Outcomes were specific to the age at which the injury occurred, suggesting that functioning is differently affected by pediatric mTBI during vulnerable periods. Groups did not differ in white matter microstructure.

Empirical assessment of the assumptions of ComBat with diffusion tensor imaging

Purpose

Diffusion tensor imaging (DTI) is a magnetic resonance imaging technique that provides unique information about white matter microstructure in the brain but is susceptible to confounding effects introduced by scanner or acquisition differences. ComBat is a leading approach for addressing these site biases. However, despite its frequent use for harmonization, ComBat's robustness toward site dissimilarities and overall cohort size have not yet been evaluated in terms of DTI.

Approach

As a baseline, we match N = 358 participants from two sites to create a "silver standard" that simulates a cohort for multi-site harmonization. Across sites, we harmonize mean fractional anisotropy and mean diffusivity, calculated using participant DTI data, for the regions of interest defined by the JHU EVE-Type III atlas. We bootstrap 10 iterations at 19 levels of total sample size, 10 levels of sample size imbalance between sites, and 6 levels of mean age difference between sites to quantify (i) β AGE , the linear regression coefficient of the relationship between FA and age; (ii) γ ^ s f * , the ComBat-estimated site-shift; and (iii) δ ^ s f * , the ComBat-estimated site-scaling. We characterize the reliability of ComBat by evaluating the root mean squared error in these three metrics and examine if there is a correlation between the reliability of ComBat and a violation of assumptions.

Results

ComBat remains well behaved for β AGE when N > 162 and when the mean age difference is less than 4 years. The assumptions of the ComBat model regarding the normality of residual distributions are not violated as the model becomes unstable.

Conclusion

Prior to harmonization of DTI data with ComBat, the input cohort should be examined for size and covariate distributions of each site. Direct assessment of residual distributions is less informative on stability than bootstrap analysis. We caution use ComBat of in situations that do not conform to the above thresholds.

Harmonized diffusion MRI data and white matter measures from the Adolescent Brain Cognitive Development Study

The Adolescent Brain Cognitive Development (ABCD) Study® has collected data from over 10,000 children across 21 sites, providing insights into adolescent brain development. However, site-specific scanner variability has made it challenging to use diffusion MRI (dMRI) data from this study. To address this, a dataset of harmonized and processed ABCD dMRI data (from release 3) has been created, comprising quality-controlled imaging data from 9,345 subjects, focusing exclusively on the baseline session, i.e., the first time point of the study. This resource required substantial computational time (approx. 50,000 CPU hours) for harmonization, whole-brain tractography, and white matter parcellation. The dataset includes harmonized dMRI data, 800 white matter clusters, 73 anatomically labeled white matter tracts in full and low resolution, and 804 different dMRI-derived measures per subject (72.3 TB total size). Accessible via the NIMH Data Archive, it offers a large-scale dMRI dataset for studying structural connectivity in child and adolescent neurodevelopment. Additionally, several post-harmonization experiments were conducted to demonstrate the success of the harmonization process on the ABCD dataset.

Tractography-based automated identification of the retinogeniculate visual pathway with novel microstructure-informed supervised contrastive learning

The retinogeniculate visual pathway (RGVP) is responsible for carrying visual information from the retina to the lateral geniculate nucleus. Identification and visualization of the RGVP are important in studying the anatomy of the visual system and can inform the treatment of related brain diseases. Diffusion MRI (dMRI) tractography is an advanced imaging method that uniquely enables in vivo mapping of the 3D trajectory of the RGVP. Currently, identification of the RGVP from tractography data relies on expert (manual) selection of tractography streamlines, which is time-consuming, has high clinical and expert labor costs, and is affected by inter-observer variability. In this paper, we present a novel deep learning framework, DeepRGVP , to enable fast and accurate identification of the RGVP from dMRI tractography data. We design a novel microstructure-informed supervised contrastive learning method that leverages both streamline label and tissue microstructure information to determine positive and negative pairs. We propose a simple and successful streamline-level data augmentation method to address highly imbalanced training data, where the number of RGVP streamlines is much lower than that of non-RGVP streamlines. We perform comparisons with several state-of-the-art deep learning methods that were designed for tractography parcellation, and we show superior RGVP identification results using DeepRGVP. In addition, we demonstrate a good generalizability of DeepRGVP to dMRI tractography data from neurosurgical patients with pituitary tumors and we show DeepRGVP can successfully identify RGVPs despite the effect of lesions affecting the RGVPs. Overall, our study shows the high potential of using deep learning to automatically identify the RGVP.

CROSS-AGE AND CROSS-SITE DOMAIN SHIFT IMPACTS ON DEEP LEARNING-BASED WHITE MATTER FIBER ESTIMATION IN NEWBORN AND BABY BRAINS

Deep learning models have shown great promise in estimating tissue microstructure from limited diffusion magnetic resonance imaging data. However, these models face domain shift challenges when test and train data are from different scanners and protocols, or when the models are applied to data with inherent variations such as the developing brains of infants and children scanned at various ages. Several techniques have been proposed to address some of these challenges, such as data harmonization or domain adaptation in the adult brain. However, those techniques remain unexplored for the estimation of fiber orientation distribution functions in the rapidly developing brains of infants. In this work, we extensively investigate the age effect and domain shift within and across two different cohorts of 201 newborns and 165 babies using the Method of Moments and fine-tuning strategies. Our results show that reduced variations in the microstructural development of babies in comparison to newborns directly impact the deep learning models' cross-age performance. We also demonstrate that a small number of target domain samples can significantly mitigate domain shift problems.

Prediction of misfolded proteins spreading in Alzheimer's disease using machine learning and spreading models

The pervasive impact of Alzheimer's disease on aging society represents one of the main challenges at this time. Current investigations highlight 2 specific misfolded proteins in its development: Amyloid-$\beta$ and tau. Previous studies focused on spreading for misfolded proteins exploited simulations, which required several parameters to be empirically estimated. Here, we provide an alternative view based on 2 machine learning approaches which we compare with known simulation models. The first approach applies an autoregressive model constrained by structural connectivity, while the second is based on graph convolutional networks. The aim is to predict concentrations of Amyloid-$\beta$ 2 yr after a provided baseline. We also evaluate its real-world effectiveness and suitability by providing a web service for physicians and researchers. In experiments, the autoregressive model generally outperformed state-of-the-art models resulting in lower prediction errors. While it is important to note that a comprehensive prognostic plan cannot solely rely on amyloid beta concentrations, their prediction, achieved by the discussed approaches, can be valuable for planning therapies and other cures, especially when dealing with asymptomatic patients for whom novel therapies could prove effective.

Structural impairment in superficial and deep white matter in schizophrenia

OBJECTIVE

Although disconnectivity among brain regions has been one of the main hypotheses for schizophrenia, the superficial white matter (SWM) has received less attention in schizophrenia research than the deep white matter (DWM) owing to the challenge of consistent reconstruction across subjects.

METHODS

We obtained the diffusion magnetic resonance imaging (dMRI) data of 223 healthy controls and 143 patients with schizophrenia. After harmonising the raw dMRIs from three different studies, we performed whole-brain two-tensor tractography and fibre clustering on the tractography data. We compared the fractional anisotropy (FA) of white matter tracts between healthy controls and patients with schizophrenia. Spearman's rho was adopted for the associations with clinical symptoms measured by the Positive and Negative Syndrome Scale (PANSS). The Bonferroni correction was used to adjust multiple testing.

RESULTS

Among the 33 DWM and 8 SWM tracts, patients with schizophrenia had a lower FA in 14 DWM and 4 SWM tracts than healthy controls, with small effect sizes. In the patient group, the FA deviations of the corticospinal and superficial-occipital tracts were negatively correlated with the PANSS negative score; however, this correlation was not evident after adjusting for multiple testing.

CONCLUSION

We observed the structural impairments of both the DWM and SWM tracts in patients with schizophrenia. The SWM could be a potential target of interest in future research on neural biomarkers for schizophrenia.

MidRISH: Unbiased harmonization of rotationally invariant harmonics of the diffusion signal

Objective

Data harmonization is necessary for removing confounding effects in multi-site diffusion image analysis. One such harmonization method, LinearRISH, scales rotationally invariant spherical harmonic (RISH) features from one site ("target") to the second ("reference") to reduce confounding scanner effects. However, reference and target site designations are not arbitrary and resultant diffusion metrics (fractional anisotropy, mean diffusivity) are biased by this choice. In this work we propose MidRISH: rather than scaling reference RISH features to target RISH features, we project both sites to a mid-space.

Methods

We validate MidRISH with the following experiments: harmonizing scanner differences from 37 matched patients free of cognitive impairment, and harmonizing acquisition and study differences on 117 matched patients free of cognitive impairment.

Conclusion

MidRISH reduces bias of reference selection while preserving harmonization efficacy of LinearRISH.

Significance

Users should be cautious when performing LinearRISH harmonization. To select a reference site is to choose diffusion metric effect-size. Our proposed method eliminates the bias-inducing site selection step.

Normative values of the topological metrics of the structural connectome: A multi-site reproducibility study across the Italian Neuroscience network

PURPOSE

The use of topological metrics to derive quantitative descriptors from structural connectomes is receiving increasing attention but deserves specific studies to investigate their reproducibility and variability in the clinical context. This work exploits the harmonization of diffusion-weighted acquisition for neuroimaging data performed by the Italian Neuroscience and Neurorehabilitation Network initiative to obtain normative values of topological metrics and to investigate their reproducibility and variability across centers.

METHODS

Different topological metrics, at global and local level, were calculated on multishell diffusion-weighted data acquired at high-field (e.g. 3 T) Magnetic Resonance Imaging scanners in 13 different centers, following the harmonization of the acquisition protocol, on young and healthy adults. A "traveling brains" dataset acquired on a subgroup of subjects at 3 different centers was also analyzed as reference data. All data were processed following a common processing pipeline that includes data pre-processing, tractography, generation of structural connectomes and calculation of graph-based metrics. The results were evaluated both with statistical analysis of variability and consistency among sites with the traveling brains range. In addition, inter-site reproducibility was assessed in terms of intra-class correlation variability.

RESULTS

The results show an inter-center and inter-subject variability of <10%, except for "clustering coefficient" (variability of 30%). Statistical analysis identifies significant differences among sites, as expected given the wide range of scanners' hardware.

CONCLUSIONS

The results show low variability of connectivity topological metrics across sites running a harmonised protocol.

RELIEF: A structured multivariate approach for removal of latent inter-scanner effects

Combining data collected from multiple study sites is becoming common and is advantageous to researchers to increase the generalizability and replicability of scientific discoveries. However, at the same time, unwanted inter-scanner biases are commonly observed across neuroimaging data collected from multiple study sites or scanners, rendering difficulties in integrating such data to obtain reliable findings. While several methods for handling such unwanted variations have been proposed, most of them use univariate approaches that could be too simple to capture all sources of scanner-specific variations. To address these challenges, we propose a novel multivariate harmonization method called RELIEF (REmoval of Latent Inter-scanner Effects through Factorization) for estimating and removing both explicit and latent scanner effects. Our method is the first approach to introduce the simultaneous dimension reduction and factorization of interlinked matrices to a data harmonization context, which provides a new direction in methodological research for correcting inter-scanner biases. Analyzing diffusion tensor imaging (DTI) data from the Social Processes Initiative in Neurobiology of the Schizophrenia (SPINS) study and conducting extensive simulation studies, we show that RELIEF outperforms existing harmonization methods in mitigating inter-scanner biases and retaining biological associations of interest to increase statistical power. RELIEF is publicly available as an R package.

Image harmonization: A review of statistical and deep learning methods for removing batch effects and evaluation metrics for effective harmonization

Magnetic resonance imaging and computed tomography from multiple batches (e.g. sites, scanners, datasets, etc.) are increasingly used alongside complex downstream analyses to obtain new insights into the human brain. However, significant confounding due to batch-related technical variation, called batch effects, is present in this data; direct application of downstream analyses to the data may lead to biased results. Image harmonization methods seek to remove these batch effects and enable increased generalizability and reproducibility of downstream results. In this review, we describe and categorize current approaches in statistical and deep learning harmonization methods. We also describe current evaluation metrics used to assess harmonization methods and provide a standardized framework to evaluate newly-proposed methods for effective harmonization and preservation of biological information. Finally, we provide recommendations to end-users to advocate for more effective use of current methods and to methodologists to direct future efforts and accelerate development of the field.

The organization of frontostriatal brain wiring in non-affective early psychosis compared with healthy subjects using a novel diffusion imaging fiber cluster analysis

BACKGROUND

Alterations in brain connectivity may underlie neuropsychiatric conditions such as schizophrenia. We here assessed the degree of convergence of frontostriatal fiber projections in 56 young adult healthy controls (HCs) and 108 matched Early Psychosis-Non-Affective patients (EP-NAs) using our novel fiber cluster analysis of whole brain diffusion magnetic resonance imaging tractography.

METHODS

Using whole brain tractography and our fiber clustering methodology on harmonized diffusion magnetic resonance imaging data from the Human Connectome Project for Early Psychosis we identified 17 white matter fiber clusters that connect frontal cortex (FCtx) and caudate (Cd) per hemisphere in each group. To quantify the degree of convergence and, hence, topographical relationship of these fiber clusters, we measured the inter-cluster mean distances between the endpoints of the fiber clusters at the level of the FCtx and of the Cd, respectively.

RESULTS

We found (1) in both groups, bilaterally, a non-linear relationship, yielding convex curves, between FCtx and Cd distances for FCtx-Cd connecting fiber clusters, driven by a cluster projecting from inferior frontal gyrus; however, in the right hemisphere, the convex curve was more flattened in EP-NAs; (2) that cluster pairs in the right (p = 0.03), but not left (p = 0.13), hemisphere were significantly more convergent in HCs vs EP-NAs; (3) in both groups, bilaterally, similar clusters projected significantly convergently to the Cd; and, (4) a significant group by fiber cluster pair interaction for 2 right hemisphere fiber clusters (numbers 5, 11; p = .00023; p = .00023) originating in selective PFC subregions.

CONCLUSIONS

In both groups, we found the FCtx-Cd wiring pattern deviated from a strictly topographic relationship and that similar clusters projected significantly more convergently to the Cd. Interestingly, we also found a significantly more convergent pattern of connectivity in HCs in the right hemisphere and that 2 clusters from PFC subregions in the right hemisphere significantly differed in their pattern of connectivity between groups.

Harmonized diffusion MRI data and white matter measures from the Adolescent Brain Cognitive Development Study

The Adolescent Brain Cognitive Development (ABCD) study has collected data from over 10,000 children across 21 sites, providing valuable insights into adolescent brain development. However, site-specific scanner variability has made it challenging to use diffusion MRI (dMRI) data from this study. To address this, a database of harmonized and processed ABCD dMRI data has been created, comprising quality-controlled imaging data from 9345 subjects. This resource required significant computational effort, taking ~50,000 CPU hours to harmonize the data, perform white matter parcellation, and run whole brain tractography. The database includes harmonized dMRI data, 800 white matter clusters, 73 anatomically labeled white matter tracts both in full-resolution (for analysis) and low-resolution (for visualization), and 804 different dMRI-derived measures per subject. It is available via the NIMH Data Archive and offers tremendous potential for scientific discoveries in structural connectivity studies of neurodevelopment in children and adolescents. Additionally, several post-harmonization experiments were conducted to demonstrate the success of the harmonization process on the ABCD dataset.

Characterization of the extracellular free water signal in schizophrenia using multi-site diffusion MRI harmonization

Studies applying Free Water Imaging have consistently reported significant global increases in extracellular free water (FW) in populations of individuals with early psychosis. However, these published studies focused on homogenous clinical participant groups (e.g., only first episode or chronic), thereby limiting our understanding of the time course of free water elevations across illness stages. Moreover, the relationship between FW and duration of illness has yet to be directly tested. Leveraging our multi-site diffusion magnetic resonance imaging(dMRI) harmonization approach, we analyzed dMRI scans collected by 12 international sites from 441 healthy controls and 434 individuals diagnosed with schizophrenia-spectrum disorders at different illness stages and ages (15-58 years). We characterized the pattern of age-related FW changes by assessing whole brain white matter in individuals with schizophrenia and healthy controls. In individuals with schizophrenia, average whole brain FW was higher than in controls across all ages, with the greatest FW values observed from 15 to 23 years (effect size range = [0.70-0.87]). Following this peak, FW exhibited a monotonic decrease until reaching a minima at the age of 39 years. After 39 years, an attenuated monotonic increase in FW was observed, but with markedly smaller effect sizes when compared to younger patients (effect size range = [0.32-0.43]). Importantly, FW was found to be negatively associated with duration of illness in schizophrenia (p = 0.006), independent of the effects of other clinical and demographic data. In summary, our study finds in a large, age-diverse sample that participants with schizophrenia with a shorter duration of illness showed higher FW values compared to participants with more prolonged illness. Our findings provide further evidence that elevations in the FW are present in individuals with schizophrenia, with the greatest differences in the FW being observed in those at the early stages of the disorder, which might suggest acute extracellular processes.

Individualised profiling of white matter organisation in moderate-to-severe traumatic brain injury patients

BACKGROUND AND PURPOSE

Approximately 65% of moderate-to-severe traumatic brain injury (m-sTBI) patients present with poor long-term behavioural outcomes, which can significantly impair activities of daily living. Numerous diffusion-weighted MRI studies have linked these poor outcomes to decreased white matter integrity of several commissural tracts, association fibres and projection fibres in the brain. However, most studies have focused on group-based analyses, which are unable to deal with the substantial between-patient heterogeneity in m-sTBI. As a result, there is increasing interest and need in conducting individualised neuroimaging analyses.

MATERIALS AND METHODS

Here, we generated a detailed subject-specific characterisation of microstructural organisation of white matter tracts in 5 chronic patients with m-sTBI (29 - 49y, 2 females), presented as a proof-of-concept. We developed an imaging analysis framework using fixel-based analysis and TractLearn to determine whether the values of fibre density of white matter tracts at the individual patient level deviate from the healthy control group (n = 12, 8F, M_age = 35.7y, age range 25 - 64y).

RESULTS

Our individualised analysis revealed unique white matter profiles, confirming the heterogenous nature of m-sTBI and the need of individualised profiles to properly characterise the extent of injury. Future studies incorporating clinical data, as well as utilising larger reference samples and examining the test-retest reliability of the fixel-wise metrics are warranted.

CONCLUSIONS

Individualised profiles may assist clinicians in tracking recovery and planning personalised training programs for chronic m-sTBI patients, which is necessary to achieve optimal behavioural outcomes and improved quality of life.

Diffusion-Weighted Imaging in Mild Traumatic Brain Injury: A Systematic Review of the Literature

There is evidence that diffusion-weighted imaging (DWI) is able to detect tissue alterations following mild traumatic brain injury (mTBI) that may not be observed on conventional neuroimaging; however, findings are often inconsistent between studies. This systematic review assesses patterns of differences in DWI metrics between those with and without a history of mTBI. A PubMed literature search was performed using relevant indexing terms for articles published prior to May 14, 2020. Findings were limited to human studies using DWI in mTBI. Articles were excluded if they were not full-length, did not contain original data, if they were case studies, pertained to military populations, had inadequate injury severity classification, or did not report post-injury interval. Findings were reported independently for four subgroups: acute/subacute pediatric mTBI, acute/subacute adult mTBI, chronic adult mTBI, and sport-related concussion, and all DWI acquisition and analysis methods used were included. Patterns of findings between studies were reported, along with strengths and weaknesses of the current state of the literature. Although heterogeneity of sample characteristics and study methods limited the consistency of findings, alterations in DWI metrics were most commonly reported in the corpus callosum, corona radiata, internal capsule, and long association pathways. Many acute/subacute pediatric studies reported higher FA and lower ADC or MD in various regions. In contrast, acute/subacute adult studies most commonly indicate lower FA within the context of higher MD and RD. In the chronic phase of recovery, FA may remain low, possibly indicating overall demyelination or Wallerian degeneration over time. Longitudinal studies, though limited, generally indicate at least a partial normalization of DWI metrics over time, which is often associated with functional improvement. We conclude that DWI is able to detect structural mTBI-related abnormalities that may persist over time, although future DWI research will benefit from larger samples, improved data analysis methods, standardized reporting, and increasing transparency.

The Role of Neuroimaging in Evolving TBI Research and Clinical Practice

Neuroimaging technologies such as computed tomography (CT) and magnetic resonance imaging (MRI) have been widely adopted in the clinical diagnosis and management of traumatic brain injury (TBI), particularly at the more acute and severe levels of injury. Additionally, a number of advanced applications of MRI have been employed in TBI-related clinical research with great promise, and researchers have used these techniques to better understand underlying mechanisms, progression of secondary injury and tissue perturbation over time, and relation of focal and diffuse injury to later outcome. However, the acquisition and analysis time, the cost of these and other imaging modalities, and the need for specialized expertise have represented historical barriers in extending these tools in clinical practice. While group studies are important in detecting patterns, heterogeneity among patient presentation and limited sample sizes from which to compare individual level data to well-developed normative data have also played a role in the limited translatability of imaging to wider clinical application. Fortunately, the field of TBI has benefitted from increased public and scientific awareness of the prevalence and impact of TBI, particularly in head injury related to recent military conflicts and sport-related concussion. This awareness parallels an increase in federal funding in the United States and other countries allocated to investigation in these areas. In this article we summarize funding and publication trends since the mainstream adoption of imaging in TBI to elucidate evolving trends and priorities in the application of different techniques and patient populations. We also review recent and ongoing efforts to advance the field through promoting reproducibility, data sharing, big data analytic methods, and team science. Finally, we discuss international collaborative efforts to combine and harmonize neuroimaging, cognitive, and clinical data, both prospectively and retrospectively. Each of these represent unique, but related, efforts that facilitate closing gaps between the use of advanced imaging solely as a research tool and the use of it in clinical diagnosis, prognosis, and treatment planning and monitoring.

White matter microstructure shows sex differences in late childhood: Evidence from 6797 children

Sex differences in white matter microstructure have been robustly demonstrated in the adult brain using both conventional and advanced diffusion-weighted magnetic resonance imaging approaches. However, sex differences in white matter microstructure prior to adulthood remain poorly understood; previous developmental work focused on conventional microstructure metrics and yielded mixed results. Here, we rigorously characterized sex differences in white matter microstructure among over 6000 children from the Adolescent Brain Cognitive Development study who were between 9 and 10 years old. Microstructure was quantified using both the conventional model-diffusion tensor imaging (DTI)-and an advanced model, restriction spectrum imaging (RSI). DTI metrics included fractional anisotropy (FA) and mean, axial, and radial diffusivity (MD, AD, RD). RSI metrics included normalized isotropic, directional, and total intracellular diffusion (N0, ND, NT). We found significant and replicable sex differences in DTI or RSI microstructure metrics in every white matter region examined across the brain. Sex differences in FA were regionally specific. Across white matter regions, boys exhibited greater MD, AD, and RD than girls, on average. Girls displayed increased N0, ND, and NT compared to boys, on average, suggesting greater cell and neurite density in girls. Together, these robust and replicable findings provide an important foundation for understanding sex differences in health and disease.

Functional magnetic resonance imaging in migraine: A systematic review

BACKGROUND

Migraine is a highly prevalent primary headache disorder. Despite a high burden of disease, key disease mechanisms are not entirely understood. Functional magnetic resonance imaging is an imaging method using the blood-oxygen-level-dependent signal, which has been increasingly used in migraine research over recent years. This systematic review summarizes recent findings employing functional magnetic resonance imaging for the investigation of migraine.

METHODS

We conducted a systematic search and selection of functional magnetic resonance imaging applications in migraine from April 2014 to December 2021 (PubMed and references of identified articles according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines). Methodological details and main findings were extracted and synthesized.

RESULTS

Out of 224 articles identified, 114 were included after selection. Repeatedly emerging structures of interest included the insula, brainstem, limbic system, hypothalamus, thalamus, and functional networks. Assessment of functional brain changes in response to treatment is emerging, and machine learning has been used to investigate potential functional magnetic resonance imaging-based markers of migraine.

CONCLUSIONS

A wide variety of functional magnetic resonance imaging-based metrics were found altered across the brain for heterogeneous migraine cohorts, partially correlating with clinical parameters and supporting the concept to conceive migraine as a brain state. However, a majority of findings from previous studies have not been replicated, and studies varied considerably regarding image acquisition and analyses techniques. Thus, while functional magnetic resonance imaging appears to have the potential to advance our understanding of migraine pathophysiology, replication of findings in large representative datasets and precise, standardized reporting of clinical data would likely benefit the field and further increase the value of observations.

Posttraumatic Stress and Traumatic Brain Injury: Cognition, Behavior, and Neuroimaging Markers in Vietnam Veterans

BACKGROUND

Posttraumatic stress disorder (PTSD) and traumatic brain injury (TBI) are common in Veterans and linked to behavioral disturbances, increased risk of cognitive decline, and Alzheimer's disease.

OBJECTIVE

We studied the synergistic effects of PTSD and TBI on behavioral, cognitive, and neuroimaging measures in Vietnam war Veterans.

METHODS

Data were acquired at baseline and after about one-year from male Veterans categorized into: PTSD, TBI, PTSD+TBI, and Veteran controls without PTSD or TBI. We applied manual tractography to examine white matter microstructure of three fiber tracts: uncinate fasciculus (N = 91), cingulum (N = 87), and inferior longitudinal fasciculus (N = 95). ANCOVAs were used to compare Veterans' baseline behavioral and cognitive functioning (N = 285), white matter microstructure, amyloid-β (N = 230), and tau PET (N = 120). Additional ANCOVAs examined scores' differences from baseline to follow-up.

RESULTS

Veterans with PTSD and PTSD+TBI, but not Veterans with TBI only, exhibited poorer behavioral and cognitive functioning at baseline than controls. The groups did not differ in baseline white matter, amyloid-β, or tau, nor in behavioral and cognitive functioning, and tau accumulation change. Progression of white matter abnormalities of the uncinate fasciculus in Veterans with PTSD compared to controls was observed; analyses in TBI and PTSD+TBI were not run due to insufficient sample size.

CONCLUSIONS

PTSD and PTSD+TBI negatively affect behavioral and cognitive functioning, while TBI does not contribute independently. Whether progressive decline in uncinate fasciculus microstructure in Veterans with PTSD might account for cognitive decline should be further studied. Findings did not support an association between PTSD, TBI, and Alzheimer's disease pathology based on amyloid and tau PET.

Challenges for machine learning in clinical translation of big data imaging studies

Combining deep learning image analysis methods and large-scale imaging datasets offers many opportunities to neuroscience imaging and epidemiology. However, despite these opportunities and the success of deep learning when applied to a range of neuroimaging tasks and domains, significant barriers continue to limit the impact of large-scale datasets and analysis tools. Here, we examine the main challenges and the approaches that have been explored to overcome them. We focus on issues relating to data availability, interpretability, evaluation, and logistical challenges and discuss the problems that still need to be tackled to enable the success of "big data" deep learning approaches beyond research.

Improved sensitivity and precision in multicentre diffusion MRI network analysis using thresholding and harmonization

PURPOSE

To investigate if network thresholding and raw data harmonization improve consistency of diffusion MRI (dMRI)-based brain networks while also increasing precision and sensitivity to detect disease effects in multicentre datasets.

METHODS

Brain networks were reconstructed from dMRI of five samples with cerebral small vessel disease (SVD; 629 patients, 166 controls), as a clinically relevant exemplar condition for studies on network integrity. We evaluated consistency of network architecture in age-matched controls, by calculating cross-site differences in connection probability and fractional anisotropy (FA). Subsequently we evaluated precision and sensitivity to disease effects by identifying connections with low FA in sporadic SVD patients relative to controls, using more severely affected patients with a pure form of genetically defined SVD as reference.

RESULTS

In controls, thresholding and harmonization improved consistency of network architecture, minimizing cross-site differences in connection probability and FA. In patients relative to controls, thresholding improved precision to detect disrupted connections by removing false positive connections (precision, before: 0.09-0.19; after: 0.38-0.70). Before harmonization, sensitivity was low within individual sites, with few connections surviving multiple testing correction (k = 0-25 connections). Harmonization and pooling improved sensitivity (k = 38), while also achieving higher precision when combined with thresholding (0.97).

CONCLUSION

We demonstrated that network consistency, precision and sensitivity to detect disease effects in SVD are improved by thresholding and harmonization. We recommend introducing these techniques to leverage large existing multicentre datasets to better understand the impact of disease on brain networks.

How Machine Learning is Powering Neuroimaging to Improve Brain Health

This report presents an overview of how machine learning is rapidly advancing clinical translational imaging in ways that will aid in the early detection, prediction, and treatment of diseases that threaten brain health. Towards this goal, we aresharing the information presented at a symposium, "Neuroimaging Indicators of Brain Structure and Function - Closing the Gap Between Research and Clinical Application", co-hosted by the McCance Center for Brain Health at Mass General Hospital and the MIT HST Neuroimaging Training Program on February 12, 2021. The symposium focused on the potential for machine learning approaches, applied to increasingly large-scale neuroimaging datasets, to transform healthcare delivery and change the trajectory of brain health by addressing brain care earlier in the lifespan. While not exhaustive, this overview uniquely addresses many of the technical challenges from image formation, to analysis and visualization, to synthesis and incorporation into the clinical workflow. Some of the ethical challenges inherent to this work are also explored, as are some of the regulatory requirements for implementation. We seek to educate, motivate, and inspire graduate students, postdoctoral fellows, and early career investigators to contribute to a future where neuroimaging meaningfully contributes to the maintenance of brain health.

Diffusion Tensor Imaging Reveals Elevated Diffusivity of White Matter Microstructure that Is Independently Associated with Long-Term Outcome after Mild Traumatic Brain Injury: A TRACK-TBI Study

Diffusion tensor imaging (DTI) literature on single-center studies contains conflicting results regarding acute effects of mild traumatic brain injury (mTBI) on white matter (WM) microstructure and the prognostic significance. This larger-scale multi-center DTI study aimed to determine how acute mTBI affects WM microstructure over time and how early WM changes affect long-term outcome. From Transforming Research and Clinical Knowledge in Traumatic Brain Injury (TRACK-TBI), a cohort study at 11 United States level 1 trauma centers, a total of 391 patients with acute mTBI ages 17 to 60 years were included and studied at two weeks and six months post-injury. Demographically matched friends or family of the participants were the control group (n = 148). Axial diffusivity (AD), fractional anisotropy (FA), mean diffusivity (MD), and radial diffusivity (RD) were the measures of WM microstructure. The primary outcome was the Glasgow Outcome Scale Extended (GOSE) score of injury-related functional limitations across broad life domains at six months post-injury. The AD, MD, and RD were higher and FA was lower in mTBI versus friend control (FC) at both two weeks and six months post-injury throughout most major WM tracts of the cerebral hemispheres. In the mTBI group, AD and, to a lesser extent, MD decreased in WM from two weeks to six months post-injury. At two weeks post-injury, global WM AD and MD were both independently associated with six-month incomplete recovery (GOSE <8 vs = 8) even after accounting for demographic, clinical, and other imaging factors. DTI provides reliable imaging biomarkers of dynamic WM microstructural changes after mTBI that have utility for patient selection and treatment response in clinical trials. Continued technological advances in the sensitivity, specificity, and precision of diffusion magnetic resonance imaging hold promise for routine clinical application in mTBI.

Cognitive deficits, clinical variables, and white matter microstructure in schizophrenia: a multisite harmonization study

Cognitive deficits are among the best predictors of real-world functioning in schizophrenia. However, our understanding of how cognitive deficits relate to neuropathology and clinical presentation over the disease lifespan is limited. Here, we combine multi-site, harmonized cognitive, imaging, demographic, and clinical data from over 900 individuals to characterize a) cognitive deficits across the schizophrenia lifespan and b) the association between cognitive deficits, clinical presentation, and white matter (WM) microstructure. Multimodal harmonization was accomplished using T-scores for cognitive data, previously reported standardization methods for demographic and clinical data, and an established harmonization method for imaging data. We applied t-tests and correlation analysis to describe cognitive deficits in individuals with schizophrenia. We then calculated whole-brain WM fractional anisotropy (FA) and utilized regression-mediation analyses to model the association between diagnosis, FA, and cognitive deficits. We observed pronounced cognitive deficits in individuals with schizophrenia (p < 0.006), associated with more positive symptoms and medication dosage. Regression-mediation analyses showed that WM microstructure mediated the association between schizophrenia and language/processing speed/working memory/non-verbal memory. In addition, processing speed mediated the influence of diagnosis and WM microstructure on the other cognitive domains. Our study highlights the critical role of cognitive deficits in schizophrenia. We further show that WM is crucial when trying to understand the role of cognitive deficits, given that it explains the association between schizophrenia and cognitive deficits (directly and via processing speed).

Personalized dMRI Harmonization on Cortical Surface

The inter-site variability of diffusion magnetic resonance imaging (dMRI) hinders the aggregation of dMRI data from multiple centers. This necessitates dMRI harmonization for removing non-biological site-effects. Recently, the emergence of high-resolution dMRI data across various connectome imaging studies allows the large-scale analysis of cortical micro-structure. Existing harmonization methods, however, perform poorly in the harmonization of dMRI data in cortical areas because they rely on image registration methods to factor out anatomical variations, which have known difficulty in aligning cortical folding patterns. To overcome this fundamental challenge in dMRI harmonization, we propose a framework of personalized dMRI harmonization on the cortical surface to improve the dMRI harmonization of gray matter by adaptively estimating the inter-site harmonization mappings. In our experiments, we demonstrate the effectiveness of the proposed method by applying it to harmonize dMRI across the Human Connectome Project (HCP) and the Lifespan Human Connectome Projects in Development (HCPD) studies and achieved much better performance in comparison with conventional methods based on image registration.

Advanced diffusion MRI and image texture analysis detect widespread brain structural differences between relapsing-remitting and secondary progressive multiple sclerosis

Introduction

Disease development in multiple sclerosis (MS) causes dramatic structural changes, but the exact changing patterns are unclear. Our objective is to investigate the differences in brain structure locally and spatially between relapsing-remitting MS (RRMS) and its advanced form, secondary progressive MS (SPMS), through advanced analysis of diffusion magnetic resonance imaging (MRI) and image texture.

Methods

A total of 20 patients with RRMS and nine patients with SPMS from two datasets underwent 3T anatomical and diffusion tensor imaging (DTI). The DTI was harmonized, augmented, and then modeled, which generated six voxel- and sub-voxel-scale measures. Texture analysis focused on T2 and FLAIR MRI, which produced two phase-based measures, namely, phase congruency and weighted mean phase. Data analysis was 3-fold, i.e., histogram analysis of whole-brain normal appearing white matter (NAWM); region of interest (ROI) analysis of NAWM and lesions within three critical white matter tracts, namely, corpus callosum, corticospinal tract, and optic radiation; and along-tract statistics. Furthermore, by calculating the z-score of core-rim pathology within lesions based on diffusion measures, we developed a novel method to define chronic active lesions and compared them between cohorts.

Results

Histogram features from diffusion and all but one texture measure differentiated between RRMS and SPMS. Within-tract ROI analysis detected cohort differences in both NAWM and lesions of the corpus callosum body in three measures of neurite orientation and anisotropy. Along-tract statistics detected cohort differences from multiple measures, particularly lesion extent, which increased significantly in SPMS in posterior corpus callosum and optic radiations. The number of chronic active lesions were also significantly higher (by 5-20% over z-scores 0.5 and 1.0) in SPMS than RRMS based on diffusion anisotropy, neurite content, and diameter.

Conclusion

Advanced diffusion MRI and texture analysis may be promising approaches for thorough understanding of brain structural changes from RRMS to SPMS, thereby providing new insight into disease development mechanisms in MS.

Reduced asymmetry of the hand knob area and decreased sensorimotor u-fiber connectivity in middle-aged adults with autism

Individuals with autism spectrum disorder (ASD) frequently present with impairments in motor skills (e.g., limb coordination, handwriting and balance), which are observed across the lifespan but remain largely untreated. Many adults with ASD may thus experience adverse motor outcomes in aging, when physical decline naturally occurs. The 'hand knob' of the sensorimotor cortex is an area that is critical for motor control of the fingers and hands. However, this region has received little attention in ASD research, especially in adults after midlife. The hand knob area of the precentral (PrC_hand) and postcentral (PoC_hand) gyri was semi-manually delineated in 49 right-handed adults (25 ASD, 24 typical comparison [TC] participants, aged 41-70 years). Using multimodal (T1-weighted, diffusion-weighted, and resting-state functional) MRI, we examined the morphology, ipsilateral connectivity and laterality of these regions. We also explored correlations between hand knob measures with motor skills and autism symptoms, and between structural and functional connectivity measures. Bayesian analyses indicated moderate evidence of group effects with greater right PrC_hand volume and reduced leftward laterality of PrC_hand and PoC_hand volume in the ASD relative to TC group. Furthermore, the right PoC-PrC_hand u-fibers showed increased mean diffusivity in the ASD group. In the ASD group, right u-fiber volume positively correlated with corresponding functional connectivity but did not survive multiple comparisons correction. Correlations of hand knob measures and behavior were observed in the ASD group but did not survive multiple comparisons correction. Our findings suggest that morphological laterality and u-fiber connectivity of the sensorimotor network, putatively involved in hand motor/premotor function, may be diminished in middle-aged adults with ASD, perhaps rendering them more vulnerable to motor decline in old age. The altered morphology may relate to atypical functional motor asymmetries found in ASD earlier in life, possibly reflecting altered functional asymmetries over time.

Cross-site harmonization of multi-shell diffusion MRI measures based on rotational invariant spherical harmonics (RISH)

Quantification methods based on the acquisition of diffusion magnetic resonance imaging (dMRI) with multiple diffusion weightings (e.g., multi-shell) are becoming increasingly applied to study the in-vivo brain. Compared to single-shell data for diffusion tensor imaging (DTI), multi-shell data allows to apply more complex models such as diffusion kurtosis imaging (DKI), which attempts to capture both diffusion hindrance and restriction effects, or biophysical models such as NODDI, which attempt to increase specificity by separating biophysical components. Because of the strong dependence of the dMRI signal on the measurement hardware, DKI and NODDI metrics show scanner and site differences, much like other dMRI metrics. These effects limit the implementation of multi-shell approaches in multicenter studies, which are needed to collect large sample sizes for robust analyses. Recently, a post-processing technique based on rotation invariant spherical harmonics (RISH) features was introduced to mitigate cross-scanner differences in DTI metrics. Unlike statistical harmonization methods, which require repeated application to every dMRI metric of choice, RISH harmonization is applied once on the raw data, and can be followed by any analysis. RISH features harmonization has been tested on DTI features but not its generalizability to harmonize multi-shell dMRI. In this work, we investigated whether performing the RISH features harmonization of multi-shell dMRI data removes cross-site differences in DKI and NODDI metrics while retaining longitudinal effects. To this end, 46 subjects underwent a longitudinal (up to 3 time points) two-shell dMRI protocol at 3 imaging sites. DKI and NODDI metrics were derived before and after harmonization and compared both at the whole brain level and at the voxel level. Then, the harmonization effects on cross-sectional and on longitudinal group differences were evaluated. RISH features averaged for each of the 3 sites exhibited prominent between-site differences in the frontal and posterior part of the brain. Statistically significant differences in fractional anisotropy, mean diffusivity and mean kurtosis were observed both at the whole brain and voxel level between all the acquisition sites before harmonization, but not after. The RISH method also proved effective to harmonize NODDI metrics, particularly in white matter. The RISH based harmonization maintained the magnitude and variance of longitudinal changes as compared to the non-harmonized data of all considered metrics. In conclusion, the application of RISH feature based harmonization to multi-shell dMRI data can be used to remove cross-site differences in DKI metrics and NODDI analyses, while retaining inherent relations between longitudinal acquisitions.

Data harmonisation for information fusion in digital healthcare: A state-of-the-art systematic review, meta-analysis and future research directions

Removing the bias and variance of multicentre data has always been a challenge in large scale digital healthcare studies, which requires the ability to integrate clinical features extracted from data acquired by different scanners and protocols to improve stability and robustness. Previous studies have described various computational approaches to fuse single modality multicentre datasets. However, these surveys rarely focused on evaluation metrics and lacked a checklist for computational data harmonisation studies. In this systematic review, we summarise the computational data harmonisation approaches for multi-modality data in the digital healthcare field, including harmonisation strategies and evaluation metrics based on different theories. In addition, a comprehensive checklist that summarises common practices for data harmonisation studies is proposed to guide researchers to report their research findings more effectively. Last but not least, flowcharts presenting possible ways for methodology and metric selection are proposed and the limitations of different methods have been surveyed for future research.

Sensitivity of Diffusion MRI to White Matter Pathology: Influence of Diffusion Protocol, Magnetic Field Strength, and Processing Pipeline in Systemic Lupus Erythematosus

There are many ways to acquire and process diffusion MRI (dMRI) data for group studies, but it is unknown which maximizes the sensitivity to white matter (WM) pathology. Inspired by this question, we analyzed data acquired for diffusion tensor imaging (DTI) and diffusion kurtosis imaging (DKI) at 3T (3T-DTI and 3T-DKI) and DTI at 7T in patients with systemic lupus erythematosus (SLE) and healthy controls (HC). Parameter estimates in 72 WM tracts were obtained using TractSeg. The impact on the sensitivity to WM pathology was evaluated for the diffusion protocol, the magnetic field strength, and the processing pipeline. Sensitivity was quantified in terms of Cohen's d for group comparison. Results showed that the choice of diffusion protocol had the largest impact on the effect size. The effect size in fractional anisotropy (FA) across all WM tracts was 0.26 higher when derived by DTI than by DKI and 0.20 higher in 3T compared with 7T. The difference due to the diffusion protocol was larger than the difference due to magnetic field strength for the majority of diffusion parameters. In contrast, the difference between including or excluding different processing steps was near negligible, except for the correction of distortions from eddy currents and motion which had a clearly positive impact. For example, effect sizes increased on average by 0.07 by including motion and eddy correction for FA derived from 3T-DTI. Effect sizes were slightly reduced by the incorporation of denoising and Gibbs-ringing removal (on average by 0.011 and 0.005, respectively). Smoothing prior to diffusion model fitting generally reduced effect sizes. In summary, 3T-DTI in combination with eddy current and motion correction yielded the highest sensitivity to WM pathology in patients with SLE. However, our results also indicated that the 3T-DKI and 7T-DTI protocols used here may be adjusted to increase effect sizes.

From the Group to the Individual in Schizophrenia Spectrum Disorders: Biomarkers of Social Cognitive Impairments and Therapeutic Translation

People with schizophrenia spectrum disorders (SSDs) often experience persistent social cognitive impairments, associated with poor functional outcome. There are currently no approved treatment options for these debilitating symptoms, highlighting the need for novel therapeutic strategies. Work to date has elucidated differential social processes and underlying neural circuitry affected in SSDs, which may be amenable to modulation using neurostimulation. Further, advances in functional connectivity mapping and electric field modeling may be used to identify individualized treatment targets to maximize the impact of brain stimulation on social cognitive networks. Here, we review literature supporting a roadmap for translating functional connectivity biomarker discovery to individualized treatment development for social cognitive impairments in SSDs. First, we outline the relevance of social cognitive impairments in SSDs. We review machine learning approaches for dimensional brain-behavior biomarker discovery, emphasizing the importance of individual differences. We synthesize research showing that brain stimulation techniques, such as repetitive transcranial magnetic stimulation, can be used to target relevant networks. Further, functional connectivity-based individualized targeting may enhance treatment response. We then outline recent approaches to account for neuroanatomical variability and optimize coil positioning to individually maximize target engagement. Overall, the synthesized literature provides support for the utility and feasibility of this translational approach to precision treatment. The proposed roadmap to translate biomarkers of social cognitive impairments to individualized treatment is currently under evaluation in precision-guided trials. Such a translational approach may also be applicable across conditions and generalizable for the development of individualized neurostimulation targeting other behavioral deficits.

The SACT Template: A Human Brain Diffusion Tensor Template for School-age Children

School-age children are in a specific development stage corresponding to juvenility, when the white matter of the brain experiences ongoing maturation. Diffusion-weighted magnetic resonance imaging (DWI), especially diffusion tensor imaging (DTI), is extensively used to characterize the maturation by assessing white matter properties in vivo. In the analysis of DWI data, spatial normalization is crucial for conducting inter-subject analyses or linking the individual space with the reference space. Using tensor-based registration with an appropriate diffusion tensor template presents high accuracy regarding spatial normalization. However, there is a lack of a standardized diffusion tensor template dedicated to school-age children with ongoing brain development. Here, we established the school-age children diffusion tensor (SACT) template by optimizing tensor reorientation on high-quality DTI data from a large sample of cognitively normal participants aged 6-12 years. With an age-balanced design, the SACT template represented the entire age range well by showing high similarity to the age-specific templates. Compared with the tensor template of adults, the SACT template revealed significantly higher spatial normalization accuracy and inter-subject coherence upon evaluation of subjects in two different datasets of school-age children. A practical application regarding the age associations with the normalized DTI-derived data was conducted to further compare the SACT template and the adult template. Although similar spatial patterns were found, the SACT template showed significant effects on the distributions of the statistical results, which may be related to the performance of spatial normalization. Looking forward, the SACT template could contribute to future studies of white matter development in both healthy and clinical populations. The SACT template is publicly available now ( https://figshare.com/articles/dataset/SACT_template/14071283 ).

White matter impairments in patients with schizophrenia: A multisite diffusion MRI study

There is a lack of convincing and replicative findings regarding white matter abnormalities in schizophrenia. Several multisite diffusion magnetic resonance imaging (dMRI) studies have been conducted to increase statistical power and reveal subtle white matter changes. Data pooling methods are crucial in joint analysis to compensate for the use of different scanners and image acquisition parameters. A harmonization method using raw dMRI data was developed to overcome the limited generalizability of previous data pooling methods. We obtained dMRI data of 242 healthy controls and 190 patients with schizophrenia from four different study sites. After applying the harmonization method to the raw dMRI data, a two-tensor whole-brain tractography was performed, and diffusion measures were compared between the two groups. The correlation of fractional anisotropy (FA) with the positive and negative symptoms was evaluated, and the interaction effect of diagnosis-by-age, age-squared, and sex was examined. The following white matter tracts showed significant group differences in the FA: the right superior longitudinal fascicle (SLF), the left-to-right lateral orbitofrontal commissural tract, pars orbitalis (pOr-pOr) commissural tract, and pars triangularis (pTr-pTr) commissural tract. The FA of the right SLF and pTr-pTr commissural tract were significantly associated with the Positive and Negative Syndrome Scale (PANSS) positive and negative scores. No significant interaction effect was observed. These findings add to the evidence on structural brain abnormalities in schizophrenia and can aid in obtaining a better understanding of the biological foundations of schizophrenia.

A multi-scanner neuroimaging data harmonization using RAVEL and ComBat

Modern neuroimaging studies frequently combine data collected from multiple scanners and experimental conditions. Such data often contain substantial technical variability associated with image intensity scale (image intensity scales are not the same in different images) and scanner effects (images obtained from different scanners contain substantial technical biases). Here we evaluate and compare results of data analysis methods without any data transformation (RAW), with intensity normalization using RAVEL, with regional harmonization methods using ComBat, and a combination of RAVEL and ComBat. Methods are evaluated on a unique sample of 16 study participants who were scanned on both 1.5T and 3T scanners a few months apart. Neuroradiological evaluation was conducted for 7 different regions of interest (ROI's) pertinent to Alzheimer's disease (AD). Cortical measures and results indicate that: (1) RAVEL substantially improved the reproducibility of image intensities; (2) ComBat is preferred over RAVEL and the RAVEL-ComBat combination in terms of regional level harmonization due to more consistent harmonization across subjects and image-derived measures; (3) RAVEL and ComBat substantially reduced bias compared to analysis of RAW images, but RAVEL also resulted in larger variance; and (4) the larger root mean square deviation (RMSD) of RAVEL compared to ComBat is due mainly to its larger variance.