Lesion Symptom Mapping Guide

Lesion-Symptom Mapping Guide

Purpose

This skill encodes expert methodological knowledge for lesion-symptom mapping in clinical and cognitive neuroscience. A competent programmer without neuropsychology and neuroimaging training will get this wrong because:

• Lesion distributions are not random. Stroke lesions cluster in the middle cerebral artery (MCA) territory, creating systematic collinearity between brain regions. Standard voxelwise methods confound truly causal regions with regions that are simply co-damaged (Sperber, 2020).
• Lesion volume is a massive confound. Larger lesions produce worse behavioral deficits simply because more tissue is damaged. Any analysis that does not control for total lesion volume will attribute behavioral deficits to whichever regions are most often part of large lesions (DeMarco & Turkeltaub, 2018).
• White matter disconnection matters as much as grey matter damage. A focal lesion can disrupt distant regions via white matter pathway disruption (diaschisis). VLSM misses this entirely; disconnection analysis is needed (Foulon et al., 2018).
• Standard multiple comparison correction is insufficient. Permutation-based FWE correction is required because voxelwise tests are massively non-independent (lesion voxels are spatially correlated), violating assumptions of FDR and parametric corrections (Kimberg et al., 2007).
• Small samples produce false localizations. With fewer than 50 patients, VLSM lacks power and produces unreliable maps that do not replicate (Sperber, 2020).

When to Use This Skill

• Planning a voxel-based lesion-symptom mapping (VLSM) study
• Choosing between VLSM, multivariate, and disconnection-based approaches
• Setting statistical thresholds and correction methods for lesion analyses
• Evaluating whether a published lesion study adequately controlled for confounds
• Implementing lesion segmentation and registration pipelines
• Conducting network-based lesion mapping with normative connectome data

Do NOT use this skill for:

• Radiological lesion diagnosis (requires clinical neuroradiology)
• General fMRI analysis (see fmri-glm-analysis-guide)
• White matter tractography methodology (upstream of this skill)

Research Planning Protocol

Before executing the domain-specific steps below, you MUST:

1. State the research question -- What specific question is this analysis/paradigm addressing?
2. Justify the method choice -- Why is this approach appropriate? What alternatives were considered?
3. Declare expected outcomes -- What results would support vs. refute the hypothesis?
4. Note assumptions and limitations -- What does this method assume? Where could it mislead?
5. Present the plan to the user and WAIT for confirmation before proceeding.

For detailed methodology guidance, see the research-literacy skill.

⚠️ Verification Notice

This skill was generated by AI from academic literature. All parameters, thresholds, and citations require independent verification before use in research. If you find errors, please open an issue.

Method Selection Decision Tree

What is your research question?
|
+-- "Which brain voxels are associated with a behavioral deficit?"
| |
| +-- N >= 50 patients, continuous outcome
| | --> VLSM (mass-univariate)
| |
| +-- N >= 50 patients, binary outcome
| | --> VLSM with Brunner-Munzel or chi-square
| |
| +-- N >= 100 patients, distributed representations expected
| | --> SVR-LSM (multivariate)
| |
| +-- N < 50 patients
| --> Underpowered for VLSM; consider ROI-based approach
| or case-series descriptive analysis
|
+-- "Which white matter pathways mediate the deficit?"
| --> Disconnection analysis (BCBToolkit, Disconnectome)
| Can supplement VLSM or replace it when tracts are the question
|
+-- "Which brain networks, when lesioned, produce this symptom?"
 --> Lesion network mapping (normative connectome)
 Maps lesion location to network disruption

Lesion Segmentation

Methods

Domain judgment: Manual tracing remains the gold standard. Semi-automated methods should always be visually inspected and manually corrected. Never trust a fully automated segmentation without visual QC of every patient (Brett et al., 2001).

Registration to Standard Space

1. Cost function masking (CRITICAL): When registering a lesioned brain to MNI template, the lesion MUST be masked out of the cost function. Without masking, the registration algorithm warps healthy tissue to fill the lesion, distorting the spatial normalization (Brett et al., 2001).

2. Recommended pipeline:

• Segment lesion on native T1 (or FLAIR)
• Create binary lesion mask
• Register T1 to MNI using nonlinear registration (e.g., ANTs SyN) with lesion mask as cost function mask
• Apply the same warp to the lesion mask
• Verify registration quality visually

3. Enantiomorphic normalization: For large lesions, replace the lesioned hemisphere with a flipped version of the intact hemisphere before registration. This improves normalization quality for large lesions (Nachev et al., 2008).

Voxel-Based Lesion-Symptom Mapping (VLSM)

Prerequisites and Sample Size

Domain judgment: VLSM with N < 50 is severely underpowered. Sperber (2020) showed that with N = 30, VLSM produces maps that fail to replicate and have unacceptably high false discovery rates. If N < 50, consider ROI-based analyses with a priori regions or descriptive lesion overlap approaches.

Statistical Tests

Multiple Comparison Correction

Domain judgment: Permutation testing is strongly preferred because lesion maps violate the independence assumptions of FDR and parametric corrections. The spatial correlation structure of lesions means neighboring voxels are highly non-independent. Permutation testing implicitly accounts for this correlation structure (Kimberg et al., 2007).

Controlling for Lesion Volume

Lesion volume MUST be controlled. Methods (DeMarco & Turkeltaub, 2018):

1. Direct regression: Include total lesion volume as a covariate in the voxelwise regression model
2. Behavioral residualization: Regress behavior on lesion volume first; use residuals as the dependent variable in VLSM
3. Both yield similar results, but direct regression is preferred for interpretability (DeMarco & Turkeltaub, 2018)

CRITICAL: Failing to control for lesion volume is the single most common error in VLSM studies. Large lesions damage more regions and produce worse deficits, creating a spurious correlation between any frequently-damaged voxel and behavioral impairment (DeMarco & Turkeltaub, 2018).

Multivariate Lesion-Symptom Mapping

SVR-LSM (Zhang et al., 2014)

Support vector regression-based lesion-symptom mapping considers the full lesion pattern simultaneously, addressing the collinearity problem of mass-univariate VLSM.

Advantages over VLSM (Zhang et al., 2014):

• Considers spatial covariance of lesion damage (voxels are analyzed jointly, not independently)
• Better handles the collinearity created by vascular territories
• Can detect distributed patterns

Limitations:

• Computationally expensive (hours per analysis with permutation testing)
• Requires larger samples (N >= 80--100 recommended) for stable SVR weights
• Interpretation of SVR weight maps is less straightforward than VLSM t-maps

Machine Learning-Based Mapping (MLBM)

Other multivariate approaches (random forests, LASSO regression) can be applied:

• LASSO is useful for variable selection among brain regions
• Random forests handle nonlinear relationships
• All require permutation-based significance testing
• Cross-validation (leave-one-out or k-fold) is mandatory

Disconnection Analysis

Rationale

A focal lesion disrupts not only the damaged tissue but also white matter pathways passing through the lesion, disconnecting distant brain regions (Foulon et al., 2018). VLSM maps only the lesion site; disconnection analysis maps the affected structural connections.

Methods

BCBToolkit Pipeline (Foulon et al., 2018)

1. Register lesion mask to MNI space
2. Use the normative tractography atlas (from 170 healthy controls) to identify tracts passing through each lesion
3. Generate a disconnection probability map: for each voxel outside the lesion, the probability that it is disconnected
4. Use disconnection maps (instead of or in addition to lesion masks) as predictors in VLSM-like analyses

Advantages Over VLSM

• Captures remote effects (diaschisis) that VLSM misses entirely
• Accounts for the fact that two lesions in the same voxel can disconnect different pathways depending on their extent
• Particularly important for white matter lesions, where the damaged tissue itself has no gray matter function

Network Lesion Mapping

Lesion Network Mapping (Boes et al., 2015; Fox, 2018)

Maps a lesion to the brain-wide functional network it disrupts, using normative resting-state fMRI data:

1. Use the lesion as a seed region in a normative resting-state fMRI connectome (e.g., 1000-subject dataset)
2. Compute the functional connectivity map of the lesion location in healthy brains
3. The resulting map shows which brain regions are functionally connected to the lesion site
4. Across patients, identify the network that is commonly disrupted

Key Considerations

• Normative connectome: Results depend on the quality and size of the normative dataset. Larger datasets (N >= 500 healthy controls) provide more reliable connectivity estimates (Fox, 2018).
• Seed definition: Use the entire lesion as a seed, not just the peak voxel. Consider weighting by lesion probability if lesion masks are probabilistic.
• Validation: Compare network maps to known functional anatomy. The disrupted network should be biologically plausible for the observed deficit.
• Limitations: Assumes that healthy-brain connectivity predicts the effect of a lesion. Does not account for reorganization or compensation (Fox, 2018).