Surface based analysis of diffusion orientation for identifying architectonic domains in the in vivo human cortex

Romer-EPTI: rotating-view motion-robust super-resolution EPTI for SNR-efficient distortion-free in-vivo mesoscale dMRI and microstructure imaging

Purpose

To overcome the major challenges in dMRI acquisition, including low SNR, distortion/blurring, and motion vulnerability.

Methods

A novel Romer-EPTI technique is developed to provide distortion-free dMRI with significant SNR gain, high motion-robustness, sharp spatial resolution, and simultaneous multi-TE imaging. It introduces a ROtating-view Motion-robust supEr-Resolution technique (Romer) combined with a distortion/blurring-free EPTI encoding. Romer enhances SNR by a simultaneous multi-thick-slice acquisition with rotating-view encoding, while providing high motion-robustness through a motion-aware super-resolution reconstruction, which also incorporates slice-profile and real-value diffusion, to resolve high-isotropic-resolution volumes. The in-plane encoding is performed using distortion/blurring-free EPTI, which further improves effective spatial resolution and motion robustness by preventing not only T2/T2*-blurring but also additional blurring resulting from combining encoded volumes with inconsistent geometries caused by dynamic distortions. Self-navigation was incorporated to enable efficient phase correction. Additional developments include strategies to address slab-boundary artifacts, achieve minimal TE for SNR gain at 7T, and achieve high robustness to strong phase variations at high b-values.

Results

Using Romer-EPTI, we demonstrate distortion-free whole-brain mesoscale in-vivo dMRI at both 3T (500-μm-iso) and 7T (485-μm-iso) for the first time, with high SNR efficiency (e.g., 25 × ), and high image quality free from distortion and slab-boundary artifacts with minimal blurring. Motion experiments demonstrate Romer-EPTI's high motion-robustness and ability to recover sharp images in the presence of motion. Romer-EPTI also demonstrates significant SNR gain and robustness in high b-value (b=5000s/mm2) and time-dependent dMRI.

Conclusion

Romer-EPTI significantly improves SNR, motion-robustness, and image quality, providing a highly efficient acquisition for high-resolution dMRI and microstructure imaging.

Diffusion and structural MRI as potential biomarkers in people with Parkinson's disease and cognitive impairment

OBJECTIVE

To explore the neuroimage change in Parkinson's disease (PD) patients with cognitive impairments, this study investigated the correlation between plasma biomarkers and morphological brain changes in patients with normal cognition and mild cognitive impairment. The objective was to identify the potential target deposition regions of the plasma biomarkers and to search for the relevant early neuroimaging biomarkers on the basis of different cognitive domains.

METHODS

Structural brain MRI and diffusion weighted images were analyzed from 49 eligible PD participants (male/female: 27/22; mean age: 73.4 ± 8.5 years) from a retrospective analysis. Plasma levels of α-synuclein, amyloid beta peptide, and total tau were collected. A comprehensive neuropsychological assessment of the general and specific cognitive domains was performed. Difference between PD patients with normal cognition and impairment was examined. Regression analysis was performed to evaluate the correlation between image-derived index and plasma biomarkers or neuropsychological assessments.

RESULTS

Significant correlation was found between plasma Aβ-42 level and fractional anisotropy of the middle occipital, angular, and middle temporal gyri of the left brain, as well as plasma T-tau level and the surface area of the isthmus or the average thickness of the posterior part of right cingulate gyrus. Visuospatial and executive function is positively correlated with axial diffusivity in bilateral cingulate gyri.

CONCLUSION

In nondemented PD patients, the target regions for plasma deposition might be located in the cingulate, middle occipital, angular, and middle temporal gyri. Changes from multiple brain regions can be correlated to the performance of different cognitive domains.

CLINICAL RELEVANCE STATEMENT

Cognitive impairment in Parkinson's disease is primarily linked to biomarkers associated with Alzheimer's disease rather than those related to Parkinson's disease and resembles the frontal variant of Alzheimer's disease, which may guide management strategies for cognitive impairment in Parkinson's disease.

KEY POINTS

• Fractional anisotropy, surface area, and thickness in the cingulate, middle occipital, angular, and middle temporal gyri can be significantly correlated with plasma Aβ-42 and T-tau level. • Axial diffusivity in the cingulate gyri was correlated with visuospatial and executive function. • The pattern of cognitive impairment in Parkinson's disease can be similar to the frontal variant than typical Alzheimer's disease.

Spatiotemporal patterns of cortical microstructural maturation in children and adolescents with diffusion MRI

Neocortical maturation is a dynamic process that proceeds in a hierarchical manner; however, the spatiotemporal organization of cortical microstructure with diffusion MRI has yet to be fully defined. This study characterized cortical microstructural maturation using diffusion MRI (fwe-diffusion tensor imaging [DTI] and neurite orientation dispersion and density imaging [NODDI] multicompartment modeling) in a cohort of 637 children and adolescents between 8 and 21 years of age. We found spatially heterogeneous developmental patterns broadly demarcated into functional domains where NODDI metrics increased, and fwe-DTI metrics decreased with age. By applying nonlinear growth models in a vertex-wise analysis, we observed a general posterior-to-anterior pattern of maturation, where the fwe-DTI measures mean diffusivity and radial diffusivity reached peak maturation earlier than the NODDI metrics neurite density index. Using non-negative matrix factorization, we found occipito-parietal cortical regions that correspond to lower order sensory domains mature earlier than fronto-temporal higher order association domains. Our findings corroborate previous histological and neuroimaging studies that show spatially varying patterns of cortical maturation that may reflect unique developmental processes of cytoarchitectonically determined regional patterns of change.

Mapping the individual human cortex using multidimensional MRI and unsupervised learning

Human evolution has seen the development of higher-order cognitive and social capabilities in conjunction with the unique laminar cytoarchitecture of the human cortex. Moreover, early-life cortical maldevelopment has been associated with various neurodevelopmental diseases. Despite these connections, there is currently no noninvasive technique available for imaging the detailed cortical laminar structure. This study aims to address this scientific and clinical gap by introducing an approach for imaging human cortical lamina. This method combines diffusion-relaxation multidimensional MRI with a tailored unsupervised machine learning approach that introduces enhanced microstructural sensitivity. This new imaging method simultaneously encodes the microstructure, the local chemical composition and importantly their correlation within complex and heterogenous tissue. To validate our approach, we compared the intra-cortical layers obtained using our ex vivo MRI-based method with those derived from Nissl staining of postmortem human brain specimens. The integration of unsupervised learning with diffusion-relaxation correlation MRI generated maps that demonstrate sensitivity to areal differences in cytoarchitectonic features observed in histology. Significantly, our observations revealed layer-specific diffusion-relaxation signatures, showing reductions in both relaxation times and diffusivities at the deeper cortical levels. These findings suggest a radial decrease in myelin content and changes in cell size and anisotropy, reflecting variations in both cytoarchitecture and myeloarchitecture. Additionally, we demonstrated that 1D relaxation and high-order diffusion MRI scalar indices, even when aggregated and used jointly in a multimodal fashion, cannot disentangle the cortical layers. Looking ahead, our technique holds the potential to open new avenues of research in human neurodevelopment and the vast array of disorders caused by disruptions in neurodevelopment.

Sampling strategies and integrated reconstruction for reducing distortion and boundary slice aliasing in high-resolution 3D diffusion MRI

PURPOSE

To develop a new method for high-fidelity, high-resolution 3D multi-slab diffusion MRI with minimal distortion and boundary slice aliasing.

METHODS

Our method modifies 3D multi-slab imaging to integrate blip-reversed acquisitions for distortion correction and oversampling in the slice direction (kz ) for reducing boundary slice aliasing. Our aim is to achieve robust acceleration to keep the scan time the same as conventional 3D multi-slab acquisitions, in which data are acquired with a single direction of blip traversal and without kz -oversampling. We employ a two-stage reconstruction. In the first stage, the blip-up/down images are respectively reconstructed and analyzed to produce a field map for each diffusion direction. In the second stage, the blip-reversed data and the field map are incorporated into a joint reconstruction to produce images that are corrected for distortion and boundary slice aliasing.

RESULTS

We conducted experiments at 7T in six healthy subjects. Stage 1 reconstruction produces images from highly under-sampled data (R = 7.2) with sufficient quality to provide accurate field map estimation. Stage 2 joint reconstruction substantially reduces distortion artifacts with comparable quality to fully-sampled blip-reversed results (2.4× scan time). Whole-brain in-vivo results acquired at 1.22 mm and 1.05 mm isotropic resolutions demonstrate improved anatomical fidelity compared to conventional 3D multi-slab imaging. Data demonstrate good reliability and reproducibility of the proposed method over multiple subjects.

CONCLUSION

The proposed acquisition and reconstruction framework provide major reductions in distortion and boundary slice aliasing for 3D multi-slab diffusion MRI without increasing the scan time, which can potentially produce high-quality, high-resolution diffusion MRI.

Three-dimensional diffusion MRI using simultaneous multislab with blipped-CAIPI in a 4D k-space framework

PURPOSE

To develop an efficient simultaneous multislab imaging method with blipped-controlled aliasing in parallel imaging (blipped-SMSlab) in a 4D k-space framework, and to demonstrate its efficacy in high-resolution diffusion MRI (dMRI).

THEORY AND METHODS

First, the SMSlab 4D k-space signal expression is formulated, and the phase interferences from intraslab and interslab encodings on the same physical z-axis are analyzed. Then, the blipped-SMSlab dMRI sequence is designed, with blipped-controlled aliasing in parallel imaging (blipped-CAIPI) gradients for interslab encoding, and a 2D multiband accelerated navigator for inter-kz-shot phase correction. Third, strategies are developed to remove the phase interferences, by RF phase modulation and/or phase correction during reconstruction, thus decoupling intraslab and interslab encodings that are otherwise entangled. In vivo experiments are performed to validate the blipped-SMSlab method and preliminarily evaluate its performance in high-resolution dMRI compared with traditional 2D imaging.

RESULTS

In the 4D k-space framework, interslab and intraslab phase interferences of blipped-SMSlab are successfully removed using the proposed strategies. Compared with non-CAIPI sampling, the blipped-SMSlab acquisition reduces the g-factor and g-factor-related SNR penalty by about 12%. In addition, in vivo experiments show the SNR advantage of blipped-SMSlab dMRI over traditional 2D dMRI for 1.3-mm and 1.0-mm isotropic resolution imaging with matched acquisition time.

CONCLUSION

Removing interslab and intraslab phase interferences enables SMSlab dMRI with blipped-CAIPI in a 4D k-space framework. The proposed blipped-SMSlab dMRI is demonstrated to be more SNR-efficient than 2D dMRI and thus capable of high-quality, high-resolution fiber orientation detection.

Mapping tissue microstructure across the human brain on a clinical scanner with soma and neurite density image metrics

Soma and neurite density image (SANDI) is an advanced diffusion magnetic resonance imaging biophysical signal model devised to probe in vivo microstructural information in the gray matter (GM). This model requires acquisitions that include b values that are at least six times higher than those used in clinical practice. Such high b values are required to disentangle the signal contribution of water diffusing in soma from that diffusing in neurites and extracellular space, while keeping the diffusion time as short as possible to minimize potential bias due to water exchange. These requirements have limited the use of SANDI only to preclinical or cutting-edge human scanners. Here, we investigate the potential impact of neglecting water exchange in the SANDI model and present a 10-min acquisition protocol that enables to characterize both GM and white matter (WM) on 3 T scanners. We implemented analytical simulations to (i) evaluate the stability of the fitting of SANDI parameters when diminishing the number of shells; (ii) estimate the bias due to potential exchange between neurites and extracellular space in such reduced acquisition scheme, comparing it with the bias due to experimental noise. Then, we demonstrated the feasibility and assessed the repeatability and reproducibility of our approach by computing microstructural metrics of SANDI with AMICO toolbox and other state-of-the-art models on five healthy subjects. Finally, we applied our protocol to five multiple sclerosis patients. Results suggest that SANDI is a practical method to characterize WM and GM tissues in vivo on performant clinical scanners.

The acute effects of repetitive transcranial magnetic stimulation on laminar diffusion anisotropy of neocortical gray matter

Repetitive transcranial magnetic stimulation (rTMS) is increasingly used to treat neuropsychiatric disorders. Inhibitory and excitatory regimens have been both adopted but the exact mechanism of action remains unclear, and investigating their differential effects on laminar diffusion profiles of neocortex may add important evidence. Twenty healthy participants were randomly assigned to receive a low-frequency/inhibitory or high-frequency/excitatory rTMS targeting the left dorsolateral prefrontal cortex (DLPFC). With the brand-new submillimeter diffusion tensor imaging of whole brain and specialized surface-based laminar analysis, fractional anisotropy (FA) and mean diffusion (MD) profiles of cortical layers at different cortical depths were characterized before/after rTMS. Inhibitory and excitatory rTMS both showed impacts on diffusion metrics of somatosensory, limbic, and sensory regions, but different patterns of changes were observed-increased FA with inhibitory rTMS, whereas decreased FA with excitatory rTMS. More importantly, laminar analysis indicated laminar specificity of changes in somatosensory regions during different rTMS patterns-inhibitory rTMS affected the superficial layers contralateral to the DLPFC, while excitatory rTMS led to changes in the intermediate/deep layers bilateral to the DLPFC. These findings provide novel insights into acute neurobiological effects on diffusion profiles of rTMS that may add critical evidence relevant to different protocols of rTMS on neocortex.

High-fidelity mesoscale in-vivo diffusion MRI through gSlider-BUDA and circular EPI with S-LORAKS reconstruction

PURPOSE

To develop a high-fidelity diffusion MRI acquisition and reconstruction framework with reduced echo-train-length for less T2* image blurring compared to typical highly accelerated echo-planar imaging (EPI) acquisitions at sub-millimeter isotropic resolution.

METHODS

We first proposed a circular-EPI trajectory with partial Fourier sampling on both the readout and phase-encoding directions to minimize the echo-train-length and echo time. We then utilized this trajectory in an interleaved two-shot EPI acquisition with reversed phase-encoding polarity, to aid in the correction of off-resonance-induced image distortions and provide complementary k-space coverage in the missing partial Fourier regions. Using model-based reconstruction with structured low-rank constraint and smooth phase prior, we corrected the shot-to-shot phase variations across the two shots and recover the missing k-space data. Finally, we combined the proposed acquisition/reconstruction framework with an SNR-efficient RF-encoded simultaneous multi-slab technique, termed gSlider, to achieve high-fidelity 720 µm and 500 µm isotropic resolution in-vivo diffusion MRI.

RESULTS

Both simulation and in-vivo results demonstrate the effectiveness of the proposed acquisition and reconstruction framework to provide distortion-corrected diffusion imaging at the mesoscale with markedly reduced T2*-blurring. The in-vivo results of 720 µm and 500 µm datasets show high-fidelity diffusion images with reduced image blurring and echo time using the proposed approaches.

CONCLUSIONS

The proposed method provides high-quality distortion-corrected diffusion-weighted images with ∼40% reduction in the echo-train-length and T2* blurring at 500µm-isotropic-resolution compared to standard multi-shot EPI.

Multi-tensor diffusion abnormalities of gray matter in an animal model of cortical dysplasia

Focal cortical dysplasias are a type of malformations of cortical development that are a common cause of drug-resistant focal epilepsy. Surgical treatment is a viable option for some of these patients, with their outcome being highly related to complete surgical resection of lesions visible in magnetic resonance imaging (MRI). However, subtle lesions often go undetected on conventional imaging. Several methods to analyze MRI have been proposed, with the common goal of rendering subtle cortical lesions visible. However, most image-processing methods are targeted to detect the macroscopic characteristics of cortical dysplasias, which do not always correspond to the microstructural disarrangement of these cortical malformations. Quantitative analysis of diffusion-weighted MRI (dMRI) enables the inference of tissue characteristics, and novel methods provide valuable microstructural features of complex tissue, including gray matter. We investigated the ability of advanced dMRI descriptors to detect diffusion abnormalities in an animal model of cortical dysplasia. For this purpose, we induced cortical dysplasia in 18 animals that were scanned at 30 postnatal days (along with 19 control animals). We obtained multi-shell dMRI, to which we fitted single and multi-tensor representations. Quantitative dMRI parameters derived from these methods were queried using a curvilinear coordinate system to sample the cortical mantle, providing inter-subject anatomical correspondence. We found region- and layer-specific diffusion abnormalities in experimental animals. Moreover, we were able to distinguish diffusion abnormalities related to altered intra-cortical tangential fibers from those associated with radial cortical fibers. Histological examinations revealed myelo-architectural abnormalities that explain the alterations observed through dMRI. The methods for dMRI acquisition and analysis used here are available in clinical settings and our work shows their clinical relevance to detect subtle cortical dysplasias through analysis of their microstructural properties.

Diffusion MRI data analysis assisted by deep learning synthesized anatomical images (DeepAnat)

Diffusion MRI is a useful neuroimaging tool for non-invasive mapping of human brain microstructure and structural connections. The analysis of diffusion MRI data often requires brain segmentation, including volumetric segmentation and cerebral cortical surfaces, from additional high-resolution T1-weighted (T1w) anatomical MRI data, which may be unacquired, corrupted by subject motion or hardware failure, or cannot be accurately co-registered to the diffusion data that are not corrected for susceptibility-induced geometric distortion. To address these challenges, this study proposes to synthesize high-quality T1w anatomical images directly from diffusion data using convolutional neural networks (CNNs) (entitled "DeepAnat"), including a U-Net and a hybrid generative adversarial network (GAN), and perform brain segmentation on synthesized T1w images or assist the co-registration using synthesized T1w images. The quantitative and systematic evaluations using data of 60 young subjects provided by the Human Connectome Project (HCP) show that the synthesized T1w images and results for brain segmentation and comprehensive diffusion analysis tasks are highly similar to those from native T1w data. The brain segmentation accuracy is slightly higher for the U-Net than the GAN. The efficacy of DeepAnat is further validated on a larger dataset of 300 more elderly subjects provided by the UK Biobank. Moreover, the U-Nets trained and validated on the HCP and UK Biobank data are shown to be highly generalizable to the diffusion data from Massachusetts General Hospital Connectome Diffusion Microstructure Dataset (MGH CDMD) acquired with different hardware systems and imaging protocols and therefore can be used directly without retraining or with fine-tuning for further improved performance. Finally, it is quantitatively demonstrated that the alignment between native T1w images and diffusion images uncorrected for geometric distortion assisted by synthesized T1w images substantially improves upon that by directly co-registering the diffusion and T1w images using the data of 20 subjects from MGH CDMD. In summary, our study demonstrates the benefits and practical feasibility of DeepAnat for assisting various diffusion MRI data analyses and supports its use in neuroscientific applications.

Column-based cortical depth analysis of the diffusion anisotropy and radiality in submillimeter whole-brain diffusion tensor imaging of the human cortical gray matter in vivo

High-resolution diffusion tensor imaging (DTI) can noninvasively probe the microstructure of cortical gray matter in vivo. In this study, 0.9-mm isotropic whole-brain DTI data were acquired in healthy subjects with an efficient multi-band multi-shot echo-planar imaging sequence. A column-based analysis that samples the fractional anisotropy (FA) and radiality index (RI) along radially oriented cortical columns was then performed to quantitatively analyze the FA and RI dependence on the cortical depth, cortical region, cortical curvature, and cortical thickness across the whole brain, which has not been simultaneously and systematically investigated in previous studies. The results showed characteristic FA and RI vs. cortical depth profiles, with an FA local maximum and minimum (or two inflection points) and a single RI maximum at intermediate cortical depths in most cortical regions, except for the postcentral gyrus where no FA peaks and a lower RI were observed. These results were consistent between repeated scans from the same subjects and across different subjects. They were also dependent on the cortical curvature and cortical thickness in that the characteristic FA and RI peaks were more pronounced i) at the banks than at the crown of gyri or at the fundus of sulci and ii) as the cortical thickness increases. This methodology can help characterize variations in microstructure along the cortical depth and across the whole brain in vivo, potentially providing quantitative biomarkers for neurological disorders.

Spatiotemporal patterns of cortical microstructural maturation in children and adolescents with diffusion MRI

Neocortical maturation is a dynamic process that proceeds in a hierarchical manner; however, the spatiotemporal organization of cortical microstructure with diffusion MRI has yet to be fully defined. This study characterized cortical microstructural maturation using diffusion MRI (fwe-DTI and NODDI multi-compartment modeling) in a cohort of 637 children and adolescents between 8 and 21 years of age. We found spatially heterogeneous developmental patterns broadly demarcated into functional domains where NODDI metrics increased and fwe-DTI metrics decreased with age. Using non-negative matrix factorization, we found cortical regions that correspond to lower-order sensory regions mature earlier than higher-order association regions. Our findings corroborate previous histological and neuroimaging studies that show spatially-varying patterns of cortical maturation that may reflect unique developmental processes of cytoarchitectonically-determined regional patterns of change.

COnstrained Reference frame diffusion TEnsor Correlation Spectroscopic (CORTECS) MRI: A practical framework for high-resolution diffusion tensor distribution imaging

High-resolution imaging studies have consistently shown that in cortical tissue water diffuses preferentially along radial and tangential orientations with respect to the cortical surface, in agreement with histology. These dominant orientations do not change significantly even if the relative contributions from microscopic water pools to the net voxel signal vary across experiments that use different diffusion times, b-values, TEs, and TRs. With this in mind, we propose a practical new framework for imaging non-parametric diffusion tensor distributions (DTDs) by constraining the microscopic diffusion tensors of the DTD to be diagonalized using the same orthonormal reference frame of the mesoscopic voxel. In each voxel, the constrained DTD (cDTD) is completely determined by the correlation spectrum of the microscopic principal diffusivities associated with the axes of the voxel reference frame. Consequently, all cDTDs are inherently limited to the domain of positive definite tensors and can be reconstructed efficiently using Inverse Laplace Transform methods. Moreover, the cDTD reconstruction can be performed using only data acquired efficiently with single diffusion encoding, although it also supports datasets with multiple diffusion encoding. In tissues with a well-defined architecture, such as the cortex, we can further constrain the cDTD to contain only cylindrically symmetric diffusion tensors and measure the 2D correlation spectra of principal diffusivities along the radial and tangential orientation with respect to the cortical surface. To demonstrate this framework, we perform numerical simulations and analyze high-resolution dMRI data from a fixed macaque monkey brain. We estimate 2D cDTDs in the cortex and derive, in each voxel, the marginal distributions of the microscopic principal diffusivities, the corresponding distributions of the microscopic fractional anisotropies and mean diffusivities along with their 2D correlation spectra to quantify the cDTD shape-size characteristics. Signal components corresponding to specific bands in these cDTD-derived spectra show high specificity to cortical laminar structures observed with histology. Our framework drastically simplifies the measurement of non-parametric DTDs in high-resolution datasets with mesoscopic voxel sizes much smaller than the radius of curvature of the underlying anatomy, e.g., cortical surface, and can be applied retrospectively to analyze existing diffusion MRI data from fixed cortical tissues.

High-resolution cortical MAP-MRI reveals areal borders and laminar substructures observed with histological staining

The variations in cellular composition and tissue architecture measured with histology provide the biological basis for partitioning the brain into distinct cytoarchitectonic areas and for characterizing neuropathological tissue alterations. Clearly, there is an urgent need to develop whole-brain neuroradiological methods that can assess cortical cyto- and myeloarchitectonic features non-invasively. Mean apparent propagator (MAP) MRI is a clinically feasible diffusion MRI method that quantifies efficiently and comprehensively the net microscopic displacements of water molecules diffusing in tissues. We investigate the sensitivity of high-resolution MAP-MRI to detecting areal and laminar variations in cortical cytoarchitecture and compare our results with observations from corresponding histological sections in the entire brain of a rhesus macaque monkey. High-resolution images of MAP-derived parameters, in particular the propagator anisotropy (PA), non-gaussianity (NG), and the return-to-axis probability (RTAP) reveal cortical area-specific lamination patterns in good agreement with the corresponding histological stained sections. In a few regions, the MAP parameters provide superior contrast to the five histological stains used in this study, delineating more clearly boundaries and transition regions between cortical areas and laminar substructures. Throughout the cortex, various MAP parameters can be used to delineate transition regions between specific cortical areas observed with histology and to refine areal boundaries estimated using atlas registration-based cortical parcellation. Using surface-based analysis of MAP parameters we quantify the cortical depth dependence of diffusion propagators in multiple regions-of-interest in a consistent and rigorous manner that is largely independent of the cortical folding geometry. The ability to assess cortical cytoarchitectonic features efficiently and non-invasively, its clinical feasibility, and translatability make high-resolution MAP-MRI a promising 3D imaging tool for studying whole-brain cortical organization, characterizing abnormal cortical development, improving early diagnosis of neurodegenerative diseases, identifying targets for biopsies, and complementing neuropathological investigations.

Comparison of individualized behavioral predictions across anatomical, diffusion and functional connectivity MRI

A fundamental goal across the neurosciences is the characterization of relationships linking brain anatomy, functioning, and behavior. Although various MRI modalities have been developed to probe these relationships, direct comparisons of their ability to predict behavior have been lacking. Here, we compared the ability of anatomical T1, diffusion and functional MRI (fMRI) to predict behavior at an individual level. Cortical thickness, area and volume were extracted from anatomical T1 images. Diffusion Tensor Imaging (DTI) and approximate Neurite Orientation Dispersion and Density Imaging (NODDI) models were fitted to the diffusion images. The resulting metrics were projected to the Tract-Based Spatial Statistics (TBSS) skeleton. We also ran probabilistic tractography for the diffusion images, from which we extracted the stream count, average stream length, and the average of each DTI and NODDI metric across tracts connecting each pair of brain regions. Functional connectivity (FC) was extracted from both task and resting-state fMRI. Individualized prediction of a wide range of behavioral measures were performed using kernel ridge regression, linear ridge regression and elastic net regression. Consistency of the results were investigated with the Human Connectome Project (HCP) and Adolescent Brain Cognitive Development (ABCD) datasets. In both datasets, FC-based models gave the best prediction performance, regardless of regression model or behavioral measure. This was especially true for the cognitive component. Furthermore, all modalities were able to predict cognition better than other behavioral components. Combining all modalities improved prediction of cognition, but not other behavioral components. Finally, across all behaviors, combining resting and task FC yielded prediction performance similar to combining all modalities. Overall, our study suggests that in the case of healthy children and young adults, behaviorally-relevant information in T1 and diffusion features might reflect a subset of the variance captured by FC.

Inter- and intra-individual variation in brain structural-cognition relationships in aging

The sources of inter- and intra-individual variability in age-related cognitive decline remain poorly understood. We examined the association between 20-year trajectories of cognitive decline and multimodal brain structure and morphology in older age. We used the Whitehall II Study, an extensively characterised cohort with 3T brain magnetic resonance images acquired at older age (mean age = 69.52 ± 4.9) and 5 repeated cognitive performance assessments between mid-life (mean age = 53.2 ±4.9 years) and late-life (mean age = 67.7 ± 4.9). Using non-negative matrix factorization, we identified 10 brain components integrating cortical thickness, surface area, fractional anisotropy, and mean and radial diffusivities. We observed two latent variables describing distinct brain-cognition associations. The first describes variations in 5 structural components associated with low mid-life performance across multiple cognitive domains, decline in reasoning, but maintenance of fluency abilities. The second describes variations in 6 structural components associated with low mid-life performance in fluency and memory, but retention of multiple abilities. Expression of latent variables predicts future cognition 3.2 years later (mean age = 70.87 ± 4.9). This data-driven approach highlights brain-cognition relationships wherein individuals degrees of cognitive decline and maintenance across diverse cognitive functions are both positively and negatively associated with markers of cortical structure.

High B-value diffusion tensor imaging for early detection of hippocampal microstructural alteration in a mouse model of multiple sclerosis

Several studies have highlighted the value of diffusion tensor imaging (DTI) with strong diffusion weighting to reveal white matter microstructural lesions, but data in gray matter (GM) remains scarce. Herein, the effects of b-values combined with different numbers of diffusion-encoding directions (NDIRs) on DTI metrics to capture the normal hippocampal microstructure and its early alterations were investigated in a mouse model of multiple sclerosis (experimental autoimmune encephalomyelitis [EAE]). Two initial DTI datasets (B2700-43Dir acquired with b = 2700 s.mm⁻² and NDIR = 43; B1000-22Dir acquired with b = 1000 s.mm⁻² and NDIR = 22) were collected from 18 normal and 18 EAE mice at 4.7 T. Three additional datasets (B2700-22Dir, B2700-12Dir and B1000-12Dir) were extracted from the initial datasets. In healthy mice, we found a significant influence of b-values and NDIR on all DTI metrics. Confronting unsupervised hippocampal layers classification to the true anatomical classification highlighted the remarkable discrimination of the molecular layer with B2700-43Dir compared with the other datasets. Only DTI from the B2700 datasets captured the dendritic loss occurring in the molecular layer of EAE mice. Our findings stress the needs for both high b-values and sufficient NDIR to achieve a GM DTI with more biologically meaningful correlations, though DTI-metrics should be interpreted with caution in these settings.

In vivo irreversible and reversible transverse relaxation rates in human cerebral cortex via line scans at 7 T with 250 micron resolution perpendicular to the cortical surface

Understanding how and why MR signals and their associated relaxation rates vary with cortical depth could ultimately enable the noninvasive investigation of the laminar architecture of cerebral cortex in the living human brain. However, cortical gray matter is typically only a few millimeters thick, making it challenging to sample many cortical depths with the voxel sizes commonly used in MRI studies. Line-scan techniques provide a way to overcome this challenge and here we implemented a novel line-scan GESSE pulse sequence that allowed us to measure irreversible and reversible transverse relaxation rates-R₂ and R₂´, respectively-with extremely high resolution (250 μm) in the radial direction, perpendicular to the cortical surface. Eight healthy human subjects were scanned at 7 T using this sequence, with primary visual cortex (V1) targeted in three subjects and primary motor (M1) and somatosensory cortex (S1) targeted in the other five. In all three cortical areas, a peak in R₂ values near the central depths was seen consistently across subjects-an observation that has not been made before, to our knowledge. On the other hand, no consistent pattern was apparent for R₂´ values as a function of cortical depth. The intracortical R₂ peak reported here is unlikely to be explained by myelin content or by deoxyhemoglobin in the microvasculature; however, this peak is in accord with the laminar distribution of non-heme iron in these cortical areas, known from prior histology studies. Obtaining information about tissue microstructure via measurements of transverse relaxation (and other quantitative MR contrast mechanisms) at the extremely high radial resolutions achievable through the use of line-scan techniques could therefore bring us closer to being able to perform "in vivo histology" of the cerebral cortex.

High-resolution neural network-driven mapping of multiple diffusion metrics leveraging asymmetries in the balanced steady-state free precession frequency profile

We propose to utilize the rich information content about microstructural tissue properties entangled in asymmetric balanced steady-state free precession (bSSFP) profiles to estimate multiple diffusion metrics simultaneously by neural network (NN) parameter quantification. A 12-point bSSFP phase-cycling scheme with high-resolution whole-brain coverage is employed at 3 and 9.4 T for NN input. Low-resolution target diffusion data are derived based on diffusion-weighted spin-echo echo-planar-imaging (SE-EPI) scans, that is, mean, axial, and radial diffusivity (MD, AD, and RD), fractional anisotropy (FA), as well as the spherical coordinates (azimuth Φ and inclination ϴ) of the principal diffusion eigenvector. A feedforward NN is trained with incorporated probabilistic uncertainty estimation. The NN predictions yielded highly reliable results in white matter (WM) and gray matter structures for MD. The quantification of FA, AD, and RD was overall in good agreement with the reference but the dependence of these parameters on WM anisotropy was somewhat biased (e.g. in corpus callosum). The inclination ϴ was well predicted for anisotropic WM structures, while the azimuth Φ was overall poorly predicted. The findings were highly consistent across both field strengths. Application of the optimized NN to high-resolution input data provided whole-brain maps with rich structural details. In conclusion, the proposed NN-driven approach showed potential to provide distortion-free high-resolution whole-brain maps of multiple diffusion metrics at high to ultrahigh field strengths in clinically relevant scan times.

Brainstem Diffusion Tensor Tractography and Clinical Applications in Pain

The brainstem is one of the most vulnerable brain structures in many neurological conditions, such as pain, sleep problems, autonomic dysfunctions, and neurodegenerative disorders. Diffusion tensor imaging and tractography provide structural details and quantitative measures of brainstem fiber pathways. Until recently, diffusion tensor tractographic studies have mainly focused on whole-brain MRI acquisition. Due to the brainstem's spatial localization, size, and tissue characteristics, and limits of imaging techniques, brainstem diffusion MRI poses particular challenges in tractography. We provide a brief overview on recent advances in diffusion tensor tractography in revealing human pathways connecting the brainstem to the subcortical regions (e.g., basal ganglia, mesolimbic, basal forebrain), and cortical regions. Each of these pathways contains different distributions of fiber tracts from known neurotransmitter-specific nuclei in the brainstem. We compare the brainstem tractographic approaches in literature and our in-lab developed automated brainstem tractography in terms of atlas building, technical advantages, and neuroanatomical implications on neurotransmitter systems. Lastly, we summarize recent investigations of using brainstem tractography as a promising tool in association with pain.

In vivo Correlation Tensor MRI reveals microscopic kurtosis in the human brain on a clinical 3T scanner

Diffusion MRI (dMRI) has become one of the most important imaging modalities for noninvasively probing tissue microstructure. Diffusional Kurtosis MRI (DKI) quantifies the degree of non-gaussian diffusion, which in turn has been shown to increase sensitivity towards, e.g., disease and orientation mapping in neural tissue. However, the specificity of DKI is limited as different sources can contribute to the total intravoxel diffusional kurtosis, including: variance in diffusion tensor magnitudes (K_iso), variance due to diffusion anisotropy (K_aniso), and microscopic kurtosis (μK) related to restricted diffusion, microstructural disorder, and/or exchange. Interestingly, μK is typically ignored in diffusion MRI signal modeling as it is assumed to be negligible in neural tissues. However, recently, Correlation Tensor MRI (CTI) based on Double-Diffusion-Encoding (DDE) was introduced for kurtosis source separation, revealing non negligible μK in preclinical imaging. Here, we implemented CTI for the first time on a clinical 3T scanner and investigated the sources of total kurtosis in healthy subjects. A robust framework for kurtosis source separation in humans is introduced, followed by estimation of μK (and the other kurtosis sources) in the healthy brain. Using this clinical CTI approach, we find that μK significantly contributes to total diffusional kurtosis both in gray and white matter tissue but, as expected, not in the ventricles. The first μK maps of the human brain are presented, revealing that the spatial distribution of μK provides a unique source of contrast, appearing different from isotropic and anisotropic kurtosis counterparts. Moreover, group average templates of these kurtosis sources have been generated for the first time, which corroborated our findings at the underlying individual-level maps. We further show that the common practice of ignoring μK and assuming the multiple gaussian component approximation for kurtosis source estimation introduces significant bias in the estimation of other kurtosis sources and, perhaps even worse, compromises their interpretation. Finally, a twofold acceleration of CTI is discussed in the context of potential future clinical applications. We conclude that CTI has much potential for future in vivo microstructural characterizations in healthy and pathological tissue.

SDnDTI: Self-supervised deep learning-based denoising for diffusion tensor MRI

Diffusion tensor magnetic resonance imaging (DTI) is a widely adopted neuroimaging method for the in vivo mapping of brain tissue microstructure and white matter tracts. Nonetheless, the noise in the diffusion-weighted images (DWIs) decreases the accuracy and precision of DTI derived microstructural parameters and leads to prolonged acquisition time for achieving improved signal-to-noise ratio (SNR). Deep learning-based image denoising using convolutional neural networks (CNNs) has superior performance but often requires additional high-SNR data for supervising the training of CNNs, which reduces the feasibility of supervised learning-based denoising in practice. In this work, we develop a self-supervised deep learning-based method entitled "SDnDTI" for denoising DTI data, which does not require additional high-SNR data for training. Specifically, SDnDTI divides multi-directional DTI data into many subsets of six DWI volumes and transforms DWIs from each subset to along the same diffusion-encoding directions through the diffusion tensor model, generating multiple repetitions of DWIs with identical image contrasts but different noise observations. SDnDTI removes noise by first denoising each repetition of DWIs using a deep 3-dimensional CNN with the average of all repetitions with higher SNR as the training target, following the same approach as normal supervised learning based denoising methods, and then averaging CNN-denoised images for achieving higher SNR. The denoising efficacy of SDnDTI is demonstrated in terms of the similarity of output images and resultant DTI metrics compared to the ground truth generated using substantially more DWI volumes on two datasets with different spatial resolutions, b-values and numbers of input DWI volumes provided by the Human Connectome Project (HCP) and the Lifespan HCP in Aging. The SDnDTI results preserve image sharpness and textural details and substantially improve upon those from the raw data. The results of SDnDTI are comparable to those from supervised learning-based denoising and outperform those from state-of-the-art conventional denoising algorithms including BM4D, AONLM and MPPCA. By leveraging domain knowledge of diffusion MRI physics, SDnDTI makes it easier to use CNN-based denoising methods in practice and has the potential to benefit a wider range of research and clinical applications that require accelerated DTI acquisition and high-quality DTI data for mapping of tissue microstructure, fiber tracts and structural connectivity in the living human brain.

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 dMRI 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 dMRI and where the field is headed in the coming years.

Value of ultra-high field MRI in patients with suspected focal epilepsy and negative 3 T MRI (EpiUltraStudy): protocol for a prospective, longitudinal therapeutic study

PURPOSE

Resective epilepsy surgery is a well-established, evidence-based treatment option in patients with drug-resistant focal epilepsy. A major predictive factor of good surgical outcome is visualization and delineation of a potential epileptogenic lesion by MRI. However, frequently, these lesions are subtle and may escape detection by conventional MRI (≤ 3 T).

METHODS

We present the EpiUltraStudy protocol to address the hypothesis that application of ultra-high field (UHF) MRI increases the rate of detection of structural lesions and functional brain aberrances in patients with drug-resistant focal epilepsy who are candidates for resective epilepsy surgery. Additionally, therapeutic gain will be addressed, testing whether increased lesion detection and tailored resections result in higher rates of seizure freedom 1 year after epilepsy surgery. Sixty patients enroll the study according to the following inclusion criteria: aged ≥ 12 years, diagnosed with drug-resistant focal epilepsy with a suspected epileptogenic focus, negative conventional 3 T MRI during pre-surgical work-up.

RESULTS

All patients will be evaluated by 7 T MRI; ten patients will undergo an additional 9.4 T MRI exam. Images will be evaluated independently by two neuroradiologists and a neurologist or neurosurgeon. Clinical and UHF MRI will be discussed in the multidisciplinary epilepsy surgery conference. Demographic and epilepsy characteristics, along with postoperative seizure outcome and histopathological evaluation, will be recorded.

CONCLUSION

This protocol was reviewed and approved by the local Institutional Review Board and complies with the Declaration of Helsinki and principles of Good Clinical Practice. Results will be submitted to international peer-reviewed journals and presented at international conferences.

TRIAL REGISTRATION NUMBER

www.trialregister.nl : NTR7536.

Common principles in the lateralisation of auditory cortex structure and function for vocal communication in primates and rodents

This review summarises recent findings on the lateralisation of communicative sound processing in the auditory cortex (AC) of humans, non-human primates, and rodents. Functional imaging in humans has demonstrated a left hemispheric preference for some acoustic features of speech, but it is unclear to which degree this is caused by bottom-up acoustic feature selectivity or top-down modulation from language areas. Although non-human primates show a less pronounced functional lateralisation in AC, the properties of AC fields and behavioral asymmetries are qualitatively similar. Rodent studies demonstrate microstructural circuits that might underlie bottom-up acoustic feature selectivity in both hemispheres. Functionally, the left AC in the mouse appears to be specifically tuned to communication calls, whereas the right AC may have a more 'generalist' role. Rodents also show anatomical AC lateralisation, such as differences in size and connectivity. Several of these functional and anatomical characteristics are also lateralized in human AC. Thus, complex vocal communication processing shares common features among rodents and primates. We argue that a synthesis of results from humans, non-human primates, and rodents is necessary to identify the neural circuitry of vocal communication processing. However, data from different species and methods are often difficult to compare. Recent advances may enable better integration of methods across species. Efforts to standardise data formats and analysis tools would benefit comparative research and enable synergies between psychological and biological research in the area of vocal communication processing.

High-fidelity fast volumetric brain MRI using synergistic wave-controlled aliasing in parallel imaging and a hybrid denoising generative adversarial network (HDnGAN)

PURPOSE

The goal of this study is to leverage an advanced fast imaging technique, wave-controlled aliasing in parallel imaging (Wave-CAIPI), and a generative adversarial network (GAN) for denoising to achieve accelerated high-quality high-signal-to-noise-ratio (SNR) volumetric MRI.

METHODS

Three-dimensional (3D) T₂ -weighted fluid-attenuated inversion recovery (FLAIR) image data were acquired on 33 multiple sclerosis (MS) patients using a prototype Wave-CAIPI sequence (acceleration factor R = 3×2, 2.75 minutes) and a standard T₂ -SPACE FLAIR sequence (R = 2, 7.25 minutes). A hybrid denoising GAN entitled "HDnGAN" consisting of a 3D generator and a 2D discriminator was proposed to denoise highly accelerated Wave-CAIPI images. HDnGAN benefits from the improved image synthesis performance provided by the 3D generator and increased training samples from a limited number of patients for training the 2D discriminator. HDnGAN was trained and validated on data from 25 MS patients with the standard FLAIR images as the target and evaluated on data from 8 MS patients not seen during training. HDnGAN was compared to other denoising methods including AONLM, BM4D, MU-Net, and 3D GAN in qualitative and quantitative analysis of output images using the mean squared error (MSE) and VGG perceptual loss compared to standard FLAIR images, and a reader assessment by two neuroradiologists regarding sharpness, SNR, lesion conspicuity, and overall quality. Finally, the performance of these denoising methods was compared at higher noise levels using simulated data with added Rician noise.

RESULTS

HDnGAN effectively denoised low-SNR Wave-CAIPI images with sharpness and rich textural details, which could be adjusted by controlling the contribution of the adversarial loss to the total loss when training the generator. Quantitatively, HDnGAN (λ = 10^-3 ) achieved low MSE and the lowest VGG perceptual loss. The reader study showed that HDnGAN (λ = 10^-3 ) significantly improved the SNR of Wave-CAIPI images (P<0.001), outperformed AONLM (P = 0.015), BM4D (P<0.001), MU-Net (P<0.001) and 3D GAN (λ = 10^-3 ) (P<0.001) regarding image sharpness, and outperformed MU-Net (P<0.001) and 3D GAN (λ = 10^-3 ) (P = 0.001) regarding lesion conspicuity. The overall quality score of HDnGAN (λ = 10^-3 ) (4.25±0.43) was significantly higher than those from Wave-CAIPI (3.69±0.46, P = 0.003), BM4D (3.50±0.71, P = 0.001), MU-Net (3.25±0.75, P<0.001), and 3D GAN (λ = 10^-3 ) (3.50±0.50, P<0.001), with no significant difference compared to standard FLAIR images (4.38±0.48, P = 0.333). The advantages of HDnGAN over other methods were more obvious at higher noise levels.

CONCLUSION

HDnGAN provides robust and feasible denoising while preserving rich textural detail in empirical volumetric MRI data. Our study using empirical patient data and systematic evaluation supports the use of HDnGAN in combination with modern fast imaging techniques such as Wave-CAIPI to achieve high-fidelity fast volumetric MRI and represents an important step to the clinical translation of GANs. This article is protected by copyright. All rights reserved.

Microstructural Properties of Human Brain Revealed by Fractional Anisotropy Can Predict the After-Effect of Intermittent Theta Burst Stimulation

Intermittent theta burst stimulation (iTBS) delivered by transcranial magnetic stimulation (TMS) produces a long-term potentiation-like after-effect useful for investigations of cortical function and of potential therapeutic value. However, the iTBS after-effect over the primary motor cortex (M1) as measured by changes in motor evoked potential (MEP) amplitude exhibits a largely unexplained variability across individuals. Here, we present evidence that individual differences in white matter (WM) and gray matter (GM) microstructural properties revealed by fractional anisotropy (FA) predict the magnitude of the iTBS-induced after-effect over M1. The MEP amplitude change in the early phase (5–10 min post-iTBS) was associated with FA values in WM tracts such as right superior longitudinal fasciculus and corpus callosum. By contrast, the MEP amplitude change in the late phase (15–30 min post-iTBS) was associated with FA in GM, primarily in right frontal cortex. These results suggest that the microstructural properties of regions connected directly or indirectly to the target region (M1) are crucial determinants of the iTBS after-effect. FA values indicative of these microstructural differences can predict the potential effectiveness of repetitive TMS for both investigational use and clinical application.

Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome

The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.

Tractography-Pathology Correlations in Traumatic Brain Injury: A TRACK-TBI Study

Diffusion tractography magnetic resonance imaging (MRI) can infer changes in network connectivity in patients with traumatic brain injury (TBI), but the pathological substrates of disconnected tracts have not been well defined because of a lack of high-resolution imaging with histopathological validation. We developed an ex vivo MRI protocol to analyze tract terminations at 750-μm isotropic resolution, followed by histopathological evaluation of white matter pathology, and applied these methods to a 60-year-old man who died 26 days after TBI. Analysis of 74 cerebral hemispheric white matter regions revealed a heterogeneous distribution of tract disruptions. Associated histopathology identified variable white matter injury with patchy deposition of amyloid precursor protein (APP), loss of neurofilament-positive axonal processes, myelin dissolution, astrogliosis, microgliosis, and perivascular hemosiderin-laden macrophages. Multiple linear regression revealed that tract disruption strongly correlated with the density of APP-positive axonal swellings and neurofilament loss. Ex vivo diffusion MRI can detect tract disruptions in the human brain that reflect axonal injury.

Improving distortion correction for isotropic high-resolution 3D diffusion MRI by optimizing Jacobian modulation

PURPOSE

To improve distortion correction for isotropic high-resolution whole-brain 3D diffusion MRI when in a time-saving acquisition scenario.

THEORY AND METHODS

Data were acquired using simultaneous multi-slab (SMSlab) acquisitions, with a b = 0 image pair encoded by reversed polarity gradients (RPG) for phase encoding (PE) and diffusion weighted images encoded by a single PE direction. Eddy current-induced distortions were corrected first. During the following susceptibility distortion correction, image deformation was first corrected by the field map estimated from the b = 0 image pair. Intensity variation was subsequently corrected by Jacobian modulation. Two Jacobian modulation methods were compared. They calculated the Jacobian modulation map from the field map, or from the deformation corrected b = 0 image pair, termed as J_Field and J_RPG , respectively. A modified version of the J_RPG method, with proper smoothing, was further proposed for improved correction performance, termed as J_RPG-smooth .

RESULTS

Compared to J_Field modulation, less remaining distortions are observed when using the J_RPG and J_RPG-smooth methods, especially in areas with large B0 field inhomogeneity. The original J_RPG method causes signal-to-noise ratio (SNR) deficiency problem, which manifests as degraded SNR of the diffusion weighted images, while the J_RPG-smooth method maintains the original image SNR. Less estimation errors of diffusion metrics are observed when using the J_RPG-smooth method.

CONCLUSION

This study improves the distortion correction for isotropic high-resolution whole-brain 3D diffusion MRI by optimizing Jacobian modulation. The optimized method outperforms the conventional J_Field method regarding intensity variation correction and accuracy of diffusion metrics estimation, and outperforms the original J_RPG method regarding SNR performance.

High-resolution whole-brain diffusion MRI at 3T using simultaneous multi-slab (SMSlab) acquisition

High-resolution diffusion MRI (dMRI) is a crucial tool in neuroscience studies to detect fine fiber structure, depict complex fiber architecture and analyze cortical anisotropy. However, high-resolution dMRI is limited by its intrinsically low SNR due to diffusion attenuation. A series of techniques have been proposed to improve the SNR performance, but most of them are at the cost of long scan time, which in turn sacrifice the SNR efficiency, especially for large FOV imaging, such as whole-brain imaging. Recently, a combination of 3D multi-slab acquisition and simultaneous multi-slice (SMS) excitation, namely simultaneous multi-slab (SMSlab), has been demonstrated to have potential for high-resolution diffusion imaging with high SNR and SNR efficiency. In our previous work, we have proposed a 3D Fourier encoding and reconstruction framework for SMSlab acquisition. In this study, we extend this 3D k-space framework to diffusion imaging, by developing a novel navigator acquisition strategy and exploring a k-space-based phase correction method. In vivo brain data are acquired using the proposed SMSlab method and compared with a series of different acquisitions, including the traditional 3D multi-slab, 2D SMS and 2D single-shot EPI (ss-EPI) acquisitions. The results demonstrate that SMSlab has a better SNR performance compared with 3D multi-slab and 2D SMS. The detection capacity of fine fiber structures is improved using SMSlab, compared with the low-resolution diffusion imaging using conventional 2D ss-EPI.

In vivo human whole-brain Connectom diffusion MRI dataset at 760 µm isotropic resolution

We present a whole-brain in vivo diffusion MRI (dMRI) dataset acquired at 760 μm isotropic resolution and sampled at 1260 q-space points across 9 two-hour sessions on a single healthy participant. The creation of this benchmark dataset is possible through the synergistic use of advanced acquisition hardware and software including the high-gradient-strength Connectom scanner, a custom-built 64-channel phased-array coil, a personalized motion-robust head stabilizer, a recently developed SNR-efficient dMRI acquisition method, and parallel imaging reconstruction with advanced ghost reduction algorithm. With its unprecedented resolution, SNR and image quality, we envision that this dataset will have a broad range of investigational, educational, and clinical applications that will advance the understanding of human brain structures and connectivity. This comprehensive dataset can also be used as a test bed for new modeling, sub-sampling strategies, denoising and processing algorithms, potentially providing a common testing platform for further development of in vivo high resolution dMRI techniques. Whole brain anatomical T1-weighted and T2-weighted images at submillimeter scale along with field maps are also made available.

Automated cerebral cortex segmentation based solely on diffusion tensor imaging for investigating cortical anisotropy

To extract Diffusion Tensor Imaging (DTI) parameters from the human cortex, the inner and outer boundaries of the cortex are usually defined on 3D-T1-weighted images and then applied to the co-registered DTI. However, this analysis requires the acquisition of an additional high-resolution structural image that may not be practical in various imaging studies. Here an automatic cortical boundary segmentation method was developed to work directly only on the native DTI images by using fractional anisotropy (FA) maps and mean diffusion weighted images (DWI), the latter with acceptable gray-white matter image contrast. This new method was compared to the conventional cortical segmentations generated from high-resolution T1 structural images in 5 participants. In addition, the proposed method was applied to 15 healthy young adults (10 cross-sectional, 5 test-retest) to measure FA, MD, and radiality of the primary eigenvector across the cortex on whole-brain 1.5 mm isotropic images acquired in 3.5 min at 3T. The proposed method generated reasonable segmentations of the cortical boundaries for all individuals and large proportions of the proposed method segmentations (more than 85%) were within ±1 mm from those generated with the conventional approach on higher resolution T1 structural images. Both FA (~0.15) and MD (~0.77 × 10^-3 mm²/s) extracted halfway between the cortical boundaries were relatively stable across the cortex, although focal regions such as the posterior bank of the central sulcus, anterior insula, and medial temporal lobe showed higher FA. The primary eigenvectors were primarily oriented radially to the middle cortical surface, but there were tangential orientations in the sulcal fundi as well as in the posterior bank of the central sulcus. The proposed method demonstrates the feasibility and accuracy of cortical analysis in native DTI space while avoiding the acquisition of other imaging contrasts like 3D T1-weighted scans.

Compartmental diffusion and microstructural properties of human brain gray and white matter studied with double diffusion encoding magnetic resonance spectroscopy of metabolites and water

Double diffusion encoding (DDE) of the water signal offers a unique ability to separate the effect of microscopic anisotropic diffusion in structural units of tissue from the overall macroscopic orientational distribution of cells. However, the specificity in detected microscopic anisotropy is limited as the signal is averaged over different cell types and across tissue compartments. Performing side-by-side water and metabolite DDE spectroscopic (DDES) experiments provides complementary measures from which intracellular and extracellular microscopic fractional anisotropies (μFA) and diffusivities can be estimated. Metabolites are largely confined to the intracellular space and therefore provide a benchmark for intracellular μFA and diffusivities of specific cell types. By contrast, water DDES measurements allow examination of the separate contributions to water μFA and diffusivity from the intra- and extracellular spaces, by using a wide range of b values to gradually eliminate the extracellular contribution. Here, we aimed to estimate tissue and compartment specific human brain microstructure by combining water and metabolites DDES experiments. We performed our DDES measurements in two brain regions that contain widely different amounts of white matter (WM) and gray matter (GM): parietal white matter (PWM) and occipital gray matter (OGM) in a total of 20 healthy volunteers at 7 Tesla. Metabolite DDES measurements were performed at b = 7199 s/mm², while water DDES measurements were performed with a range of b values from 918 to 7199 s/mm². The experimental framework we employed here resulted in a set of insights pertaining to the morphology of the intracellular and extracellular spaces in both gray and white matter. Results of the metabolite DDES experiments in both PWM and OGM suggest a highly anisotropic intracellular space within neurons and glia, with the possible exception of gray matter glia. The water μFA obtained from the DDES results at high b values in both regions converged with that of the metabolite DDES, suggesting that the signal from the extracellular space is indeed effectively suppressed at the highest b value. The μFA measured in the OGM significantly decreased at lower b values, suggesting a considerably lower anisotropy of the extracellular space in GM compared to WM. In PWM, the water μFA remained high even at the lowest b value, indicating a high degree of organization in the interstitial space in WM. Tortuosity values in the cytoplasm for water and tNAA, obtained with correlation analysis of microscopic parallel diffusivity with respect to GM/WM tissue fraction in the volume of interest, are remarkably similar for both molecules, while exhibiting a clear difference between gray and white matter, suggesting a more crowded cytoplasm and more complex cytomorphology of neuronal cell bodies and dendrites in GM than those found in long-range axons in WM.

Brain tissues have single-voxel signatures in multi-spectral MRI

Since the seminal works by Brodmann and contemporaries, it is well-known that different brain regions exhibit unique cytoarchitectonic and myeloarchitectonic features. Transferring the approach of classifying brain tissues - and other tissues - based on their intrinsic features to the realm of magnetic resonance (MR) is a longstanding endeavor. In the 1990s, atlas-based segmentation replaced earlier multi-spectral classification approaches because of the large overlap between the class distributions. Here, we explored the feasibility of performing global brain classification based on intrinsic MR features, and used several technological advances: ultra-high field MRI, q-space trajectory diffusion imaging revealing voxel-intrinsic diffusion properties, chemical exchange saturation transfer and semi-solid magnetization transfer imaging as a marker of myelination and neurochemistry, and current neural network architectures to analyze the data. In particular, we used the raw image data as well to increase the number of input features. We found that a global brain classification of roughly 97 brain regions was feasible with gross classification accuracy of 60%; and that mapping from voxel-intrinsic MR data to the brain region to which the data belongs is possible. This indicates the presence of unique MR signals of different brain regions, similar to their cytoarchitectonic and myeloarchitectonic fingerprints.

Multiple b-values improve discrimination of cortical gray matter regions using diffusion MRI: an experimental validation with a data-driven approach

Objective

To investigate whether varied or repeated b-values provide better diffusion MRI data for discriminating cortical areas with a data-driven approach.

Methods

Data were acquired from three volunteers at 1.5T with b-values of 800, 1400, 2000 s/mm² along 64 diffusion-encoding directions. The diffusion signal was sampled from gray matter in seven regions of interest (ROIs). Rotational invariants of the local diffusion profile were extracted as features that characterize local tissue properties. Random forest classification experiments assessed whether classification accuracy improved when data with multiple b-values were used over repeated acquisition of the same (1400 s/mm²) b-value to compare all possible pairs of the seven ROIs. Three data sets from the Human Connectome Project were subjected to similar processing and analysis pipelines in eight ROIs.

Results

Three different b-values showed an average improvement in correct classification rates of 5.6% and 4.6%, respectively, in the local and HCP data over repeated measurements of the same b-value. The improvement in correct classification rate reached as high as 16% for individual binary classification experiments between two ROIs. Often using only two of the available three b-values were adequate to make such an improvement in classification rates.

Conclusion

Acquisitions with varying b-values are more suitable for discriminating cortical areas.

Distortion-free, high-isotropic-resolution diffusion MRI with gSlider BUDA-EPI and multicoil dynamic B0 shimming

PURPOSE

We combine SNR-efficient acquisition and model-based reconstruction strategies with newly available hardware instrumentation to achieve distortion-free in vivo diffusion MRI of the brain at submillimeter-isotropic resolution with high fidelity and sensitivity on a clinical 3T scanner.

METHODS

We propose blip-up/down acquisition (BUDA) for multishot EPI using interleaved blip-up/blip-down phase encoding and incorporate B₀ forward-modeling into structured low-rank reconstruction to enable distortion-free and navigator-free diffusion MRI. We further combine BUDA-EPI with an SNR-efficient simultaneous multislab acquisition (generalized slice-dithered enhanced resolution ["gSlider"]), to achieve high-isotropic-resolution diffusion MRI. To validate gSlider BUDA-EPI, whole-brain diffusion data at 860-μm and 780-μm data sets were acquired. Finally, to improve the conditioning and minimize noise penalty in BUDA reconstruction at very high resolutions where B₀ inhomogeneity can have a detrimental effect, the level of B₀ inhomogeneity was reduced by incorporating slab-by-slab dynamic shimming with a 32-channel AC/DC coil into the acquisition. Whole-brain 600-μm diffusion data were then acquired with this combined approach of gSlider BUDA-EPI with dynamic shimming.

RESULTS

The results of 860-μm and 780-μm datasets show high geometry fidelity with gSlider BUDA-EPI. With dynamic shimming, the BUDA reconstruction's noise penalty was further alleviated. This enables whole-brain 600-μm isotropic resolution diffusion imaging with high image quality.

CONCLUSIONS

The gSlider BUDA-EPI method enables high-quality, distortion-free diffusion imaging across the whole brain at submillimeter resolution, where the use of multicoil dynamic B₀ shimming further improves reconstruction performance, which can be particularly useful at very high resolutions.

Improved cortical surface reconstruction using sub-millimeter resolution MPRAGE by image denoising

Automatic cerebral cortical surface reconstruction is a useful tool for cortical anatomy quantification, analysis and visualization. Recently, the Human Connectome Project and several studies have shown the advantages of using T₁-weighted magnetic resonance (MR) images with sub-millimeter isotropic spatial resolution instead of the standard 1-mm isotropic resolution for improved accuracy of cortical surface positioning and thickness estimation. Nonetheless, sub-millimeter resolution images are noisy by nature and require averaging multiple repetitions to increase the signal-to-noise ratio for precisely delineating the cortical boundary. The prolonged acquisition time and potential motion artifacts pose significant barriers to the wide adoption of cortical surface reconstruction at sub-millimeter resolution for a broad range of neuroscientific and clinical applications. We address this challenge by evaluating the cortical surface reconstruction resulting from denoised single-repetition sub-millimeter T₁-weighted images. We systematically characterized the effects of image denoising on empirical data acquired at 0.6 mm isotropic resolution using three classical denoising methods, including denoising convolutional neural network (DnCNN), block-matching and 4-dimensional filtering (BM4D) and adaptive optimized non-local means (AONLM). The denoised single-repetition images were found to be highly similar to 6-repetition averaged images, with a low whole-brain averaged mean absolute difference of ~0.016, high whole-brain averaged peak signal-to-noise ratio of ~33.5 dB and structural similarity index of ~0.92, and minimal gray matter–white matter contrast loss (2% to 9%). The whole-brain mean absolute discrepancies in gray matter–white matter surface placement, gray matter–cerebrospinal fluid surface placement and cortical thickness estimation were lower than 165 μm, 155 μm and 145 μm—sufficiently accurate for most applications. These discrepancies were approximately one third to half of those from 1-mm isotropic resolution data. The denoising performance was equivalent to averaging ~2.5 repetitions of the data in terms of image similarity, and 1.6–2.2 repetitions in terms of the cortical surface placement accuracy. The scan-rescan variability of the cortical surface positioning and thickness estimation was lower than 170 μm. Our unique dataset and systematic characterization support the use of denoising methods for improved cortical surface reconstruction at sub-millimeter resolution.

Ultra-high field magnetic resonance imaging in human epilepsy: A systematic review

RATIONALE

Resective epilepsy surgery is an evidence-based curative treatment option for patients with drug-resistant focal epilepsy. The major preoperative predictor of a good surgical outcome is detection of an epileptogenic lesion by magnetic resonance imaging (MRI). Application of ultra-high field (UHF) MRI, i.e. field strengths ≥ 7 Tesla (T), may increase the sensitivity to detect such a lesion.

METHODS

A keyword search strategy was submitted to Pubmed, EMBASE, Cochrane Database and clinicaltrials.gov to select studies on UHF MRI in patients with epilepsy. Follow-up study selection and data extraction were performed following PRISMA guidelines. We focused on I) diagnostic gain of UHF- over conventional MRI, II) concordance of MRI-detected lesion, seizure onset zone and surgical decision-making, and III) postoperative histopathological diagnosis and seizure outcome.

RESULTS

Sixteen observational cohort studies, all using 7T MRI were included. Diagnostic gain of 7T over conventional MRI ranged from 8% to 67%, with a pooled gain of 31%. Novel techniques to visualize pathological processes in epilepsy and lesion detection are discussed. Seizure freedom was achieved in 73% of operated patients; no seizure outcome comparison was made between 7T MRI positive, 7T negative and 3T positive patients. 7T could influence surgical decision-making, with high concordance of lesion and seizure onset zone. Focal cortical dysplasia (54%), hippocampal sclerosis (12%) and gliosis (8.1%) were the most frequently diagnosed histopathological entities.

SIGNIFICANCE

UHF MRI increases, yet variably, the sensitivity to detect an epileptogenic lesion, showing potential for use in clinical practice. It remains to be established whether this results in improved seizure outcome after surgical treatment. Prospective studies with larger cohorts of epilepsy patients, uniform scan and sequence protocols, and innovative post-processing technology are equally important as further increasing field strengths. Besides technical ameliorations, improved correlation of imaging features with clinical semiology, histopathology and clinical outcome has to be established.

Diffusion MRI revealed altered inter-hippocampal projections in the mouse brain after intrauterine inflammation

Diffusion MRI (dMRI) is commonly used to map large axonal pathways in the white matter. Recent technical advances have also enabled dMRI to resolve the small and complex axonal and dendritic projections in the gray matter. This study investigated whether high-resolution dMRI can resolve the hippocampal neuronal projections and detect abnormal connections due to neurological injury. We performed 3D high spatial and angular resolution dMRI of the mouse brains of the offspring survivors from a model of intrauterine (UI) inflammation, who had known functional deficiency in the hippocampus. We used a novel hippocampal connection mapping method to quantify the intra- and inter-hippocampal projections among 34 automatically segmented hippocampal sub-regions. The results demonstrated wide-spread intra-hippocampal projections, but rather specific intra-hippocampal projections that primarily connected through the CA3 region. Compared with the control group (n = 9), UI-injured mice (n = 11) exhibited significantly reduced inter-hippocampal projection strength (p < 0.01), which correlated well with the neurobehavioral assessments (R² = 0.47). Furthermore, using a whole-brain fixel-based analysis, we identified reduced fiber-density in the CA3 and the ventral hippocampal commissure of the UI-injured mice, which may explain the reduced inter-hippocampal projections. Histological findings also indicated reduced commissural fibers due to the UI-injury. Our study suggested that the dMRI-based connectivity mapping technique can potentially characterize abnormal hippocampal projections in neurological disorders.

Improving in vivo human cerebral cortical surface reconstruction using data-driven super-resolution

Accurate and automated reconstruction of the in vivo human cerebral cortical surface from anatomical magnetic resonance (MR) images facilitates the quantitative analysis of cortical structure. Anatomical MR images with sub-millimeter isotropic spatial resolution improve the accuracy of cortical surface and thickness estimation compared to the standard 1-millimeter isotropic resolution. Nonetheless, sub-millimeter resolution acquisitions require averaging multiple repetitions to achieve sufficient signal-to-noise ratio and are therefore long and potentially vulnerable to subject motion. We address this challenge by synthesizing sub-millimeter resolution images from standard 1-millimeter isotropic resolution images using a data-driven supervised machine learning-based super-resolution approach achieved via a deep convolutional neural network. We systematically characterize our approach using a large-scale simulated dataset and demonstrate its efficacy in empirical data. The super-resolution data provide improved cortical surfaces similar to those obtained from native sub-millimeter resolution data. The whole-brain mean absolute discrepancy in cortical surface positioning and thickness estimation is below 100 μm at the single-subject level and below 50 μm at the group level for the simulated data, and below 200 μm at the single-subject level and below 100 μm at the group level for the empirical data, making the accuracy of cortical surfaces derived from super-resolution sufficient for most applications.

Diffusion kurtosis imaging of gray matter in young adults with autism spectrum disorder

Prior ex vivo histological postmortem studies of autism spectrum disorder (ASD) have shown gray matter microstructural abnormalities, however, in vivo examination of gray matter microstructure in ASD has remained scarce due to the relative lack of non-invasive methods to assess it. The aim of this work was to evaluate the feasibility of employing diffusional kurtosis imaging (DKI) to describe gray matter abnormalities in ASD in vivo. DKI data were examined for 16 male participants with a diagnosis of ASD and IQ>80 and 17 age- and IQ-matched male typically developing (TD) young adults 18-25 years old. Mean (MK), axial (AK), radial (RK) kurtosis and mean diffusivity (MD) metrics were calculated for lobar and sub-lobar regions of interest. Significantly decreased MK, RK, and MD were found in ASD compared to TD participants in the frontal and temporal lobes and several sub-lobar regions previously associated with ASD pathology. In ASD participants, decreased kurtosis in gray matter ROIs correlated with increased repetitive and restricted behaviors and poor social interaction symptoms. Decreased kurtosis in ASD may reflect a pathology associated with a less restrictive microstructural environment such as decreased neuronal density and size, atypically sized cortical columns, or limited dendritic arborizations.

Atypical measures of diffusion at the gray-white matter boundary in autism spectrum disorder in adulthood

Autism spectrum disorder (ASD) is a highly complex neurodevelopmental condition that is accompanied by neuroanatomical differences on the macroscopic and microscopic level. Findings from histological, genetic, and more recently in vivo neuroimaging studies converge in suggesting that neuroanatomical abnormalities, specifically around the gray‐white matter (GWM) boundary, represent a crucial feature of ASD. However, no research has yet characterized the GWM boundary in ASD based on measures of diffusion. Here, we registered diffusion tensor imaging data to the structural T1‐weighted images of 92 adults with ASD and 92 matched neurotypical controls in order to examine between‐group differences and group‐by‐sex interactions in fractional anisotropy and mean diffusivity sampled at the GWM boundary, and at different sampling depths within the superficial white and into the gray matter. As hypothesized, we observed atypical diffusion at and around the GWM boundary in ASD, with between‐group differences and group‐by‐sex interactions depending on tissue class and sampling depth. Furthermore, we identified that altered diffusion at the GWM boundary partially (i.e., ~50%) overlapped with atypical gray‐white matter tissue contrast in ASD. Our study thus replicates and extends previous work highlighting the GWM boundary as a crucial target of neuropathology in ASD, and guides future work elucidating etiological mechanisms.

Detection of gray matter microstructural changes in Alzheimer's disease continuum using fiber orientation

BACKGROUND

This study aimed to investigate feasible gray matter microstructural biomarkers with high sensitivity for early Alzheimer's disease (AD) detection. We propose a diffusion tensor imaging (DTI) measure, "radiality", as an early AD biomarker. It is the dot product of the normal vector of the cortical surface and primary diffusion direction, which reflects the fiber orientation within the cortical column.

METHODS

We analyzed neuroimages from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database, including images from 78 cognitively normal (CN), 50 early mild cognitive impairment (EMCI), 34 late mild cognitive impairment (LMCI), and 39 AD patients. We then evaluated the cortical thickness (CTh), mean diffusivity (MD), which are conventional AD magnetic resonance imaging (MRI) biomarkers, and the amount of accumulated amyloid and tau using positron emission tomography (PET). Radiality was projected on the gray matter surface to compare and validate the changes with different stages alongside other neuroimage biomarkers.

RESULTS

The results revealed decreased radiality primarily in the entorhinal, insula, frontal, and temporal cortex with further progression of disease. In particular, radiality could delineate the difference between the CN and EMCI groups, while the other biomarkers could not. We examined the relationship between radiality and other biomarkers to validate its pathological evidence in AD. Overall, radiality showed a high association with conventional biomarkers. Additional ROI analysis revealed the dynamics of AD-related changes as stages onward.

CONCLUSION

Radiality in cortical gray matter showed AD-specific changes and relevance with other conventional AD biomarkers with high sensitivity. Moreover, radiality could identify the group differences seen in EMCI, representative of changes in early AD, which supports its superiority in early diagnosis compared to that possible with conventional biomarkers. We provide evidence of structural changes with cognitive impairment and suggest radiality as a sensitive biomarker for identifying early AD.

DeepDTI: High-fidelity six-direction diffusion tensor imaging using deep learning

Diffusion tensor magnetic resonance imaging (DTI) is unsurpassed in its ability to map tissue microstructure and structural connectivity in the living human brain. Nonetheless, the angular sampling requirement for DTI leads to long scan times and poses a critical barrier to performing high-quality DTI in routine clinical practice and large-scale research studies. In this work we present a new processing framework for DTI entitled DeepDTI that minimizes the data requirement of DTI to six diffusion-weighted images (DWIs) required by conventional voxel-wise fitting methods for deriving the six unique unknowns in a diffusion tensor using data-driven supervised deep learning. DeepDTI maps the input non-diffusion-weighted (b ​= ​0) image and six DWI volumes sampled along optimized diffusion-encoding directions, along with T1-weighted and T2-weighted image volumes, to the residuals between the input and high-quality output b = 0 image and DWI volumes using a 10-layer three-dimensional convolutional neural network (CNN). The inputs and outputs of DeepDTI are uniquely formulated, which not only enables residual learning to boost CNN performance but also enables tensor fitting of resultant high-quality DWIs to generate orientational DTI metrics for tractography. The very deep CNN used by DeepDTI leverages the redundancy in local and non-local spatial information and across diffusion-encoding directions and image contrasts in the data. The performance of DeepDTI was systematically quantified in terms of the quality of the output images, DTI metrics, DTI-based tractography and tract-specific analysis results. We demonstrate rotationally-invariant and robust estimation of DTI metrics from DeepDTI that are comparable to those obtained with two b ​= ​0 images and 21 DWIs for the primary eigenvector derived from DTI and two b ​= ​0 images and 26-30 DWIs for various scalar metrics derived from DTI, achieving 3.3-4.6 × ​acceleration, and twice as good as those of a state-of-the-art denoising algorithm at the group level. The twenty major white-matter tracts can be accurately identified from the tractography of DeepDTI results. The mean distance between the core of the major white-matter tracts identified from DeepDTI results and those from the ground-truth results using 18 ​b ​= ​0 images and 90 DWIs measures around 1-1.5 ​mm. DeepDTI leverages domain knowledge of diffusion MRI physics and power of deep learning to render DTI, DTI-based tractography, major white-matter tracts identification and tract-specific analysis more feasible for a wider range of neuroscientific and clinical studies.

Improving a probabilistic cytoarchitectonic atlas of auditory cortex using a novel method for inter-individual alignment

The human superior temporal plane, the site of the auditory cortex, displays high inter-individual macro-anatomical variation. This questions the validity of curvature-based alignment (CBA) methods for in vivo imaging data. Here, we have addressed this issue by developing CBA+, which is a cortical surface registration method that uses prior macro-anatomical knowledge. We validate this method by using cytoarchitectonic areas on 10 individual brains (which we make publicly available). Compared to volumetric and standard surface registration, CBA+ results in a more accurate cytoarchitectonic auditory atlas. The improved correspondence of micro-anatomy following the improved alignment of macro-anatomy validates the superiority of CBA+ compared to CBA. In addition, we use CBA+ to align in vivo and postmortem data. This allows projection of functional and anatomical information collected in vivo onto the cytoarchitectonic areas, which has the potential to contribute to the ongoing debate on the parcellation of the human auditory cortex.

Probing in vivo cortical myeloarchitecture in humans via line-scan diffusion acquisitions at 7 T with 250-500 micron radial resolution

PURPOSE

The goal of this study was to measure diffusion signals within the cerebral cortex using the line-scan technique to achieve extremely high resolution in the radial direction (ie, perpendicular to the cortical surface) and to demonstrate the utility of these measurements for investigating laminar architecture in the living human brain.

METHODS

Line-scan diffusion data with 250-500 micron radial resolution were acquired at 7 T on 8 healthy volunteers, with each line prescribed perpendicularly to primary somatosensory cortex (S1) and primary motor cortex (M1). Apparent diffusion coefficients, fractional anisotropy values, and radiality indices were measured as a function of cortical depth.

RESULTS

In the deep layers of S1, we found evidence for high anisotropy and predominantly tangential diffusion, with low anisotropy observed in superficial S1. In M1, moderate anisotropy and predominantly radial diffusion was seen at almost all cortical depths. These patterns were consistent across subjects and were conspicuous without averaging data across different locations on the cortical sheet.

CONCLUSION

Our results are in accord with the myeloarchitecture of S1 and M1, known from prior histology studies: in S1, dense bands of tangential myelinated fibers run through the deep layers but not the superficial ones, and in M1, radial myelinated fibers are prominent at most cortical depths. This work therefore provides support for the idea that high-resolution diffusion signals, measured with the line-scan technique and receiving a boost in SNR at 7 T, may serve as a sensitive probe of in vivo laminar architecture.