@brainpilot/skills 0.0.6

package/skills/10_Clinical_Neuropsychology/lesion-symptom-mapping-guide/SKILL.md

@@ -0,0 +1,367 @@
+---
+name: "lesion-symptom-mapping-guide"
+description: "Advises on lesion-symptom mapping methods: VLSM, disconnection analysis, network lesion mapping"
+domain: "clinical-neuroscience"
+version: "1.0.0"
+authors:
+ - "Claude Code Agent"
+papers:
+ - "Bates et al., 2003"
+ - "Kimberg et al., 2007"
+ - "Foulon et al., 2018"
+ - "Zhang et al., 2014"
+ - "Boes et al., 2015"
+ - "Sperber, 2020"
+dependencies:
+ required:
+ - research-literacy
+ recommended:
+ - fmri-preprocessing-pipeline-guide
+review_status: "ai-generated"
+---
+
+# 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](https://github.com/HaoxuanLiTHUAI/awesome_cognitive_and_neuroscience_skills/issues).
+
+## 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
+
+| Method | Description | Accuracy | Time per patient | Source |
+|---|---|---|---|---|
+| **Manual tracing** | Expert traces lesion on each slice | Gold standard | 30--60 min | Brett et al., 2001 |
+| **Semi-automated (lesion_gnb)** | Gaussian naive Bayes on FLAIR | Good for chronic WM lesions | 5--15 min | Pustina et al., 2016 |
+| **LINDA** | Random forest on T1 | Good for chronic stroke | 5--10 min | Pustina et al., 2016 |
+| **U-net / deep learning** | Trained CNNs | Approaching manual accuracy | < 1 min | Kamnitsas et al., 2017 |
+
+> **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
+
+| Requirement | Minimum | Recommended | Source |
+|---|---|---|---|
+| Sample size | N >= **50** | N >= **100** | Sperber, 2020; Kimberg et al., 2007 |
+| Lesion overlap per voxel | >= **10% of sample** (or N >= 10) | >= 15% | Kimberg et al., 2007 |
+| Behavioral measure | Continuous preferred | -- | Bates et al., 2003 |
+
+> **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
+
+| Test | When to Use | Advantages | Source |
+|---|---|---|---|
+| **t-test** | Continuous behavior, binary lesion status per voxel | Simple, widely used | Bates et al., 2003 |
+| **Brunner-Munzel** | Non-normal behavioral data, unequal variances | Robust to non-normality and unequal group sizes | Rorden et al., 2007 |
+| **Regression** | Controlling for covariates (age, lesion volume) | Flexible; includes confound control | Sperber, 2020 |
+| **Liebermeister** | Binary behavioral outcome (impaired/spared) | For binary classification | Rorden et al., 2007 |
+
+### Multiple Comparison Correction
+
+| Method | Description | Recommended? | Source |
+|---|---|---|---|
+| **Permutation-based FWE** | Permute behavioral scores 5000+ times; threshold at 5th percentile of max statistic | **YES -- gold standard** | Kimberg et al., 2007 |
+| **FDR** (Benjamini-Hochberg) | Controls false discovery rate | Acceptable alternative but assumes independence | Kimberg et al., 2007 |
+| **Bonferroni** | Divide alpha by number of voxels | **Too conservative**; almost never detects effects | Expert consensus |
+| **Uncorrected** | No correction | **NEVER for publication** | Expert consensus |
+
+> **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.
+
+| Parameter | Recommended Value | Source |
+|---|---|---|
+| Kernel | **Linear** | Zhang et al., 2014 |
+| C parameter | **30** (default for LSM) | Zhang et al., 2014 |
+| Feature reduction | Remove voxels with < 10% lesion overlap | Zhang et al., 2014 |
+| Statistical inference | **Permutation testing** (5000+ permutations) | Zhang et al., 2014 |
+
+**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
+
+| Tool | Approach | Data Required | Source |
+|---|---|---|---|
+| **BCBToolkit** | Maps lesion to disconnected tracts using normative tractography atlas | Lesion mask in MNI space | Foulon et al., 2018 |
+| **Disconnectome maps** | Pre-computed: for each brain voxel, which tracts are disconnected | Lesion mask in MNI space | Thiebaut de Schotten et al., 2015 |
+| **Individual tractography** | DTI/DWI tractography in each patient | Patient DWI data | Expert consensus |
+
+### 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).
+
+---
+
+## Common Confounds
+
+### 1. Lesion Volume Correlation
+
+Larger lesions produce worse behavioral deficits and damage more voxels. Without controlling for lesion volume, VLSM maps reflect voxels that are part of large lesions rather than voxels critical for the behavior (DeMarco & Turkeltaub, 2018).
+
+### 2. Non-Random Lesion Distribution (MCA Territory Bias)
+
+Stroke lesions are not uniformly distributed across the brain. The MCA territory (lateral frontal, temporal, parietal, insular cortex) is disproportionately affected. This means:
+- High statistical power in MCA territory, low power elsewhere
+- Voxels outside MCA territory may be critical but undetectable
+- Collinearity between MCA-territory voxels inflates false positives for non-critical regions (Sperber, 2020)
+
+### 3. Time Post-Onset
+
+| Phase | Time | Concern | Source |
+|---|---|---|---|
+| Acute (< 2 weeks) | Edema, diaschisis, penumbra | Lesion extent overestimated; behavior worst | Karnath et al., 2004 |
+| Subacute (2 weeks -- 3 months) | Recovery, reorganization | Lesion stabilizing; behavior improving | Expert consensus |
+| Chronic (> 3 months) | Stable lesion | **Preferred for VLSM**; most stable brain-behavior relationship | Sperber, 2020 |
+
+> **Domain judgment**: Chronic-phase data (> 3 months post-onset) is strongly preferred for VLSM because both lesion extent and behavioral deficits have stabilized. Acute-phase data confounds true lesion effects with transient diaschisis and edema (Karnath et al., 2004; Sperber, 2020).
+
+### 4. Covariates
+
+Always consider controlling for:
+- **Age**: Older patients have worse outcomes independent of lesion
+- **Education**: Affects cognitive test performance
+- **Time post-onset**: If sample includes mixed acute/chronic patients
+- **Handedness**: Affects lateralization of language
+- **Lesion hemisphere**: If analyzing bilateral samples, hemisphere effects must be modeled
+
+---
+
+## Software
+
+| Software | Methods | Language | Source |
+|---|---|---|---|
+| **NiiStat** | VLSM, ROI analysis | MATLAB | Rorden et al., 2007 |
+| **VLSM2** | VLSM with permutation testing | MATLAB | Bates et al., 2003 |
+| **SVR-LSM toolbox** | Multivariate SVR-LSM | MATLAB | Zhang et al., 2014 |
+| **BCBToolkit** | Disconnection analysis, disconnectome | GUI/Python | Foulon et al., 2018 |
+| **LESYMAP** | VLSM, SVR-LSM, SCCAN | R | Pustina et al., 2018 |
+| **ANTs** | Registration with cost function masking | C++/Python | Avants et al., 2011 |
+| **FSL** | Registration, lesion masking | Python/C++ | Jenkinson et al., 2012 |
+
+---
+
+## Common Pitfalls
+
+### 1. Not Controlling for Lesion Volume
+
+The most frequent and most damaging error. Always include lesion volume as a covariate or use residualized behavioral scores (DeMarco & Turkeltaub, 2018).
+
+### 2. Too Small a Sample
+
+VLSM with N < 50 produces unreliable maps. With N = 30, false positive rates can exceed 50% in some simulations (Sperber, 2020). If your sample is small, use ROI-based approaches or descriptive methods.
+
+### 3. Not Using Cost Function Masking During Registration
+
+Registering lesioned brains to template without masking the lesion distorts the normalization, warping healthy tissue into the lesion cavity and misaligning the rest of the brain (Brett et al., 2001).
+
+### 4. Using Uncorrected or Bonferroni Correction
+
+Uncorrected thresholds produce massive false positives. Bonferroni is too conservative due to spatial correlation. Use permutation-based FWE (Kimberg et al., 2007).
+
+### 5. Ignoring the Vascular Architecture
+
+Interpreting a VLSM map as showing "the region responsible for function X" ignores that vascular territory collinearity may have driven the result. Consider supplementing with disconnection analysis or using multivariate methods (Sperber, 2020).
+
+### 6. Mixing Acute and Chronic Patients Without Controlling for Time
+
+Acute patients have larger effective lesions (edema) and worse behavior, confounding time-post-onset with lesion severity. Analyze chronic patients separately or include time as a covariate (Karnath et al., 2004).
+
+---
+
+## Minimum Reporting Checklist
+
+Based on Sperber (2020), DeMarco & Turkeltaub (2018), and Kimberg et al. (2007):
+
+- [ ] Sample size and patient demographics (age, sex, education, handedness)
+- [ ] Etiology (ischemic stroke, hemorrhagic, tumor, etc.) and affected hemisphere(s)
+- [ ] Time post-onset for each patient (or range)
+- [ ] Lesion segmentation method (manual, semi-automated, automated) and software
+- [ ] Registration method and whether cost function masking was used
+- [ ] Minimum lesion overlap threshold for voxel inclusion
+- [ ] Statistical test used (t-test, Brunner-Munzel, regression, SVR)
+- [ ] Multiple comparison correction method and parameters (number of permutations)
+- [ ] Whether and how lesion volume was controlled
+- [ ] Behavioral measure(s) and their psychometric properties
+- [ ] Software and version for all analysis steps
+- [ ] For disconnection analysis: normative tractography dataset and pipeline details
+
+---
+
+## Key References
+
+- Bates, E., Wilson, S. M., Saygin, A. P., Dick, F., Sereno, M. I., Knight, R. T., & Dronkers, N. F. (2003). Voxel-based lesion-symptom mapping. *Nature Neuroscience*, 6(5), 448--450.
+- Boes, A. D., Prasad, S., Liu, H., Liu, Q., Pascual-Leone, A., Caviness, V. S., & Fox, M. D. (2015). Network localization of neurological symptoms from focal brain lesions. *Brain*, 138(10), 3061--3075.
+- Brett, M., Leff, A. P., Rorden, C., & Ashburner, J. (2001). Spatial normalization of brain images with focal lesions using cost function masking. *NeuroImage*, 14(2), 486--500.
+- DeMarco, A. T., & Turkeltaub, P. E. (2018). A multivariate lesion symptom mapping toolbox and examination of lesion-volume biases and correction methods in lesion-symptom mapping. *Human Brain Mapping*, 39(11), 4169--4182.
+- Foulon, C., Cerliani, L., Kinkingnehun, S., Levy, R., Rosso, C., Urbanski, M., Volle, E., & Thiebaut de Schotten, M. (2018). Advanced lesion symptom mapping analyses and implementation as BCBtoolkit. *GigaScience*, 7(3), giy004.
+- Fox, M. D. (2018). Mapping symptoms to brain networks with the human connectome. *New England Journal of Medicine*, 379(23), 2237--2245.
+- Kamnitsas, K., Ledig, C., Newcombe, V. F. J., Simpson, J. P., Kane, A. D., Menon, D. K., Rueckert, D., & Glocker, B. (2017). Efficient multi-scale 3D CNN with fully connected CRF for accurate brain lesion segmentation. *Medical Image Analysis*, 36, 61--78.
+- Karnath, H.-O., Fruhmann Berger, M., Kuker, W., & Rorden, C. (2004). The anatomy of spatial neglect based on voxelwise statistical analysis: A study of 140 patients. *Cerebral Cortex*, 14(10), 1164--1172.
+- Kimberg, D. Y., Coslett, H. B., & Schwartz, M. F. (2007). Power in voxel-based lesion-symptom mapping. *Journal of Cognitive Neuroscience*, 19(7), 1067--1080.
+- Nachev, P., Coulthard, E., Jager, H. R., Kennard, C., & Husain, M. (2008). Enantiomorphic normalization of focally lesioned brains. *NeuroImage*, 39(3), 1215--1226.
+- Pustina, D., Coslett, H. B., Turkeltaub, P. E., Tustison, N., Schwartz, M. F., & Avants, B. (2016). Automated segmentation of chronic stroke lesions using LINDA: Lesion identification with neighborhood data analysis. *Human Brain Mapping*, 37(4), 1405--1421.
+- Rorden, C., Karnath, H.-O., & Bonilha, L. (2007). Improving lesion-symptom mapping. *Journal of Cognitive Neuroscience*, 19(7), 1081--1088.
+- Sperber, C. (2020). Rethinking causality and data complexity in brain lesion-behaviour inference and its implications for lesion-behaviour modelling. *Cortex*, 126, 49--62.
+- Zhang, Y., Kimberg, D. Y., Coslett, H. B., Schwartz, M. F., & Wang, Z. (2014). Multivariate lesion-symptom mapping using support vector regression. *Human Brain Mapping*, 35(12), 5861--5876.
+
+See `references/vlsm-pipeline.md` for step-by-step VLSM analysis workflow.
+See `references/disconnection-guide.md` for detailed BCBToolkit and disconnection analysis procedures.

package/skills/10_Clinical_Neuropsychology/lesion-symptom-mapping-guide/references/disconnection-guide.md

@@ -0,0 +1,152 @@
+# Disconnection Analysis and Network Lesion Mapping Guide
+
+This reference document supplements `SKILL.md` with detailed procedures for disconnection analysis using BCBToolkit and network lesion mapping approaches.
+
+## Disconnection Analysis with BCBToolkit
+
+### Overview
+
+BCBToolkit (Foulon et al., 2018) maps focal brain lesions to the white matter pathways they disrupt, using a normative tractography atlas derived from 170 healthy control diffusion MRI datasets.
+
+### Step-by-Step Procedure
+
+#### Step 1: Prepare Lesion Masks
+
+- Lesion masks must be in MNI152 space (1mm or 2mm resolution)
+- Binary format (0 = healthy, 1 = lesioned)
+- One NIfTI file per patient
+
+#### Step 2: Run Disconnection Mapping
+
+For each patient lesion mask:
+1. The toolbox identifies which tractography streamlines pass through the lesion
+2. For each voxel outside the lesion, compute the probability of disconnection:
+ ```
+ P(disconnection at voxel v) = N(streamlines through lesion AND through v) / N(total streamlines through v)
+ ```
+3. Output: a whole-brain disconnection probability map per patient
+
+#### Step 3: Threshold Disconnection Maps
+
+- Apply a disconnection probability threshold (typically **> 0.5** or **> 0.7**)
+- Lower thresholds are more sensitive but include more false positives
+- Higher thresholds are more specific but may miss partially disconnected regions
+
+#### Step 4: Statistical Analysis
+
+Use disconnection probability maps in the same statistical framework as VLSM:
+- **AnaCOM**: Anatomical connectivity-based mapping (included in BCBToolkit)
+- **Voxelwise regression**: Disconnection probability as predictor of behavior
+- **Tract-specific analysis**: Sum disconnection probability within specific tracts (e.g., arcuate fasciculus, corticospinal tract)
+
+### Tract-Specific Analysis
+
+Major tracts to consider in lesion studies:
+
+| Tract | Function | Common Lesion Association | Source |
+|---|---|---|---|
+| Arcuate fasciculus | Language (phonological) | Conduction aphasia, repetition deficits | Catani & Mesulam, 2008 |
+| Inferior fronto-occipital fasciculus | Semantic processing | Semantic deficits | Duffau et al., 2005 |
+| Superior longitudinal fasciculus (II, III) | Spatial attention, praxis | Neglect, apraxia | Thiebaut de Schotten et al., 2011 |
+| Uncinate fasciculus | Ventral semantic pathway | Naming deficits (controversial) | Papagno et al., 2011 |
+| Corticospinal tract | Motor function | Hemiparesis | Zhu et al., 2010 |
+| Inferior longitudinal fasciculus | Visual processing | Visual recognition deficits | Catani et al., 2003 |
+| Cingulum | Memory, emotional processing | Memory deficits | Bubb et al., 2018 |
+| Corpus callosum | Interhemispheric transfer | Disconnection syndromes | Gazzaniga, 2005 |
+
+### Comparing VLSM and Disconnection Results
+
+Run both VLSM and disconnection analysis on the same dataset to compare:
+
+1. **Overlap**: Regions identified by both methods are strong candidates for critical involvement
+2. **VLSM-only regions**: May reflect direct grey matter damage effects
+3. **Disconnection-only regions**: May reflect remote effects via white matter disruption
+4. **Statistical comparison**: Test whether disconnection probability explains additional behavioral variance beyond lesion location (using hierarchical regression)
+
+## Network Lesion Mapping
+
+### Lesion Network Mapping Procedure (Boes et al., 2015; Fox, 2018)
+
+#### Step 1: Obtain Normative Connectome Data
+
+- Use a large resting-state fMRI dataset (N >= 500 healthy controls)
+- Common choices:
+ - GSP 1000 connectome (Harvard/MGH)
+ - Human Connectome Project (HCP)
+ - UK Biobank
+- Preprocess normative data with standard resting-state pipeline (motion correction, bandpass filtering 0.01--0.1 Hz, nuisance regression)
+
+#### Step 2: Compute Lesion-Seeded Connectivity
+
+For each patient:
+1. Use the lesion mask as a seed region
+2. Extract the mean BOLD time series from the seed location using the normative data
+3. Compute whole-brain connectivity (Pearson correlation) from the seed
+4. Fisher z-transform the correlation map
+
+#### Step 3: Identify Common Network
+
+Across all patients:
+1. Threshold individual connectivity maps (e.g., at z > 0.2)
+2. Compute the overlap: which regions are functionally connected to ALL (or most) lesion sites
+3. Alternatively, correlate patient-specific connectivity patterns with behavioral scores
+
+#### Step 4: Validate
+
+- Compare the identified network to known functional anatomy
+- Test whether network disruption score predicts behavioral deficit severity
+- Cross-validate using leave-one-out or split-half approaches
+
+### Network Lesion Mapping Considerations
+
+| Consideration | Details | Source |
+|---|---|---|
+| Normative data quality | Larger, better-preprocessed datasets yield more reliable maps | Fox, 2018 |
+| Lesion size effects | Large lesions produce more diffuse connectivity maps; consider size-matched controls | Fox, 2018 |
+| Compensation/reorganization | Normative connectivity does not account for post-lesion plasticity | Fox, 2018 |
+| Symptom specificity | Different symptoms from lesions in the same location may map to different networks | Boes et al., 2015 |
+| Laterality | Ensure normative data includes adequate coverage of both hemispheres | Expert consensus |
+
+### When to Use Each Approach
+
+| Research Question | Recommended Method | Rationale |
+|---|---|---|
+| Which voxels are associated with deficit? | VLSM | Direct, established, widely accepted |
+| Which tracts mediate the deficit? | Disconnection analysis | Captures white matter effects |
+| Which brain network is disrupted? | Lesion network mapping | Captures distributed network effects |
+| What explains the most variance? | Combine all three | Compare models; use hierarchical regression |
+
+## Integration of Multiple Approaches
+
+### Recommended Multi-Method Pipeline
+
+1. **Start with VLSM**: Establish voxel-level associations (requires N >= 50)
+2. **Add disconnection analysis**: Identify tract-level contributions (BCBToolkit)
+3. **Add network mapping**: Identify network-level disruption (if normative connectome available)
+4. **Compare**: Use hierarchical regression to test which level (voxel, tract, network) best explains behavioral variance
+5. **Report all**: Provide maps at all three levels for comprehensive brain-behavior inference
+
+### Statistical Framework for Model Comparison
+
+```
+Model 1: Behavior ~ Lesion volume (baseline)
+Model 2: Behavior ~ Lesion volume + VLSM ROI damage
+Model 3: Behavior ~ Lesion volume + Tract disconnection
+Model 4: Behavior ~ Lesion volume + Network disruption
+Model 5: Behavior ~ Lesion volume + All levels
+```
+
+Compare using AIC/BIC or likelihood ratio tests to determine which level provides the best explanation.
+
+## References
+
+- Boes, A. D., et al. (2015). Network localization of neurological symptoms from focal brain lesions. *Brain*, 138(10), 3061--3075.
+- Bubb, E. J., Metzler-Baddeley, C., & Aggleton, J. P. (2018). The cingulum bundle: Anatomy, function, and dysfunction. *Neuroscience & Biobehavioral Reviews*, 92, 104--127.
+- Catani, M., & Mesulam, M. (2008). The arcuate fasciculus and the disconnection theme in language and aphasia. *Cortex*, 44(8), 953--961.
+- Catani, M., et al. (2003). Occipito-temporal connections in the human brain. *Brain*, 126(9), 2093--2107.
+- Duffau, H., et al. (2005). New insights into the anatomo-functional connectivity of the semantic system. *NeuroReport*, 16(13), 1461--1465.
+- Foulon, C., et al. (2018). Advanced lesion symptom mapping analyses and implementation as BCBtoolkit. *GigaScience*, 7(3), giy004.
+- Fox, M. D. (2018). Mapping symptoms to brain networks with the human connectome. *New England Journal of Medicine*, 379(23), 2237--2245.
+- Gazzaniga, M. S. (2005). Forty-five years of split-brain research and still going strong. *Nature Reviews Neuroscience*, 6(8), 653--659.
+- Thiebaut de Schotten, M., et al. (2011). A lateralized brain network for visuospatial attention. *Nature Neuroscience*, 14(10), 1245--1246.
+- Thiebaut de Schotten, M., et al. (2015). Brain disconnections link structural connectivity with function and behaviour. *Nature Communications*, 6, 7462.