Use of multi-perturbation Shapley analysis in lesion studies of functional networks: The case of upper limb paresis

Understanding the impact of variation in lesion topography on the expression of functional impairments following stroke is important, as it may pave the way to modeling structure-function relations in statistical terms while pointing to constraints for adaptive remapping and functional recovery. Multi-perturbation Shapley-value analysis (MSA) is a relatively novel game-theoretical approach for multivariate lesion-symptom mapping. In this methodological paper, we provide a comprehensive explanation of MSA. We use synthetic data to assess the method's accuracy and perform parameter optimization. We then demonstrate its application using a cohort of 107 first-event subacute stroke patients, assessed for upper limb (UL) motor impairment (Fugl-Meyer Assessment scale). Under the conditions tested, MSA could correctly detect simulated ground-truth lesion-symptom relationships with a sensitivity of 75% and specificity of similar to 90%. For real behavioral data, MSA disclosed a strong hemispheric effect in the relative contribution of specific regions-of-interest (ROIs): poststroke UL motor function was mostly contributed by damage to ROIs associated with movement planning (supplementary motor cortex and superior frontal gyrus) following left-hemispheric damage (LHD) and by ROIs associated with movement execution (primary motor and somatosensory cortices and the ventral brainstem) following right-hemispheric damage (RHD). Residual UL motor ability following LHD was found to depend on a wider array of brain structures compared to the residual motor ability of RHD patients. The results demonstrate that MSA can provide a unique insight into the relative importance of different hubs in neural networks, which is difficult to obtain using standard univariate methods.