A four-dimensional computational model of dynamic contrast-enhanced magnetic resonance imaging measurement of subtle blood-brain barrier leakage.

Dynamic contrast-enhanced MRI (DCE-MRI) is increasingly used to quantify and map the spatial distribution of blood-brain barrier (BBB) leakage in neurodegenerative disease, including cerebral small vessel disease and dementia. However, the subtle nature of leakage and resulting small signal changes make quantification challenging. While simplified one-dimensional simulations have probed the impact of noise, scanner drift, and model assumptions, the impact of spatio-temporal effects such as gross motion, k-space sampling and motion artefacts on parametric leakage maps has been overlooked. Moreover, evidence on which to base the design of imaging protocols is lacking due to practical difficulties and the lack of a reference method. To address these problems, we present an open-source computational model of the DCE-MRI acquisition process for generating four dimensional Digital Reference Objects (DROs), using a high-resolution brain atlas and incorporating realistic patient motion, extra-cerebral signals, noise and k-space sampling. Simulations using the DROs demonstrated a dominant influence of spatio-temporal effects on both the visual appearance of parameter maps and on measured tissue leakage rates. The computational model permits greater understanding of the sensitivity and limitations of subtle BBB leakage measurement and provides a non-invasive means of testing and optimising imaging protocols for future studies.

Assessment of perivascular space filtering methods using a three-dimensional computational model.

Growing interest surrounds the assessment of perivascular spaces (PVS) on magnetic resonance imaging (MRI) and their validation as a clinical biomarker of adverse brain health. Nonetheless, the limits of validity of current state-of-the-art segmentation methods are still unclear. Here, we propose an open-source three-dimensional computational framework comprising 3D digital reference objects and evaluate the performance of three PVS filtering methods under various spatiotemporal imaging considerations (including sampling, motion artefacts, and Rician noise). Specifically, we study the performance of the Frangi, Jerman and RORPO filters in enhancing PVS-like structures to facilitate segmentation. Our findings were three-fold. First, as long as voxels are isotropic, RORPO outperforms the other two filters, regardless of imaging quality. Unlike the Frangi and Jerman filters, RORPO's performance does not deteriorate as PVS volume increases. Second, the performance of all "vesselness" filters is heavily influenced by imaging quality, with sampling and motion artefacts being the most damaging for these types of analyses. Third, none of the filters can distinguish PVS from other hyperintense structures (e.g. white matter hyperintensities, stroke lesions, or lacunes) effectively, the area under precision-recall curve dropped substantially (Frangi: from 94.21 [IQR 91.60, 96.16] to 43.76 [IQR 25.19, 63.38]; Jerman: from 94.51 [IQR 91.90, 95.37] to 58.00 [IQR 35.68, 64.87]; RORPO: from 98.72 [IQR 95.37, 98.96] to 71.87 [IQR 57.21, 76.63] without and with other hyperintense structures, respectively). The use of our computational model enables comparing segmentation methods and identifying their advantages and disadvantages, thereby providing means for testing and optimising pipelines for ongoing and future studies.