Neuroplasticity and immune system are related to altered grey matter networks: a cohort study in sporadic Alzheimer's disease.
The 2 matches
- [1] § Methods › Single-subject GM networks ↔ Extract_individual_GM_networks_v20131106/batch_extract_networks20131106.m, lines 1–66 · score 0.83 · GM network, standard deviation, clustering, MATLAB, angles, segmentation
- [2] § Methods › MRI acquisition and preprocessing ↔ Extract_individual_GM_networks_v20131106/batch_extract_networks20131106.m, lines 69–120 · score 0.60 · isotropic voxels, SPM, resliced, volumes, segment, scans
Paper
Loaded from Europe PMC by your browser, not stored by OSCR: doi.org · Europe PMC
The paper is loaded when this pane is shown.
The authors' code
MATLAB · 208 lines · 8.4 KB · MIT · 2 matches
- % Batch script: Extracts networks from grey matter segmentations.
- % Network analyses is done with other scripts.
- %
- % Depends on:
- % - the matlab nifti toolbox, download from: www/mathworks.com/matlabcentral/fileexchange/8797-tools-for-nifti-and-analyze-image
- % - the output created with the 'create_cube_template.m script'.
- % - SPM5 or SPM8 imcalc function
- % - the functions that live in the "ALl_network_scripts" dir:
- % - determine_rois_with_minimum_nz : Determines which template extracts networks of minimum size.
- % - create_rrois : Permutes all the values in the scan.
- % - fast_cross_correlation: Computes the correlation matrices
- %
- % PLEASE ADJUST the following lines:
- % - line 43: add the path to the directory where you want to store the results
- % - line 46: add the full path to the directory were All_network_scripts/functions/ live.
- % - line 49: add the full path + file name of the txt file that contains the dir+names of all the scans (1 line per scan)
- % - line 52: add the full path of the directory were the resliced images need to be written to.
- % - line 55: add the full path of the directory where the spm canonical brain lives.
- % - line 58: add the full path of the directory where /bl_ind directory lives, this directory contains the cube templates (output from 'create_cube_template.m')
- %
- %
- %
- % Output (stored in the subject directory unless specified otherwise):
- % - iso2mm_s (subject number) .img : resliced images (this is stored in designated directory)
- % - creates a directory for each subject (1:n).
- % - Sa = 3D volume that contains the MRI data in a matrix form
- % - Va = header file that contains the scan info
- % - nz = number of cubes (i.e., size of the network)
- % - nancount = number of cubes that have standard deviation of 0. (these are excluded from further calculations)
- % - off_set = the template used to extract the minimum number of cubes.
- % - bind.mat (in single precision): this file contains the indices of the cubes that contain grey matter values
- % - rois.mat and rrois.mat (in single precision): these files contain the cubes
- % - lookup.mat and rlookup.mat : provide the link between the cubes from rois and rrois and the original scans (bind). These files are needed later when results are written in images (e.g., the degree or clustering values of the cubes).
- % - rotcorr.mat and rrcorr.mat (in single precision) are the correlation matrices maximised for with reflection over all angles and rotation with multiples of 45degrees.
- % - th = threshold to binarise the matrices: corresponds to the correlation value where in the random image
- % - fp = percentage of positive random correlations
- % - sp = sparsity of the matrix with threshold p_corrected = 0.05
- %
- % Author: Betty Tijms, 2011, version 12.08.2013
- % -------------------------------------------------------------------------------------------------------------------------------
- %cd /usr/local/MATLAB/R2011a/bin/
- matlab -nodesktop
- % --- ADJUST FOLLOWING ----%
- % Go to the directory where you want to store your output (i.e., networks)
- result_dir ='/ENTER_PATH_NAME_HERE/results';
- cd result_dir
- % FILL IN THE FILE NAME with full path that contains the text file with all the grey matter segmentations
- [CN_a1]=textread('/ENTER_PATH_NAME_HERE/grey_matter_segmentation_name_file.txt','%s');
- % FILL IN THE FULL PATH of the directory where the resliced images should live.
- reslice_dir = '/ENTER_PATH_NAME_HERE/reslice/';
- % Add the path to the Nifti toolbox (see read me file for download link)
- addpath /ENTER_PATH_NAME_HERE/ENTER_NAME_OF_NIFTI_TOOLBOX_HERE/
- % PATH where the canonical SPM image "avg305T1.nii´ can be found --> this is used to reslice the images
- P1= '/ENTER_PATH_TO_SPM/spm8/canonical/avg305T1.nii';
- % Please provide path to the folder Extract_individual_GM_networks
- all_network_scripts_path ='/ENTER_PATH_HERE/Extract_individual_GM_networks_v20131106/';
- % Please provide path to the folder bl_ind
- bl_dir = '/home/jagust/share/connect/networks_struct/Extract_individual_GM_networks_v20131106/bl_ind/';
- %% -----Finished with the adjustments ---%
- % add path to the functions in All_network_scripts folder
- addpath(strcat(all_network_scripts_path,'functions'))
- % Get content of all_dirs.txt files ( this file contains all the dirs of the cube template for all possible off set values)
- CN_a2 =textread(strcat(bl_dir, 'all_dirs.txt'),'%s');
- % end of loop
- numimages=size(CN_a1,1);
- % number of dimensions
- n=3;
- s=n^3;
- % Loop through all the grey matter segmentations listed in CN_a1 and extract network. Takes ~25 minutes per scan.
- for im=1:numimages
- % Get this scan and do the following loop to threshold and then reslice it
- this_scan=char(CN_a1(im));
- % make the directory for this subject
- mkdir(strcat('s', int2str(im),'/data/rotation/'))
- mkdir(strcat('s', int2str(im),'/images/'))
- % Go to this directory
- cd(strcat('s', int2str(im)))
- % which directory are we now?
- t_result_dir = strcat(result_dir, '/s',int2str(im));
- % %% Next reslice the scan to 2x2x2 mm isotropic voxels to reduce amount of data.
- P = strvcat(P1, this_scan);
- Q = strcat(reslice_dir,'/iso2mm_s', num2str(im), '.nii');
- f = 'i2';
- flags = {[],[],[],[]};
- Q = spm_imcalc_ui(P,Q,f,flags);
- %
- % %close the figure window
- close
- % %% Now extract Sa and Va --> Va contains all the info from the .hrd file, Sa is the actual image, it contains the grey matter intensity values.
- Va=spm_vol(strcat(reslice_dir, '/iso2mm_s', num2str(im), '.nii'));
- % Va=spm_vol(this_scan);
- Sa=spm_read_vols(Va);
- % Save them in the current directory
- save Va.mat Va
- save Sa.mat Sa
- % clear variables that aren't needed anymore
- clear this_scan P Q f flags
- %% Get the off_set bl_ind (this corresponds to the indices that make up the 3d cubes) with the minimum number if cubes (min. nz).
- [nz, nan_count, off_set] = determine_rois_with_minimumNZ(CN_a2, bl_dir, n, Sa, t_result_dir);
- % Store nz.m : the number of cubes, this is the size of the network
- save data/nz.mat nz
- % nan_count = number of cubes that have a variance of 0, these are excluded because cannot compute correlation coefficient for these cubes.
- save data/nan_count.mat nan_count
- % off_set indicates which specific template was used to get the cubes.
- save data/off_set.mat off_set
- % Get bind and store bind too. Bind contains the indices of the cubes, so we can efficiently do computations later. It is a long vector of which every 27 consecutive voxels are a cube.
- %tt=strcat(bl_dir, off_set, '/data/bind.m');
- load data/bind.m
- % Convert to single, for memory reasons
- bind=single(bind);
- delete data/bind.m
- save data/bind.mat bind
- % Now create:
- % - the rois: This is a matrix of which each column corresponds to a cube, we will compute use this to compute the correlations
- % - lookup table --> to lookup the cubes that correspond to the correlations.
- lb=length(bind); % End point for loop
- col=1; % iteration counter
- rois=zeros(s,nz, 'single'); % Create variable to store the rois (i.e., cubes)
- lookup=zeros(s,nz, 'single'); % Create variable to store the lookup table --> this table links the cubes to the indices in the original iso2mm image.
- % The next loop gos through bind for indices of voxels that belong to each ROI (i.e. cube)
- % lb is the last voxel that belongs to a roi
- %lookup is lookup table to go from corr index to bind to Sa
- for i=1:s:lb
- rois(:,col)=Sa(bind(i:(i+(s-1))));
- lookup(:,col)=i:(i+(s-1));
- col=col+1;
- end
- % Save the rois and lookup table
- save data/rois.mat rois
- % look up table --> to go from corr index to bind to Sa
- save data/lookup.mat lookup
- %clear unneeded variables
- clear lookup lb
- %% Randomise ROIS: Create a 'random brain' to estimate the threshold with for later stages
- [rrois, rlookup]= create_rrois (rois, n, Sa, Va, off_set, bind, bl_dir, nz);
- % save rrois and rlookup
- save data/rrois.mat rrois
- save data/rlookup.mat rlookup
- % Remove all variables that take space and aren't needed anymore
- clear bind lookup rlookup Sa Va
- % Set the following variable to 2 when correlation is maximised for rotation with angle multiples of 45degrees.
- forty=2;
- [rotcorr, rrcorr] = fast_cross_correlation(rois, rrois, n,forty);
- % Save it
- save data/rotation/rotcorr.mat rotcorr
- save data/rotation/rrcorr.mat rrcorr
- %% Now get the threshold and save it
- [th, fp, sp] = auto_threshold(rotcorr, rrcorr, nz);
- % Add this threshold to all_th
- %all_th(im,:)=[th,fp,sp];
- %save this and get later
- save data/th.mat th
- save data/fp.mat fp
- save data/sp.mat sp
- % Clear all variables
- clear rotcorr rrcorr th fp sp
- cd ..
- im
- end
batch_extract_networks20131106.m at commit 183d3c1, under MIT · at the source
Overview
- Alzheimer Center Amsterdam, Neurology, Vrije Universiteit Amsterdam, Amsterdam UMC Location VUmc, 1081HV Amsterdam, The Netherlands
- Amsterdam Neuroscience, Neurodegeneration, 1081HV Amsterdam, The Netherlands
- Department of Neurology, Mayo Clinic, Rochester, MN 55905, USA
- Department of Laboratory Medicine, Neurochemistry Lab, Amsterdam UMC, Vrije Universiteit Amsterdam, 1081HV Amsterdam, The Netherlands
- Department of Radiology and Nuclear Medicine, Amsterdam UMC, Vrije Universiteit, 1081HV Amsterdam, The Netherlands
- Queen Square Institute of Neurology and Centre for Medical Image Computing, University College London, WC1N 3BG London, UK
Abstract
Grey matter network topology is altered in Alzheimer’s disease and these alterations are related to cognitive decline. Understanding the biological underpinnings of loss of brain connectivity may provide insights into mechanisms related to developing Alzheimer’s dementia (i.e. dementia A+). We investigated which biological processes as measured in CSF proteomics were associated with loss of brain connections across the Alzheimer’s disease continuum. We included 347 individuals with abnormal CSF amyloid [mean age ± standard deviation (SD) 66 ± 8; 98 cognitively unimpaired—A+, 88 mild cognitive impairment—A+, 161 dementia A+] and 146 cognitively unimpaired individuals with normal CSF amyloid (mean age ± SD 62 ± 8) and available T1w MRI-scans and CSF proteomic data (3097 proteins using tandem mass tag spectrometry) from the Amsterdam Dementia Cohort. We used an automated pipeline to construct grey matter networks from 3D-T1 sequences and for each network, calculated the small-worldness coefficient, which we previously found to be robustly related to cognitive decline. Linear models were applied to test associations between CSF protein levels and connectivity measures using an interaction term for clinical stage while controlling for connectivity density, age and sex. We validated our results in data from the Alzheimer’s disease Neuroimaging Initiative (ADNI). Pathway enrichment analysis was performed for proteins associated with loss of brain connectivity (P < 0.05) using the Gene Ontology database. Individuals across the Alzheimer’s disease continuum had lower small-worldness coefficients compared with controls (ANOVA P < 0.001). In amyloid positive individuals, higher levels of 222 proteins and lower levels of 482 proteins were associated with lower small-worldness coefficients and were enriched for innate immune system and neuroplasticity pathways, respectively. Stratified for disease stage, most protein associations with lower small-worldness coefficients were found in mild cognitive impairment A+ (n = 527 proteins) and dementia A+ (n = 799 proteins) with considerable overlap (n = 239 proteins). Proteins in these stages were enriched for complement activation and synaptic integrity. In cognitive unimpairment A+, we found proteins enriched for processes involved in apoptosis. We did not find any enriched biological processes in controls. Repeating analyses in ADNI indicated that similar biological processes were associated with altered grey matter network connectivity. Higher CSF levels of proteins involved in immune responses and lower levels of proteins related to neuroplasticity were associated with lower small-worldness coefficients across the Alzheimer’s disease continuum. This suggests that preserving cognitive function in the presence of amyloid and prevention of dementia A+ may require therapies that strengthen synapses and targets the innate immune system in addition to amyloid and tau.
Reproduced under the paper's license (CC BY), from the paper cited above.
Repositories
Its files are read in the Code ↔ Paper reader above, with 2 matches between paragraphs and lines of code.
bettytijms/Single_Subject_Grey_Matter_Networks
183d3c1d7eaddfb62c056b4b32e9cc1b4ee93e8f, 4 June 2020Availability: 1 check, the latest on 27 September 2026: the link answers
- 27 September 2026: the link answers
15 files
- Extract_individual_GM_ne
tworks_v20131106/ — MATLAB, 208 lines, 2 matchesbatch_extract_networks20 131106.m - Extract_individual_GM_ne
tworks_v20131106/ — MATLAB, 119 linescreate_cube_templates.m - Extract_individual_GM_ne
tworks_v20131106/ — MATLAB, 38 linesfunctions/ auto_threshold.m - Extract_individual_GM_ne
tworks_v20131106/ — MATLAB, 48 linesfunctions/ create_rrois.m - Extract_individual_GM_ne
tworks_v20131106/ — MATLAB, 131 linesfunctions/ determine_rois_with_mini mumNZ.m - Extract_individual_GM_ne
tworks_v20131106/ — MATLAB, 106 linesfunctions/ fast_cross_correlation.m - Extract_individual_GM_ne
tworks_v20131106/ — MATLAB, 119 linesfunctions/ rotation_xyz.m - Extract_individual_GM_ne
tworks_v20150902/ — MATLAB, 119 linescreate_cube_templates.m - Extract_individual_GM_ne
tworks_v20150902/ — MATLAB, 38 linesfunctions/ auto_threshold.m - Extract_individual_GM_ne
tworks_v20150902/ — MATLAB, 48 linesfunctions/ create_rrois.m - Extract_individual_GM_ne
tworks_v20150902/ — MATLAB, 131 linesfunctions/ determine_rois_with_mini mumNZ.m - Extract_individual_GM_ne
tworks_v20150902/ — MATLAB, 106 linesfunctions/ fast_cross_correlation.m - Extract_individual_GM_ne
tworks_v20150902/ — MATLAB, 119 linesfunctions/ rotation_xyz.m - LICENSE.md — License, 24 lines
- README.md — Text, 47 lines
sites.google.com/site/bctnet
Availability: 1 check, the latest on 27 September 2026: the link answers (HTTP 200)
- 27 September 2026: the link answers (HTTP 200)
Tracing map
Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.
What the map holds:
- 2 repositories of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
- 13 scripts, each with its path and the digest of its content;
- 2 matches between paragraphs of the paper and lines of the code (method lexical-v1);
- neither the text of the paper nor the code itself.
Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.
Data
No dataset and no data link were found in the paper.
Data availability
The data from the ADC cohort that support the findings of this study are available from the corresponding author, upon reasonable request. The data from the ADNI cohort can be requested online at https://
Reproduced under the paper's license (CC BY), from the paper cited above.
Versions
The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.
Version 1, 27 September 2026: the first record
Recorded: type, language, journal, volume, issue, pages, dates, 9 authors, 3 keywords, 2 funders, 83 references.
Cite
This paper
de Leeuw, D. M., Duits, F. H., Dicks, E., Vromen, E. M., Teunissen, C. E., Barkhof, F., van der Flier, W. M., Visser, P. J., & Tijms, B. M. (2026). Neuroplasticity and immune system are related to altered grey matter networks: a cohort study in sporadic Alzheimer's disease. Brain communications, 8(4), fcag257. https://
BibTeX
@article{deleeuw2026neur
author = {de Leeuw, Diederick M and Duits, Flora H and Dicks, Ellen and Vromen, Eleonora M and Teunissen, Charlotte E and Barkhof, Frederik and van der Flier, Wiesje M and Visser, Pieter Jelle and Tijms, Betty M},
title = {{Neuroplasticity and immune system are related to altered grey matter networks: a cohort study in sporadic Alzheimer's disease}},
journal = {Brain communications},
year = {2026},
month = jul,
volume = {8},
number = {4},
pages = {fcag257},
publisher = {Oxford University Press},
issn = {2632-1297},
doi = {10.1093/
url = {https://
pmid = {42465730},
pmcid = {PMC13373788}
}
RIS
TY - JOUR
AU - de Leeuw, Diederick M
AU - Duits, Flora H
AU - Dicks, Ellen
AU - Vromen, Eleonora M
AU - Teunissen, Charlotte E
AU - Barkhof, Frederik
AU - van der Flier, Wiesje M
AU - Visser, Pieter Jelle
AU - Tijms, Betty M
TI - Neuroplasticity and immune system are related to altered grey matter networks: a cohort study in sporadic Alzheimer's disease
T2 - Brain communications
J2 - Brain Commun
PY - 2026
DA - 2026/
VL - 8
IS - 4
SP - fcag257
SN - 2632-1297
PB - Oxford University Press
DO - 10.1093/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1093/
"type": "article-journal",
"title": "Neuroplasticity and immune system are related to altered grey matter networks: a cohort study in sporadic Alzheimer's disease",
"container-title": "Brain communications",
"author": [
{
"family": "de Leeuw",
"given": "Diederick M"
},
{
"family": "Duits",
"given": "Flora H"
},
{
"family": "Dicks",
"given": "Ellen"
},
{
"family": "Vromen",
"given": "Eleonora M"
},
{
"family": "Teunissen",
"given": "Charlotte E"
},
{
"family": "Barkhof",
"given": "Frederik"
},
{
"family": "van der Flier",
"given": "Wiesje M"
},
{
"family": "Visser",
"given": "Pieter Jelle"
},
{
"family": "Tijms",
"given": "Betty M"
}
],
"container-title-short":
"volume": "8",
"issue": "4",
"page": "fcag257",
"DOI": "10.1093/
"PMID": "42465730",
"PMCID": "PMC13373788",
"ISSN": "2632-1297",
"publisher": "Oxford University Press",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
7,
8
]
]
}
}
The tracing map gets a citation of its own once an author has validated it and it has a DOI.
Similar papers
The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.
- [1] doi:10.1093/braincomms/fcag196 [code]
- Music processing in behavioural variant frontotemporal dementia and Alzheimer's disease: a functional MRI study.Journal: Brain communicationsIn common: SPM, Alzheimer's / dementia, 3 references, 2 authors
- [2] doi:10.64898/2026.06.03.26354793 [code]
- Bayesian networks to estimate prognosis in vascular cognitive impairment and small vessel disease: integrated analyses of interdependent contributors to multiple outcomesJournal: medRxiv (preprint)In common: Alzheimer's / dementia, 2 authors
- [3] doi:10.1093/braincomms/fcag341
- Distinct and shared multimodal neuroimaging patterns in clinical frontotemporal dementia syndromes.Journal: Brain communicationsIn common: Alzheimer's / dementia, 5 references
- [4] doi:10.1093/braincomms/fcag176 [code]
- Tau topography subtypes account for clinical heterogeneity and longitudinal trajectories in early-onset Alzheimer's disease.Journal: Brain communicationsIn common: SPM, Alzheimer's / dementia, 4 references
- [5] doi:10.1038/s43856-026-01707-2 [code]
- Decreased amyloid-related structure-function coupling in preclinical Alzheimer's disease.Journal: Communications medicineIn common: SPM, Alzheimer's / dementia, author Frederik Barkhof
- [6] doi:10.1038/s41598-026-51725-2 [code]
- AI-powered evaluation of dementia severity based on clinical data and visual scoring systems (MTA, ERICA, GCA) from MRI.Journal: Scientific reportsIn common: Alzheimer's / dementia, 4 references
- [7] doi:10.1016/j.ynirp.2026.100351
- 3D quantitative MRI for cognitive performance: T1 relaxation times of the right putamen and volumes of the left hippocampus as key biomarkers.Journal: Neuroimage. ReportsIn common: Alzheimer's / dementia, 3 references
- [8] doi:10.1038/s41598-026-56387-8 [code]
- Alterations in topological and dynamical parameters correlate with disease biomarkers and neuropsychological scores in prodromic stages of dementia.Journal: Scientific reportsIn common: Alzheimer's / dementia, 3 references
- [9] doi:10.1002/alz.71609
- Cortical gray-white matter contrast alterations precede amyloid-β positivity and macrostructural changes in older adults without dementia.Journal: Alzheimer's & dementia : the journal of the Alzheimer's AssociationIn common: Alzheimer's / dementia, author Frederik Barkhof
- [10] doi:10.1093/braincomms/fcag169 [code]
- Clinicoanatomic localization of iron-rich gliosis in aphasic presentations of globular glial tauopathy.Journal: Brain communicationsIn common: Alzheimer's / dementia, 3 references
Contribute
The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.
Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.
Claim this paper
Correct its record
Say what each link of this record is, remove the ones that are not the paper's, add the ones that are missing. The correction becomes a new version of the record, in its Versions section.
Validate its tracing map
You validate the map as this page shows it: 2 repositories of the authors' code, each at its verified commit and with its license, 13 scripts, and 2 matches between paragraphs and code (see the Code and Map sections). It then receives a DOI on Zenodo, with you (your ORCID iD) and OSCR as its creators; the code itself is not deposited.
The map's fingerprint: sha256:ef8f6d9841488fe0…
Add the badge to its README
The badge links the code to this page. Copy one of these into the README of the paper's code: only you decide where it goes, and nothing is changed for you.
Markdown
[, paste the snippet at the top, then “Commit changes…” and, to review it first, “Create a new branch and start a pull request”. You open the pull request; OSCR asks for no permission.
Request its removal
To ask OSCR to remove this record, the copies of its authors' scripts or its tracing map, use the removal request page: signed in, you say who you are, what to remove and why, then review and confirm the request. Published rules decide every request (how).
Discussion, reproductions, activity
Discussion: questions and error reports about this paper and its code, from signed-in readers and its authors. It opens with sign-in.
Reproductions: reports from readers who ran the authors' code: what they reproduced, with which environment, commit and data. It opens with sign-in.
Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.
