Using GPUs to accelerate computational diffusion MRI: From microstructure estimation to tractography and connectomes

Clinical applicability of automated tractography for stroke rehabilitation: Z-score conversion of fractional anisotropy

[Purpose] To expand the applicability of diffusion-tensor tractography fractional anisotropy for stroke rehabilitation, this study aimed to provide references for representative neural tracts from non-lesioned hemispheres. Therefore, we applied the assessment of neural integrity to representative stroke patients using Z-score conversion. [Participants and Methods] Fractional anisotropy values were assessed in neural tracts, including the corticospinal tract, inferior fronto-occipital fasciculus, uncinate fasciculus, and anterior thalamic radiation, of stroke patients receiving acute care. [Results] Data were collected from 60 patients for the non-lesioned right hemisphere and 68 patients for the non-lesioned left hemisphere. Mean fractional anisotropy values in the corticospinal tract and inferior fronto-occipital fasciculus were notably elevated, reaching approximately 0.6 and 0.5, respectively. The mean fractional anisotropy values for other neural tracts were approximately 0.4, and, the overall standard deviations were approximately 0.04. In two typical stroke patients assessed using Z-scores, the scores in the corticospinal tract corresponded to the severity of the hemiparesis. The scores in the anterior thalamic radiation and inferior fronto-occipital fasciculus were associated with more significant brain dysfunction, including inattention and aphasia. [Conclusion] In this study, the Z-score findings related to stroke symptoms align with those reported in the literature, indicating the appropriateness of the methodology used and its potential in future applications.

Predicting 2-year neurodevelopmental outcomes in preterm infants using multimodal structural brain magnetic resonance imaging with local connectivity

The neurodevelopmental outcomes of preterm infants can be stratified based on the level of prematurity. We explored brain structural networks in extremely preterm (EP; < 28 weeks of gestation) and very-to-late (V-LP; ≥ 28 and < 37 weeks of gestation) preterm infants at term-equivalent age to predict 2-year neurodevelopmental outcomes. Using MRI and diffusion MRI on 62 EP and 131 V-LP infants, we built a multimodal feature set for volumetric and structural network analysis. We employed linear and nonlinear machine learning models to predict the Bayley Scales of Infant and Toddler Development, Third Edition (BSID-III) scores, assessing predictive accuracy and feature importance. Our findings revealed that models incorporating local connectivity features demonstrated high predictive performance for BSID-III subsets in preterm infants. Specifically, for cognitive scores in preterm (variance explained, 17%) and V-LP infants (variance explained, 17%), and for motor scores in EP infants (variance explained, 15%), models with local connectivity features outperformed others. Additionally, a model using only local connectivity features effectively predicted language scores in preterm infants (variance explained, 15%). This study underscores the value of multimodal feature sets, particularly local connectivity, in predicting neurodevelopmental outcomes, highlighting the utility of machine learning in understanding microstructural changes and their implications for early intervention.

Spatial and temporal pattern of structure-function coupling of human brain connectome with development

Brain structural circuitry shapes a richly patterned functional synchronization, supporting for complex cognitive and behavioural abilities. However, how coupling of structural connectome (SC) and functional connectome (FC) develops and its relationships with cognitive functions and transcriptomic architecture remain unclear. We used multimodal magnetic resonance imaging data from 439 participants aged 5.7 to 21.9 years to predict functional connectivity by incorporating intracortical and extracortical structural connectivity, characterizing SC-FC coupling. Our findings revealed that SC-FC coupling was strongest in the visual and somatomotor networks, consistent with evolutionary expansion, myelin content, and functional principal gradient. As development progressed, SC-FC coupling exhibited heterogeneous alterations dominated by an increase in cortical regions, broadly distributed across the somatomotor, frontoparietal, dorsal attention, and default mode networks. Moreover, we discovered that SC-FC coupling significantly predicted individual variability in general intelligence, mainly influencing frontoparietal and default mode networks. Finally, our results demonstrated that the heterogeneous development of SC-FC coupling is positively associated with genes in oligodendrocyte-related pathways and negatively associated with astrocyte-related genes. This study offers insight into the maturational principles of SC-FC coupling in typical development.

Feature similarity gradients detect alterations in the neonatal cortex associated with preterm birth

The early life environment programmes cortical architecture and cognition across the life course. A measure of cortical organisation that integrates information from multimodal MRI and is unbound by arbitrary parcellations has proven elusive, which hampers efforts to uncover the perinatal origins of cortical health. Here, we use the Vogt-Bailey index to provide a fine-grained description of regional homogeneities and sharp variations in cortical microstructure based on feature gradients, and we investigate the impact of being born preterm on cortical development at term-equivalent age. Compared with term-born controls, preterm infants have a homogeneous microstructure in temporal and occipital lobes, and the medial parietal, cingulate, and frontal cortices, compared with term infants. These observations replicated across two independent datasets and were robust to differences that remain in the data after matching samples and alignment of processing and quality control strategies. We conclude that cortical microstructural architecture is altered in preterm infants in a spatially distributed rather than localised fashion.

GPU-Based Parallel Processing Techniques for Enhanced Brain Magnetic Resonance Imaging Analysis: A Review of Recent Advances

The approach of using more than one processor to compute in order to overcome the complexity of different medical imaging methods that make up an overall job is known as GPU (graphic processing unit)-based parallel processing. It is extremely important for several medical imaging techniques such as image classification, object detection, image segmentation, registration, and content-based image retrieval, since the GPU-based parallel processing approach allows for time-efficient computation by a software, allowing multiple computations to be completed at once. On the other hand, a non-invasive imaging technology that may depict the shape of an anatomy and the biological advancements of the human body is known as magnetic resonance imaging (MRI). Implementing GPU-based parallel processing approaches in brain MRI analysis with medical imaging techniques might be helpful in achieving immediate and timely image capture. Therefore, this extended review (the extension of the IWBBIO2023 conference paper) offers a thorough overview of the literature with an emphasis on the expanding use of GPU-based parallel processing methods for the medical analysis of brain MRIs with the imaging techniques mentioned above, given the need for quicker computation to acquire early and real-time feedback in medicine. Between 2019 and 2023, we examined the articles in the literature matrix that include the tasks, techniques, MRI sequences, and processing results. As a result, the methods discussed in this review demonstrate the advancements achieved until now in minimizing computing runtime as well as the obstacles and problems still to be solved in the future.

On the stability of canonical correlation analysis and partial least squares with application to brain-behavior associations

Associations between datasets can be discovered through multivariate methods like Canonical Correlation Analysis (CCA) or Partial Least Squares (PLS). A requisite property for interpretability and generalizability of CCA/PLS associations is stability of their feature patterns. However, stability of CCA/PLS in high-dimensional datasets is questionable, as found in empirical characterizations. To study these issues systematically, we developed a generative modeling framework to simulate synthetic datasets. We found that when sample size is relatively small, but comparable to typical studies, CCA/PLS associations are highly unstable and inaccurate; both in their magnitude and importantly in the feature pattern underlying the association. We confirmed these trends across two neuroimaging modalities and in independent datasets with n ≈ 1000 and n = 20,000, and found that only the latter comprised sufficient observations for stable mappings between imaging-derived and behavioral features. We further developed a power calculator to provide sample sizes required for stability and reliability of multivariate analyses. Collectively, we characterize how to limit detrimental effects of overfitting on CCA/PLS stability, and provide recommendations for future studies.

Age of onset of obsessive-compulsive disorder differentially affects white matter microstructure

Previous diffusion MRI studies have reported mixed findings on white matter microstructure alterations in obsessive-compulsive disorder (OCD), likely due to variation in demographic and clinical characteristics, scanning methods, and underpowered samples. The OCD global study was created across five international sites to overcome these challenges by harmonizing data collection to identify consistent brain signatures of OCD that are reproducible and generalizable. Single-shell diffusion measures (e.g., fractional anisotropy), multi-shell Neurite Orientation Dispersion and Density Imaging (NODDI) and fixel-based measures, were extracted from skeletonized white matter tracts in 260 medication-free adults with OCD and 252 healthy controls. We additionally performed structural connectome analysis. We compared cases with controls and cases with early (<18) versus late (18+) OCD onset using mixed-model and Bayesian multilevel analysis. Compared with healthy controls, adult OCD individuals showed higher fiber density in the sagittal stratum (B[SE] = 0.10[0.05], P = 0.04) and credible evidence for higher fiber density in several other tracts. When comparing early (n = 145) and late-onset (n = 114) cases, converging evidence showed lower integrity of the posterior thalamic radiation -particularly radial diffusivity (B[SE] = 0.28[0.12], P = 0.03)-and lower global efficiency of the structural connectome (B[SE] = 15.3[6.6], P = 0.03) in late-onset cases. Post-hoc analyses indicated divergent direction of effects of the two OCD groups compared to healthy controls. Age of OCD onset differentially affects the integrity of thalamo-parietal/occipital tracts and the efficiency of the structural brain network. These results lend further support for the role of the thalamus and its afferent fibers and visual attentional processes in the pathophysiology of OCD.

Applicability of multiple quantitative magnetic resonance methods in genetic brain white matter disorders

BACKGROUND AND PURPOSE

Magnetic resonance imaging (MRI) measures of tissue microstructure are important for monitoring brain white matter (WM) disorders like leukodystrophies and multiple sclerosis. They should be sensitive to underlying pathological changes. Three whole-brain isotropic quantitative methods were applied and compared within a cohort of controls and leukodystrophy patients: two novel myelin water imaging (MWI) techniques (multi-compartment relaxometry diffusion-informed MWI: MCR-DIMWI, and multi-echo T2 relaxation imaging with compressed sensing: METRICS) and neurite orientation dispersion and density imaging (NODDI).

METHODS

For 9 patients with different leukodystrophies (age range 0.4-62.4 years) and 15 control subjects (2.3-61.3 years), T1-weighted MRI, fluid-attenuated inversion recovery, multi-echo gradient echo with variable flip angles, METRICS, and multi-shell diffusion-weighted imaging were acquired on 3 Tesla. MCR-DIMWI, METRICS, NODDI, and quality control measures were extracted to evaluate differences between patients and controls in WM and deep gray matter (GM) regions of interest (ROIs). Pearson correlations, effect size calculations, and multi-level analyses were performed.

RESULTS

MCR-DIMWI and METRICS-derived myelin water fractions (MWFs) were lower and relaxation times were higher in patients than in controls. Effect sizes of MWF values and relaxation times were large for both techniques. Differences between patients and controls were more pronounced in WM ROIs than in deep GM. MCR-DIMWI-MWFs were more homogeneous within ROIs and more bilaterally symmetrical than METRICS-MWFs. The neurite density index was more sensitive in detecting differences between patients and controls than fractional anisotropy. Most measures obtained from MCR-DIMWI, METRICS, NODDI, and diffusion tensor imaging correlated strongly with each other.

CONCLUSION

This proof-of-concept study shows that MCR-DIMWI, METRICS, and NODDI are sensitive techniques to detect changes in tissue microstructure in WM disorders.

Longitudinal development of the human white matter structural connectome and its association with brain transcriptomic and cellular architecture

From childhood to adolescence, the spatiotemporal development pattern of the human brain white matter connectome and its underlying transcriptomic and cellular mechanisms remain largely unknown. With a longitudinal diffusion MRI cohort of 604 participants, we map the developmental trajectory of the white matter connectome from global to regional levels and identify that most brain network properties followed a linear developmental trajectory. Importantly, connectome-transcriptomic analysis reveals that the spatial development pattern of white matter connectome is potentially regulated by the transcriptomic architecture, with positively correlated genes involve in ion transport- and development-related pathways expressed in excitatory and inhibitory neurons, and negatively correlated genes enriches in synapse- and development-related pathways expressed in astrocytes, inhibitory neurons and microglia. Additionally, the macroscale developmental pattern is also associated with myelin content and thicknesses of specific laminas. These findings offer insights into the underlying genetics and neural mechanisms of macroscale white matter connectome development from childhood to adolescence.

Applicability of fractional anisotropy from standardized automated tractography for outcome prediction of patients after stroke

[Purpose] Diffusion-tensor fractional anisotropy has been used for outcome prediction in stroke patients. We assessed the clinical applicability of the two major fractional anisotropy methodologies-fractional anisotropy derived from segmentation maps in the standard brain (region of interest) and fractional anisotropy derived from standardized automated tractography-in relation to outcomes. [Participants and Methods] The study design was a retrospective survey of medical records collected from October 2021 to September 2022. Diffusion-tensor imaging was conducted in the second week after stroke onset. Outcomes were assessed using the total score of the motor component of the Stroke Impairment Assessment Set (null to full, 0 to 25). Correlations between fractional anisotropy and the outcomes were then assessed. [Results] Fourteen patients with hemorrhagic stroke were sampled. The fractional anisotropy from standardized automated tractography of the corticospinal tract on the lesion side (mean ± standard deviation, 0.403 ± 0.070) was significantly and tightly correlated (r=0.813) with the outcomes (13.4 ± 9.2), whereas the fractional anisotropy from a region of interest set in the cerebral peduncle on the lesion side (0.548 ± 0.064) was not significantly correlated with the outcomes (r=0.507). [Conclusion] The findings suggest that fractional anisotropy derived from standardized automated tractography can be more applicable to outcome prediction than that derived from a region of interest defined in the standard brain.

Automated Tractography for the Assessment of Aphasia in Acute Care Stroke Rehabilitation: A Case Series

Background

Aphasia is a common disorder among stroke patients. Assessment of aphasia is essential for scheduling appropriate rehabilitative treatment. Although this is conventionally accomplished using neuropsychological test batteries, these tests are not always accessible because of attention and/or consciousness disturbances during acute care. To overcome this issue, we have introduced a newly developed automated tractography known as XTRACT.

Cases

Diffusion-tensor images were acquired from three patients on days 10-14. Brain images were processed by XTRACT, which automatically extracts neural tracts using standardized protocols. Fractional anisotropy (FA) values were then bilaterally evaluated in the following neural tracts associated with aphasia: arcuate fasciculus, inferior fronto-occipital fasciculus, middle longitudinal fasciculus, inferior longitudinal fasciculus, and uncinate fasciculus. Case 1 had word-finding difficulty on admission. FA values in the lesioned left hemisphere were not decreased in all tracts and this patient fully recovered during acute care. Case 2 had reduced spontaneous speech and a low FA value in the left arcuate fasciculus. Rehabilitative treatment was scheduled to improve the verbal output of sentences and word recall. Case 3 could not complete the conventional aphasia test battery because of attention disturbance. He had low FA values in all tracts in the left hemisphere. Rehabilitative treatment was designed to focus on both speaking and auditory comprehension.

Discussion

Automated tractography enables quantitative assessment of the neural damage associated with aphasia, even in patients with attention and/or consciousness disturbances. This modality can aid in the assessment of aphasia and allows the planning of appropriate rehabilitative treatment.

Non-reversibility outperforms functional connectivity in characterisation of brain states in MEG data

Characterising brain states during tasks is common practice for many neuroscientific experiments using electrophysiological modalities such as electroencephalography (EEG) and magnetoencephalography (MEG). Brain states are often described in terms of oscillatory power and correlated brain activity, i.e. functional connectivity. It is, however, not unusual to observe weak task induced functional connectivity alterations in the presence of strong task induced power modulations using classical time-frequency representation of the data. Here, we propose that non-reversibility, or the temporal asymmetry in functional interactions, may be more sensitive to characterise task induced brain states than functional connectivity. As a second step, we explore causal mechanisms of non-reversibility in MEG data using whole brain computational models. We include working memory, motor, language tasks and resting-state data from participants of the Human Connectome Project (HCP). Non-reversibility is derived from the lagged amplitude envelope correlation (LAEC), and is based on asymmetry of the forward and reversed cross-correlations of the amplitude envelopes. Using random forests, we find that non-reversibility outperforms functional connectivity in the identification of task induced brain states. Non-reversibility shows especially better sensitivity to capture bottom-up gamma induced brain states across all tasks, but also alpha band associated brain states. Using whole brain computational models we find that asymmetry in the effective connectivity and axonal conduction delays play a major role in shaping non-reversibility across the brain. Our work paves the way for better sensitivity in characterising brain states during both bottom-up as well as top-down modulation in future neuroscientific experiments.

An open resource combining multi-contrast MRI and microscopy in the macaque brain

Understanding brain structure and function often requires combining data across different modalities and scales to link microscale cellular structures to macroscale features of whole brain organisation. Here we introduce the BigMac dataset, a resource combining in vivo MRI, extensive postmortem MRI and multi-contrast microscopy for multimodal characterisation of a single whole macaque brain. The data spans modalities (MRI and microscopy), tissue states (in vivo and postmortem), and four orders of spatial magnitude, from microscopy images with micrometre or sub-micrometre resolution, to MRI signals on the order of millimetres. Crucially, the MRI and microscopy images are carefully co-registered together to facilitate quantitative multimodal analyses. Here we detail the acquisition, curation, and first release of the data, that together make BigMac a unique, openly-disseminated resource available to researchers worldwide. Further, we demonstrate example analyses and opportunities afforded by the data, including improvement of connectivity estimates from ultra-high angular resolution diffusion MRI, neuroanatomical insight provided by polarised light imaging and myelin-stained histology, and the joint analysis of MRI and microscopy data for reconstruction of the microscopy-inspired connectome. All data and code are made openly available.

Microstructural Properties of the Cerebellar Peduncles in Children with Developmental Language Disorder

Children with developmental language disorder (DLD) struggle to learn their native language for no apparent reason. While research on the neurobiological underpinnings of the disorder has focused on the role of cortico-striatal systems, little is known about the role of the cerebellum in DLD. Cortico-cerebellar circuits might be involved in the disorder as they contribute to complex sensorimotor skill learning, including the acquisition of spoken language. Here, we used diffusion-weighted imaging data from 77 typically developing and 54 children with DLD and performed probabilistic tractography to identify the cerebellum's white matter tracts: the inferior, middle, and superior cerebellar peduncles. Children with DLD showed lower fractional anisotropy (FA) in the inferior cerebellar peduncles (ICP), fiber tracts that carry motor and sensory input via the inferior olive to the cerebellum. Lower FA in DLD was driven by lower axial diffusivity. Probing this further with more sophisticated modeling of diffusion data, we found higher orientation dispersion but no difference in neurite density in the ICP of DLD. Reduced FA is therefore unlikely to be reflecting microstructural differences in myelination in this tract, rather the organization of axons in these pathways is disrupted. ICP microstructure was not associated with language or motor coordination performance in our sample. We also found no differences in the middle and superior peduncles, the main pathways connecting the cerebellum with the cortex. To conclude, it is not cortico-cerebellar but atypical olivocerebellar white matter connections that characterize DLD and suggest the involvement of the olivocerebellar system in speech acquisition and development.

Characterization of Atypical Corticospinal Tract Microstructure and Hand Impairments in Early-Onset Hemiplegic Cerebral Palsy: Preliminary Findings

Unilateral brain injuries occurring before at or shortly after full-term can result in hemiplegic cerebral palsy (HCP). HCP affects one side of the body and can be characterized in the hand with measures of weakness and a loss of independent hand control resulting in mirror movements. Hand impairment severity is extremely heterogeneous across individuals with HCP and the neural basis for this variability is unclear. We used diffusion MRI and tractography to investigate the relationship between structural morphology of the supraspinal corticospinal tract (CST) and the severity of two typical hand impairments experienced by individuals with HCP, grasp weakness and mirror movements. Results from nine children with HCP and eight children with typical development show that there is a significant hemispheric association between CST microstructure and hand impairment severity that may be explained by atypical development and fiber distribution of motor pathways. Further analysis in the non-lesioned (dominant) hemisphere shows significant differences for CST termination in the cortex between participants with HCP and those with typical development. These findings suggest that structural disparities at the cellular level in the seemingly unaffected hemisphere after early unilateral brain injury may be the cause of heterogeneous hand impairments seen in this population.Clinical Relevance- Quantitative measurement of the variability in hand function in individuals with HCP is necessary to represent the distinct impairments experienced by each person. Further understanding of the structural neural morphology underlying distal upper extremity motor deficits after early unilateral brain injury will help lead to the development of more specific targeted interventions that increase functional outcomes.

A framework for focal and connectomic mapping of transiently disrupted brain function

The distributed nature of the neural substrate, and the difficulty of establishing necessity from correlative data, combine to render the mapping of brain function a far harder task than it seems. Methods capable of combining connective anatomical information with focal disruption of function are needed to disambiguate local from global neural dependence, and critical from merely coincidental activity. Here we present a comprehensive framework for focal and connective spatial inference based on sparse disruptive data, and demonstrate its application in the context of transient direct electrical stimulation of the human medial frontal wall during the pre-surgical evaluation of patients with focal epilepsy. Our framework formalizes voxel-wise mass-univariate inference on sparsely sampled data within the statistical parametric mapping framework, encompassing the analysis of distributed maps defined by any criterion of connectivity. Applied to the medial frontal wall, this transient dysconnectome approach reveals marked discrepancies between local and distributed associations of major categories of motor and sensory behaviour, revealing differentiation by remote connectivity to which purely local analysis is blind. Our framework enables disruptive mapping of the human brain based on sparsely sampled data with minimal spatial assumptions, good statistical efficiency, flexible model formulation, and explicit comparison of local and distributed effects.

A spatiotemporal complexity architecture of human brain activity

The human brain operates in large-scale functional networks. These networks are an expression of temporally correlated activity across brain regions, but how global network properties relate to the neural dynamics of individual regions remains incompletely understood. Here, we show that the brain's network architecture is tightly linked to critical episodes of neural regularity, visible as spontaneous "complexity drops" in functional magnetic resonance imaging signals. These episodes closely explain functional connectivity strength between regions, subserve the propagation of neural activity patterns, and reflect interindividual differences in age and behavior. Furthermore, complexity drops define neural activity states that dynamically shape the connectivity strength, topological configuration, and hierarchy of brain networks and comprehensively explain known structure-function relationships within the brain. These findings delineate a principled complexity architecture of neural activity-a human "complexome" that underpins the brain's functional network organization.

Knee Cartilage Defect Assessment by Graph Representation and Surface Convolution

Knee osteoarthritis (OA) is the most common osteoarthritis and a leading cause of disability. Cartilage defects are regarded as major manifestations of knee OA, which are visible by magnetic resonance imaging (MRI). Thus early detection and assessment for knee cartilage defects are important for protecting patients from knee OA. In this way, many attempts have been made on knee cartilage defect assessment by applying convolutional neural networks (CNNs) to knee MRI. However, the physiologic characteristics of the cartilage may hinder such efforts: the cartilage is a thin curved layer, implying that only a small portion of voxels in knee MRI can contribute to the cartilage defect assessment; heterogeneous scanning protocols further challenge the feasibility of the CNNs in clinical practice; the CNN-based knee cartilage evaluation results lack interpretability. To address these challenges, we model the cartilages structure and appearance from knee MRI into a graph representation, which is capable of handling highly diverse clinical data. Then, guided by the cartilage graph representation, we design a non-Euclidean deep learning network with the self-attention mechanism, to extract cartilage features in the local and global, and to derive the final assessment with a visualized result. Our comprehensive experiments show that the proposed method yields superior performance in knee cartilage defect assessment, plus its convenient 3D visualization for interpretability.

Applicability of automated tractography during acute care stroke rehabilitation

[Purpose] To assess the clinical applicability of a novel automated tractography tool named XTRACT during acute stroke rehabilitation. [Participants and Methods] Three patients with left hemisphere stroke were sampled. Diffusion tensor images were acquired on the second week, and automated tractography was then applied. Tractography images and fractional anisotropy (FA) values in the corticospinal tract (CST) and arcuate fasciculus (AF) were assessed in relation to hemiparesis and aphasia. [Results] Patient 1 was nearly asymptomatic; FA in the left CST was 0.610 and that in the AF was 0.509. Patient 2 had severe hemiparesis and mild motor aphasia. Tractography images of the CST and AF were blurred; FA in the left CST was 0.295 and that in the AF was 0.304. Patient 3 showed no hemiparesis or aphasia at initial assessment. Tractography image of the CST was intact but that of the AF was less clear; FA in the left CST was 0.586 and that in the AF was 0.338. Considering the less clear images of the AF and lower FA value in Patients 2 and 3, further examinations for aphasia were performed, which revealed agraphia. [Conclusion] Visualization and quantification of neural fibers using automated tractography promoted planning acute care rehabilitative treatment in patients with stroke.

White matter alterations associated with lifetime and current depression in adolescents: Evidence for cingulum disruptions

INTRODUCTION

Compared to research on adults with depression, relatively little work has examined white matter microstructure differences in depression arising earlier in life. Here we tested hypotheses about disruptions to white matter structure in adolescents with current and past depression, with an a priori focus on the cingulum bundles, uncinate fasciculi, corpus collosum, and superior longitudinal fasciculus.

METHODS

One hundred thirty-one children from the Preschool Depression Study were assessed using a Human Connectome Project style diffusion imaging sequence which was processed with HCP pipelines and TRACULA to generate estimates of fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD) and radial diffusivity (RD).

RESULTS

We found that reduced FA, reduced AD, and increased RD in the dorsal cingulum bundle were associated with a lifetime diagnosis of major depression and greater cumulative and current depression severity. Reduced FA, reduced AD, and increased RD in the ventral cingulum were associated with greater cumulative depression severity.

CONCLUSION

These findings support the emergence of white matter differences detected in adolescence associated with earlier life and concurrent depression. They also highlight the importance of connections of the cingulate to other brain regions in association with depression, potentially relevant to understanding emotion dysregulation and functional connectivity differences in depression.

Early structural connectivity within the sensorimotor network: Deviations related to prematurity and association to neurodevelopmental outcome

Consisting of distributed and interconnected structures that interact through cortico-cortical connections and cortico-subcortical loops, the sensorimotor (SM) network undergoes rapid maturation during the perinatal period and is thus particularly vulnerable to preterm birth. However, the impact of prematurity on the development and integrity of the emerging SM connections and their relationship to later motor and global impairments are still poorly understood. In this study we aimed to explore to which extent the early microstructural maturation of SM white matter (WM) connections at term-equivalent age (TEA) is modulated by prematurity and related with neurodevelopmental outcome at 18 months corrected age. We analyzed 118 diffusion MRI datasets from the developing Human Connectome Project (dHCP) database: 59 preterm (PT) low-risk infants scanned near TEA and a control group of full-term (FT) neonates paired for age at MRI and sex. We delineated WM connections between the primary SM cortices (S1, M1 and paracentral region) and subcortical structures using probabilistic tractography, and evaluated their microstructure with diffusion tensor imaging (DTI) and neurite orientation dispersion and density imaging (NODDI) models. To go beyond tract-specific univariate analyses, we computed a maturational distance related to prematurity based on the multi-parametric Mahalanobis distance of each PT infant relative to the FT group. Our results confirmed the presence of microstructural differences in SM tracts between PT and FT infants, with effects increasing with lower gestational age at birth. Maturational distance analyses highlighted that prematurity has a differential effect on SM tracts with higher distances and thus impact on (i) cortico-cortical than cortico-subcortical connections; (ii) projections involving S1 than M1 and paracentral region; and (iii) the most rostral cortico-subcortical tracts, involving the lenticular nucleus. These different alterations at TEA suggested that vulnerability follows a specific pattern coherent with the established WM caudo-rostral progression of maturation. Finally, we highlighted some relationships between NODDI-derived maturational distances of specific tracts and fine motor and cognitive outcomes at 18 months. As a whole, our results expand understanding of the significant impact of premature birth and early alterations on the emerging SM network even in low-risk infants, with possible relationship with neurodevelopmental outcomes. This encourages further exploration of these potential neuroimaging markers for prediction of neurodevelopmental disorders, with special interest for subtle neuromotor impairments frequently observed in preterm-born children.

Estimating axial diffusivity in the NODDI model

To estimate microstructure-related parameters from diffusion MRI data, biophysical models make strong, simplifying assumptions about the underlying tissue. The extent to which many of these assumptions are valid remains an open research question. This study was inspired by the disparity between the estimated intra-axonal axial diffusivity from literature and that typically assumed by the Neurite Orientation Dispersion and Density Imaging (NODDI) model (d_∥=1.7μm²/ms). We first demonstrate how changing the assumed axial diffusivity results in considerably different NODDI parameter estimates. Second, we illustrate the ability to estimate axial diffusivity as a free parameter of the model using high b-value data and an adapted NODDI framework. Using both simulated and in vivo data we investigate the impact of fitting to either real-valued or magnitude data, with Gaussian and Rician noise characteristics respectively, and what happens if we get the noise assumptions wrong in this high b-value and thus low SNR regime. Our results from real-valued human data estimate intra-axonal axial diffusivities of ∼2-2.5μm²/ms, in line with current literature. Crucially, our results demonstrate the importance of accounting for both a rectified noise floor and/or a signal offset to avoid biased parameter estimates when dealing with low SNR data.

Diffusion magnetic resonance imaging connectome features are predictive of functional lateralization of semantic processing in the anterior temporal lobes

Assessment of regional language lateralization is crucial in many scenarios, but not all populations are suited for its evaluation via task-functional magnetic resonance imaging (fMRI). In this study, the utility of structural connectome features for the classification of language lateralization in the anterior temporal lobes (ATLs) was investigated. Laterality indices for semantic processing in the ATL were computed from task-fMRI in 1038 subjects from the Human Connectome Project who were labeled as stronger rightward lateralized (RL) or stronger leftward to bilaterally lateralized (LL) in a data-driven approach. Data of unrelated subjects (n = 432) were used for further analyses. Structural connectomes were generated from diffusion-MRI tractography, and graph theoretical metrics (node degree, betweenness centrality) were computed. A neural network (NN) and a random forest (RF) classifier were trained on these metrics to classify subjects as RL or LL. After classification, comparisons of network measures were conducted via permutation testing. Degree-based classifiers produced significant above-chance predictions both during cross-validation (NN: AUC-ROC[CI] = 0.68[0.64-0.73], accuracy[CI] = 68.34%[63-73.2%]; RF: AUC-ROC[CI] = 0.7[0.66-0.73], accuracy[CI] = 64.81%[60.9-68.5]) and testing (NN: AUC-ROC[CI] = 0.69[0.53-0.84], accuracy[CI] = 68.09[53.2-80.9]; RF: AUC-ROC[CI] = 0.68[0.53-0.84], accuracy[CI] = 68.09[55.3-80.9]). Comparison of network metrics revealed small effects of increased node degree within the right posterior middle temporal gyrus (pMTG) in subjects with RL, while degree was decreased in the right posterior cingulate cortex (PCC). Above-chance predictions of functional language lateralization in the ATL are possible based on diffusion-MRI connectomes alone. Increased degree within the right pMTG as a right-sided homologue of a known semantic hub, and decreased hubness of the right PCC may form a structural basis for rightward-lateralized semantic processing.

The structural brain network topology of episodic memory

Episodic memory is supported by a distributed network of brain regions, and this complex network of regions does not operate in isolation. To date, neuroscience research in this area has typically focused on the activation levels in specific regions or pairwise connectivity between such regions. However, research has yet to investigate how the complex interactions of structural brain networks influence episodic memory abilities. We applied graph theory methods to diffusion-based anatomical networks in order to examine the structural architecture of the medial temporal lobe needed to support effective episodic memory functioning. We examined the relationship between performance on tests of verbal and non-verbal episodic memory with node strength, which indexes how well connected a brain region is in the network. Findings mapped onto the Posterior Medial memory system, subserved by the parahippocampal cortex and overlapped with findings of previous studies of episodic memory employing different methodologies. This expands our current understanding by providing independent evidence for the importance of identified regions and suggesting the particular manner in which these regions support episodic memory.

Laterality and Sex Differences of Human Lateral Habenula Afferent and Efferent Fiber Tracts

Introduction

The lateral habenula (LHb) is an epithalamic nucleus associated with negative valence and affective disorders. It receives input via the stria medullaris (SM) and sends output via the fasciculus retroflexus (FR). Here, we use tractography to reconstruct and characterize this pathway.

Methods

Multi-shell human diffusion magnetic resonance imaging (dMRI) data was obtained from the human connectome project (HCP) (n = 20, 10 males) and from healthy controls (n = 10, 6 males) scanned at our institution. We generated LHb afferents and efferents using probabilistic tractography by selecting the pallidum as the seed region and the ventral tegmental area as the output target.

Results

We were able to reconstruct the intended streamlines in all individuals from the HCP dataset and our dataset. Our technique also aided in identification of the LHb. In right-handed individuals, the streamlines were significantly more numerous in the left hemisphere (mean ratio 1.59 ± 0.09, p = 0.04). In left-handed individuals, there was no hemispheric asymmetry on average (mean ratio 1.00 ± 0.09, p = 1.0). Additionally, these streamlines were significantly more numerous in females than in males (619.9 ± 159.7 vs. 225.9 ± 66.03, p = 0.04).

Conclusion

We developed a method to reconstruct the SM and FR without manual identification of the LHb. This technique enables targeting of these fiber tracts as well as the LHb. Furthermore, we have demonstrated that there are sex and hemispheric differences in streamline number. These findings may have therapeutic implications and warrant further investigation.

Abnormal thalamocortical connectivity of preterm infants with elevated thyroid stimulating hormone identified with diffusion tensor imaging

While thyroid disturbances during perinatal and postnatal periods in preterm infants with congenital hypothyroidism reportedly disrupt neuronal development, no study has considered the effect of thyroid disturbances in premature infants with subclinical hypothyroidism with elevations of thyroid stimulating hormone. We aimed to identify altered fiber integrity from the thalamus to cortices in preterm infants with subclinical hypothyroidism. All preterm infants born were categorized according to thyroid stimulating hormone levels through serial thyroid function tests (36 preterm controls and 29 preterm infants with subclinical hypothyroidism). Diffusion tensor images were acquired to determine differences in thalamocortical fiber lengths between the groups, and cerebral asymmetries were investigated to observe neurodevelopmental changes. Thalamocortical fiber lengths in the subclinical hypothyroidism group were significantly reduced in the bilateral superior temporal gyrus, heschl’s gyrus, lingual gyrus, and calcarine cortex (all p < 0.05). According to the asymmetric value in the orbitofrontal regions, there is a left dominance in the subclinical hypothyroidism group contrary to the controls (p = 0.012), and that of the cuneus areas showed significant decreases in the subclinical hypothyroidism group (p = 0.035). These findings could reflect altered neurodevelopment, which could help treatment plans using biomarkers for subclinical hypothyroidism.

Combining Neuroimaging and Omics Datasets for Disease Classification Using Graph Neural Networks

Both neuroimaging and genomics datasets are often gathered for the detection of neurodegenerative diseases. Huge dimensionalities of neuroimaging data as well as omics data pose tremendous challenge for methods integrating multiple modalities. There are few existing solutions that can combine both multi-modal imaging and multi-omics datasets to derive neurological insights. We propose a deep neural network architecture that combines both structural and functional connectome data with multi-omics data for disease classification. A graph convolution layer is used to model functional magnetic resonance imaging (fMRI) and diffusion tensor imaging (DTI) data simultaneously to learn compact representations of the connectome. A separate set of graph convolution layers are then used to model multi-omics datasets, expressed in the form of population graphs, and combine them with latent representations of the connectome. An attention mechanism is used to fuse these outputs and provide insights on which omics data contributed most to the model's classification decision. We demonstrate our methods for Parkinson's disease (PD) classification by using datasets from the Parkinson's Progression Markers Initiative (PPMI). PD has been shown to be associated with changes in the human connectome and it is also known to be influenced by genetic factors. We combine DTI and fMRI data with multi-omics data from RNA Expression, Single Nucleotide Polymorphism (SNP), DNA Methylation and non-coding RNA experiments. A Matthew Correlation Coefficient of greater than 0.8 over many combinations of multi-modal imaging data and multi-omics data was achieved with our proposed architecture. To address the paucity of paired multi-modal imaging data and the problem of imbalanced data in the PPMI dataset, we compared the use of oversampling against using CycleGAN on structural and functional connectomes to generate missing imaging modalities. Furthermore, we performed ablation studies that offer insights into the importance of each imaging and omics modality for the prediction of PD. Analysis of the generated attention matrices revealed that DNA Methylation and SNP data were the most important omics modalities out of all the omics datasets considered. Our work motivates further research into imaging genetics and the creation of more multi-modal imaging and multi-omics datasets to study PD and other complex neurodegenerative diseases.

Autologous cellular therapy for cerebral palsy: a randomized, crossover trial

We examined an autologous mononuclear-cell-therapy-based approach to treat cerebral palsy using autologous umbilical cord blood or bone-marrow-derived mononuclear cells. The primary objective was to determine if autologous cells are safe to administer in children with cerebral palsy. The secondary objectives were to determine if there was improvement in motor function of patients 12 months after infusion using the Gross Motor Function Measure and to evaluate impact of treatment on corticospinal tract microstructure as determined by radial diffusivity measurement. This Phase 1/2a trial was a randomized, blinded, placebo-controlled, crossover study in children aged 2-10 years of age with cerebral palsy enrolled between November 2013 and November 2016. Participants were randomized to 2:1 treatment:placebo. Treatment was either autologous bone-marrow-derived mononuclear cells or autologous umbilical cord blood. All participants who enrolled and completed their baseline visit planned to return for follow-up visits at 6 months, 12 months and 24 months after the baseline visit. At the 12-month post-treatment visit, participants who originally received the placebo received either bone-marrow-derived mononuclear cell or umbilical cord blood treatment. Twenty participants were included; 7 initially randomized to placebo, and 13 randomized to treatment. Five participants randomized to placebo received bone-marrow-derived mononuclear cells, and 2 received umbilical cord blood at the 12-month visit. None of the participants experienced adverse events related to the stem cell infusion. Cell infusion at the doses used in our study did not dramatically alter motor function. We observed concordant bilateral changes in radial diffusivity in 10 of 15 cases where each corticospinal tract could be reconstructed in each hemisphere. In 60% of these cases (6/10), concordant decreases in bilateral corticospinal tract radial diffusivity occurred post-treatment. In addition, 100% of unilateral corticospinal tract cases (3/3) exhibited decreased corticospinal tract radial diffusivity post-treatment. In our discordant cases (n = 5), directionality of changes in corticospinal tract radial diffusivity appeared to coincide with handedness. There was a significant improvement in corticospinal tract radial diffusivity that appears related to handedness. Connectivity strength increased in either or both pathways (corticio-striatal and thalamo-cortical) in each participant at 12 months post-treatment. These data suggest that both stem cell infusions are safe. There may be an improvement in myelination in some groups of patients that correlate with small improvements in the Gross Motor Function Measure scales. A larger autologous cord blood trial is impractical at current rates of blood banking. Either increased private banking or matched units would be required to perform a larger-scale trial.

Methodological evaluation of individual cognitive prediction based on the brain white matter structural connectome

An emerging trend is to use regression‐based machine learning approaches to predict cognitive functions at the individual level from neuroimaging data. However, individual prediction models are inherently influenced by the vast options for network construction and model selection in machine learning pipelines. In particular, the brain white matter (WM) structural connectome lacks a systematic evaluation of the effects of different options in the pipeline on predictive performance. Here, we focused on the methodological evaluation of brain structural connectome‐based predictions. For network construction, we considered two parcellation schemes for defining nodes and seven strategies for defining edges. For the regression algorithms, we used eight regression models. Four cognitive domains and brain age were targeted as predictive tasks based on two independent datasets (Beijing Aging Brain Rejuvenation Initiative [BABRI]: 633 healthy older adults; Human Connectome Projects in Aging [HCP‐A]: 560 healthy older adults). Based on the results, the WM structural connectome provided a satisfying predictive ability for individual age and cognitive functions, especially for executive function and attention. Second, different parcellation schemes induce a significant difference in predictive performance. Third, prediction results from different data sets showed that dMRI with distinct acquisition parameters may plausibly result in a preference for proper fiber reconstruction algorithms and different weighting options. Finally, deep learning and Elastic‐Net models are more accurate and robust in connectome‐based predictions. Together, significant effects of different options in WM network construction and regression algorithms on the predictive performances are identified in this study, which may provide important references and guidelines to select suitable options for future studies in this field.

Comparison of CPU and GPU bayesian estimates of fibre orientations from diffusion MRI

Background

The correct estimation of fibre orientations is a crucial step for reconstructing human brain tracts. Bayesian Estimation of Diffusion Parameters Obtained using Sampling Techniques (bedpostx) is able to estimate several fibre orientations and their diffusion parameters per voxel using Markov Chain Monte Carlo (MCMC) in a whole brain diffusion MRI data, and it is capable of running on GPUs, achieving speed-up of over 100 times compared to CPUs. However, few studies have looked at whether the results from the CPU and GPU algorithms differ. In this study, we compared CPU and GPU bedpostx outputs by running multiple trials of both algorithms on the same whole brain diffusion data and compared each distribution of output using Kolmogorov-Smirnov tests.

Results

We show that distributions of fibre fraction parameters and principal diffusion direction angles from bedpostx and bedpostx_gpu display few statistically significant differences in shape and are localized sparsely throughout the whole brain. Average output differences are small in magnitude compared to underlying uncertainty.

Conclusions

Despite small amount of differences in output between CPU and GPU bedpostx algorithms, results are comparable given the difference in operation order and library usage between CPU and GPU bedpostx.

Canonical Sentence Processing and the Inferior Frontal Cortex: Is There a Connection?

The role of left inferior frontal cortex (LIFC) in canonical sentence comprehension is controversial. Many studies have found involvement of LIFC in sentence production or complex sentence comprehension, but negative or mixed results are often found in comprehension of simple or canonical sentences. We used voxel-, region-, and connectivity-based lesion symptom mapping (VLSM, RLSM, CLSM) in left-hemisphere chronic stroke survivors to investigate canonical sentence comprehension while controlling for lexical-semantic, executive, and phonological processes. We investigated how damage and disrupted white matter connectivity of LIFC and two other language-related regions, the left anterior temporal lobe (LATL) and posterior temporal-inferior parietal area (LpT-iP), affected sentence comprehension. VLSM and RLSM revealed that LIFC damage was not associated with canonical sentence comprehension measured by a sensibility judgment task. LIFC damage was associated instead with impairments in a lexical semantic similarity judgment task with high semantic/executive demands. Damage to the LpT-iP, specifically posterior middle temporal gyrus (pMTG), predicted worse sentence comprehension after controlling for visual lexical access, semantic knowledge, and auditory-verbal short-term memory (STM), but not auditory single-word comprehension, suggesting pMTG is vital for auditory language comprehension. CLSM revealed that disruption of left-lateralized white-matter connections from LIFC to LATL and LpT-iP was associated with worse sentence comprehension, controlling for performance in tasks related to lexical access, auditory word comprehension, and auditory-verbal STM. However, the LIFC connections were accounted for by the lexical semantic similarity judgment task, which had high semantic/executive demands. This suggests that LIFC connectivity is relevant to canonical sentence comprehension when task-related semantic/executive demands are high.

A populational connection distribution map for the whole brain white matter reveals ordered cortical wiring in the space of white matter

The white matter in the brain is composed of neural fibers that wire the cortical and subcortical brain systems. Considering the complexity in terms of interconnections of many neural populations within the narrow space surrounded by the folding walls of the gray matter, the brain requires a certain way of structured wiring. To explore the three-dimensional organization of wiring in an accessible manner, we focused on voxel-level wiring patterns in the white matter according to cortical distributions in which each voxel mediates the wiring. We constructed a voxel-wise connection distribution map from the whole white matter voxels to 360 cortical regions using probabilistic tractography of the 100 diffusion imaging data in the Human Connectome Project. We then explored the spatial organization of the fiber bundles at the white matter voxels in terms of the maximal and second maximal cortical connection labels and the local gradient and entropy of cortical connection density using the population connection distribution map. We presented dominant cortical connection labels, local gradient, and connection entropy for the most representative white matter regions, including the internal capsule, external capsule, corpus callosum, cingulum bundle, and uncinate fascicles, most of which were introduced in the current study. Those major tracts showed a gradient organization of connection distributions for individual voxels. This organized pattern, as reflected in the whole brain connection distribution map, could be established to optimize wiring in the entire brain within the physical space of the white matter.

Right fronto-parietal networks mediate the neurocognitive benefits of enriched environments

Exposure to enriched environments throughout a lifetime, providing so-called reserve, protects against cognitive decline in later years. It has been hypothesized that high levels of alertness necessitated by enriched environments might strengthen the right fronto-parietal networks to facilitate this neurocognitive resilience. We have previously shown that enriched environments offset age-related deficits in selective attention by preserving grey matter within right fronto-parietal regions. Here, using neurite orientation dispersion and density imaging, we examined the relationship between enriched environments, microstructural properties of fronto-parietal white matter association pathways (three branches of the superior longitudinal fasciculus), structural brain health (atrophy), and attention (alertness, orienting and executive control) in a group of older adults. We show that exposure to enriched environments is associated with a lower orientation dispersion index within the right superior longitudinal fasciculus 1 which in turn mediates the relationship between enriched environments and alertness, as well as grey and white matter atrophy. This suggests that enriched environments may induce white matter plasticity (and prevent age-related dispersion of axons) within the right fronto-parietal networks to facilitate the preservation of neurocognitive health in later years.

DNA methylation in relation to gestational age and brain dysmaturation in preterm infants

Preterm birth is associated with dysconnectivity of structural brain networks and is a leading cause of neurocognitive impairment in childhood. Variation in DNA methylation is associated with early exposure to extrauterine life but there has been little research exploring its relationship with brain development. Using genome-wide DNA methylation data from the saliva of 258 neonates, we investigated the impact of gestational age on the methylome and performed functional analysis to identify enriched gene sets from probes that contributed to differentially methylated probes or regions. We tested the hypothesis that variation in DNA methylation could underpin the association between low gestational age at birth and atypical brain development by linking differentially methylated probes with measures of white matter connectivity derived from diffusion MRI metrics: peak width skeletonized mean diffusivity, peak width skeletonized fractional anisotropy and peak width skeletonized neurite density index. Gestational age at birth was associated with widespread differential methylation at term equivalent age, with genome-wide significant associations observed for 8870 CpG probes (P < 3.6 × 10-8) and 1767 differentially methylated regions. Functional analysis identified 14 enriched gene ontology terms pertaining to cell-cell contacts and cell-extracellular matrix contacts. Principal component analysis of probes with genome-wide significance revealed a first principal component that explained 23.5% of the variance in DNA methylation, and this was negatively associated with gestational age at birth. The first principal component was associated with peak width of skeletonized mean diffusivity (β = 0.349, P = 8.37 × 10-10) and peak width skeletonized neurite density index (β = 0.364, P = 4.15 × 10-5), but not with peak width skeletonized fraction anisotropy (β = -0.035, P = 0.510); these relationships mirrored the imaging metrics' associations with gestational age at birth. Low gestational age at birth has a profound and widely distributed effect on the neonatal saliva methylome that is apparent at term equivalent age. Enriched gene ontology terms related to cell-cell contacts reveal pathways that could mediate the effect of early life environmental exposures on development. Finally, associations between differential DNA methylation and image markers of white matter tract microstructure suggest that variation in DNA methylation may provide a link between preterm birth and the dysconnectivity of developing brain networks that characterizes atypical brain development in preterm infants.

Tract-specific statistics based on diffusion-weighted probabilistic tractography

Diffusion-weighted neuroimaging approaches provide rich evidence for estimating the structural integrity of white matter in vivo, but typically do not assess white matter integrity for connections between two specific regions of the brain. Here, we present a method for deriving tract-specific diffusion statistics, based upon predefined regions of interest. Our approach derives a population distribution using probabilistic tractography, based on the Nathan Kline Institute (NKI) Enhanced Rockland sample. We determine the most likely geometry of a path between two regions and express this as a spatial distribution. We then estimate the average orientation of streamlines traversing this path, at discrete distances along its trajectory, and the fraction of diffusion directed along this orientation for each participant. The resulting participant-wise metrics (tract-specific anisotropy; TSA) can then be used for statistical analysis on any comparable population. Based on this method, we report both negative and positive associations between age and TSA for two networks derived from published meta-analytic studies (the “default mode” and “what-where” networks), along with more moderate sex differences and age-by-sex interactions. The proposed method can be applied to any arbitrary set of brain regions, to estimate both the spatial trajectory and DWI-based anisotropy specific to those regions.

Andrew Reid et al. use publicly available data to present a method for deriving tract-specific statistics based on diffusion-weighted MRI, based upon arbitrarily-defined regions of interest. Their approach enables them to report both negative and positive associations between age and tract-specific anisotropy along with more moderate sex differences and age-by-sex interactions.

Electrophysiological abnormalities as indicators of early-stage pathology in Primary Progressive Aphasia (PPA): A case study in semantic variant PPA

Language induced and spontaneous oscillatory activity was measured using MEG in a patient with the semantic variant of Primary Progressive Aphasia (svPPA) and 15 healthy controls.The patient showed oscillatory slowing in the left anterior temporal lobe (ATL) that extended into non-atrophied brain tissue in left and right frontal areas. The white matter connections were reduced to the left and right ATL and left frontal regions, exhibiting electrophysiological abnormalities. Altered diffusion metrics in all four language tracts, indicted compromised white matter integrity. Task-related and spontaneous oscillatory abnormalities can indicate early neurodegeneration in svPPA, providing promising targets for intervention.

DIFFnet: Diffusion Parameter Mapping Network Generalized for Input Diffusion Gradient Schemes and b-Value

In MRI, deep neural networks have been proposed to reconstruct diffusion model parameters. However, the inputs of the networks were designed for a specific diffusion gradient scheme (i.e., diffusion gradient directions and numbers) and a specific b-value that are the same as the training data. In this study, a new deep neural network, referred to as DIFFnet, is developed to function as a generalized reconstruction tool of the diffusion-weighted signals for various gradient schemes and b-values. For generalization, diffusion signals are normalized in a q-space and then projected and quantized, producing a matrix (Qmatrix) as an input for the network. To demonstrate the validity of this approach, DIFFnet is evaluated for diffusion tensor imaging (DIFFnet_DTI) and for neurite orientation dispersion and density imaging (DIFFnet_NODDI). In each model, two datasets with different gradient schemes and b-values are tested. The results demonstrate accurate reconstruction of the diffusion parameters at substantially reduced processing time (approximately 8.7 times and 2240 times faster processing time than conventional methods in DTI and NODDI, respectively; less than 4% mean normalized root-mean-square errors (NRMSE) in DTI and less than 8% in NODDI). The generalization capability of the networks was further validated using reduced numbers of diffusion signals from the datasets and a public dataset from Human Connection Project. Different from previously proposed deep neural networks, DIFFnet does not require any specific gradient scheme and b-value for its input. As a result, it can be adopted as an online reconstruction tool for various complex diffusion imaging.

Coupling cognitive and brainstem dysfunction in multiple sclerosis-related chronic neuropathic limb pain

Chronic pain in multiple sclerosis is common and difficult to treat. Its mechanisms remain incompletely understood. Dysfunction of the descending pain modulatory system is known to contribute to human chronic pain conditions. However, it is not clear how alterations in executive function influence this network, despite healthy volunteer studies linking function of the descending pain modulatory system, to cognition. In adults with multiple sclerosis-associated chronic neuropathic limb pain, compared to those without pain, we hypothesized altered functional connectivity of the descending pain modulatory system, coupled to executive dysfunction. Specifically we hypothesized reduced mental flexibility, because of potential importance in stimulus reappraisal. To investigate these hypotheses, we conducted a case-control cross-sectional study of 47 adults with relapsing remitting multiple sclerosis (31 with chronic neuropathic limb pain, 16 without pain), employing clinical, neuropsychological, structural, and functional MRI measures. We measured brain lesions and atrophy affecting descending pain modulatory system structures. Both cognitive and affective dysfunctions were confirmed in the chronic neuropathic limb pain group, including reduced mental flexibility (Delis Kaplan Executive Function System card sorting tests P < 0.001). Functional connectivity of rostral anterior cingulate and ventrolateral periaqueductal gray, key structures of the descending pain modulatory system, was significantly lower in the group experiencing chronic neuropathic pain. There was no significant between-group difference in whole-brain grey matter or lesion volumes, nor lesion volume affecting white matter tracts between rostral anterior cingulate and periaqueductal gray. Brainstem-specific lesion volume was higher in the chronic neuropathic limb pain group (P = 0.0017). Differential functional connectivity remained after correction for brainstem-specific lesion volume. Gabapentinoid medications were more frequently used in the chronic pain group. We describe executive dysfunction in people with multiple sclerosis affected by chronic neuropathic pain, along with functional and structural MRI evidence compatible with dysfunction of the descending pain modulatory system. These findings extend understanding of close inter-relationships between cognition, function of the descending pain modulatory system, and chronic pain, both in multiple sclerosis and more generally in human chronic pain conditions. These findings could support application of pharmacological and cognitive interventions in chronic neuropathic pain associated with multiple sclerosis.

Differential structure-function network coupling in the inattentive and combined types of attention deficit hyperactivity disorder

The heterogeneous presentation of inattentive and hyperactive-impulsive core symptoms in attention deficit hyperactivity disorder (ADHD) warrants further investigation into brain network connectivity as a basis for subtype divisions in this prevalent disorder. With diffusion and resting-state functional magnetic resonance imaging data from the Healthy Brain Network database, we analyzed both structural and functional network efficiency and structure-functional network (SC-FC) coupling at the default mode (DMN), executive control (ECN), and salience (SAN) intrinsic networks in 201 children diagnosed with the inattentive subtype (ADHD-I), the combined subtype (ADHD-C), and typically developing children (TDC) to characterize ADHD symptoms relative to TDC and to test differences between ADHD subtypes. Relative to TDC, children with ADHD had lower structural connectivity and network efficiency in the DMN, without significant group differences in functional networks. Children with ADHD-C had higher SC-FC coupling, a finding consistent with diminished cognitive flexibility, for all subnetworks compared to TDC. The ADHD-C group also demonstrated increased SC-FC coupling in the DMN compared to the ADHD-I group. The correlation between SC-FC coupling and hyperactivity scores was negative in the ADHD-I, but not in the ADHD-C group. The current study suggests that ADHD-C and ADHD-I may differ with respect to their underlying neuronal connectivity and that the added dimensionality of hyperactivity may not explain this distinction.

MASiVar: Multisite, multiscanner, and multisubject acquisitions for studying variability in diffusion weighted MRI

PURPOSE

Diffusion-weighted imaging allows investigators to identify structural, microstructural, and connectivity-based differences between subjects, but variability due to session and scanner biases is a challenge.

METHODS

To investigate DWI variability, we present MASiVar, a multisite data set consisting of 319 diffusion scans acquired at 3 T from b = 1000 to 3000 s/mm2 across 14 healthy adults, 83 healthy children (5 to 8 years), three sites, and four scanners as a publicly available, preprocessed, and de-identified data set. With the adult data, we demonstrate the capacity of MASiVar to simultaneously quantify the intrasession, intersession, interscanner, and intersubject variability of four common DWI processing approaches: (1) a tensor signal representation, (2) a multi-compartment neurite orientation dispersion and density model, (3) white-matter bundle segmentation, and (4) structural connectomics. Respectively, we evaluate region-wise fractional anisotropy, mean diffusivity, and principal eigenvector; region-wise CSF volume fraction, intracellular volume fraction, and orientation dispersion index; bundle-wise shape, volume, fractional anisotropy, and length; and whole connectome correlation and maximized modularity, global efficiency, and characteristic path length.

RESULTS

We plot the variability in these measures at each level and find that it consistently increases with intrasession to intersession to interscanner to intersubject effects across all processing approaches and that sometimes interscanner variability can approach intersubject variability.

CONCLUSIONS

This study demonstrates the potential of MASiVar to more globally investigate DWI variability across multiple levels and processing approaches simultaneously and suggests harmonization between scanners for multisite analyses should be considered before inference of group differences on subjects.

A critical review of connectome validation studies

Diffusion MRI tractography is the most widely used macroscale method for mapping connectomes in vivo. However, tractography is prone to various errors and biases, and thus tractography-derived connectomes require careful validation. Here, we critically review studies that have developed or utilized phantoms and tracer maps to validate tractography-derived connectomes, either quantitatively or qualitatively. We identify key factors impacting connectome reconstruction accuracy, including streamline seeding, propagation and filtering methods, and consider the strengths and limitations of state-of-the-art connectome phantoms and associated validation studies. These studies demonstrate the inherent limitations of current fiber orientation models and tractography algorithms and their impact on connectome reconstruction accuracy. Reconstructing connectomes with both high sensitivity and high specificity is challenging, given that some tractography methods can generate an abundance of spurious connections, while others can overlook genuine fiber bundles. We argue that streamline filtering can minimize spurious connections and potentially improve the biological plausibility of connectomes derived from tractography. We find that algorithmic choices such as the tractography seeding methodology, angular threshold, and streamline propagation method can substantially impact connectome reconstruction accuracy. Hence, careful application of tractography is necessary to reconstruct accurate connectomes. Improvements in diffusion MRI acquisition techniques will not necessarily overcome current tractography limitations without accompanying modeling and algorithmic advances.

The Human Connectome Project: A retrospective

The Human Connectome Project (HCP) was launched in 2010 as an ambitious effort to accelerate advances in human neuroimaging, particularly for measures of brain connectivity; apply these advances to study a large number of healthy young adults; and freely share the data and tools with the scientific community. NIH awarded grants to two consortia; this retrospective focuses on the "WU-Minn-Ox" HCP consortium centered at Washington University, the University of Minnesota, and University of Oxford. In just over 6 years, the WU-Minn-Ox consortium succeeded in its core objectives by: 1) improving MR scanner hardware, pulse sequence design, and image reconstruction methods, 2) acquiring and analyzing multimodal MRI and MEG data of unprecedented quality together with behavioral measures from more than 1100 HCP participants, and 3) freely sharing the data (via the ConnectomeDB database) and associated analysis and visualization tools. To date, more than 27 Petabytes of data have been shared, and 1538 papers acknowledging HCP data use have been published. The "HCP-style" neuroimaging paradigm has emerged as a set of best-practice strategies for optimizing data acquisition and analysis. This article reviews the history of the HCP, including comments on key events and decisions associated with major project components. We discuss several scientific advances using HCP data, including improved cortical parcellations, analyses of connectivity based on functional and diffusion MRI, and analyses of brain-behavior relationships. We also touch upon our efforts to develop and share a variety of associated data processing and analysis tools along with detailed documentation, tutorials, and an educational course to train the next generation of neuroimagers. We conclude with a look forward at opportunities and challenges facing the human neuroimaging field from the perspective of the HCP consortium.

Identifying brain regions supporting amygdalar functionality: Application of a novel graph theory technique

Effective amygdalar functionality depends on the concerted activity of a complex network of regions. Thus, the role of the amygdala cannot be fully understood without identifying the set of brain structures that allow the processes performed by the amygdala to emerge. However, this identification has yet to occur, hampering our ability to understand both normative and pathological processes that rely on the amygdala. We developed and applied novel graph theory methods to diffusion-based anatomical networks in a large sample (n = 1,052, 54.28% female, mean age=28.75) to identify nodes that critically support amygdalar interactions with the larger brain network. We examined three graph properties, each indexing a different emergent aspect of amygdalar network communication: current-flow betweenness centrality (amygdalar influence on information flowing between other pairs of nodes), node communicability (clarity of communication between the amygdala and other nodes), and subgraph centrality (amygdalar influence over local network processing). Findings demonstrate that each of these aspects of amygdalar communication is associated with separable sets of regions and, in some cases, these sets map onto previously identified sub-circuits. For example, betweenness and communicability were each associated with different sub-circuits that have been identified in previous work as supporting distinct aspects of memory-guided behavior. Other regions identified span basic (e.g., visual cortex) to higher-order (e.g., insula) sensory processing and executive functions (e.g., dorsolateral prefrontal cortex). Present findings expand our current understanding of amygdalar function by showing that there is no single 'amygdala network', but rather multiple networks, each supporting different modes of amygdalar interaction with the larger brain network. Additionally, our novel method allowed for the identification of how such regions support the amygdala, which has not been previously explored.

Diffusion MRI data, sulcal anatomy, and tractography for eight species from the Primate Brain Bank

Large-scale comparative neuroscience requires data from many species and, ideally, at multiple levels of description. Here, we contribute to this endeavor by presenting diffusion and structural MRI data from eight primate species that have not or rarely been described in the literature. The selected samples from the Primate Brain Bank cover a prosimian, New and Old World monkeys, and a great ape. We present preliminary labelling of the cortical sulci and tractography of the optic radiation, dorsal part of the cingulum bundle, and dorsal parietal-frontal and ventral temporal-frontal longitudinal white matter tracts. Both dorsal and ventral association fiber systems could be observed in all samples, with the dorsal tracts occupying much less relative volume in the prosimian than in other species. We discuss the results in the context of known primate specializations and present hypotheses for further research. All data and results presented here are available online as a resource for the scientific community.

Fiber Clustering Acceleration With a Modified Kmeans++ Algorithm Using Data Parallelism

Fiber clustering methods are typically used in brain research to study the organization of white matter bundles from large diffusion MRI tractography datasets. These methods enable exploratory bundle inspection using visualization and other methods that require identifying brain white matter structures in individuals or a population. Some applications, such as real-time visualization and inter-subject clustering, need fast and high-quality intra-subject clustering algorithms. This work proposes a parallel algorithm using a General Purpose Graphics Processing Unit (GPGPU) for fiber clustering based on the FFClust algorithm. The proposed GPGPU implementation exploits data parallelism using both multicore and GPU fine-grained parallelism present in commodity architectures, including current laptops and desktop computers. Our approach implements all FFClust steps in parallel, improving execution times in all of them. In addition, our parallel approach includes a parallel Kmeans++ algorithm implementation and defines a new variant of Kmeans++ to reduce the impact of choosing outliers as initial centroids. The results show that our approach provides clustering quality results very similar to FFClust, and it requires an execution time of 3.5 s for processing about a million fibers, achieving a speedup of 11.5 times compared to FFClust.

Ventralis intermedius nucleus anatomical variability assessment by MRI structural connectivity

The ventralis intermedius nucleus (Vim) is centrally placed in the dentato-thalamo-cortical pathway (DTCp) and is a key surgical target in the treatment of severe medically refractory tremor. It is not visible on conventional MRI sequences; consequently, stereotactic targeting currently relies on atlas-based coordinates. This fails to capture individual anatomical variability, which may lead to poor long-term clinical efficacy. Probabilistic tractography, combined with known anatomical connectivity, enables localisation of thalamic nuclei at an individual subject level. There are, however, a number of confounds associated with this technique that may influence results. Here we focused on an established method, using probabilistic tractography to reconstruct the DTCp, to identify the connectivity-defined Vim (cd-Vim) in vivo. Using 100 healthy individuals from the Human Connectome Project, our aim was to quantify cd-Vim variability across this population, measure the discrepancy with atlas-defined Vim (ad-Vim), and assess the influence of potential methodological confounds. We found no significant effect of any of the confounds. The mean cd-Vim coordinate was located within 1.88 mm (left) and 2.12 mm (right) of the average midpoint and 3.98 mm (left) and 5.41 mm (right) from the ad-Vim coordinates. cd-Vim location was more variable on the right, which reflects hemispheric asymmetries in the probabilistic DTC reconstructed. The method was reproducible, with no significant cd-Vim location differences in a separate test-retest cohort. The superior cerebellar peduncle was identified as a potential source of artificial variance. This work demonstrates significant individual anatomical variability of the cd-Vim that atlas-based coordinate targeting fails to capture. This variability was not related to any methodological confound tested. Lateralisation of cerebellar functions, such as speech, may contribute to the observed asymmetry. Tractography-based methods seem sensitive to individual anatomical variability that is missed by conventional neurosurgical targeting; these findings may form the basis for translational tools to improve efficacy and reduce side-effects of thalamic surgery for tremor.

The autonomic brain: Multi-dimensional generative hierarchical modelling of the autonomic connectome

The autonomic nervous system governs the body's multifaceted internal adaptation to diverse changes in the external environment, a role more complex than is accessible to the methods-and data scales-hitherto used to illuminate its operation. Here we apply generative graphical modelling to large-scale multimodal neuroimaging data encompassing normal and abnormal states to derive a comprehensive hierarchical representation of the autonomic brain. We demonstrate that whereas conventional structural and functional maps identify regions jointly modulated by parasympathetic and sympathetic systems, only graphical analysis discriminates between them, revealing the cardinal roles of the autonomic system to be mediated by high-level distributed interactions. We provide a novel representation of the autonomic system-a multidimensional, generative network-that renders its richness tractable within future models of its function in health and disease.

Minimal specifications for non-human primate MRI: Challenges in standardizing and harmonizing data collection

Recent methodological advances in MRI have enabled substantial growth in neuroimaging studies of non-human primates (NHPs), while open data-sharing through the PRIME-DE initiative has increased the availability of NHP MRI data and the need for robust multi-subject multi-center analyses. Streamlined acquisition and analysis protocols would accelerate and improve these efforts. However, consensus on minimal standards for data acquisition protocols and analysis pipelines for NHP imaging remains to be established, particularly for multi-center studies. Here, we draw parallels between NHP and human neuroimaging and provide minimal guidelines for harmonizing and standardizing data acquisition. We advocate robust translation of widely used open-access toolkits that are well established for analyzing human data. We also encourage the use of validated, automated pre-processing tools for analyzing NHP data sets. These guidelines aim to refine methodological and analytical strategies for small and large-scale NHP neuroimaging data. This will improve reproducibility of results, and accelerate the convergence between NHP and human neuroimaging strategies which will ultimately benefit fundamental and translational brain science.

Empty Sella Syndrome as a Window Into the Neuroprotective Effects of Prolactin

Background: The goal of this study was to relate diffusion MR measures of white matter integrity of the retinofugal visual pathway with prolactin levels in a patient with downward herniation of the optic chiasm secondary to medical treatment of a prolactinoma.

Methods: A 36-year-old woman with a prolactinoma presented with progressive bilateral visual field defects 9 years after initial diagnosis and medical treatment. She was diagnosed with empty-sella syndrome and instructed to stop cabergoline. Hormone testing was conducted in tandem with routine clinical evaluations over 1 year and the patient was followed with diffusion magnetic resonance imaging (dMRI), optical coherence tomography (OCT), and automated perimetry at three time points. Five healthy controls underwent a complementary battery of clinical and neuroimaging tests at a single time point.

Results: Shortly after discontinuing cabergoline, diffusion metrics in the optic tracts were within the range of values observed in healthy controls. However, following a brief period where the patient resumed cabergoline (of her own volition), there was a decrease in serum prolactin with a corresponding decrease in visual ability and increase in radial diffusivity (p < 0.001). Those measures again returned to their baseline ranges after discontinuing cabergoline a second time.

Conclusions: These results demonstrate the sensitivity of dMRI to detect rapid and functionally significant microstructural changes in white matter tracts secondary to alterations in serum prolactin levels. The inverse relations between prolactin and measures of white matter integrity and visual function are consistent with the hypothesis that prolactin can play a neuroprotective role in the injured nervous system.

Supplementary and Premotor Aspects of the Corticospinal Tract Show Links with Restricted and Repetitive Behaviors in Middle-Aged Adults with Autism Spectrum Disorder

Individuals with autism spectrum disorder (ASD) show motor impairment into adulthood and risk decline during aging, but little is known about brain changes in aging adults with ASD. Few studies of ASD have directly examined the corticospinal tract (CST)-the major descending pathway in the brain responsible for voluntary motor behavior-outside its primary motor (M1) connections. In 26 middle-aged adults with ASD and 26 age-matched typical comparison participants, we used diffusion imaging to examine the microstructure and volume of CST projections from M1, dorsal premotor (PMd), supplementary motor area (SMA), and primary somatosensory (S1) cortices with respect to age. We also examined relationships between each CST sub-tract (-cst), motor skills, and autism symptoms. We detected no significant group or age-related differences in tracts extending from M1 or other areas. However, sub-tracts of the CST extending from secondary (but not primary) motor areas were associated with core autism traits. Increased microstructural integrity of left PMd-cst and SMA-cst were associated with less-severe restricted and repetitive behaviors (RRB) in the ASD group. These findings suggest that secondary motor cortical areas, known to be involved in selecting motor programs, may be implicated in cognitive motor processes underlying RRB in ASD.