Search

Diffusion imaging is a critical component in the pursuit of developing a better understanding of the human brain. Recent technical advances promise enabling the advancement in the quality of data that can be obtained. In this review the context for different approaches relative to the Human Connectome Project are compared. Significant new gains are anticipated from the use of high-performance head gradients. These gains can be particularly large when the high-performance gradients are employed together with ultrahigh magnetic fields. Transmit array designs are critical in realizing high accelerations in diffusion-weighted (d)MRI acquisitions, while maintaining large field of view (FOV) coverage, and several techniques for optimal signal-encoding are now available. Reconstruction and processing pipelines that precisely disentangle the acquired neuroanatomical information are established and provide the foundation for the application of deep learning in the advancement of dMRI for complex tissues. Level of Evidence: 3 Technical Efficacy Stage: Stage 3.

In a companion paper by Cohen-Adad et al. we introduce the spine generic quantitative MRI protocol that provides valuable metrics for assessing spinal cord macrostructural and microstructural integrity. This protocol was used to acquire a single subject dataset across 19 centers and a multi-subject dataset across 42 centers (for a total of 260 participants), spanning the three main MRI manufacturers: GE, Philips and Siemens. Both datasets are publicly available via git-annex. Data were analysed using the Spinal Cord Toolbox to produce normative values as well as inter/intra-site and inter/intra-manufacturer statistics. Reproducibility for the spine generic protocol was high across sites and manufacturers, with an average inter-site coefficient of variation of less than 5% for all the metrics. Full documentation and results can be found at https://spine-generic.rtfd.io/. The datasets and analysis pipeline will help pave the way towards accessible and reproducible quantitative MRI in the spinal cord., Measurement(s)spinal cordTechnology Type(s)magnetic resonance imagingFactor Type(s)manufacturer • siteSample Characteristic - OrganismHomo sapiensSample Characteristic - LocationCanada • Switzerland • Australia • United States of America • United Kingdom • Germany • French Republic • Czech Republic • Italy • Japan • Kingdom of Spain • China Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.14052269

Mild traumatic brain injury (mTBI) is a significant cause of disability, especially when symptoms become chronic. This chronicity is often linked to oculomotor dysfunction (OMD). To our knowledge, this is the first prospective study to localize aberrations in brain function between mTBI cohorts, by comparing patients with mTBI with OMD with an mTBI control group without OMD, using task and resting-state functional magnetic resonance imaging (fMRI). Ten subjects with mTBI who had OMD (OMD group) were compared with nine subjects with mTBI who had no findings of OMD (control group). These groups were determined by a developmental optometrist using objective testing for OMD. The (convergence) task fMRI data demonstrated significantly decreased brain activity, measured as decreases in the blood oxygen level dependent (BOLD) signal, in the OMD group compared with the control group in three brain regions: the left posterior lingual gyrus, the bilateral anterior lingual gyrus and cuneus, and the parahippocampal gyrus. When doing a seed-based resting state fMRI analysis in the lingual/parahippocampal region, a large cluster covering the left middle frontal gyrus and the dorsolateral pre-frontal cortex (Brodmann areas 9 and 10), with decreased functional correlation in the OMD group, was identified. Together these observations provide evidence for neural networks of interactions involving the control of eye movement for visual processing, reading comprehension, spatial localization and navigation, and spatial working memory that appear to be decreased in mTBI patients with OMD compared with mTBI patients without OMD. The clinical symptomatology associated with post-traumatic OMD correlates well with these MRI findings.

BACKGROUND AND PURPOSE Non-myelopathic degenerative cervical spinal cord compression (NMDC) frequently occurs throughout aging and may progress to potentially irreversible degenerative cervical myelopathy (DCM). Whereas standard clinical magnetic resonance imaging (MRI) and electrophysiological measures assess compression severity and neurological dysfunction, respectively, underlying microstructural deficits still have to be established in NMDC and DCM patients. The study aims to establish tract-specific diffusion MRI markers of electrophysiological deficits to predict the progression of asymptomatic NMDC to symptomatic DCM. METHODS High-resolution 3 T diffusion MRI was acquired for 103 NMDC and 21 DCM patients compared to 60 healthy controls to reveal diffusion alterations and relationships between tract-specific diffusion metrics and corresponding electrophysiological measures and compression severity. Relationship between the degree of DCM disability, assessed by the modified Japanese Orthopaedic Association scale, and tract-specific microstructural changes in DCM patients was also explored. RESULTS The study identified diffusion-derived abnormalities in the gray matter, dorsal and lateral tracts congruent with trans-synaptic degeneration and demyelination in chronic degenerative spinal cord compression with more profound alterations in DCM than NMDC. Diffusion metrics were affected in the C3-6 area as well as above the compression level at C3 with more profound rostral deficits in DCM than NMDC. Alterations in lateral motor and dorsal sensory tracts correlated with motor and sensory evoked potentials, respectively, whereas electromyography outcomes corresponded with gray matter microstructure. DCM disability corresponded with microstructure alteration in lateral columns. CONCLUSIONS Outcomes imply the necessity of high-resolution tract-specific diffusion MRI for monitoring degenerative spinal pathology in longitudinal studies.

Spinocerebellar ataxia type 1 is a progressive neurodegenerative, movement disorder. With potential therapies on the horizon, it is critical to identify biomarkers that (i) differentiate between unaffected and spinocerebellar ataxia Type 1-affected individuals; (ii) track disease progression; and (iii) are directly related to clinical changes of the patient. Magnetic resonance imaging of volumetric changes in the brain may be a suitable source of biomarkers for spinocerebellar ataxia Type 1. In a previous report on a longitudinal study of patients with spinocerebellar ataxia Type 1, we evaluated the volume and magnetic resonance spectroscopy measures of the cerebellum and pons, showing pontine volume and pontine N-acetylaspartate-to-myo-inositol ratio were sensitive to change over time. As a follow-up, the current study conducts a whole brain exploration of volumetric MRI measures with the aim to identify biomarkers for spinocerebellar ataxia Type 1 progression. We adapted a joint label fusion approach using multiple, automatically generated, morphologically matched atlases to label brain regions including cerebellar sub-regions. We adjusted regional volumes by total intracranial volume allowing for linear and power-law relationships. We then utilized Bonferroni corrected linear mixed effects models to (i) determine group differences in regional brain volume and (ii) identify change within affected patients only. We then evaluated the rate of change within each brain region to identify areas that changed most rapidly. Lastly, we used a penalized, linear mixed effects model to determine the strongest brain predictors of motor outcomes. Decrease in pontine volume and accelerating decrease in putamen volume: (i) reliably differentiated spinocerebellar ataxia Type 1-affected and -unaffected individuals; (ii) were observable in affected individuals without referencing an unaffected comparison group; (iii) were detectable within ∼6–9 months; and (iv) were associated with increased disease burden. In conclusion, volumetric change in the pons and putamen may provide powerful biomarkers to track disease progression in spinocerebellar ataxia Type 1. The methods employed here are readily translatable to current clinical settings, providing a framework for study and usage of volumetric neuroimaging biomarkers for clinical trials., There is an urgent need for biomarkers that track neurodegenerative disease progression. Pontine and putamen volume decreases are specific to affected individuals, track disease progression within 6–9 months and are directly related to clinical changes in individuals affected by spinocerebellar ataxia Type 1, making them potential biomarkers for clinical trials., Graphical Abstract

The robustness of brain structural networks, estimated from diffusion MRI data, may be relevant to cognition. We investigate whether measures of network robustness, such as Ollivier-Ricci curvature, can explain cognitive impairment in multiple sclerosis (MS). We assessed whether local (i.e., cortical area) and/or global (i.e., whole brain) robustness, differs between cognitively impaired (MSCI) and non-impaired (MSNI) MS patients. Fifty patients, with Expanded Disability Status Scale mean (m): 3.2, disease duration m: 12 years, and age m: 40 years, were enrolled. Cognitive impairment scores were estimated from the Minimal Assessment of Cognitive Function in Multiple Sclerosis. Images were obtained in a 3T MRI using a diffusion protocol with a 2 min acquisition time. Brain structural networks were created using 333 cortical areas. Local and global robustness was estimated for each individual, and comparisons were performed between MSCI and MSNI patients. 31 MSCI and 10 MSNI patients were included in the analyses. Brain structural network robustness and centrality showed significant correlations with cognitive impairment. Measures of network robustness and centrality identified specific cortical areas relevant to MS-related cognitive impairment. These measures can be obtained on clinical scanners and are succinct yet accurate potential biomarkers of cognitive impairment.

Diffusion magnetic resonance imaging (dMRI) proved promising in patients with non-myelopathic degenerative cervical cord compression (NMDCCC), i.e., without clinically manifested myelopathy. Aim of the study is to present a fast multi-shell HARDI-ZOOMit dMRI protocol and validate its usability to detect microstructural myelopathy in NMDCCC patients. In 7 young healthy volunteers, 13 age-comparable healthy controls, 18 patients with mild NMDCCC and 15 patients with severe NMDCCC, the protocol provided higher signal-to-noise ratio, enhanced visualization of white/gray matter structures in microstructural maps, improved dMRI metric reproducibility, preserved sensitivity (SE = 87.88%) and increased specificity (SP = 92.31%) of control-patient group differences when compared to DTI-RESOLVE protocol (SE = 87.88%, SP = 76.92%). Of the 56 tested microstructural parameters, HARDI-ZOOMit yielded significant patient-control differences in 19 parameters, whereas in DTI-RESOLVE data, differences were observed in 10 parameters, with mostly lower robustness. Novel marker the white-gray matter diffusivity gradient demonstrated the highest separation. HARDI-ZOOMit protocol detected larger number of crossing fibers (5–15% of voxels) with physiologically plausible orientations than DTI-RESOLVE protocol (0–8% of voxels). Crossings were detected in areas of dorsal horns and anterior white commissure. HARDI-ZOOMit protocol proved to be a sensitive and practical tool for clinical quantitative spinal cord imaging.

Development of imaging biomarkers for rare neurodegenerative diseases such as spinocerebellar ataxia (SCA) is important to non-invasively track progression of disease pathology and monitor response to interventions. Diffusion MRI (dMRI) has been shown to identify cross-sectional degeneration of white matter (WM) microstructure and connectivity between healthy controls and patients with SCAs, using various analysis methods. In this paper, we present dMRI data in SCAs type 1, 2, 3, and 6 and matched controls, including longitudinal acquisitions at 12–24-month intervals in a subset of the cohort, with up to 5 visits. The SCA1 cohort also contained 3 premanifest patients at baseline, with 2 showing ataxia symptoms at the time of the follow-up scans. We focused on two aspects: first, multimodal evaluation of the dMRI data in a cross-sectional approach, and second, longitudinal trends in dMRI data in SCAs. Three different pipelines were used to perform cross-sectional analyses in WM: region of interest (ROI), tract-based spatial statistics (TBSS), and fixel-based analysis (FBA). We further analyzed longitudinal changes in dMRI metrics throughout the brain using ROI-based analysis. Both ROI and TBSS analyses identified higher mean (MD), axial (AD), and radial (RD) diffusivity and lower fractional anisotropy (FA) in the cerebellum for all SCAs compared to controls, as well as some cerebral alterations in SCA1, 2, and 3. FBA showed lower fiber density (FD) and fiber crossing (FC) regions similar to those identified by ROI and TBSS analyses. FBA also highlighted corticospinal tract (CST) abnormalities, which was not detected by the other two pipelines. Longitudinal ROI-based analysis showed significant increase in AD in the middle cerebellar peduncle (MCP) for patients with SCA1, suggesting that the MCP may be a good candidate region to monitor disease progression. The patient who remained symptom-free throughout the study displayed no microstructural abnormalities. On the other hand, the two patients who were at the premanifest stage at baseline, and showed ataxia symptoms in their follow-up visits, displayed AD values in the MCP that were already in the range of symptomatic patients with SCA1 at their baseline visit, demonstrating that microstructural abnormalities are detectable prior to the onset of ataxia.

The success of deep brain stimulation (DBS) surgeries for the treatment of movement disorders relies on the accurate placement of an electrode within the motor portion of subcortical brain targets. However, the high number of electrodes requiring relocation indicates that today's methods do not ensure sufficient accuracy for all patients. Here, with the goal of aiding DBS targeting, we use 7 Tesla (T) MRI data to identify the functional territories and parcellate the globus pallidus pars interna (GPi) into motor, associative and limbic regions in individual subjects. 7 T MRI scans were performed in seventeen patients (prior to DBS surgery) and one healthy control. Tractography-based parcellation of each patient's GPi was performed. The cortex was divided into four masks representing motor, limbic, associative and “other” regions. Given that no direct connections between the GPi and the cortex have been shown to exist, the parcellation was carried out in two steps: 1) The thalamus was parcellated based on the cortical targets, 2) The GPi was parcellated using the thalamus parcels derived from step 1. Reproducibility, via repeated scans of a healthy subject, and validity of the findings, using different anatomical pathways for parcellation, were assessed. Lastly, post-operative imaging data was used to validate and determine the clinical relevance of the parcellation. The organization of the functional territories of the GPi observed in our individual patient population agrees with that previously reported in the literature: the motor territory was located posterolaterally, followed anteriorly by the associative region, and further antero-ventrally by the limbic territory. While this organizational pattern was observed across patients, there was considerable variability among patients. The organization of the functional territories of the GPi was remarkably reproducible in intra-subject scans. Furthermore, the organizational pattern was observed consistently by performing the parcellation of the GPi via the thalamus and via a different pathway, going through the striatum. Finally, the active therapeutic contact of the DBS electrode, identified with a combination of post-operative imaging and post-surgery DBS programming, overlapped with the high-probability “motor” region of the GPi as defined by imaging-based methods. The consistency, validity, and clinical relevance of our findings have the potential for improving DBS targeting, by increasing patient-specific knowledge of subregions of the GPi to be targeted or avoided, at the stage of surgical planning, and later, at the stage when stimulation is adjusted.

Diffusion MRI of the cortical gray matter is challenging because the micro-environment probed by water molecules is much more complex than within the white matter. High spatial and angular resolutions are therefore necessary to uncover anisotropic diffusion patterns and laminar structures, which provide complementary (e.g. to anatomical and functional MRI) microstructural information about the cortex architectonic. Several ex-vivo and in-vivo MRI studies have recently addressed this question, however predominantly with an emphasis on specific cortical areas. There is currently no whole brain in-vivo data leveraging multi-shell diffusion MRI acquisition at high spatial resolution, and depth dependent analysis, to characterize the complex organization of cortical fibers. Here, we present unique in-vivo human 7T diffusion MRI data, and a dedicated cortical depth dependent analysis pipeline. We leverage the high spatial (1.05 mm isotropic) and angular (198 diffusion gradient directions) resolution of this whole brain dataset to improve cortical fiber orientations mapping, and study neurites (axons and/or dendrites) trajectories across cortical depths. Tangential fibers in superficial cortical depths and crossing fiber configurations in deep cortical depths are identified. Fibers gradually inserting into the gyral walls are visualized, which contributes to mitigating the gyral bias effect. Quantitative radiality maps and histograms in individual subjects and cortex-based aligned datasets further support our results.

PURPOSE: To develop a fast and automated volume-of-interest (VOI) prescription pipeline (AutoVOI) for single-voxel MR spectroscopy (MRS) that removes the need for manual VOI placement, allows flexible VOI planning in any brain region, and enables high inter- and intra-subject consistency of VOI prescription. METHODS: AutoVOI was designed to transfer pre-defined VOIs from an atlas to the 3D anatomical data of the subject during the scan. The AutoVOI pipeline was optimized for consistency in VOI placement (precision), enhanced coverage of the targeted tissue (accuracy) and fast computation speed. The tool was evaluated against manual VOI placement using existing T(1)-weighted datasets and corresponding VOI prescriptions. Finally, it was implemented on two scanner platforms to acquire MRS data from clinically-relevant VOIs that span the cerebrum, cerebellum and the brainstem. RESULTS: The AutoVOI pipeline includes skull stripping, non-linear registration of the atlas to the subject’s brain, and computation of the VOI coordinates and angulations using a minimum oriented bounding box algorithm. When compared against manual prescription, AutoVOI showed higher intra- and inter-subject spatial consistency, as quantified by generalized Dice coefficients (GDC), lower intra- and inter-subject variability in tissue composition (gray matter, white matter and cerebrospinal fluid) and higher or equal accuracy, as quantified by GDC of prescribed VOI with targeted tissues. High quality spectra were obtained on Siemens and Philips 3T systems from 6 automatically prescribed VOIs by the tool. CONCLUSION: Robust automatic VOI prescription is feasible and can help facilitate clinical adoption of MRS by avoiding operator dependence of manual selection.

PURPOSE: Investigating the utility of RF parallel transmission (pTx) for Human Connectome Project (HCP)-style whole-brain diffusion MRI (dMRI) data at 7 Tesla (7T). METHODS: Healthy subjects were scanned in pTx and single-transmit (1Tx) modes. Multiband (MB), single-spoke pTx pulses were designed to image sagittal slices. HCP-style dMRI data (i.e., 1.05-mm resolutions, MB2, b-values=1000/2000 s/mm(2), 286 images and 40-minute scan) and data with higher accelerations (MB3 and MB4) were acquired with pTx. RESULTS: pTx significantly improved flip-angle detected signal uniformity across the brain, yielding ~19% increase in temporal signal-to-noise ratio (tSNR) averaged over the brain relative to 1Tx. This allowed significantly enhanced estimation of multiple fiber orientations (with ~21% decrease in dispersion) in HCP-style 7T dMRI datasets. Additionally, pTx pulses achieved substantially lower power deposition, permitting higher accelerations, enabling collection of the same data in 2/3 and 1/2 the scan time or of more data in the same scan time. CONCLUSION: pTx provides a solution to two major limitations for slice-accelerated high-resolution whole-brain dMRI at 7T; it improves flip-angle uniformity, and enables higher slice acceleration relative to current state-of-the-art. As such, pTx provides significant advantages for rapid acquisition of high-quality, high-resolution truly whole brain dMRI data.

We present a sparse Bayesian unmixing algorithm BusineX: Bayesian Unmixing for Sparse Inference-based Estimation of Fiber Crossings (X), for estimation of white matter fiber parameters from compressed (under-sampled) diffusion MRI (dMRI) data. BusineX combines compressive sensing with linear unmixing and introduces sparsity to the previously proposed multiresolution data fusion algorithm RubiX, resulting in a method for improved reconstruction, especially from data with lower number of diffusion gradients. We formulate the estimation of fiber parameters as a sparse signal recovery problem and propose a linear unmixing framework with sparse Bayesian learning for the recovery of sparse signals, the fiber orientations and volume fractions. The data is modeled using a parametric spherical deconvolution approach and represented using a dictionary created with the exponential decay components along different possible diffusion directions. Volume fractions of fibers along these directions define the dictionary weights. The proposed sparse inference, which is based on the dictionary representation, considers the sparsity of fiber populations and exploits the spatial redundancy in data representation, thereby facilitating inference from under-sampled q-space. The algorithm improves parameter estimation from dMRI through data-dependent local learning of hyperparameters, at each voxel and for each possible fiber orientation, that moderate the strength of priors governing the parameter variances. Experimental results on synthetic and in-vivo data show improved accuracy with a lower uncertainty in fiber parameter estimates. BusineX resolves a higher number of second and third fiber crossings. For under-sampled data, the algorithm is also shown to produce more reliable estimates.

BACKGROUND: Spinocerebellar ataxia type 1 (SCA1) causes progressive degeneration of the cerebellum and brainstem. Volumetric magnetic resonance imaging (MRI) was shown to be more sensitive to disease progression than the most sensitive clinical measure, the Scale for the Assessment and Rating of Ataxia (SARA), in longitudinal studies, and magnetic resonance spectroscopy (MRS) was shown to detect neurochemical abnormalities with high sensitivity cross‐sectionally in SCA1. OBJECTIVES: The objectives of this study were to compare the sensitivities to change of volumetric MRI, MRS, and SARA in a 3‐year longitudinal study in SCA1. METHODS: A total of 16 early‐to‐moderate stage patients with SCA1 (SARA 0‐14) and 21 matched healthy participants were scanned up to 3 times with 1.5‐year intervals. Ataxia severity was assessed with SARA. T(1)‐weighted images and magnetic resonance spectra from the cerebellar vermis, cerebellar white matter, and pons were acquired at 3T. RESULTS: The pontine total N‐acetylaspartate‐to‐myo‐inositol ratio was the most sensitive MRS measure to change (−3.9 ± 4.6%/yr in SCA1 vs. −0.3 ± 3.5%/yr in controls; P

International audience; Computational methods are crucial for the analysis of diffusion magnetic resonance imaging (MRI) of the brain. Computational diffusion MRI can provide rich information at many size scales, including local microstructure measures such as diffusion anisotropies or apparent axon diameters, whole‐brain connectivity information that describes the brain's wiring diagram and population‐based studies in health and disease. Many of the diffusion MRI analyses performed today were not possible five, ten or twenty years ago, due to the requirements for large amounts of computer memory or processor time. In addition, mathematical frameworks had to be developed or adapted from other fields to create new ways to analyze diffusion MRI data. The purpose of this review is to highlight recent computational and statistical advances in diffusion MRI and to put these advances into context by comparison with the more traditional computational methods that are in popular clinical and scientific use. We aim to provide a high‐level overview of interest to diffusion MRI researchers, with a more in‐depth treatment to illustrate selected computational advances.

Adolescence is a period of rapid brain maturation and a critical period for the onset of Major Depressive Disorder (MDD) that usually leads to serious outcomes such as suicide. Although changes in anatomical connectivity in MDD have been reported, changes in network topology for MDD remain unclear. Additionally, whether the changes are the same for adolescent MDD and adult MDD remains unclear as well. This paper explores anatomical features including: a) anatomical connectivity defined by diffusion tensor imaging measurements between a pair of brain regions, and b) topological measurements from anatomical networks, and apply machine learning approaches to identify responsive biomarkers distinguishing MDD patients from healthy subjects. In addition to statistical tests, univariate classifiers are designed to evaluate the discriminating power of features. Furthermore, multivariate classifiers are trained for distinguishing healthy subjects from MDD patients. The best classifier achieves an accuracy of 76.56%, 81.08% sensitivity, 70.37% specificity and 78.95% precision for 64 subjects (37 MDD and 27 matched healthy control). The selected features include: 1) betweenness centrality of the right lingual gyrus of the ADC network at 12% sparsity, 2) participation coefficient of the right pars opercularis of the AD network at 16% sparsity, 3) participation coefficient of the left insular cortex of the MD network at 21% sparsity, and 4) participation coefficient of the right lateral orbitofrontal cortex in the ADC network at 10% sparsity. These features reflect changes in the topological structure of the brain anatomical network in MDD.

The adolescent major depressive disorder is one of the top 10 debilitating psychiatric illnesses and the effectiveness of current treatment methods are constrained by the limited understanding of biological causes. In this paper, we use diffusion tensor imaging to explore changes in anatomical connectivity between the MDD group (n=37) and control group (n=27). Furthermore, along-track analysis is performed to identify locations of alterations along connections with significant connectivity change. For the connection between the hippocampus and the insular cortex in the right hemisphere, decreased connectivity in axial diffusivity (AD), mean diffusivity (MD), and radial diffusivity (RD) was discovered. Additionally, the number of tracks for the same connection is increased for the MDD group. Moreover, for the connection between the parahippocampal gyrus and the insular cortex in the left hemisphere, increased connectivity is observed in fractional anisotropy (FA), AD, MD, and RD. Furthermore, the locations of significant alterations are identified to be between 65% to 100% from the insular cortex to the hippocampus in the right hemisphere and at the 80% location from the insular cortex to the parahippocampal gyrus in the left hemisphere. The significant and consistent white matter changes at the hippocampus end of the insula-hippocampus connection suggest potential correlations to the previously reported grey matter shrinkage and functional abnormalities.

Due to the high resistance (35%) to the current treatment methods in adolescent Major Depressive Disorder (MDD) and its tragic outcomes, the discovery of treatmentrelated responders is critical to developing effective treatments. In this paper, the permutation test is performed to identify statistically significant changes in anatomical characteristics during pairwise comparisons among the control group (n=27), treated MDD group (n=37), and untreated MDD group (n=15). The anatomical characteristics include: 1) anatomical connectivity defined using DTI metrics between a pair of brain regions, and 2) topological measurements of anatomical networks. With the Bonferroni correction for multiple-comparison, significant alterations in community structure and local topology were identified as the p-value < 5%, which include: 1) a reduced nodal centrality (degree and strength) on right hippocampus for treated compared to untreated group, 2) an elevated clustering coefficient and local efficiency on right lateral orbitofrontal cortex for untreated compared to the combination of control and treated groups, 3) an increased participation coefficient for untreated patients on left insula cortex in the meandiffusivity network compared to the combination of control and treated groups, and 4) a degraded module degree z-score on right caudate nucleus for all the patients compared to the control group. Two connections, hippocampus-insula in the right hemisphere and parahippocampal-insula in the left hemisphere, were found significantly altered in TR, AD, and FA due to MDD.

Adolescent Major Depressive Disorder (MDD) is a common and serious mental illness that could lead to tragic outcomes including chronic adult disability and suicide. In this paper, we explore anatomical features and apply machine learning approaches to identify responsive biomarkers distinguishing MDD patients from healthy subjects. The features of interest include metrics in two categories: a) anatomical connectivity defined by diffusion tensor imaging measurements between a pair of brain regions, and b) topological measurements from anatomical networks. A combination of p-value based filtering and minimum redundancy maximum relevance method is performed to select features for optimal classification accuracy. A leave-one-out cross-validation method is used for the classification performance evaluation. The proposed methodology achieves an improved accuracy of 78%, 90.39% sensitivity, and 79.66% precision for 79 subjects. The most distinguishing features are the betweenness centrality of the right lingual gyrus of the ADC network at 12% sparsity, the participation coefficient of the right lateral occipital sulcus of the ADC network at 22% sparsity, the participation coefficient of the right pars opercularis of the AD network at 16% sparsity, and the participation coefficient of the right lateral orbitofrontal cortex in the ADC network at 10% sparsity. Those network measures reflect the change of connectivity between the regions and their associated anatomical subnetworks.

Identification of the treatment-related responders for adolescent Major Depressive Disorder (MDD) is urgently needed to develop effective treatments. In this paper, machine learning based classifiers are used to reveal anatomical features as responders for distinguishing MDD patients who have received treatment from those who never received any treatment. The features are drawn from two sets of measurements: 1) anatomical connectivity defined by diffusion tensor imaging measurements between a pair of brain regions, and 2) topological measurements from anatomical networks. Feature selection was performed based on p-value and minimum redundancy maximum relevance (mRMR) method to achieve improved classification accuracy. The classification performance is evaluated with a leave-one-out cross-validation method using 37 treated and 15 untreated subjects. The proposed methodology achieves 73% accuracy, 100% specificity, and 100% precision for 52 subjects. The most distinguishing features are the strength of the right hippocampus of the mean diffusivity (MD) network at 18% density and of the track-count (TR) network, the participation coefficient of the left middle temporal gyrus of the radial diffusivity (RD) network at 20% density, the axial diffusivity (AD) connectivity between right middle temporal gyrus and right supramarginal gyrus, the betweenness centrality of the right hippocampus of the TR network at 11% density, the apparent diffusion coefficient (ADC) connectivity between the left pars opercularis and the left rostral anterior cingulate cortex, the clustering coefficient of the middle anterior corpus callosum of the TR network at 11% density, and the AD connectivity between the left pars opercularis and the left rostral anterior cingulate cortex.

Diffusion Magnetic Resonance Imaging (dMRI) is a widely-utilized method for assessment of microstructural properties in the central nervous system i.e., the brain and spinal cord (SC). In the SC, almost all previous human studies utilized Diffusion Tensor Imaging (DTI), which cannot accurately model areas where white matter (WM) pathways cross or diverge. While High Angular Diffusion Resolution Imaging (HARDI) can overcome some of these limitations, longer acquisition times critically limit its applicability to clinical human studies. In addition, previous human HARDI studies have used limited spatial resolution, with typically a few slices and voxel size ~1 × 1 × 5 mm3 being acquired in tens of minutes. Thus, we have optimized a novel fast HARDI protocol that allows collecting dMRI data at high angular and spatial resolutions in clinically-feasible time. Our data was acquired, using a 3T Siemens Prisma scanner, in less than 9 min. It has a total of 75 diffusion-weighted volumes and high spatial resolution of 0.67 × 0.67 × 3 mm3 (after interpolation in Fourier space) covering the cervical segments C4–C6. Our preliminary results demonstrate applicability of our technique in healthy individuals with good correspondence between low fractional anisotropy (FA) gray matter areas from the dMRI scans, and the same regions delineated on T2-weighted MR images with spatial resolution of 0.35 × 0.35 × 2.5 mm3. Our data also allows the detection of crossing fibers that were previously shown in vivo only in animal studies.

We propose a sparse Bayesian learning algorithm for improved estimation of white matter fiber parameters from compressed (under-sampled q-space) multi-shell diffusion MRI data. The multi-shell data is represented in a dictionary form using a non-monoexponential decay model of diffusion, based on continuous gamma distribution of diffusivities. The fiber volume fractions with predefined orientations, which are the unknown parameters, form the dictionary weights. These unknown parameters are estimated with a linear un-mixing framework, using a sparse Bayesian learning algorithm. A localized learning of hyperparameters at each voxel and for each possible fiber orientations improves the parameter estimation. Our experiments using synthetic data from the ISBI 2012 HARDI reconstruction challenge and in-vivo data from the Human Connectome Project demonstrate the improvements.

Although brain functionality is often remarkably robust to lesions and other insults, it may be fragile when these take place in specific locations. Previous attempts to quantify robustness and fragility sought to understand how the functional connectivity of brain networks is affected by structural changes, using either model-based predictions or empirical studies of the effects of lesions. We advance a geometric viewpoint relying on a notion of network curvature, the so-called Ollivier-Ricci curvature. This approach has been proposed to assess financial market robustness and to differentiate biological networks of cancer cells from healthy ones. Here, we apply curvature-based measures to brain structural networks to identify robust and fragile brain regions in healthy subjects. We show that curvature can also be used to track changes in brain connectivity related to age and autism spectrum disorder (ASD), and we obtain results that are in agreement with previous MRI studies., The brain can often continue to function despite lesions in many areas, but damage to particular locations may have serious effects. Here, the authors use the concept of Ollivier-Ricci curvature to investigate the robustness of brain networks.

Diffusion MRI (dMRI) reveals microstructural features of the brain white matter by quantifying the anisotropic diffusion of water molecules within axonal bundles. Yet, identifying features such as axonal orientation dispersion, density, diameter, etc., in complex white matter fiber configurations (e.g. crossings) has proved challenging. Besides optimized data acquisition and advanced biophysical models, computational procedures to fit such models to the data are critical. However, these procedures have been largely overlooked by the dMRI microstructure community and new, more versatile, approaches are needed to solve complex biophysical model fitting problems. Existing methods are limited to models assuming single fiber orientation, relevant to limited brain areas like the corpus callosum, or multiple orientations but without the ability to extract detailed microstructural features. Here, we introduce a new and versatile optimization technique (MIX), which enables microstructure imaging of crossing white matter fibers. We provide a MATLAB implementation of MIX, and demonstrate its applicability to general microstructure models in fiber crossings using synthetic as well as ex-vivo and in-vivo brain data.

PURPOSE: In diffusion MRI, a technique known as diffusion spectrum imaging reconstructs the propagator with a discrete Fourier transform, from a Cartesian sampling of the diffusion signal. Alternatively, it is possible to directly reconstruct the orientation distribution function in q-ball imaging, providing so-called high angular resolution diffusion imaging. In between these two techniques, acquisitions on several spheres in q-space offer an interesting trade-off between the angular resolution and the radial information gathered in diffusion MRI. A careful design is central in the success of multishell acquisition and reconstruction techniques. METHODS: The design of acquisition in multishell is still an open and active field of research, however. In this work, we provide a general method to design multishell acquisition with uniform angular coverage. This method is based on a generalization of electrostatic repulsion to multishell. RESULTS: We evaluate the impact of our method using simulations, on the angular resolution in one and two bundles of fiber configurations. Compared to more commonly used radial sampling, we show that our method improves the angular resolution, as well as fiber crossing discrimination. DISCUSSION: We propose a novel method to design sampling schemes with optimal angular coverage and show the positive impact on angular resolution in diffusion MRI.

The acquisition of high angular resolution diffusion MRI is particularly long and subject motion can become an issue. The orientation distribution function (ODF) can be reconstructed online incrementally from diffusion-weighted MRI with a Kalman filtering framework. This online reconstruction provides real-time feedback throughout the acquisition process. In this article, the Kalman filter is first adapted to the reconstruction of the ODF in constant solid angle. Then, a method called STAR (STatistical Analysis of Residuals) is presented and applied to the online detection of motion in high angular resolution diffusion images. Compared to existing techniques, this method is image based and is built on top of a Kalman filter. Therefore, it introduces no additional scan time and does not require additional hardware. The performance of STAR is tested on simulated and real data and compared to the classical generalized likelihood ratio test. Successful detection of small motion is reported (rotation under 2°) with no delay and robustness to noise., National Institutes of Health (U.S.) (NIH grant Grant P41 RR008079), National Institutes of Health (U.S.) (NIH grant P41 EB015894), National Institutes of Health (U.S.) (NIH grant P30 NS057091), National Institutes of Health (U.S.) (Human Connectome Project U54 MH091657), United States. Air Force Office of Scientific Research (NSSEFF), National Science Foundation (U.S.), United States. Army Research Office, United States. Defense Advanced Research Projects Agency, United States. National Geospatial-Intelligence Agency, France. Agence nationale de la recherche (ANR NucleiPark), Institut national de recherche en informatique et en automatique (France), National Institutes of Health (U.S.) (Human Connectome Project, Grant R01 EB008432)

We introduce a novel data fitting procedure of multi compartment models of the brain white matter for diffusion MRI (dMRI) data. These biophysical models aim to characterize important micro-structure quantities like axonal radius, density and orientations. In order to describe the underlying tissue properties, a variety of models for intra-/extra-axonal diffusion signals have been proposed. Combinations of these analytic models are used to predict the diffusion MRI signal in multi-compartment settings. However, parameter estimation from these multi-compartment models is an ill-posed problem. Consequently, many existing fitting algorithms either rely on an initial grid search to find a good start point, or have strong assumptions like single fiber orientation to estimate some of these parameters from simpler models like the diffusion tensor (DT). In both cases, there is a trade-off between computational complexity and accuracy of the estimated parameters. Here, we describe a novel algorithm based on the separation of the Nonlinear Least Squares (NLLS) fitting problem, via Variable Projection Method , to search for non-linearly and linearly entering parameters independently. We use stochastic global search algorithms to find a global minimum, while estimating non-linearly entering parameters. The approach is independent of any starting point, and does not rely on estimates from simpler models. We show that the suggested algorithm is faster than algorithms involving grid search, and its greater accuracy and robustness are demonstrated on synthetic as well as ex-/in-vivo data.

We address the problem of the segmentation of cerebral white matter structures from diffusion tensor images (DTI). A DTI produces, from a set of diffusion-weighted MR images, tensor-valued images where each voxel is assigned with a 3 times 3 symmetric, positive-definite matrix. This second order tensor is simply the covariance matrix of a local Gaussian process, with zero-mean, modeling the average motion of water molecules. As we will show in this paper, the definition of a dissimilarity measure and statistics between such quantities is a nontrivial task which must be tackled carefully. We claim and demonstrate that, by using the theoretically well-founded differential geometrical properties of the manifold of multivariate normal distributions, it is possible to improve the quality of the segmentation results obtained with other dissimilarity measures such as the Euclidean distance or the Kullback-Leibler divergence. The main goal of this paper is to prove that the choice of the probability metric, i.e., the dissimilarity measure, has a deep impact on the tensor statistics and, hence, on the achieved results. We introduce a variational formulation, in the level-set framework, to estimate the optimal segmentation of a DTI according to the following hypothesis: Diffusion tensors exhibit a Gaussian distribution in the different partitions. We must also respect the geometric constraints imposed by the interfaces existing among the cerebral structures and detected by the gradient of the DTI. We show how to express all the statistical quantities for the different probability metrics. We validate and compare the results obtained on various synthetic data-sets, a biological rat spinal cord phantom and human brain DTIs

Registration of diffusion-weighted magnetic resonance images (DW-MRIs) is a key step for population studies, or construction of brain atlases, among other important tasks. Given the high dimensionality of the data, registration is usually performed by relying on scalar representative images, such as the fractional anisotropy (FA) and non diffusion-weighted (b0) images, thereby ignoring much of the directional information conveyed by DW-MR datasets itself. Alternatively, model-based registration algorithms have been proposed to exploit information on the preferred fiber orientation(s) at each voxel. Models such as the diffusion tensor or orientation distribution function (ODF) have been used for this purpose. Tensor-based registration methods rely on a model that does not completely capture the information contained in DW-MRIs, and largely depends on the accurate estimation of tensors. ODF-based approaches are more recent and computationally challenging, but also better describe complex fiber configurations thereby potentially improving the accuracy of DW-MRI registration. A new algorithm based on angular interpolation of the diffusion-weighted volumes was proposed for affine registration, and does not rely on any specific local diffusion model. In this work, we first extensively compare the performance of registration algorithms based on (i) angular interpolation, (ii) non diffusion-weighted scalar volume (b0), and (iii) diffusion tensor image (DTI). Moreover, we generalize the concept of angular interpolation (AI) to non-linear image registration, and implement it in the FMRIB Software Library (FSL). We demonstrate that AI registration of DW-MRIs is a powerful alternative to volume and tensor-based approaches. In particular, we show that AI improves the registration accuracy in many cases over existing state-of-the-art algorithms, while providing registered raw DW-MRI data, which can be used for any subsequent analysis.

Basal ganglia circuits are affected in neurological disorders such as Parkinson's disease (PD), essential tremor, dystonia and Tourette syndrome. Understanding the structural and functional connectivity of these circuits is critical for elucidating the mechanisms of the movement and neuropsychiatric disorders, and is vital for developing new therapeutic strategies such as deep brain stimulation (DBS). Knowledge about the connectivity of the human basal ganglia and thalamus has rapidly evolved over recent years through non-invasive imaging techniques, but has remained incomplete because of insufficient resolution and sensitivity of these techniques. Here, we present an imaging and computational protocol designed to generate a comprehensive in vivo and subject-specific, three-dimensional model of the structure and connections of the human basal ganglia. High-resolution structural and functional magnetic resonance images were acquired with a 7-Tesla magnet. Capitalizing on the enhanced signal-to-noise ratio (SNR) and enriched contrast obtained at high-field MRI, detailed structural and connectivity representations of the human basal ganglia and thalamus were achieved. This unique combination of multiple imaging modalities enabled the in-vivo visualization of the individual human basal ganglia and thalamic nuclei, the reconstruction of seven white-matter pathways and their connectivity probability that, to date, have only been reported in animal studies, histologically, or group-averaged MRI population studies. Also described are subject-specific parcellations of the basal ganglia and thalamus into sub-territories based on their distinct connectivity patterns. These anatomical connectivity findings are supported by functional connectivity data derived from resting-state functional MRI (R-fMRI). This work demonstrates new capabilities for studying basal ganglia circuitry, and opens new avenues of investigation into the movement and neuropsychiatric disorders, in individual human subjects.

International audience; The orientation distribution function (ODF) can be reconstructed online incrementally from diffusion-weighted MRI with a Kalman filtering framework. This online reconstruction can provide real-time feedback to the practitioner, especially appreciated for long acquisition protocols typical in Q-ball imaging. On top of the Kalman filter, we propose a method to evaluate online the reconstruction accuracy of the estimated ODF in constant solid angle. In addition, monitoring the residuals of the Kalman filter, we design, based on statistical tests, two algorithms for online detection of subject motion. The proposed techniques, tested on real and synthetic data under various experimental conditions, can detect rotation by angle less than 3°.

International audience; The diffusion orientation distribution function (ODF) can be reconstructed from q-ball imaging (QBI) to map the complex intravoxel structure of water diffusion. As acquisition time is particularly large for high angular resolution diffusion imaging (HARDI), fast estimation algorithms have recently been proposed, as an on-line feedback on the reconstruction accuracy. Thus the acquisition could be stopped or continued on demand. We adapt these real-time algorithms to the mathematically correct definition of ODF in constant solid angle (CSA), and develop a motion detection algorithm upon this reconstruction. Results of improved fiber crossing detection by CSA ODF are shown, and motion detection was implemented and tested in vivo.

International audience; OBJECTIVE: To determine whether there are diffusion abnormalities along the fibers connecting sensorimotor regions, including the primary sensorimotor areas and the striatum, in patients with writer's cramp using voxel-based diffusion analysis and fiber tracking. Recent studies have shown structural changes in these regions in writer's cramp. DESIGN: Patient and control group comparison. SETTING: Referral center for movement disorders. PARTICIPANTS: Twenty-six right-handed patients with writer's cramp and 26 right-handed healthy control subjects matched for sex and age. INTERVENTIONS: Clinical motor evaluations. MAIN OUTCOME MEASURES: Fractional anisotropy changes and results of fiber tracking in writer's cramp. RESULTS: Diffusion-tensor imaging revealed increased fractional anisotropy bilaterally in the white matter of the posterior limb of the internal capsule and adjacent structures in the patients with writer's cramp. Fiber tracking demonstrated that fractional anisotropy changes involve fiber tracts connecting the primary sensorimotor areas with subcortical structures. CONCLUSIONS: Diffusion abnormalities are present in fiber tracts connecting the primary sensorimotor areas with subcortical structures in writer's cramp. These abnormalities strengthen the role of the corticosubcortical pathways in the pathophysiologic mechanisms of writer's cramp.

In this article, we review recent mathematical models and computational methods for the processing of diffusion Magnetic Resonance Images, including state-of-the-art reconstruction of diffusion models, cerebral white matter connectivity analysis, and segmentation techniques. We focus on Diffusion Tensor Images (DTI) and Q-Ball Images (QBI).

We introduce an original approach for the cerebral white matter connectivity mapping from diffusion tensor imaging (DTI). Our method relies on a global modeling of the acquired magnetic resonance imaging volume as a Riemannian manifold whose metric directly derives from the diffusion tensor. These tensors will be used to measure physical three-dimensional distances between different locations of a brain diffusion tensor image. The key concept is the notion of geodesic distance that will allow us to find optimal paths in the white matter. We claim that such optimal paths are reasonable approximations of neural fiber bundles. The geodesic distance function can be seen as the solution of two theoretically equivalent but, in practice, significantly different problems in the partial differential equation framework: an initial value problem which is intrinsically dynamic, and a boundary value problem which is, on the contrary, intrinsically stationary. The two approaches have very different properties which make them more or less adequate for our problem and more or less computationally efficient. The dynamic formulation is quite easy to implement but has several practical drawbacks. On the contrary, the stationary formulation is much more tedious to implement; we will show, however, that it has many virtues which make it more suitable for our connectivity mapping problem. Finally, we will present different possible measures of connectivity, reflecting the degree of connectivity between different regions of the brain. We will illustrate these notions on synthetic and real DTI datasets.

International audience; Tensors are nowadays an increasing research domain in different areas, especially in image processing, motivated for example by diffusion tensor magnetic resonance imaging (DT-MRI). Up to now, algorithms and tools developed to deal with tensors were founded on the assumption of a matrix vector space with the constraint of remaining symmetric positive definite matrices. On the contrary, our approach is grounded on the theoretically well-founded differential geometrical properties of the space of multivariate normal distributions, where it is possible to define an affine-invariant Riemannian metric and express statistics on the manifold of symmetric positive definite matrices. In this paper, we focus on the contribution of these tools to the anisotropic filtering and regularization of tensor fields. To validate our approach we present promising results on both synthetic and real DT-MRI data.

International audience; We propose a novel, fast and robust technique for the computation of anatomical connectivity in the brain. Our approach exploits the information provided by Diffusion Tensor Magnetic Resonance Imaging (or DTI) and models the white matter by using Riemannian geometry and control theory. We show that it is possible, from a region of interest, to compute the geodesic distance to any other point and the associated optimal vector field. The latter can be used to trace shortest paths coinciding with neural fiber bundles. We also demonstrate that no explicit computation of those 3D curves is necessary to assess the degree of connectivity of the region of interest with the rest of the brain. We finally introduce a general local connectivity measure whose statistics along the optimal paths may be used to evaluate the degree of connectivity of any pair of voxels. All those quantities can be computed simultaneously in a Fast Marching framework, directly yielding the connectivity maps. Apart from being extremely fast, this method has other advantages such as the strict respect of the convoluted geometry of white matter, the fact that it is parameter-free, and its robustness to noise. We illustrate our technique by showing results on real and synthetic datasets. Our GCM (Geodesic Connectivity Mapping) algorithm is implemented in C++ and will be soon available on the web.

We survey the recent activities of the Odyssee Laboratory in the area of the application of mathematics to the design of models for studying brain anatomy and function. We start with the problem of reconstructing sources in MEG and EEG, and discuss the variational approach we have developed for solving these inverse problems. This motivates the need for geometric models of the head. We present a method for automatically and accurately extracting surface meshes of several tissues of the head from anatomical magnetic resonance (MR) images. Anatomical connectivity can be extracted from diffusion tensor magnetic resonance images but, in the current state of the technology, it must be preceded by a robust estimation and regularization stage. We discuss our work based on variational principles and show how the results can be used to track fibers in the white matter (WM) as geodesics in some Riemannian space. We then go to the statistical modeling of functional magnetic resonance imaging (fMRI) signals from the viewpoint of their decomposition in a pseudo-deterministic and stochastic part that we then use to perform clustering of voxels in a way that is inspired by the theory of support vector machines and in a way that is grounded in information theory. Multimodal image matching is discussed next in the framework of image statistics and partial differential equations (PDEs) with an eye on registering fMRI to the anatomy. The paper ends with a discussion of a new theory of random shapes that may prove useful in building anatomical and functional atlases.

Quantitative spinal cord (SC) magnetic resonance imaging (MRI) presents many challenges, including a lack of standardized imaging protocols. Here we present a prospectively harmonized quantitative MRI protocol, which we refer to as the spine generic protocol, for users of 3T MRI systems from the three main manufacturers: GE, Philips and Siemens. The protocol provides guidance for assessing SC macrostructural and microstructural integrity: T1-weighted and T2-weighted imaging for SC cross-sectional area computation, multi-echo gradient echo for gray matter cross-sectional area, and magnetization transfer and diffusion weighted imaging for assessing white matter microstructure. In a companion paper from the same authors, the spine generic protocol was used to acquire data across 42 centers in 260 healthy subjects. The key details of the spine generic protocol are also available in an open-access document that can be found at https://github.com/spine-generic/protocols . The protocol will serve as a starting point for researchers and clinicians implementing new SC imaging initiatives so that, in the future, inclusion of the SC in neuroimaging protocols will be more common. The protocol could be implemented by any trained MR technician or by a researcher/clinician familiar with MRI acquisition.

Available online 22 August 2021. White matter bundle segmentation using diffusion MRI fiber tractography has become the method of choice to identify white matter fiber pathways in vivo in human brains. However, like other analyses of complex data, there is considerable variability in segmentation protocols and techniques. This can result in different reconstructions of the same intended white matter pathways, which directly affects tractography results, quantification, and interpretation. In this study, we aim to evaluate and quantify the variability that arises from different protocols for bundle segmentation. Through an open call to users of fiber tractography, including anatomists, clinicians, and algorithm developers, 42 independent teams were given processed sets of human whole-brain streamlines and asked to segment 14 white matter fascicles on six subjects. In total, we received 57 different bundle segmentation protocols, which enabled detailed volume-based and streamline-based analyses of agreement and disagreement among protocols for each fiber pathway. Results show that even when given the exact same sets of underlying streamlines, the variability across protocols for bundle segmentation is greater than all other sources of variability in the virtual dissection process, including variability within protocols and variability across subjects. In order to foster the use of tractography bundle dissection in routine clinical settings, and as a fundamental analytical tool, future endeavors must aim to resolve and reduce this heterogeneity. Although external validation is needed to verify the anatomical accuracy of bundle dissections, reducing heterogeneity is a step towards reproducible research and may be achieved through the use of standard nomenclature and definitions of white matter bundles and well-chosen constraints and decisions in the dissection process. This work was conducted in part using the resources of the Ad- vanced Computing Center for Research and Education at Vanderbilt University, Nashville, TN. KS, BL, CH were supported by the Na- tional Institutes of Health under award numbers R01EB017230, and T32EB001628, and in part by ViSE/VICTR VR3029 and the National Center for Research Resources, Grant UL1 RR024975-01. This work was also possible thanks to the support of the Institutional Research Chair in NeuroInformatics of Universitéde Sherbrooke, NSERC and Compute Canada (MD, FR). MP received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sk ł odowska-Curie grant agreement No 754462. The Wisconsin group acknowledges the support from a core grant to the Waisman Cen- ter from the National Institute of Child Health and Human Develop- ment (IDDRC U54 HD090256). NSF OAC-1916518, NSF IIS-1912270, NSF IIS-1636893, NSF BCS-1734853, NIH NIBIB 1R01EB029272-01, and a Microsoft Faculty Fellowship to F.P. LF acknowledges the sup- port of the Cluster of Excellence Matters of Activity. Image Space Mate- rial funded by the Deutsche Forschungsgemeinschaft (DFG, German Re- search Foundation) under Germany´s Excellence Strategy –EXC 2025. SW is supported by a Medical Research Council PhD Studentship UK [MR/N013913/1]. The Nottingham group’s processing was performed using the University of Nottingham’s Augusta HPC service and the Pre- cision Imaging Beacon Cluster. JPA, MA and SMS acknowledges the support of FCT - Fundação para a Ciência e a Tecnologia within CIN- TESIS, R&D Unit (reference UID/IC/4255/2013). MM was funded by the Wellcome Trust through a Sir Henry Wellcome Postdoctoral Fellow- ship [213722/Z/18/Z]. EJC-R is supported by the Swiss National Sci- ence Foundation (SNSF, Ambizione grant PZ00P2 185814/1). CMWT is supported by a Sir Henry Wellcome Fellowship (215944/Z/19/Z) and a Veni grant from the Dutch Research Council (NWO) (17331). FC acknowledges the support of the National Health and Medical Re- search Council of Australia (APP1091593 and APP1117724) and the Australian Research Council (DP170101815). NSF OAC-1916518, NSF IIS-1912270, NSF IIS-1636893, NSF BCS-1734853, Microsoft Faculty Fellowship to F.P. D.B. was partially supported by NIH NIMH T32- MH103213 to William Hetrick (Indiana University). CL is partly sup- ported by NIH grants P41 EB027061 and P30 NS076408 “Institutional Center Cores for Advanced Neuroimaging. JYMY received positional funding from the Royal Children’s Hospital Foundation (RCH 1000). JYMY, JC, and CEK acknowledge the support of the Royal Children’s Hospital Foundation, Murdoch Children’s Research Institute, The Uni- versity of Melbourne Department of Paediatrics, and the Victorian Gov- ernment’s Operational Infrastructure Support Program. C-HY is grateful to the Ministry of Science and Technology of Taiwan (MOST 109-2222- E-182-001-MY3) for the support. LC acknowledges support from CONA- CYT and UNAM. ARM acknowledges support from CONACYT. LJO, YR, and FZ were supported by NIH P41EB015902 and R01MH119222. AJG was supported by P41EB015898. NM was supported by R01MH119222, K24MH116366, and R01MH111917. This project has received funding from the European Union’s Horizon 2020 Research and Innovation Pro- gramme under Grant Agreement No. 785907 & 945539 (HBP SGA2 & SGA3), and from the ANR IFOPASUBA- 19-CE45-0022-01. PG, CR, NL and AV were partially supported by ANID-Basal FB0008 and ANID- FONDECYT 1190701 grants. We would like to acknowledge John C Gore, Hiromasa Takemura, Anastasia Yendiki, and Riccardo Galbusera for their helplful suggestions regarding the analysis, figures, and discussions.