The genetic architecture of cortical similarity networks.
The 5 matches
- [1] § Methods › Magnetic resonance imaging (MRI) data › Structural MRI feature selection ↔ get_vertex_df.py, lines 11–57 · score 0.67 · fractional anisotropy, FA, SA, volume, MC, thickness
- [2] § Methods › Magnetic resonance imaging (MRI) data › Structural MRI feature selection ↔ MIND.py, lines 8–49 · score 0.65 · surface area, cortical thickness, Vol, MIND network, biological, SA
- [3] § Methods › Magnetic resonance imaging (MRI) data › Structural MRI feature selection ↔ MIND.py, lines 8–49 · score 0.56 · surface area, Cortical thickness, summing, volume, vertex
- [4] § Methods › Magnetic resonance imaging (MRI) data › Structural MRI feature selection ↔ MIND.py, lines 8–49 · score 0.56 · surface area, Cortical thickness, summing, volume, vertex
- [5] § Methods › Magnetic resonance imaging (MRI) data › Structural MRI feature selection ↔ get_vertex_df.py, lines 11–57 · score 0.53 · FreeSurfer, vertex, MC, CT, thickness, surface
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
Python · 203 lines · 9.4 KB · no license · 2 matches
- import sys
- import numpy as np
- import os
- import pandas as pd
- from os.path import exists
- from nibabel.freesurfer.io import read_morph_data, read_annot
- from nibabel.freesurfer.mghformat import load
- from collections import defaultdict
- from MIND_helpers import calculate_mind_network, is_outlier
- def get_vertex_df(surf_dir, features, parcellation):
- '''
- INPUT SPECIFICATIONS:
- • surf_dir (str) : This is a string the location containing all relevant directories output by FreeSurfer (i.e. label, mri, surf).
- • features (list):
- This function accepts the "features" argument as a list containing items in the following forms:
- str:
- • One of ['CT','Vol','SA','MC','SD']. In this form, the function will automatically assume that the requested features are found in the surf_dir/surf directory and correspond to the following files:
- CT: ?h.thickness
- Vol: ?h.volume
- SA: ?h.area
- MC: ?h.curv
- SD: ?h.sulc
- • You may also pass a string in the form 'thickness', 'volume', 'sulc', etc, which refer directly to default files already found in the surf_dir/surf directory for both rh and lh.
- For example, the function will interpret the entry 'thickness' to refer to the files surf_dir/surf/lh.thickness and surf_dir/surf/rh.thickness
- • Finally, you may pass a string in the form of a FULL path as follows: with a question mark '?' indicating the position specifying the hemisphere such as "full/path/to/?h.feature".
- Using this formulation will cause the command to look for files that exactly match both "full/path/to/lh.feature" and "full/path/to/rh.feature". If the left and right version of the files aren't exactly the same otherwise, this won't work.
- tuple: (path/to/lh_surface_feature, path/to/rh_surface_feature)
- • If you would rather pass other features directly into the function, you must specify the locations (using paths) of both the left and right versions of each desired feature as a tuple.
- **The files must be readable by nibabel's read_morph_data function, i.e. in FreeSurfer's surface format!***
- So for example, if you have used the provided register_and_vol2surf function to generate surface maps of fractional anisotropy in a separate folder, you could pass them as an element in the list like:
- (path/to/lh.FA.mgh, path/to/rh.FA.mgh)
- A valid list of feature values combining these different input types would therefore be: ['CT','SD',(path/to/lh_feature1, path/to/rh_feature1), (path/to/lh_feature2, path/to/rh_feature2)]
- • parcellation (str): This is a string the location containing parcellation scheme to be used. The files 'lh.' + parcellation + '.annot' and 'rh.' + parcellation + '.annot' must exist inside the surf_dir/label directory.
- '''
- #specify data locations
- surfer_location = surf_dir + '/'
- #Check inputs!
- if (exists(surfer_location + '/label/lh.' + parcellation + '.annot') == False) or (exists(surfer_location + '/label/rh.' + parcellation + '.annot') == False):
- raise Exception('Parcellation files not found.')
- all_shorthand_features = ['CT','Vol','SA','MC','SD']
- all_shorthand_features_dict = dict(zip(all_shorthand_features, ['thickness','volume','area','curv','sulc']))
- lh_feature_locs = []
- rh_feature_locs = []
- #Check feature inputs, store location of files
- for feature in features:
- if feature in all_shorthand_features:
- lh_loc = surfer_location + 'surf/lh.' + all_shorthand_features_dict[feature]
- rh_loc = surfer_location + 'surf/rh.' + all_shorthand_features_dict[feature]
- if (exists(lh_loc) == False) or (exists(rh_loc) == False):
- raise Exception('Feature for input "' + feature +'" not found.')
- else:
- lh_feature_locs.append(lh_loc)
- rh_feature_locs.append(rh_loc)
- elif type(feature) is str:
- if len(feature.split('/')) == 1:
- lh_loc = surfer_location + 'surf/lh.' + feature
- rh_loc = surfer_location + 'surf/rh.' + feature
- if (exists(lh_loc) == False) or (exists(rh_loc) == False):
- raise Exception('Feature for input "' + feature +'" not found.')
- else:
- lh_feature_locs.append(lh_loc)
- rh_feature_locs.append(rh_loc)
- else:
- lh_loc = feature.split('/')
- lh_loc[-1] = 'l' + lh_loc[-1][1:]
- lh_loc = '/'.join(lh_loc)
- rh_loc = feature.split('/')
- rh_loc[-1] = 'r' + rh_loc[-1][1:]
- rh_loc = '/'.join(rh_loc)
- if (exists(lh_loc) == False) or (exists(rh_loc) == False):
- raise Exception('Feature for input "' + feature +'" not found.')
- else:
- lh_feature_locs.append(lh_loc)
- rh_feature_locs.append(rh_loc)
- elif type(feature) is tuple:
- lh_loc = feature[0]
- rh_loc = feature[1]
- if (exists(lh_loc) == False) or (exists(rh_loc) == False):
- raise Exception('Feature for input "' + feature[0] +'" or "' + feature[1] +'" not found.')
- else:
- lh_feature_locs.append(lh_loc)
- rh_feature_locs.append(rh_loc)
- else:
- raise Exception('Unrecognized format for feature input: ', feature)
- #Get annotation files
- lh_annot = read_annot(surfer_location + '/label/lh.' + parcellation + '.annot', orig_ids = True)
- rh_annot = read_annot(surfer_location + '/label/rh.' + parcellation + '.annot', orig_ids = True)
- annot_dict = {'lh':lh_annot, 'rh':rh_annot}
- '''
- The regions in the lh and rh need to be renamed and distinct.
- So, here we append lh_ or rh_ to the front of each region name and make a conversion dict.
- This will likely need to change for different processing pipelines (FS versions) and datasets etc.
- so make sure this dict is correct and looks good.
- '''
- lh_region_names = ['lh_' + str(x).split("'")[1] for x in lh_annot[2]]
- rh_region_names = ['rh_' + str(x).split("'")[1] for x in rh_annot[2]]
- lh_convert_dict = dict(zip(lh_annot[1][:,-1], lh_region_names))
- rh_convert_dict = dict(zip(rh_annot[1][:,-1], rh_region_names))
- convert_dicts = {'lh': lh_convert_dict,\
- 'rh': rh_convert_dict}
- used_labels_l = np.intersect1d(np.unique(lh_annot[0]), list(lh_convert_dict.keys()))
- used_labels_r = np.intersect1d(np.unique(rh_annot[0]), list(rh_convert_dict.keys()))
- used_labels = {'lh': used_labels_l,\
- 'rh': used_labels_r}
- used_regions_l = np.array([value for key, value in lh_convert_dict.items() if key in used_labels_l])
- used_regions_r = np.array([value for key, value in rh_convert_dict.items() if key in used_labels_r])
- combined_regions = np.hstack((used_regions_l, used_regions_r))
- unknown_regions = [x for x in combined_regions if (('?' in x) | ('unknown' in x) | ('Unknown' in x) | ('Medial_Wall' in x) | (len(x) == 3))]
- combined_regions = np.array([x for x in combined_regions if x not in unknown_regions])
- vertex_data_dict = defaultdict()
- #Now load up all the vertex-level data!
- for hemi in ['lh','rh']:
- print(hemi)
- hemi_data_dict = defaultdict()
- if hemi == 'lh':
- print('Loading left hemisphere data:')
- for i, lh_feature_loc in enumerate(lh_feature_locs):
- print(lh_feature_loc)
- #check for mgh/mgz format vs regular curv files
- if lh_feature_loc.endswith('mgh') or lh_feature_loc.endswith('mgz'):
- hemi_data_dict['Feature_' + str(i)] = load(lh_feature_loc).get_fdata().flatten()
- else:
- hemi_data_dict['Feature_' + str(i)] = read_morph_data(lh_feature_loc)
- elif hemi == 'rh':
- print('Loading right hemisphere data:')
- for i, rh_feature_loc in enumerate(rh_feature_locs):
- if rh_feature_loc.endswith('mgh') or rh_feature_loc.endswith('mgz'):
- hemi_data_dict['Feature_' + str(i)] = load(rh_feature_loc).get_fdata().flatten()
- else:
- hemi_data_dict['Feature_' + str(i)] = read_morph_data(rh_feature_loc)
- used_features = list(hemi_data_dict.keys())
- print(used_features)
- hemi_data = np.zeros((len(used_features) + 1, len(annot_dict[hemi][0])))
- hemi_data[0] = annot_dict[hemi][0]
- for i, feature in enumerate(used_features):
- print(i, feature)
- hemi_data[i + 1] = hemi_data_dict[feature]
- col_names = ['Label'] + used_features
- hemi_data = pd.DataFrame(hemi_data.T, columns = col_names)
- #Select only the vertices that map to regions.
- hemi_data = hemi_data.loc[hemi_data['Label'].isin(used_labels[hemi])]
- hemi_data["Label"] = hemi_data["Label"].map(convert_dicts[hemi])
- vertex_data_dict[hemi] = hemi_data
- vertex_data = pd.concat([vertex_data_dict['lh'], vertex_data_dict['rh']], ignore_index = True)
- #Output data
- print("features used: ")
- print(used_features)
- return vertex_data, combined_regions, used_features
get_vertex_df.py at commit 0d33445, no license · at the source
Overview
- Department of Psychiatry, University of Cambridge,Cambridge, United Kingdom
- Department of Computer Science and Technology, University of Cambridge,Cambridge, United Kingdom
- Harvard Medical School,Boston, MA US
- Department of Psychology, University of Cambridge,Cambridge, United Kingdom
- Autism Research Centre, University of Cambridge,Cambridge, United Kingdom
- Instituto de Biomedicina de Sevilla (IBiS) HUVR/CSIC/Universidad de Sevilla/CIBERSAM, ISCIII, Departamento de Fisiología Médica y Biofísica,Seville, Spain
- Department of Psychiatry, University of Pennsylvania,Philadelphia, PA USA
- Department of Child and Adolescent Psychiatry and Behavioral Science, The Children’s Hospital of Philadelphia,Philadelphia, PA USA
- Lifespan Brain Institute, The Children’s Hospital of Philadelphia,Philadelphia, PA USA
- School of Academic Psychiatry, Institute of Psychiatry, Psychology & Neuroscience, King’s College London,London, United Kingdom
- School of Biomedical Engineering & Imaging Sciences, King’s College London,London, United Kingdom
Abstract
The genetic architecture of human brain networks is central to understanding cortical organisation and evolution, the causal links between brain structure and function, and the pathogenesis of neuropsychiatric disorders. Using N > 48,000 subjects, we investigated common genetic effects on Morphometric INverse Divergence (MIND), a heritable, multi-modal structural MRI metric of inter-areal similarity and connectivity. Genetic correlations between MIND network edges were largely reducible to two gradients, each aligned with distance from one of the two phylogenetically primitive areas (paleocortex and archicortex) predicted by the dual origin theory of cortical evolution. MIND was more heritable than comparable measures of functional (f)MRI connectivity, and the paleocortically-aligned MIND gradient was genetically correlated with, and causally predictive of, fMRI connectivity. Finally, we identified genetic overlaps between MIND gradients and neuropsychiatric and biomedical traits. These results provide fresh insight into the dual origins of the cortex and their implications for brain function and health.
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 5 matches between paragraphs and lines of code.
jrs95/geni.plots
aa25803629235b681ada96741a1b0c07dc059870, 14 October 2024Availability: 1 check, the latest on 27 September 2026: the link answers
- 27 September 2026: the link answers
9 files
- R/
fig_manhattan.R , R, 1,032 lines - R/
fig_phewas.R , R, 674 lines - R/
fig_qq.R , R, 831 lines - R/
fig_region.R , R, 2,253 lines - R/
geni.plots-package.R , R, 7 lines - R/
geni_test.R , R, 221 lines - vignettes/
geni_plots.Rmd , R, 164 lines - LICENSE, License, 674 lines
- README.md, Text, 22 lines
Zenodo 7974716
Availability: 1 check, the latest on 27 September 2026: the link answers (HTTP 200)
- 27 September 2026: the link answers (HTTP 200)
7 files
- .ipynb_checkpoints/
ABCD-MIND-checkpoint.ipy , Jupyter, 429 linesnb - .ipynb_checkpoints/
ABCD-MSN-and-raw-feature , Jupyter, 668 liness-checkpoint.ipynb - MIND.py, Python, 52 lines, 1 match
- MIND_helpers.py, Python, 125 lines
- get_vertex_df.py, Python, 203 lines
- register_and_vol2surf.py
, Python, 115 lines - README.md, Text, 132 lines
isebenius/mind
0d334454eaac62e49a197801446ce7750256c882, 12 March 2024Availability: 1 check, the latest on 27 September 2026: the link answers
- 27 September 2026: the link answers
7 files
- .ipynb_checkpoints/
ABCD-MIND-checkpoint.ipy , Jupyter, 429 linesnb - .ipynb_checkpoints/
ABCD-MSN-and-raw-feature , Jupyter, 668 liness-checkpoint.ipynb - MIND.py, Python, 52 lines, 2 matches
- MIND_helpers.py, Python, 149 lines
- get_vertex_df.py, Python, 203 lines, 2 matches
- register_and_vol2surf.py
, Python, 115 lines - README.md, Text, 135 lines
Code availability
Code to calculate MIND networks can be accessed at 10.5281/
Reproduced under the paper's license (CC BY), from the paper cited above.
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:
- 3 repositories of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
- 19 scripts, each with its path and the digest of its content;
- 5 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
All summary statistics generated in this study are publicly available at 10.17863/
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, 10 authors, 2 keywords, 7 MeSH terms, 2 funders, 108 references.
Cite
This paper
Sebenius, I., Warrier, V., Bethlehem, R. A. I., Dear, R., Stauffer, E.-M., Gu, Y., Romero-Garcia, R., Seidlitz, J., Bullmore, E., & Morgan, S. (2026). The genetic architecture of cortical similarity networks. Nature communications, 17(1), 7801. https://
BibTeX
@article{sebenius2026gen
author = {Sebenius, Isaac and Warrier, Varun and Bethlehem, Richard A. I. and Dear, Richard and Stauffer, Eva-Maria and Gu, Yuanjun and Romero-Garcia, Rafael and Seidlitz, Jakob and Bullmore, Edward and Morgan, Sarah},
title = {{The genetic architecture of cortical similarity networks}},
journal = {Nature communications},
year = {2026},
month = jun,
volume = {17},
number = {1},
pages = {7801},
publisher = {Nature Publishing Group},
issn = {2041-1723},
doi = {10.1038/
url = {https://
pmid = {42323287},
pmcid = {PMC13438111}
}
RIS
TY - JOUR
AU - Sebenius, Isaac
AU - Warrier, Varun
AU - Bethlehem, Richard A. I.
AU - Dear, Richard
AU - Stauffer, Eva-Maria
AU - Gu, Yuanjun
AU - Romero-Garcia, Rafael
AU - Seidlitz, Jakob
AU - Bullmore, Edward
AU - Morgan, Sarah
TI - The genetic architecture of cortical similarity networks
T2 - Nature communications
J2 - Nat Commun
PY - 2026
DA - 2026/
VL - 17
IS - 1
SP - 7801
SN - 2041-1723
PB - Nature Publishing Group
DO - 10.1038/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1038/
"type": "article-journal",
"title": "The genetic architecture of cortical similarity networks",
"container-title": "Nature communications",
"author": [
{
"family": "Sebenius",
"given": "Isaac"
},
{
"family": "Warrier",
"given": "Varun"
},
{
"family": "Bethlehem",
"given": "Richard A. I."
},
{
"family": "Dear",
"given": "Richard"
},
{
"family": "Stauffer",
"given": "Eva-Maria"
},
{
"family": "Gu",
"given": "Yuanjun"
},
{
"family": "Romero-Garcia",
"given": "Rafael"
},
{
"family": "Seidlitz",
"given": "Jakob"
},
{
"family": "Bullmore",
"given": "Edward"
},
{
"family": "Morgan",
"given": "Sarah"
}
],
"container-title-short":
"volume": "17",
"issue": "1",
"page": "7801",
"DOI": "10.1038/
"PMID": "42323287",
"PMCID": "PMC13438111",
"ISSN": "2041-1723",
"publisher": "Nature Publishing Group",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
6,
20
]
]
}
}
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.1038/s42003-026-10276-y [code]
- The cellular correlates and adolescent reorganisation of cortical myelination networks in the common marmoset.Journal: Communications biologyIn common: Nipype, AFNI, FreeSurfer, 7 other tools, structural MRI / diffusion, 10 references, 2 authors
- [2] doi:10.64898/2026.03.06.709690
- Genetic insights on the mechanisms of human cortical foldingJournal: bioRxiv (preprint)In common: genetics / omics, 20 references, author Ed Bullmore
- [3] doi:10.1038/s41467-026-74153-2 [code]
- Regional, functional and transcriptomic decoding of multidimensional brain structure alterations in obsessive-compulsive disorder.Journal: Nature communicationsIn common: Nipype, AFNI, FreeSurfer, 7 other tools, genetics / omics, 7 references, author Isaac Sebenius
- [4] doi:10.1038/s41562-026-02486-5 [code]
- Genome-wide association studies of infant and toddler temperament in European and multi-ancestry populations.Journal: Nature human behaviourIn common: ggplot2, tidyverse, pandas, 3 other tools, genetics / omics, 9 references, author Varun Warrier
- [5] doi:10.1162/netn.a.546 [code]
- Cortical similarity networks in the rat brain: Postnatal development and sensitivity to early life stress.Journal: Network neuroscience (Cambridge, Mass.)In common: AFNI, patchwork, tidyverse, 3 other tools, structural MRI / diffusion, 6 references, author Ed Bullmore
- [6] doi:10.1038/s41593-026-02289-x [code]
- Cerebellar aging is spatially heterogeneous and supports cognitive resilience in later life.Journal: Nature neuroscienceIn common: patchwork, ggplot2, tidyverse, structural MRI / diffusion, 1 reference, 3 authors
- [7] doi:10.1038/s41467-026-73996-z [code]
- Genetic architecture of white matter microstructure captured by unsupervised deep representation learning of fractional anisotropy maps.Journal: Nature communicationsIn common: NiBabel, ggplot2, tidyverse, 4 other tools, structural MRI / diffusion, genetics / omics, 9 references
- [8] doi:10.7554/elife.103097 [code]
- Canonical neurodevelopmental trajectories of structural and functional manifolds.Journal: eLifeIn common: NiBabel, ggplot2, tidyverse, 3 other tools, structural MRI / diffusion, 5 references, author Richard Bethlehem
- [9] doi:10.1038/s42003-026-10131-0 [code]
- Shared genetic architecture between the topology of brain white matter structural connectome and fluid intelligence.Journal: Communications biologyIn common: pandas, SciPy, NumPy, structural MRI / diffusion, genetics / omics, 10 references
- [10] doi:10.1073/pnas.2609814123 [code]
- Genetic architectures of brain-related traits are shaped by strong selective constraints.Journal: Proceedings of the National Academy of Sciences of the United States of AmericaIn common: patchwork, ggplot2, tidyverse, 4 other tools, genetics / omics, 7 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: 3 repositories of the authors' code, each at its verified commit and with its license, 19 scripts, and 5 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:f08216cf0018e1ee…
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.
