Fusion in diffusion MRI for improved fibre orientation estimation: An application to the 3T and 7T data of the Human Connectome Project

Cross-species neuroanatomy in primates using tractography

Due to their integrative role in brain function, long-range white matter connections exhibit high individual variability, giving rise to personalised brain circuits. This neurovariability is more evident in the connection patterns of brain areas that have evolved more recently. Diffusion MRI tractography allows unique opportunities for comparative neuroanatomy across species to study evolution and provide unique insights into the phylogeny of brain networks, which we overview in this note, inspired by discussions at the International Society for Tractography (IST) retreat.

New method for diffusion-weighted images denoising based on patch-matching with higher-order singular value decomposition

BackgroundDiffusion-weighted imaging (DWI) is an important technique to study brain microstructure. However, diffusion-weighted (DW) images suffer from severe low signal-to-noise ratio (SNR) problem, affecting subsequent diffusion analysis.ObjectiveThe goal of this paper is to develop advanced DWI denoising technique to effectively reduce noise while improving the accuracy and reliability of subsequent diffusion model fitting and diffusion analysis, thereby facilitating the research and analysis of brain science.MethodsWe propose a new method for denoising DW images based on patch-matching with higher-order singular value decomposition (HOSVD) by combined with the variance-stabilizing transformation technique. It starts with introducing a novel non-local mean algorithm as a prefiltering stage, and then denoises the noisy data using a local HOSVD algorithm based on the HOSVD bases learned from prefiltered images.ResultsExperiments are performed on simulation, HCP and in vivo brain DWI datasets. Results show that the proposed method significantly reduces spatially invariant and variant noise, improving the most reliable diffusion analysis compared with the different denoising methods.ConclusionsThe proposed method achieves state-of-the-art performance which can improve image quality and enable accurate diffusion analysis.

Imaging the cerebral vasculature at different scales: translational tools to investigate the neurovascular interfaces

The improvements in imaging technology opened up the possibility to investigate the structure and function of cerebral vasculature and the neurovascular unit with unprecedented precision and gaining deep insights not only on the morphology of the vessels but also regarding their function and regulation related to the cerebral activity. In this review, we will dissect the different imaging capabilities regarding the cerebrovascular tree, the neurovascular unit, the haemodynamic response function, and thus, the vascular-neuronal coupling. We will discuss both clinical and preclinical setting, with a final discussion on the current scenery in cerebrovascular imaging where magnetic resonance imaging and multimodal microscopy emerge as the most potent and versatile tools, respectively, in the clinical and preclinical context.

Methodological considerations on diffusion MRI tractography in infants aged 0-2 years: a scoping review

Diffusion MRI (dMRI) enables studying the complex architectural organization of the brain's white matter (WM) through virtual reconstruction of WM fiber tracts (tractography). Despite the anticipated clinical importance of applying tractography to study structural connectivity and tract development during the critical period of rapid infant brain maturation, detailed descriptions on how to approach tractography in young infants are limited. Over the past two decades, tractography from infant dMRI has mainly been applied in research settings and focused on diffusion tensor imaging (DTI). Only few studies used techniques superior to DTI in terms of disentangling information on the brain's organizational complexity, including crossing fibers. While more advanced techniques may enhance our understanding of the intricate processes of normal and abnormal brain development and extensive knowledge has been gained from application on adult scans, their applicability in infants has remained underexplored. This may partially be due to the higher technical requirements versus the need to limit scan time in young infants. We review various previously described methodological practices for tractography in the infant brain (0-2 years-of-age) and provide recommendations to optimize advanced tractography approaches to enable more accurate reconstructions of the brain WM's complexity. IMPACT: Diffusion tensor imaging is the technique most frequently used for fiber tracking in the developing infant brain but is limited in capability to disentangle the complex white matter organization. Advanced tractography techniques allow for reconstruction of crossing fiber bundles to better reflect the brain's complex organization. Yet, they pose practical and technical challenges in the fast developing young infant's brain. Methods on how to approach advanced tractography in the young infant's brain have hardly been described. Based on a literature review, recommendations are provided to optimize tractography for the developing infant brain, aiming to advance early diagnosis and neuroprotective strategies.

A gradient of hemisphere-specific dorsal to ventral processing routes in parieto-premotor networks

Networks in the parietal and premotor cortices enable essential human abilities regarding motor processing, including attention and tool use. Even though our knowledge on its topography has steadily increased, a detailed picture of hemisphere-specific integrating pathways is still lacking. With the help of multishell diffusion magnetic resonance imaging, probabilistic tractography, and the Graph Theory Analysis, we investigated connectivity patterns between frontal premotor and posterior parietal brain areas in healthy individuals. With a two-stage node characterization approach, we defined the network role of precisely mapped cortical regions from the Julich-Brain atlas. We found evidence for a third, left-sided, medio-dorsal subpathway in a successively graded dorsal stream, referencing more specialized motor processing on the left. Supplementary motor areas had a strongly lateralized connectivity to either left dorsal or right ventral parietal domains, representing an action-attention dichotomy between hemispheres. The left sulcal parietal regions primarily coupled with areas 44 and 45, mirrored by the inferior frontal junction (IFJ) on the right, a structural lateralization we termed as "Broca's-IFJ switch." We were able to deepen knowledge on gyral and sulcal pathways as well as domain-specific contributions in parieto-premotor networks. Our study sheds new light on the complex lateralization of cortical routes for motor activity in the human brain.

Whole brain functional connectivity: Insights from next generation neural mass modelling incorporating electrical synapses

The ready availability of brain connectome data has both inspired and facilitated the modelling of whole brain activity using networks of phenomenological neural mass models that can incorporate both interaction strength and tract length between brain regions. Recently, a new class of neural mass model has been developed from an exact mean field reduction of a network of spiking cortical cell models with a biophysically realistic model of the chemical synapse. Moreover, this new population dynamics model can naturally incorporate electrical synapses. Here we demonstrate the ability of this new modelling framework, when combined with data from the Human Connectome Project, to generate patterns of functional connectivity (FC) of the type observed in both magnetoencephalography and functional magnetic resonance neuroimaging. Some limited explanatory power is obtained via an eigenmode description of frequency-specific FC patterns, obtained via a linear stability analysis of the network steady state in the neigbourhood of a Hopf bifurcation. However, direct numerical simulations show that empirical data is more faithfully recapitulated in the nonlinear regime, and exposes a key role of gap junction coupling strength in generating empirically-observed neural activity, and associated FC patterns and their evolution. Thereby, we emphasise the importance of maintaining known links with biological reality when developing multi-scale models of brain dynamics. As a tool for the study of dynamic whole brain models of the type presented here we further provide a suite of C++ codes for the efficient, and user friendly, simulation of neural mass networks with multiple delayed interactions.

A combined DTI-fMRI approach for optimizing the delineation of posteromedial versus anterolateral entorhinal cortex

In the entorhinal cortex (EC), attempts have been made to identify the human homologue regions of the medial (MEC) and lateral (LEC) subregions using either functional magnetic resonance imaging (fMRI) or diffusion tensor imaging (DTI). However, there are still discrepancies between entorhinal subdivisions depending on the choice of connectivity seed regions and the imaging modality used. While DTI can be used to follow the white matter tracts of the brain, fMRI can identify functionally connected brain regions. In this study, we used both DTI and resting-state fMRI in 103 healthy adults to investigate both structural and functional connectivity between the EC and associated cortical brain regions. Differential connectivity with these regions was then used to predict the locations of the human homologues of MEC and LEC. Our results from combining DTI and fMRI support a subdivision into posteromedial (pmEC) and anterolateral (alEC) EC and reveal a confined border between the pmEC and alEC. Furthermore, the EC subregions obtained by either imaging modality showed similar distinct whole-brain connectivity profiles. Optimizing the delineation of the human homologues of MEC and LEC with a combined, cross-validated DTI-fMRI approach allows to define a likely border between the two subdivisions and has implications for both cognitive and translational neuroscience research.

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

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

Main Paths of Brain Fibers in Diffusion Images Mixed with a Noise to Improve Performance of Tractography Algorithm-Evaluation in Phantom

Background

Some voxels may alter the tractography results due to unintentional alteration of noises and other unwanted factors.

Objective

This study aimed to investigate the effect of local phase features on tractography results providing data are mixed by a Gaussian or random distribution noise.

Material and Methods

In this simulation study, a mask was firstly designed based on the local phase features to decrease false-negative and -positive tractography results. The local phase features are calculated according to the local structures of images, which can be zero-dimensional, meaning just one point (equivalent to noise in tractography algorithm), a line (equivalent to a simple fiber), or an edge (equivalent to structures more complex than a simple fiber). A digital phantom evaluated the feasibility current model with the maximum complexities of configurations in fibers, including crossing fibers. In this paper, the diffusion images were mixed separately by a Gaussian or random distribution noise in 2 forms a zero-mean noise and a noise with a mean of data.

Results

The local mask eliminates the pixels of unfitted values with the main structures of images, due to noise or other interferer factors.

Conclusion

The local phase features of diffusion images are an innovative solution to determine principal diffusion directions.

DIMOND: DIffusion Model OptimizatioN with Deep Learning

Diffusion magnetic resonance imaging is an important tool for mapping tissue microstructure and structural connectivity non-invasively in the in vivo human brain. Numerous diffusion signal models are proposed to quantify microstructural properties. Nonetheless, accurate estimation of model parameters is computationally expensive and impeded by image noise. Supervised deep learning-based estimation approaches exhibit efficiency and superior performance but require additional training data and may be not generalizable. A new DIffusion Model OptimizatioN framework using physics-informed and self-supervised Deep learning entitled "DIMOND" is proposed to address this problem. DIMOND employs a neural network to map input image data to model parameters and optimizes the network by minimizing the difference between the input acquired data and synthetic data generated via the diffusion model parametrized by network outputs. DIMOND produces accurate diffusion tensor imaging results and is generalizable across subjects and datasets. Moreover, DIMOND outperforms conventional methods for fitting sophisticated microstructural models including the kurtosis and NODDI model. Importantly, DIMOND reduces NODDI model fitting time from hours to minutes, or seconds by leveraging transfer learning. In summary, the self-supervised manner, high efficacy, and efficiency of DIMOND increase the practical feasibility and adoption of microstructure and connectivity mapping in clinical and neuroscientific applications.

DiMANI: diffusion MRI for anatomical nuclei imaging-Application for the direct visualization of thalamic subnuclei

The thalamus is a centrally located and heterogeneous brain structure that plays a critical role in various sensory, motor, and cognitive processes. However, visualizing the individual subnuclei of the thalamus using conventional MRI techniques is challenging. This difficulty has posed obstacles in targeting specific subnuclei for clinical interventions such as deep brain stimulation (DBS). In this paper, we present DiMANI, a novel method for directly visualizing the thalamic subnuclei using diffusion MRI (dMRI). The DiMANI contrast is computed by averaging, voxelwise, diffusion-weighted volumes enabling the direct distinction of thalamic subnuclei in individuals. We evaluated the reproducibility of DiMANI through multiple approaches. First, we utilized a unique dataset comprising 8 scans of a single participant collected over a 3-year period. Secondly, we quantitatively assessed manual segmentations of thalamic subnuclei for both intra-rater and inter-rater reliability. Thirdly, we qualitatively correlated DiMANI imaging data from several patients with Essential Tremor with the localization of implanted DBS electrodes and clinical observations. Lastly, we demonstrated that DiMANI can provide similar features at 3T and 7T MRI, using varying numbers of diffusion directions. Our results establish that DiMANI is a reproducible and clinically relevant method to directly visualize thalamic subnuclei. This has significant implications for the development of new DBS targets and the optimization of DBS therapy.

versaFlow: a versatile pipeline for resolution adapted diffusion MRI processing and its application to studying the variability of the PRIME-DE database

The lack of "gold standards" in Diffusion Weighted Imaging (DWI) makes validation cumbersome. To tackle this task, studies use translational analysis where results in humans are benchmarked against findings in other species. Non-Human Primates (NHP) are particularly interesting for this, as their cytoarchitecture is closely related to humans. However, tools used for processing and analysis must be adapted and finely tuned to work well on NHP images. Here, we propose versaFlow, a modular pipeline implemented in Nextflow, designed for robustness and scalability. The pipeline is tailored to in vivo NHP DWI at any spatial resolution; it allows for maintainability and customization. Processes and workflows are implemented using cutting-edge and state-of-the-art Magnetic Resonance Imaging (MRI) processing technologies and diffusion modeling algorithms, namely Diffusion Tensor Imaging (DTI), Constrained Spherical Deconvolution (CSD), and DIstribution of Anisotropic MicrOstructural eNvironments in Diffusion-compartment imaging (DIAMOND). Using versaFlow, we provide an in-depth study of the variability of diffusion metrics computed on 32 subjects from 3 sites of the Primate Data Exchange (PRIME-DE), which contains anatomical T1-weighted (T1w) and T2-weighted (T2w) images, functional MRI (fMRI), and DWI of NHP brains. This dataset includes images acquired over a range of resolutions, using single and multi-shell gradient samplings, on multiple scanner vendors. We perform a reproducibility study of the processing of versaFlow using the Aix-Marseilles site's data, to ensure that our implementation has minimal impact on the variability observed in subsequent analyses. We report very high reproducibility for the majority of metrics; only gamma distribution parameters of DIAMOND display less reproducible behaviors, due to the absence of a mechanism to enforce a random number seed in the software we used. This should be taken into consideration when future applications are performed. We show that the PRIME-DE diffusion data exhibits a great level of variability, similar or greater than results obtained in human studies. Its usage should be done carefully to prevent instilling uncertainty in statistical analyses. This hints at a need for sufficient harmonization in acquisition protocols and for the development of robust algorithms capable of managing the variability induced in imaging due to differences in scanner models and/or vendors.

Factors predicting the transition from acute to persistent pain in people with 'sciatica': the FORECAST longitudinal prognostic factor cohort study protocol

INTRODUCTION

Sciatica is a common condition and is associated with higher levels of pain, disability, poorer quality of life, and increased use of health resources compared with low back pain alone. Although many patients recover, a third develop persistent sciatica symptoms. It remains unclear, why some patients develop persistent sciatica as none of the traditionally considered clinical parameters (eg, symptom severity, routine MRI) are consistent prognostic factors.The FORECAST study (factors predicting the transition from acute to persistent pain in people with 'sciatica') will take a different approach by exploring mechanism-based subgroups in patients with sciatica and investigate whether a mechanism-based approach can identify factors that predict pain persistence in patients with sciatica.

METHODS AND ANALYSIS

We will perform a prospective longitudinal cohort study including 180 people with acute/subacute sciatica. N=168 healthy participants will provide normative data. A detailed set of variables will be assessed within 3 months after sciatica onset. This will include self-reported sensory and psychosocial profiles, quantitative sensory testing, blood inflammatory markers and advanced neuroimaging. We will determine outcome with the Sciatica Bothersomeness Index and a Numerical Pain Rating Scale for leg pain severity at 3 and 12 months.We will use principal component analysis followed by clustering methods to identify subgroups. Univariate associations and machine learning methods optimised for high dimensional small data sets will be used to identify the most powerful predictors and model selection/accuracy.The results will provide crucial information about the pathophysiological drivers of sciatica symptoms and may identify prognostic factors of pain persistence.

ETHICS AND DISSEMINATION

The FORECAST study has received ethical approval (South Central Oxford C, 18/SC/0263). The dissemination strategy will be guided by our patient and public engagement activities and will include peer-reviewed publications, conference presentations, social media and podcasts.

TRIAL REGISTRATION NUMBER

ISRCTN18170726; Pre-results.

TrGANet: Transforming 3T to 7T dMRI using Trapezoidal Rule and Graph based Attention Modules

Diffusion MRI (dMRI) is a non-invasive tool for assessing the white matter region of the brain by approximating the fiber streamlines, structural connectivity, and estimation of microstructure. This modality can yield useful information for diagnosing several mental diseases as well as for surgical planning. The higher angular resolution diffusion imaging (HARDI) technique is helpful in obtaining more robust fiber tracts by getting a good approximation of regions where fibers cross. Moreover, HARDI is more sensitive to tissue changes and can accurately represent anatomical details in the human brain at higher magnetic strengths. In other words, magnetic strengths affect the quality of the image, and hence high magnetic strength has good tissue contrast with better spatial resolution. However, a higher magnetic strength scanner (like 7T) is costly and unaffordable to most hospitals. Hence, in this work, we have proposed a novel CNN architecture for the transformation of 3T to 7T dMRI. Additionally, we have also reconstructed the multi-shell multi-tissue fiber orientation distribution function (MSMT fODF) at 7T from single-shell 3T. The proposed architecture consists of a CNN-based ODE solver utilizing the Trapezoidal rule and graph-based attention layer alongwith L₁ and total variation loss. Finally, the model has been validated on the HCP data set quantitatively and qualitatively.

Diffusion Measures of Subcortical Structures Using High-Field MRI

The pathology of Parkinson’s disease (PD) involves the death of dopaminergic neurons in the substantia nigra (SN), which slowly influences downstream basal ganglia pathways as dopamine transport diminishes. Diffusion magnetic resonance imaging (MRI) has been used to diagnose PD by assessing white matter connectivity in some brain areas. For this study, we applied Lead-DBS to human connectome project data to automatically segment 11 subcortical structures of 49 human connectome project subjects, reducing the reliance on manual segmentation for more consistency. The Lead-connectome pipeline, which utilizes DSI Studio to generate structural connectomes from each 3T and 7T diffusion image, was applied to 3T and 7T data to investigate possible differences in diffusion measures due to different acquisition protocols. Significantly higher fractional anisotropy (FA) values were found in the 3T left SN; significantly higher MD values were found in the 3T left SN and the right amygdala, SN, and subthalamic nucleus (STN); significantly higher AD values were found in the right RN and STN; and significantly higher RD values were found in the left RN and right amygdala. We illustrate a methodology for obtaining diffusion measures of basal ganglia and basal ganglia connectivity using diffusion images, as well as show possible differences in diffusion measures that can arise due to the differences in MRI acquisitions.

Ultra-high field (10.5T) diffusion-weighted MRI of the macaque brain

Diffu0sion-weighted magnetic resonance imaging (dMRI) is a non-invasive imaging technique that provides information about the barriers to the diffusion of water molecules in tissue. In the brain, this information can be used in several important ways, including to examine tissue abnormalities associated with brain disorders and to infer anatomical connectivity and the organization of white matter bundles through the use of tractography algorithms. However, dMRI also presents certain challenges. For example, historically, the biological validation of tractography models has shown only moderate correlations with anatomical connectivity as determined through invasive tract-tracing studies. Some of the factors contributing to such issues are low spatial resolution, low signal-to-noise ratios, and long scan times required for high-quality data, along with modeling challenges like complex fiber crossing patterns. Leveraging the capabilities provided by an ultra-high field scanner combined with denoising, we have acquired whole-brain, 0.58 mm isotropic resolution dMRI with a 2D-single shot echo planar imaging sequence on a 10.5 Tesla scanner in anesthetized macaques. These data produced high-quality tractograms and maps of scalar diffusion metrics in white matter. This work demonstrates the feasibility and motivation for in-vivo dMRI studies seeking to benefit from ultra-high fields.

Accuracy and reliability of diffusion imaging models

Diffusion imaging aims to non-invasively characterize the anatomy and integrity of the brain's white matter fibers. We evaluated the accuracy and reliability of commonly used diffusion imaging methods as a function of data quantity and analysis method, using both simulations and highly sampled individual-specific data (927-1442 diffusion weighted images [DWIs] per individual). Diffusion imaging methods that allow for crossing fibers (FSL's BedpostX [BPX], DSI Studio's Constant Solid Angle Q-Ball Imaging [CSA-QBI], MRtrix3's Constrained Spherical Deconvolution [CSD]) estimated excess fibers when insufficient data were present and/or when the data did not match the model priors. To reduce such overfitting, we developed a novel Bayesian Multi-tensor Model-selection (BaMM) method and applied it to the popular ball-and-stick model used in BedpostX within the FSL software package. BaMM was robust to overfitting and showed high reliability and the relatively best crossing-fiber accuracy with increasing amounts of diffusion data. Thus, sufficient data and an overfitting resistant analysis method enhance precision diffusion imaging. For potential clinical applications of diffusion imaging, such as neurosurgical planning and deep brain stimulation (DBS), the quantities of data required to achieve diffusion imaging reliability are lower than those needed for functional MRI.

Mapping the human connectome using diffusion MRI at 300 mT/m gradient strength: Methodological advances and scientific impact

Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength diffusion MRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for diffusion MRI and where the field is headed in the coming years.

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.

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

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

Quantitative mapping of the brain's structural connectivity using diffusion MRI tractography: A review

Diffusion magnetic resonance imaging (dMRI) tractography is an advanced imaging technique that enables in vivo reconstruction of the brain's white matter connections at macro scale. It provides an important tool for quantitative mapping of the brain's structural connectivity using measures of connectivity or tissue microstructure. Over the last two decades, the study of brain connectivity using dMRI tractography has played a prominent role in the neuroimaging research landscape. In this paper, we provide a high-level overview of how tractography is used to enable quantitative analysis of the brain's structural connectivity in health and disease. We focus on two types of quantitative analyses of tractography, including: 1) tract-specific analysis that refers to research that is typically hypothesis-driven and studies particular anatomical fiber tracts, and 2) connectome-based analysis that refers to research that is more data-driven and generally studies the structural connectivity of the entire brain. We first provide a review of methodology involved in three main processing steps that are common across most approaches for quantitative analysis of tractography, including methods for tractography correction, segmentation and quantification. For each step, we aim to describe methodological choices, their popularity, and potential pros and cons. We then review studies that have used quantitative tractography approaches to study the brain's white matter, focusing on applications in neurodevelopment, aging, neurological disorders, mental disorders, and neurosurgery. We conclude that, while there have been considerable advancements in methodological technologies and breadth of applications, there nevertheless remains no consensus about the "best" methodology in quantitative analysis of tractography, and researchers should remain cautious when interpreting results in research and clinical applications.

Fiber tracing and microstructural characterization among audiovisual integration brain regions in neonates compared with young adults

Audiovisual integration has been related with cognitive-processing and behavioral advantages, as well as with various socio-cognitive disorders. While some studies have identified brain regions instantiating this ability shortly after birth, little is known about the structural pathways connecting them. The goal of the present study was to reconstruct fiber tracts linking AVI regions in the newborn in-vivo brain and assess their adult-likeness by comparing them with analogous fiber tracts of young adults. We performed probabilistic tractography and compared connective probabilities between a sample of term-born neonates (N = 311; the Developing Human Connectome Project (dHCP, http://www.developingconnectome.org) and young adults (N = 311 The Human Connectome Project; https://www.humanconnectome.org/) by means of a classification algorithm. Furthermore, we computed Dice coefficients to assess between-group spatial similarity of the reconstructed fibers and used diffusion metrics to characterize neonates' AVI brain network in terms of microstructural properties, interhemispheric differences and the association with perinatal covariates and biological sex. Overall, our results indicate that the AVI fiber bundles were successfully reconstructed in a vast majority of neonates, similarly to adults. Connective probability distributional similarities and spatial overlaps of AVI fibers between the two groups differed across the reconstructed fibers. There was a rank-order correspondence of the fibers' connective strengths across the groups. Additionally, the study revealed patterns of diffusion metrics in line with early white matter developmental trajectories and a developmental advantage for females. Altogether, these findings deliver evidence of meaningful structural connections among AVI regions in the newborn in-vivo brain.

Post-Concussive Vestibular Dysfunction Is Related to Injury to the Inferior Vestibular Nerve

Symptoms of vestibular dysfunction such as dizziness and vertigo are common after sports-related concussions (SRC) and associated with a worse outcome and a prolonged recovery. Vestibular dysfunction after SRC can be because of an impairment of the peripheral or central neural parts of the vestibular system. The aim of the present study was to establish the cause of vestibular impairment in athletes with SRC who have persisting post-concussive symptoms (PPCS). We recruited 42 participants-21 athletes with previous SRCs and PPCS ≥6 months and 21 healthy athletic age- and sex-matched controls-who underwent symptom rating, a detailed test battery of vestibular function and 7T magnetic resonance imaging with diffusion tensor imaging (DTI) and diffusion kurtosis imaging (DKI) of cerebellar white matter tracts, and T1-weighted imaging for cerebellar volumetrics. Vestibular dysfunction was observed in 13 SRC athletes and three controls (p = 0.001). Athletes with vestibular dysfunction reported more pronounced symptoms on the Dizziness Handicap Inventory (DHI; p < 0.001) and the Hospital Anxiety and Depression Scale (HADS; p < 0.001). No significant differences in DTI metrics were found, while in DKI two metrics were observed in the superior and/or inferior cerebellar tracts. Cerebellar gray and white matter volumes were similar in athletes with SRC and controls. Compared with controls, pathological video head impulse test results (vHIT; p < 0.001) and cervical vestibular evoked myogenic potentials (cVEMP; p = 0.002) were observed in athletes with SRC, indicating peripheral vestibular dysfunction and specifically suggesting injury to the inferior vestibular nerve. In athletes with persisting symptoms after SRC, vestibular dysfunction is associated with injury to the inferior vestibular nerve.

Automated hippocampal unfolding for morphometry and subfield segmentation with HippUnfold

Like neocortical structures, the archicortical hippocampus differs in its folding patterns across individuals. Here, we present an automated and robust BIDS-App, HippUnfold, for defining and indexing individual-specific hippocampal folding in MRI, analogous to popular tools used in neocortical reconstruction. Such tailoring is critical for inter-individual alignment, with topology serving as the basis for homology. This topological framework enables qualitatively new analyses of morphological and laminar structure in the hippocampus or its subfields. It is critical for refining current neuroimaging analyses at a meso- as well as micro-scale. HippUnfold uses state-of-the-art deep learning combined with previously developed topological constraints to generate uniquely folded surfaces to fit a given subject's hippocampal conformation. It is designed to work with commonly employed sub-millimetric MRI acquisitions, with possible extension to microscopic resolution. In this paper, we describe the power of HippUnfold in feature extraction, and highlight its unique value compared to several extant hippocampal subfield analysis methods.

White Matter Connectivity between Structures of the Basal Ganglia using 3T and 7T

Analysis of the basal ganglia has been important in investigating the effects of Parkinson's disease as well as treatments for Parkinson's disease. One method of analysis has been using MRI for non-invasively segmenting the basal ganglia, then investigating significant parameters that involve the basal ganglia, such as fiber orientations and positional markers for deep brain stimulation (DBS). Following enhancements to optimizations and improvements to 3T and 7T MRI acquisitions, we utilized Lead-DBS on human connectome project data to automatically segment the basal ganglia of 49 human connectome project subjects, reducing the reliance on manual segmentation for more consistency. We generated probabilistic tractography streamlines between each segmentation pair using 3T and 7T human connectome diffusion data to observe any major differences in tractography streamline patterns that can arise due to tradeoffs from different field strengths and acquisitions. Tractography streamlines generated between basal ganglia structures using 3T images showed less standard deviation in streamline count than using 7T images. Mean tractography streamline counts generated using 3T diffusion images were all higher in count than streamlines generated using 7T diffusion images. We illustrate a potential method for analyzing the structural connectivity between basal ganglia structures, as well as visualize possible differences in probabilistic tractography that can arise from different acquisition protocols.

Biophysical compartment models for single-shell diffusion MRI in the human brain: a model fitting comparison

Clinically oriented studies commonly acquire diffusion MRI (dMRI) data with a single non-zero b-value (i.e. single-shell) and diffusion weighting ofb=1000 s/mm². To produce microstructural parameter maps, the tensor model is usually used, despite known limitations. Although compartment models have demonstrated improved fits in multi-shell dMRI data, they are rarely used for single-shell parameter maps, where their effectiveness is unclear from the literature. Here, various compartment models combining isotropic balls and symmetric tensors were fitted to single-shell dMRI data to investigate model fitting optimization and extract the most information possible. Full testing was performed in 5 subjects, and 3 subjects with multi-shell data were included for comparison. The results were tested and confirmed in a further 50 subjects. The Markov chain Monte Carlo (MCMC) model fitting technique outperformed non-linear least squares. Using MCMC, the 2-fibre-orientation mono-exponential ball & stick model (BS_ME²) provided artifact-free, stable results, in little processing time. The analogous ball & zeppelin model (BZ²) also produced stable, low-noise parameter maps, though it required much greater computing resources (50 000 burn-in steps). In single-shell data, the gamma-distributed diffusivity ball & stick model (BS_GD²) underperformed relative to other models, despite being an often-used software default. It produced artifacts in the diffusivity maps even with extremely long processing times. Neither increased diffusion weighting nor a greater number of gradient orientations improved BS_GD²fits. In white matter (WM), the tensor produced the best fit as measured by Bayesian information criterion. This result contrasts with studies using multi-shell data. However, in crossing fibre regions the tensor confounded geometric effects with fractional anisotropy (FA): the planar/linear WM FA ratio was 49%, while BZ²and BS_ME²retained 76% and 83% of restricted fraction, respectively. As a result, the BZ²and BS_ME²models are strong candidates to optimize information extraction from single-shell dMRI studies.

Structural connectivity-based segmentation of the human entorhinal cortex

The medial (MEC) and lateral entorhinal cortex (LEC), widely studied in rodents, are well defined and characterized. In humans, however, the exact locations of their homologues remain uncertain. Previous functional magnetic resonance imaging (fMRI) studies have subdivided the human EC into posteromedial (pmEC) and anterolateral (alEC) parts, but uncertainty remains about the choice of imaging modality and seed regions, in particular in light of a substantial revision of the classical model of EC connectivity based on novel insights from rodent anatomy. Here, we used structural, not functional imaging, namely diffusion tensor imaging (DTI) and probabilistic tractography to segment the human EC based on differential connectivity to other brain regions known to project selectively to MEC or LEC. We defined MEC as more strongly connected with presubiculum and retrosplenial cortex (RSC), and LEC as more strongly connected with distal CA1 and proximal subiculum (dCA1pSub) and lateral orbitofrontal cortex (OFC). Although our DTI segmentation had a larger medial-lateral component than in the previous fMRI studies, our results show that the human MEC and LEC homologues have a border oriented both towards the posterior-anterior and medial-lateral axes, supporting the differentiation between pmEC and alEC.

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.

Improvement of magnetic resonance imaging using a wireless radiofrequency resonator array

In recent years, new human magnetic resonance imaging systems operating at static magnetic fields strengths of 7 Tesla or higher have become available, providing better signal sensitivity compared with lower field strengths. However, imaging human-sized objects at such high field strength and associated precession frequencies is limited due to the technical challenges associated with the wavelength effect, which substantially disturb the transmit field uniformity over the human body when conventional coils are used. Here we report a novel passive inductively-coupled radiofrequency resonator array design with a simple structure that works in conjunction with conventional coils and requires only to be tuned to the scanner's operating frequency. We show that inductive-coupling between the resonator array and the coil improves the transmit efficiency and signal sensitivity in the targeted region. The simple structure, flexibility, and cost-efficiency make the proposed array design an attractive approach for altering the transmit field distribution specially at high field systems, where the wavelength is comparable with the tissue size.

Active Cortex Tractography

Most existing diffusion tractography algorithms are affected by gyral bias, causing the termination of streamlines at gyral crowns instead of sulcal banks. In this paper, we propose a tractography technique, called active cortex tractography (ACT), to overcome gyral bias by enabling fiber streamlines to curve naturally into the cortex. We show that the cortex can play an active role in cortical tractography by providing anatomical information to overcome orientation ambiguities as the streamlines enter the superficial white matter in gyral blades and approach the cortex. This is achieved by devising a direction scouting mechanism that takes into account the white matter surface normal vectors. The scouting mechanism allows probing of directions further in space to prepare the streamlines to turn at appropriate angles. The surface normal vectors guide the streamlines to turn into the cortex, perpendicular to the white-gray matter interface. Evaluation using synthetic, macaque and human data with different streamline seeding schemes demonstrates that ACT improves cortical tractography.

Functional and diffusion MRI reveal the neurophysiological basis of neonates' noxious-stimulus evoked brain activity

Understanding the neurophysiology underlying neonatal responses to noxious stimulation is central to improving early life pain management. In this neonatal multimodal MRI study, we use resting-state and diffusion MRI to investigate inter-individual variability in noxious-stimulus evoked brain activity. We observe that cerebral haemodynamic responses to experimental noxious stimulation can be predicted from separately acquired resting-state brain activity (n = 18). Applying this prediction model to independent Developing Human Connectome Project data (n = 215), we identify negative associations between predicted noxious-stimulus evoked responses and white matter mean diffusivity. These associations are subsequently confirmed in the original noxious stimulation paradigm dataset, validating the prediction model. Here, we observe that noxious-stimulus evoked brain activity in healthy neonates is coupled to resting-state activity and white matter microstructure, that neural features can be used to predict responses to noxious stimulation, and that the dHCP dataset could be utilised for future exploratory research of early life pain system neurophysiology.

Preliminary Comparison of Subcortical Structures in Elderly Subclinical Depression: Structural Analysis with 3T MRI

Depression in the elderly population has shown increased likelihood of neurological disorders due to structural changes in the subcortical area. However, further investigation into depression related subcortical changes is needed due to mismatches in structural analysis results between studies as well as scarcities in research regarding subcortical connectivity patterns of subclinical depression populations. This study aims to investigate structural differences in subcortical regions of aged participants with subclinical depression using 3Tesla MRI. In structural analysis, volumes of each subcortical region were measured to observe the volumetric difference and asymmetry between groups, but no significant difference was found. In addition, fractional anisotropy (FA) and apparent diffusion coefficient (ADC) did not show any significant differences between groups. Structural analysis using probabilistic tractography indicated that the connection strength between left nucleus accumbens-right hippocampus, and right thalamus-right caudate was higher in the control group than the subclinical depression group. The differences in subcortical connection strength of subclinical depression groups, have shown to correlate with emotional and cognitive disorders, such as anxiety and memory impairment. We believe that the analysis of structural differences and cross-regional network measures in subcortical structures can help identify neurophysiological changes occurring in subclinical depression.

Brain/MINDS beyond human brain MRI project: A protocol for multi-level harmonization across brain disorders throughout the lifespan

Psychiatric and neurological disorders are afflictions of the brain that can affect individuals throughout their lifespan. Many brain magnetic resonance imaging (MRI) studies have been conducted; however, imaging-based biomarkers are not yet well established for diagnostic and therapeutic use. This article describes an outline of the planned study, the Brain/MINDS Beyond human brain MRI project (BMB-HBM, FY2018 ~ FY2023), which aims to establish clinically-relevant imaging biomarkers with multi-site harmonization by collecting data from healthy traveling subjects (TS) at 13 research sites. Collection of data in psychiatric and neurological disorders across the lifespan is also scheduled at 13 sites, whereas designing measurement procedures, developing and analyzing neuroimaging protocols, and databasing are done at three research sites. A high-quality scanning protocol, Harmonization Protocol (HARP), was established for five high-quality 3 T scanners to obtain multimodal brain images including T1 and T2-weighted, resting-state and task functional and diffusion-weighted MRI. Data are preprocessed and analyzed using approaches developed by the Human Connectome Project. Preliminary results in 30 TS demonstrated cortical thickness, myelin, functional connectivity measures are comparable across 5 scanners, suggesting sensitivity to subject-specific connectome. A total of 75 TS and more than two thousand patients with various psychiatric and neurological disorders are scheduled to participate in the project, allowing a mixed model statistical harmonization. The HARP protocols are publicly available online, and all the imaging, demographic and clinical information, harmonizing database will also be made available by 2024. To the best of our knowledge, this is the first project to implement a prospective, multi-level harmonization protocol with multi-site TS data. It explores intractable brain disorders across the lifespan and may help to identify the disease-specific pathophysiology and imaging biomarkers for clinical practice.

Modelling white matter in gyral blades as a continuous vector field

Many brain imaging studies aim to measure structural connectivity with diffusion tractography. However, biases in tractography data, particularly near the boundary between white matter and cortical grey matter can limit the accuracy of such studies. When seeding from the white matter, streamlines tend to travel parallel to the convoluted cortical surface, largely avoiding sulcal fundi and terminating preferentially on gyral crowns. When seeding from the cortical grey matter, streamlines generally run near the cortical surface until reaching deep white matter. These so-called “gyral biases” limit the accuracy and effective resolution of cortical structural connectivity profiles estimated by tractography algorithms, and they do not reflect the expected distributions of axonal densities seen in invasive tracer studies or stains of myelinated fibres. We propose an algorithm that concurrently models fibre density and orientation using a divergence-free vector field within gyral blades to encourage an anatomically-justified streamline density distribution along the cortical white/grey-matter boundary while maintaining alignment with the diffusion MRI estimated fibre orientations. Using in vivo data from the Human Connectome Project, we show that this algorithm reduces tractography biases. We compare the structural connectomes to functional connectomes from resting-state fMRI, showing that our model improves cross-modal agreement. Finally, we find that after parcellation the changes in the structural connectome are very minor with slightly improved interhemispheric connections (i.e, more homotopic connectivity) and slightly worse intrahemi-spheric connections when compared to tracers.

Effects of Repetitive Transcranial Magnetic Stimulation on Working Memory Performance and Brain Structure in People With Schizophrenia Spectrum Disorders: A Double-Blind, Randomized, Sham-Controlled Trial

BACKGROUND

There are currently no approved treatments for working memory deficits in schizophrenia spectrum disorders (SSDs). The objective of the present study was to assess whether repetitive transcranial magnetic stimulation (rTMS) to the bilateral dorsolateral prefrontal cortex (DLPFC) in people with SSDs 1) improves working memory deficits and 2) changes brain structure.

METHODS

We conducted a double-blind, parallel, randomized, sham-controlled study at the Centre for Addiction and Mental Health in Toronto, Canada. We randomized 83 participants with SSDs to receive either active 20 Hz rTMS applied to the bilateral DLPFC or sham rTMS for 4 weeks. The participants also completed pre/posttreatment magnetic resonance imaging. Clinical and cognitive assessments were performed at baseline, treatment end, and 1 month later. The primary outcome was change in verbal n-back working memory performance accuracy (d-prime). The secondary outcome measures were change in DLPFC thickness and fractional anisotropy of white matter tracts connecting to the DLPFC. Prespecified exploratory outcome measures were changes in general cognition; positive, negative, and depressive symptoms.

RESULTS

Compared with sham treatment, active rTMS did not lead to significant change in working memory performance; it was associated with an increase in right DLPFC thickness but not fractional anisotropy. Prespecified exploratory analysis showed a significant decrease in depressive symptoms in the active group; the decrease in depressive symptoms was correlated with an increase in right DLPFC thickness.

CONCLUSIONS

Although rTMS applied to the bilateral DLPFC was not efficacious in treating working memory deficits in SSDs, it did increase right DLPFC thickness and decrease depressive symptoms. These findings deserve further study given the lack of efficacy of antidepressant medications in SSDs.

Structural characterization of the Extended Frontal Aslant Tract trajectory: A ML-validated laterality study in 3T and 7T

The Extended Frontal Aslant Tract (exFAT) is a recently described tractography-based extension of the Frontal Aslant Tract connecting Broca's territory to both supplementary and pre-supplementary motor areas, and more anterior prefrontal regions. In this study, we aim to characterize the microstructural properties of the exFAT trajectories as a means to perform a laterality analysis to detect interhemispheric structural differences along the tracts using the Human Connectome Project (HCP) dataset. To that end, the bilateral exFAT was reconstructed for 3T and 7T HCP acquisitions in 120 randomly selected subjects. As a complementary exploration of the exFAT anatomy, we performed a white matter dissection of the exFAT trajectory of two ex-vivo left hemispheres that provide a qualitative assessment of the tract profiles. We assessed the lateralization structural differences in the exFAT by performing: (i) a laterality comparison between the mean microstructural diffusion-derived parameters for the exFAT trajectories, (ii) a laterality comparison between the tract profiles obtained by applying the Automated Fiber Quantification (AFQ) algorithm, and (iii) a cross-validated Machine Learning (ML) classifier analysis using single and combined tract profiles parameters for single-subject classification. The mean microstructural diffusion-derived parameter comparison showed statistically significant differences in mean FA values between left and right exFATs in the 3T sample. The diffusion parameters studied with the AFQ technique suggest that the inferiormost half of the exFAT trajectory has a hemispheric-dependent fingerprint of microstructural properties, with an increased measure of tissue hindrance in the orthogonal plane and a decreased measure of orientational dispersion along the main tract direction in the left exFAT compared to the right exFAT. The classification accuracy of the ML models showed a high agreement with the magnitude of those differences.

Non-negative data-driven mapping of structural connections with application to the neonatal brain

Mapping connections in the neonatal brain can provide insight into the crucial early stages of neurodevelopment that shape brain organisation and lay the foundations for cognition and behaviour. Diffusion MRI and tractography provide unique opportunities for such explorations, through estimation of white matter bundles and brain connectivity. Atlas-based tractography protocols, i.e. a priori defined sets of masks and logical operations in a template space, have been commonly used in the adult brain to drive such explorations. However, rapid growth and maturation of the brain during early development make it challenging to ensure correspondence and validity of such atlas-based tractography approaches in the developing brain. An alternative can be provided by data-driven methods, which do not depend on predefined regions of interest. Here, we develop a novel data-driven framework to extract white matter bundles and their associated grey matter networks from neonatal tractography data, based on non-negative matrix factorisation that is inherently suited to the non-negative nature of structural connectivity data. We also develop a non-negative dual regression framework to map group-level components to individual subjects. Using in-silico simulations, we evaluate the accuracy of our approach in extracting connectivity components and compare with an alternative data-driven method, independent component analysis. We apply non-negative matrix factorisation to whole-brain connectivity obtained from publicly available datasets from the Developing Human Connectome Project, yielding grey matter components and their corresponding white matter bundles. We assess the validity and interpretability of these components against traditional tractography results and grey matter networks obtained from resting-state fMRI in the same subjects. We subsequently use them to generate a parcellation of the neonatal cortex using data from 323 new-born babies and we assess the robustness and reproducibility of this connectivity-driven parcellation.

Diffusion Imaging in the Post HCP Era

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.

Multicenter dataset of multi-shell diffusion MRI in healthy traveling adults with identical settings

Multicenter diffusion magnetic resonance imaging (MRI) has drawn great attention recently due to the expanding need for large-scale brain imaging studies, whereas the variability in MRI scanners and data acquisition tends to confound reliable individual-based analysis of diffusion measures. In addition, a growing number of multi-shell diffusion models have been shown with the potential to generate various estimates of physio-pathological information, yet their reliability and reproducibility in multicenter studies remain to be assessed. In this article, we describe a multi-shell diffusion dataset collected from three traveling subjects with identical acquisition settings in ten imaging centers. Both the scanner type and imaging protocol for anatomical and diffusion imaging were well controlled. This dataset is expected to replenish individual reproducible studies via multicenter collaboration by providing an open resource for advanced and novel microstructural and tractography modelling and quantification.

The role of node dynamics in shaping emergent functional connectivity patterns in the brain

The contribution of structural connectivity to functional brain states remains poorly understood. We present a mathematical and computational study suited to assess the structure–function issue, treating a system of Jansen–Rit neural mass nodes with heterogeneous structural connections estimated from diffusion MRI data provided by the Human Connectome Project. Via direct simulations we determine the similarity of functional (inferred from correlated activity between nodes) and structural connectivity matrices under variation of the parameters controlling single-node dynamics, highlighting a nontrivial structure–function relationship in regimes that support limit cycle oscillations. To determine their relationship, we firstly calculate network instabilities giving rise to oscillations, and the so-called ‘false bifurcations’ (for which a significant qualitative change in the orbit is observed, without a change of stability) occurring beyond this onset. We highlight that functional connectivity (FC) is inherited robustly from structure when node dynamics are poised near a Hopf bifurcation, whilst near false bifurcations, and structure only weakly influences FC. Secondly, we develop a weakly coupled oscillator description to analyse oscillatory phase-locked states and, furthermore, show how the modular structure of FC matrices can be predicted via linear stability analysis. This study thereby emphasises the substantial role that local dynamics can have in shaping large-scale functional brain states.

Patterns of oscillation across the brain arise because of structural connections between brain regions. However, the type of oscillation at a site may also play a contributory role. We focus on an idealised model of a neural mass network, coupled using estimates of structural connections obtained via tractography on Human Connectome Project MRI data. Using a mixture of computational and mathematical techniques, we show that functional connectivity is inherited most strongly from structural connectivity when the network nodes are poised at a Hopf bifurcation. However, beyond the onset of this oscillatory instability a phase-locked network state can undergo a false bifurcation, and structural connectivity only weakly influences functional connectivity. This highlights the important effect that local dynamics can have on large-scale brain states.

Denoise magnitude diffusion magnetic resonance images via variance-stabilizing transformation and optimal singular-value manipulation

Although shown to have a great utility for a wide range of neuroscientific and clinical applications, diffusion-weighted magnetic resonance imaging (dMRI) faces a major challenge of low signal-to-noise ratio (SNR), especially when pushing the spatial resolution for improved delineation of brain's fine structure or increasing the diffusion weighting for increased angular contrast or both. Here, we introduce a comprehensive denoising framework for denoising magnitude dMRI. The framework synergistically combines the variance stabilizing transform (VST) with optimal singular value manipulation. The purpose of VST is to transform the Rician data to Gaussian-like data so that an asymptotically optimal singular value manipulation strategy tailored for Gaussian data can be used. The output of the framework is the estimated underlying diffusion signal for each voxel in the image domain. The usefulness of the proposed framework for denoising magnitude dMRI is demonstrated using both simulation and real-data experiments. Our results show that the proposed denoising framework can significantly improve SNR across the entire brain, leading to substantially enhanced performances for estimating diffusion tensor related indices and for resolving crossing fibers when compared to another competing method. More encouragingly, the proposed method when used to denoise a single average of 7 ​Tesla Human Connectome Project-style diffusion acquisition provided comparable performances relative to those achievable with ten averages for resolving multiple fiber populations across the brain. As such, the proposed denoising method is expected to have a great utility for high-quality, high-resolution whole-brain dMRI, desirable for many neuroscientific and clinical applications.

A field-monitoring-based approach for correcting eddy-current-induced artifacts of up to the 2nd spatial order in human-connectome-project-style multiband diffusion MRI experiment at 7T: A pilot study

Over the recent years, significant advances in Spin-Echo (SE) Echo-Planar (EP) Diffusion MRI (dMRI) have enabled improved fiber tracking conspicuity in the human brain. At the same time, pushing the spatial resolution and using higher b-values inherently expose the acquired images to further eddy-current-induced distortion and blurring. Recently developed data-driven correction techniques, capable of significantly mitigating these defects, are included in the reconstruction pipelines developed for the Human Connectome Project (HCP) driven by the NIH BRAIN initiative. In this case, however, corrections are derived from the original diffusion-weighted (DW) magnitude images affected by distortion and blurring. Considering the complexity of k-space deviations in the presence of time varying high spatial order eddy currents, distortion and blurring may not be fully reversed when relying on magnitude DW images only. An alternative approach, consisting of iteratively reconstructing DW images based on the actual magnetic field spatiotemporal evolution measured with a magnetic field monitoring camera, has been successfully implemented at 3T in single band dMRI (Wilm et al., 2017, 2015). In this study, we aim to demonstrate the efficacy of this eddy current correction method in the challenging context of HCP-style multiband (MB ​= ​2) dMRI protocol. The magnetic field evolution was measured during the EP-dMRI readout echo train with a field monitoring camera equipped with 16 ¹⁹F NMR probes. The time variation of 0^th, 1^st and 2^nd order spherical field harmonics were used to reconstruct DW images. Individual DW images reconstructed with and without field correction were compared. The impact of eddy current correction was evaluated by comparing the corresponding direction-averaged DW images and fractional anisotropy (FA) maps. ¹⁹F field monitoring data confirmed the existence of significant field deviations induced by the diffusion-encoding gradients, with variations depending on diffusion gradient amplitude and direction. In DW images reconstructed with the field correction, residual aliasing artifacts were reduced or eliminated, and when high b-values were applied, better gray/white matter delineation and sharper gyri contours were observed, indicating reduced signal blurring. The improvement in image quality further contributed to sharper contours and better gray/white matter delineation in mean DW images and FA maps. In conclusion, we demonstrate that up-to-2^nd-order-eddy-current-induced field perturbation in multiband, in-plane accelerated HCP-style dMRI acquisition at 7T can be corrected by integrating the measured field evolution in image reconstruction.

Towards HCP-Style macaque connectomes: 24-Channel 3T multi-array coil, MRI sequences and preprocessing

Macaque monkeys are an important animal model where invasive investigations can lead to a better understanding of the cortical organization of primates including humans. However, the tools and methods for noninvasive image acquisition (e.g. MRI RF coils and pulse sequence protocols) and image data preprocessing have lagged behind those developed for humans. To resolve the structural and functional characteristics of the smaller macaque brain, high spatial, temporal, and angular resolutions combined with high signal-to-noise ratio are required to ensure good image quality. To address these challenges, we developed a macaque 24-channel receive coil for 3-T MRI with parallel imaging capabilities. This coil enables adaptation of the Human Connectome Project (HCP) image acquisition protocols to the in-vivo macaque brain. In addition, we adapted HCP preprocessing methods to the macaque brain, including spatial minimal preprocessing of structural, functional MRI (fMRI), and diffusion MRI (dMRI). The coil provides the necessary high signal-to-noise ratio and high efficiency in data acquisition, allowing four- and five-fold accelerations for dMRI and fMRI. Automated FreeSurfer segmentation of cortex, reconstruction of cortical surface, removal of artefacts and nuisance signals in fMRI, and distortion correction of dMRI all performed well, and the overall quality of basic neurobiological measures was comparable with those for the HCP. Analyses of functional connectivity in fMRI revealed high sensitivity as compared with those from publicly shared datasets. Tractography-based connectivity estimates correlated with tracer connectivity similarly to that achieved using ex-vivo dMRI. The resulting HCP-style in vivo macaque MRI data show considerable promise for analyzing cortical architecture and functional and structural connectivity using advanced methods that have previously only been available in studies of the human brain.

Rationalization of Using the MR Diffusion Imaging in B12 Deficiency

Context:

The structural imaging of brain does not demonstrate any changes in the vast majority of patients with vitamin B12 deficiency, even in the advanced stages.

Aims:

We investigated the microstructural changes in the brain with diffusion imaging among patients with biochemical evidence of B12 deficiency.

Patients and Methods:

We retrospectively analyzed all diffusion-weighted MRI images between the periods 2014–2016 who had biochemical evidence of B12. The age-sex matched controls were chosen from the group with normal B12 levels. Patients with pathological findings in conventional MRI images were excluded from the study.

Results:

About 37 patients were recruited (22 women, 15 men; mean age, 34.1 ± 9.9 years; age range). They were about thirty-four age-and sex-matched controls (with normal B12 levels), which were also included in the study. The mean apparent diffusion coefficient (ADC) value of amygdala (773.8 ± 49.9 vs. 742.2 ± 24.2, P = 0.01), hypothalamus (721.3 ± 39.2 vs. 700.2 ± 38.2, P = 0.02), striate cortex (737.6 ± 77.6 vs. 704.3 ± 58.2, P = 0.04), suprafrontal gyrus (740.7 ± 46.9 vs. 711.6 ± 40.7, P = 0.007) and medulla oblongata-olivary nucleus (787.3 ± 56.4 vs. 759.7 ± 46.2, P = 0.02) were significantly higher in B12 deficiency group compared to controls, whereas ADC values were similar at hippocampus, thalamus, insula, corpus striatum, cingulate gyrus, occipital gyrus, dentate nucleus, cerebral pedicle, tegmentum, pons, and posterior medulla oblongata.

Conclusions:

Our study indicates that a significant increase in ADC values occurs in multiple brain regions in patients with vitamin B12.

Mitigating gyral bias in cortical tractography via asymmetric fiber orientation distributions

Diffusion tractography in brain connectomics often involves tracing axonal trajectories across gray-white matter boundaries in gyral blades of complex cortical convolutions. To date, gyral bias is observed in most tractography algorithms with streamlines predominantly terminating at gyral crowns instead of sulcal banks. This work demonstrates that asymmetric fiber orientation distribution functions (AFODFs), computed via a multi-tissue global estimation framework, can mitigate the effects of gyral bias, enabling fiber streamlines at gyral blades to make sharper turns into the cortical gray matter. We use ex-vivo data of an adult rhesus macaque and in-vivo data from the Human Connectome Project (HCP) to show that the fiber streamlines given by AFODFs bend more naturally into the cortex than the conventional symmetric FODFs in typical gyral blades. We demonstrate that AFODF tractography improves cortico-cortical connectivity and provides highly consistent outcomes between two different field strengths (3T and 7T).

The role of node dynamics in shaping emergent functional connectivity patterns in the brain

The contribution of structural connectivity to functional brain states remains poorly understood. We present a mathematical and computational study suited to assess the structure-function issue, treating a system of Jansen-Rit neural mass nodes with heterogeneous structural connections estimated from diffusion MRI data provided by the Human Connectome Project. Via direct simulations we determine the similarity of functional (inferred from correlated activity between nodes) and structural connectivity matrices under variation of the parameters controlling single-node dynamics, highlighting a nontrivial structure-function relationship in regimes that support limit cycle oscillations. To determine their relationship, we firstly calculate network instabilities giving rise to oscillations, and the so-called 'false bifurcations' (for which a significant qualitative change in the orbit is observed, without a change of stability) occurring beyond this onset. We highlight that functional connectivity (FC) is inherited robustly from structure when node dynamics are poised near a Hopf bifurcation, whilst near false bifurcations, and structure only weakly influences FC. Secondly, we develop a weakly coupled oscillator description to analyse oscillatory phase-locked states and, furthermore, show how the modular structure of FC matrices can be predicted via linear stability analysis. This study thereby emphasises the substantial role that local dynamics can have in shaping large-scale functional brain states.

Next-generation neural mass and field modeling

The Wilson-Cowan population model of neural activity has greatly influenced our understanding of the mechanisms for the generation of brain rhythms and the emergence of structured brain activity. As well as the many insights that have been obtained from its mathematical analysis, it is now widely used in the computational neuroscience community for building large-scale in silico brain networks that can incorporate the increasing amount of knowledge from the Human Connectome Project. Here, we consider a neural population model in the spirit of that originally developed by Wilson and Cowan, albeit with the added advantage that it can account for the phenomena of event-related synchronization and desynchronization. This derived mean-field model provides a dynamic description for the evolution of synchrony, as measured by the Kuramoto order parameter, in a large population of quadratic integrate-and-fire model neurons. As in the original Wilson-Cowan framework, the population firing rate is at the heart of our new model; however, in a significant departure from the sigmoidal firing rate function approach, the population firing rate is now obtained as a real-valued function of the complex-valued population synchrony measure. To highlight the usefulness of this next-generation Wilson-Cowan style model, we deploy it in a number of neurobiological contexts, providing understanding of the changes in power spectra observed in electro- and magnetoencephalography neuroimaging studies of motor cortex during movement, insights into patterns of functional connectivity observed during rest and their disruption by transcranial magnetic stimulation, and to describe wave propagation across cortex.

Default Mode Network structural alterations in Kocher-Monro trajectory white matter transection: A 3 and 7 tesla simulation modeling approach

The Kocher-Monro trajectory to the cerebral ventricular system represents one of the most common surgical procedures in the field of neurosurgery. Several studies have analyzed the specific white matter disruption produced during this intervention, which has no reported adverse neurological outcomes. In this study, a graph-theoretical approach was applied to quantify the structural alterations in whole-brain level connectivity. To this end, 132 subjects were randomly selected from the Human Connectome Project dataset and used to create 3 independent 44 subjects groups. Two of the groups underwent a simulated left/right Kocher-Monro trajectory and the third was kept as a control group. For the right Kocher-Monro approach, the nodal analysis revealed decreased strength in the anterior cingulate gyrus of the transected hemisphere. The network-based statistic analysis revealed a set of right lateralized subnetworks with decreased connectivity strength that is consistent with a subset of the Default Mode Network, Salience Network, and Cingulo-Opercular Network. These findings could allow for a better understanding of structural alterations caused by Kocher-Monro approaches that could reveal previously undetected clinical alterations and inform the process of designing safer and less invasive cerebral ventricular approaches.

Joint modelling of diffusion MRI and microscopy

The combination of diffusion MRI (dMRI) with microscopy provides unique opportunities to study microstructural features of tissue, particularly when acquired in the same sample. Microscopy is frequently used to validate dMRI microstructure models, addressing the indirect nature of dMRI signals. Typically, these modalities are analysed separately, and microscopy is taken as a gold standard against which dMRI-derived parameters are validated. Here we propose an alternative approach in which we combine dMRI and microscopy data obtained from the same tissue sample to drive a single, joint model. This simultaneous analysis allows us to take advantage of the breadth of information provided by complementary data acquired from different modalities. By applying this framework to a spherical-deconvolution analysis, we are able to overcome a known degeneracy between fibre dispersion and radial diffusion. Spherical-deconvolution based approaches typically estimate a global fibre response function to determine the fibre orientation distribution in each voxel. However, the assumption of a 'brain-wide' fibre response function may be challenged if the diffusion characteristics of white matter vary across the brain. Using a generative joint dMRI-histology model, we demonstrate that the fibre response function is dependent on local anatomy, and that current spherical-deconvolution based models may be overestimating dispersion and underestimating the number of distinct fibre populations per voxel.