Dataset Information

Modelling white matter in gyral blades as a continuous vector field.

ABSTRACT: 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 intrahemispheric connections when compared to tracers.

SUBMITTER: Cottaar M

PROVIDER: S-EPMC7610793 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Secondary white matter degeneration is a common occurrence after ischemic stroke, as identified by Diffusion Tensor Imaging (DTI). However, despite recent advances, the time course of the process is not completely understood. The primary aim of this study was to assess secondary degeneration using an approach whereby we create a patient-specific model of damaged fibers based on the volumetric characteristics of lesions. We also examined the effects of secondary degeneration along the modelled streamlines at different distances from the primary infarction using DTI. Eleven patients who presented with upper limb motor deficits at the time of a first-ever ischemic stroke were included. They underwent scanning at weeks 6 and 29 post-stroke. The fractional anisotropy (FA), mean diffusivity (MD), primary eigenvalue (λ1), and transverse eigenvalue (λ23) were measured. Using regions of interest based on the simulation output, the differences between the modelled fibers and matched contralateral areas were analyzed. The longitudinal change between the two time points and across five distances from the primary lesion was also assessed using the ratios of diffusion quantities (rFA, rMD, rλ1, and rλ23) between the ipsilesional and contralesional hemisphere. At week 6 post-stroke, significantly decreased λ1 was found along the ipsilesional corticospinal tract (CST) with a trend towards lower FA, reduced MD and λ23. At week 29 post-stroke, significantly decreased FA was shown relative to the non-lesioned side, with a trend towards lower λ1, unchanged MD, and higher λ23. Along the ipsilesional tract, the rFA diminished, whereas the rMD, rλ1, and rλ23 significantly increased over time. No significant variations in the time progressive effect with distance were demonstrated. The findings support previously described mechanisms of secondary degeneration and suggest that it spreads along the entire length of a damaged tract. Future investigations using higher-order tractography techniques can further explain the intravoxel alterations caused by ischemic injury.

Project description:IntroductionSubtle inhomogeneities in the scanner's magnetic fields (B0 and B1) alter the intensity levels of the structural magnetic resonance imaging (MRI) affecting the volumetric assessment of WMH changes. Here, we investigate the influence that (1) correcting the images for the B1 inhomogeneities (i.e. bias field correction (BFC)) and (2) selection of the WMH change assessment method can have on longitudinal analyses of WMH progression and discuss possible solutions.MethodsWe used brain structural MRI from 46 mild stroke patients scanned at stroke onset and 3 years later. We tested three BFC approaches: FSL-FAST, N4 and exponentially entropy-driven homomorphic unsharp masking (E(2)D-HUM) and analysed their effect on the measured WMH change. Separately, we tested two methods to assess WMH changes: measuring WMH volumes independently at both time points semi-automatically (MCMxxxVI) and subtracting intensity-normalised FLAIR images at both time points following image gamma correction. We then combined the BFC with the computational method that performed best across the whole sample to assess WMH changes.ResultsAnalysis of the difference in the variance-to-mean intensity ratio in normal tissue between BFC and uncorrected images and visual inspection showed that all BFC methods altered the WMH appearance and distribution, but FSL-FAST in general performed more consistently across the sample and MRI modalities. The WMH volume change over 3 years obtained with MCMxxxVI with vs. without FSL-FAST BFC did not significantly differ (medians(IQR)(with BFC) = 3.2(6.3) vs. 2.9(7.4)ml (without BFC), p = 0.5), but both differed significantly from the WMH volume change obtained from subtracting post-processed FLAIR images (without BFC)(7.6(8.2)ml, p < 0.001). This latter method considerably inflated the WMH volume change as subtle WMH at baseline that became more intense at follow-up were counted as increase in the volumetric change.ConclusionsMeasurement of WMH volume change remains challenging. Although the overall volumetric change was not significantly affected by the application of BFC, these methods distorted the image intensity distribution affecting subtle WMH. Subtracting the FLAIR images at both time points following gamma correction seems a promising technique but is adversely affected by subtle WMH. It is important to take into account not only the changes in volume but also in the signal intensity.

Project description:White matter hyperintensities (WMH), also known as white matter lesions, are localised white matter areas that appear hyperintense on MRI scans. WMH commonly occur in the ageing population, and are often associated with several factors such as cognitive disorders, cardiovascular risk factors, cerebrovascular and neurodegenerative diseases. Despite the fact that some links between lesion location and parametric factors such as age have already been established, the relationship between voxel-wise spatial distribution of lesions and these factors is not yet well understood. Hence, it would be of clinical importance to model the distribution of lesions at the population-level and quantitatively analyse the effect of various factors on the lesion distribution model. In this work we compare various methods, including our proposed method, to generate voxel-wise distributions of WMH within a population with respect to various factors. Our proposed Bayesian spline method models the spatio-temporal distribution of WMH with respect to a parametric factor of interest, in this case age, within a population. Our probabilistic model takes as input the lesion segmentation binary maps of subjects belonging to various age groups and provides a population-level parametric lesion probability map as output. We used a spline representation to ensure a degree of smoothness in space and the dimension associated with the parameter, and formulated our model using a Bayesian framework. We tested our algorithm output on simulated data and compared our results with those obtained using various existing methods with different levels of algorithmic and computational complexity. We then compared the better performing methods on a real dataset, consisting of 1000 subjects of the UK Biobank, divided in two groups based on hypertension diagnosis. Finally, we applied our method on a clinical dataset of patients with vascular disease. On simulated dataset, the results from our algorithm showed a mean square error (MSE) value of 7.27×10-5, which was lower than the MSE value reported in the literature, with the advantage of being robust and computationally efficient. In the UK Biobank data, we found that the lesion probabilities are higher for the hypertension group compared to the non-hypertension group and further verified this finding using a statistical t-test. Finally, when applying our method on patients with vascular disease, we observed that the overall probability of lesions is significantly higher in later age groups, which is in line with the current literature.

Dataset Information

Modelling white matter in gyral blades as a continuous vector field.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets