Dataset Information

Voxel-wise and spatial modelling of binary lesion masks: Comparison of methods with a realistic simulation framework.

ABSTRACT:

Objectives

White matter lesions are a very common finding on MRI in older adults and their presence increases the risk of stroke and dementia. Accurate and computationally efficient modelling methods are necessary to map the association of lesion incidence with risk factors, such as hypertension. However, there is no consensus in the brain mapping literature whether a voxel-wise modelling approach is better for binary lesion data than a more computationally intensive spatial modelling approach that accounts for voxel dependence.

Methods

We review three regression approaches for modelling binary lesion masks including mass-univariate probit regression modelling with either maximum likelihood estimates, or mean bias-reduced estimates, and spatial Bayesian modelling, where the regression coefficients have a conditional autoregressive model prior to account for local spatial dependence. We design a novel simulation framework of artificial lesion maps to compare the three alternative lesion mapping methods. The age effect on lesion probability estimated from a reference data set (13,680 individuals from the UK Biobank) is used to simulate a realistic voxel-wise distribution of lesions across age. To mimic the real features of lesion masks, we propose matching brain lesion summaries (total lesion volume, average lesion size and lesion count) across the reference data set and the simulated data sets. Thus, we allow for a fair comparison between the modelling approaches, under a realistic simulation setting.

Results

Our findings suggest that bias-reduced estimates for voxel-wise binary-response generalized linear models (GLMs) overcome the drawbacks of infinite and biased maximum likelihood estimates and scale well for large data sets because voxel-wise estimation can be performed in parallel across voxels. Contrary to the assumption of spatial dependence being key in lesion mapping, our results show that voxel-wise bias-reduction and spatial modelling result in largely similar estimates.

Conclusions

Bias-reduced estimates for voxel-wise GLMs are not only accurate but also computationally efficient, which will become increasingly important as more biobank-scale neuroimaging data sets become available.

SUBMITTER: Kindalova P

PROVIDER: S-EPMC8752964 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundVoxel-based lesion symptom mapping (VLSM) assesses the relation of lesion location at a voxel level with a specific clinical or functional outcome measure at a population level. Spatial normalization, that is, mapping the patient images into an atlas coordinate system, is an essential pre-processing step of VLSM. However, no consensus exists on the optimal registration approach to compute the transformation nor are downstream effects on VLSM statistics explored. In this work, we evaluate four registration approaches commonly used in VLSM pipelines: affine (AR), nonlinear (NLR), nonlinear with cost function masking (CFM), and enantiomorphic registration (ENR). The evaluation is based on a standard VLSM scenario: the analysis of statistical relations of brain voxels and regions in imaging data acquired early after stroke onset with follow-up modified Rankin Scale (mRS) values.Materials and methodsFluid-attenuated inversion recovery (FLAIR) MRI data from 122 acute ischemic stroke patients acquired between 2 and 3 days after stroke onset and corresponding lesion segmentations, and 30 days mRS values from a European multicenter stroke imaging study (I-KNOW) were available and used in this study. The relation of the voxel location with follow-up mRS was assessed by uni- as well as multi-variate statistical testing based on the lesion segmentations registered using the four different methods (AR, NLR, CFM, ENR; implementation based on the ANTs toolkit).ResultsThe brain areas evaluated as important for follow-up mRS were largely consistent across the registration approaches. However, NLR, CFM, and ENR led to distortions in the patient images after the corresponding nonlinear transformations were applied. In addition, local structures (for instance the lateral ventricles) and adjacent brain areas remained insufficiently aligned with corresponding atlas structures even after nonlinear registration.ConclusionsFor VLSM study designs and imaging data similar to the present work, an additional benefit of nonlinear registration variants for spatial normalization seems questionable. Related distortions in the normalized images lead to uncertainties in the VLSM analyses and may offset the theoretical benefits of nonlinear registration.