A gyral coordinate system predictive of fibre orientations

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.

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.

A model of tension-induced fiber growth predicts white matter organization during brain folding

The past decade has experienced renewed interest in the physical processes that fold the developing cerebral cortex. Biomechanical models and experiments suggest that growth of the cortex, outpacing growth of underlying subcortical tissue (prospective white matter), is sufficient to induce folding. However, current models do not explain the well-established links between white matter organization and fold morphology, nor do they consider subcortical remodeling that occurs during the period of folding. Here we propose a framework by which cortical folding may induce subcortical fiber growth and organization. Simulations incorporating stress-induced fiber elongation indicate that subcortical stresses resulting from folding are sufficient to induce stereotyped fiber organization beneath gyri and sulci. Model predictions are supported by high-resolution ex vivo diffusion tensor imaging of the developing rhesus macaque brain. Together, results provide support for the theory of cortical growth-induced folding and indicate that mechanical feedback plays a significant role in brain connectivity.

Loss and dispersion of superficial white matter in Alzheimer's disease: a diffusion MRI study

Pathological cerebral white matter changes in Alzheimer's disease have been shown using diffusion tensor imaging. Superficial white matter changes are relatively understudied despite their importance in cortico-cortical connections. Measuring superficial white matter degeneration using diffusion tensor imaging is challenging due to its complex organizational structure and proximity to the cortex. To overcome this, we investigated diffusion MRI changes in young-onset Alzheimer's disease using standard diffusion tensor imaging and Neurite Orientation Dispersion and Density Imaging to distinguish between disease-related changes that are degenerative (e.g. loss of myelinated fibres) and organizational (e.g. increased fibre dispersion). Twenty-nine young-onset Alzheimer's disease patients and 22 healthy controls had both single-shell and multi-shell diffusion MRI. We calculated fractional anisotropy, mean diffusivity, neurite density index, orientation dispersion index and tissue fraction (1-free water fraction). Diffusion metrics were sampled in 15 a priori regions of interest at four points along the cortical profile: cortical grey matter, grey/white boundary, superficial white matter (1 mm below grey/white boundary) and superficial/deeper white matter (2 mm below grey/white boundary). To estimate cross-sectional group differences, we used average marginal effects from linear mixed effect models of participants' diffusion metrics along the cortical profile. The superficial white matter of young-onset Alzheimer's disease individuals had lower neurite density index compared to controls in five regions (superior and inferior parietal, precuneus, entorhinal and parahippocampus) (all P < 0.05), and higher orientation dispersion index in three regions (fusiform, entorhinal and parahippocampus) (all P < 0.05). Young-onset Alzheimer's disease individuals had lower fractional anisotropy in the entorhinal and parahippocampus regions (both P < 0.05) and higher fractional anisotropy within the postcentral region (P < 0.05). Mean diffusivity was higher in the young-onset Alzheimer's disease group in the parahippocampal region (P < 0.05) and lower in the postcentral, precentral and superior temporal regions (all P < 0.05). In the overlying grey matter, disease-related changes were largely consistent with superficial white matter findings when using neurite density index and fractional anisotropy, but appeared at odds with orientation dispersion and mean diffusivity. Tissue fraction was significantly lower across all grey matter regions in young-onset Alzheimer's disease individuals (all P < 0.001) but group differences reduced in magnitude and coverage when moving towards the superficial white matter. These results show that microstructural changes occur within superficial white matter and along the cortical profile in individuals with young-onset Alzheimer's disease. Lower neurite density and higher orientation dispersion suggests underlying fibres undergo neurodegeneration and organizational changes, two effects previously indiscernible using standard diffusion tensor metrics in superficial white matter.

Active Cortex Tractography

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

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.

Modelling white matter in gyral blades as a continuous vector field

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

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.