OSCR

Brain functional-structural gradient coupling reflects development, behavior and genetic influences.

Code ↔ Paper

20 matches between paragraphs of the paper and lines of its authors' code, computed by the harvester (lexical-v1). Click a colored paragraph or line to see its counterpart.

The 20 matches · 3 of them tie a paragraph to a whole file, not to given lines: weak matches, whose lines are not tinted
  1. [1] § Methods › Functional and structural gradient construction ↔ brainspace/gradient/embedding.py, lines 57–198 · score 0.89 · anisotropic diffusion parameter, random walk, affinity matrices, diffusion map, operator, power
  2. [2] § Methods › Data harmonization ↔ Harmonization/main_BatchEffect.R, lines 210–267 · score 0.84 · ComBat, ABCD HCP, site harmonization, HCP YA, batch, fd
  3. [3] § Methods › Imaging processing and connectivity construction ↔ ADNI_example/adni_fmri_correction.sh, lines 50–115 · score 0.78 · phase encode, fMRI, distortion correction, Philips, echo, head
  4. [4] § Methods › Imaging processing and connectivity construction ↔ ADNI_example/preproc_step5_fmri.sh, lines 1–43 · score 0.77 · distortion correction, motion correction, MRI preprocessing, fMRI, slices, Siemens
  5. [5] § Methods › Imaging processing and connectivity construction ↔ ACT_example/dti_rpe_minimal_process.sh, lines 84–122 · score 0.73 · bias field correction, multi shell, EPI, Minimal, intensities, eddy
  6. [6] § Results › SFGC reveals regional and developmental genetic topographies ↔ Heritability-Analysis/imging_transcriptomics.py, lines 1–68 · score 0.69 · LH inferior, LH lateral, lingual, precuneus, frontal, superior
  7. [7] § Methods › Imaging processing and connectivity construction ↔ ADNI_example/adni_dwi_eddy_correction.sh, lines 1–41 · score 0.68 · head motion, fMRI, undistorted, eddy, T1w, distortions
  8. [8] § Methods › Imaging processing and connectivity construction ↔ ACT_example/dti_rpe_minimal_process.sh, lines 1–42 · score 0.66 · phase encoding, Diffusion imaging, platform, Minimal, PA, AP
  9. [9] § Methods › Heritability estimation using the ACE model ↔ Heritability-Analysis/ACE_heritability.py, the whole file · a weak match · score 0.65 · shared environment matrices, ACE model, covariates, phenotype, Yeo, genetic
  10. [10] § Methods › Imaging processing and connectivity construction ↔ pipeline_scripts/sbci_step5_atlas_connectivity.sh, the whole file · a weak match · score 0.62 · continuous SC, continuous FC, SC matrix, connections, smoothed, mesh
  11. [11] § Methods › Imaging transcriptomics of gradient coupling heritability ↔ Heritability-Analysis/imging_transcriptomics.py, lines 1–68 · score 0.62 · imaging transcriptomic, Gene expression, abagen, AHBA, heritability, atlas
  12. [12] § Methods › Heritability estimation using the ACE model ↔ Heritability-Analysis/ACE_heritability.py, the whole file · a weak match · score 0.61 · covariate matrix, ACE model, h2, Phenotypes, heritability
  13. [13] § Methods › Associations between SFGC and behavioral outcomes ↔ Association-Study/main_MLP.py, lines 226–309 · score 0.61 · Sigmoid, dropout, leaky, epochs, layers, optimize
  14. [14] § Methods › Imaging processing and connectivity construction ↔ ADNI_example/adni_fmri_correction.sh, lines 50–115 · score 0.56 · phase encoding, echo, head, TR, Siemens, space
  15. [15] § Methods › Imaging transcriptomics of gradient coupling heritability ↔ brainspace/null_models/variogram.py, lines 1620–1674 · score 0.54 · spin permutations, spatial autocorrelation, empirical, map
  16. [16] § Results › SFGC and mental health associations by age and sex ↔ Association-Study/main_KRR.py, lines 172–218 · score 0.53 · kernel ridge regression, KRR, train, predictive, fitting, modeling
  17. [17] § Methods › Model performance and feature importance ↔ Association-Study/main_KRR.py, lines 172–218 · score 0.53 · primal coefficients, KRR, kernel, predictive, fitting, model
  18. [18] § Methods › Model performance and feature importance ↔ brainspace/gradient/kernels.py, lines 44–98 · score 0.52 · correlation coefficients, kernel functions, vector, gradient
  19. [19] § Methods › Imaging transcriptomics of gradient coupling heritability ↔ Heritability-Analysis/imging_transcriptomics.py, lines 127–167 · score 0.52 · ranked gene, Scanpy, heritability, transcriptomics
  20. [20] § Methods › Imaging processing and connectivity construction ↔ ACT_example/dti_rpe_minimal_process.sh, lines 84–122 · score 0.51 · spherical deconvolution, tissue, field, map

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

Shell · 228 lines · 12 KB · no license · 3 matches

  1. #!/bin/bash
  2. #SBATCH -t 3-3:00:00
  3. #SBATCH --mem-per-cpu=32gb
  4. # Written by Andrew Jahn, University of Michigan, 02.25.2019
  5. # Updated 07.10.2020 to incorporate changes from MRtrix version 3.0.1
  6. # Based on Marlene Tahedl's BATMAN tutorial (http://www.miccai.org/edu/finalists/BATMAN_trimmed_tutorial.pdf)
  7. # The main difference between this script and the other one in this repository, is that this script assumes that your diffusion images were acquired with AP phase encoding
  8. # Thanks to John Plass and Bennet Fauber for useful comments
  9. # Adapted by Mia Anthony, March 2023
  10. # USAGE
  11. # $(basename $0) [Raw Diffusion] [RevPhaseImage] [AP bvec] [AP bval] [PA bvec] [PA bval] [Anatomical]
  12. # Arguments:
  13. # 1) The raw diffusion image;
  14. # 2) The image acquired with the reverse phase-encoding direction;
  15. # 3) The bvec file for the data acquired in the AP direction;
  16. # 4) The bval file for the data acquired in the AP direction;
  17. # 5) The bvec file for the data acquired in the PA direction;
  18. # 6) The bval file for the data acquired in the PA direction;
  19. # 7) The anatomical image"
  20. RAW_DWI=$1 # raw diffusion image /dwi/sub-0105_ses-01_dwi.nii.gz
  21. REV_PHASE=$2 # image acquired with the reverse phase-encoding direction
  22. # fmap/sub-0105_ses-01_dir-PA_dwi.nii.gz
  23. AP_BVEC=$3 # bvec file for AP dwi/bvec
  24. AP_BVAL=$4 # bval file for AP dwi/bval
  25. PA_BVEC=$5 # bvec file for PA not use
  26. PA_BVAL=$6 # bval file for PA not use
  27. ANAT=$7 # T1 anatomical
  28. BIDS_PATH=$8 # BIDS path
  29. ########################### STEP 0 ###################################
  30. # TO-DO before preprocessing (only once) #
  31. ######################################################################
  32. cd $BIDS_PATH
  33. # Create QC folder and assign the path to a variable named 'QC_DIR'
  34. QC_DIR=${BIDS_PATH}/QC/
  35. export QT_QPA_PLATFORM=offscreen
  36. export XDG_RUNTIME_DIR=${BIDS_PATH}/run
  37. export RUNLEVEL=3
  38. # export TMP_DIR=${BIDS_PATH}/tmp
  39. # Create the following folders in the QC dir to save QC images
  40. # /residual
  41. # /corrupt
  42. # /mask
  43. # /fod_voxels
  44. # /fod_overlay
  45. # /tissue_align
  46. # /gm_wm
  47. # mkdir -p $QC_DIR/residual $QC_DIR/corrupt $QC_DIR/mask $QC_DIR/fod_voxels $QC_DIR/fod_overlay $QC_DIR/tissue_align $QC_DIR/gm_wm
  48. ########################### STEP 1 ###################################
  49. # Convert data to .mif format and denoise #
  50. ######################################################################
  51. # Convert .nii to .mif format
  52. mrconvert $RAW_DWI raw_dwi.mif -fslgrad $AP_BVEC $AP_BVAL
  53. # Denoise dwi
  54. dwidenoise raw_dwi.mif dwi_den.mif -noise noise.mif
  55. # Calculate residual to check whether any region is disproportionately affected by noise.
  56. mrcalc raw_dwi.mif dwi_den.mif -subtract noise_residual.mif
  57. # ---------------------- QC steps ---------------------- #
  58. # Everything within the grey and white matter should be relatively uniform and blurry. If any clear anatomical landmarks are visible in the residual image, those parts of the brain have been corrupted by noise - increase the extent of the denoising filter from 5 (default) to a larger number, e.g. 7.
  59. # mkdir $QC_DIR/residual
  60. # mrview noise_residual.mif -capture.folder $QC_DIR/residual -capture.prefix $SBJ -capture.grab
  61. # dwidenoise raw_dwi.mif dwi_den.mif -extent 7 -noise noise_7extent.nii.gz
  62. # Determine whether Gibbs denoising is needed. Check diffusion data for ringing artifacts before and after to determine whether Gibbs denoising improved the data. If the data looks worse or the same, then skip Gibbs denoising.
  63. # mrdegibbs raw_dwi.mif dwi_den.mif
  64. mrdegibbs raw_dwi.mif dwi_den.mif -force
  65. # mrview raw_dwi.mif dwi_den.mif -capture.folder $QC_DIR/gibbs -capture.prefix $SBJ -capture.grab
  66. ##### Check for gibbs ringing before continuing to process
  67. # -------------------------------------------------------- #
  68. # Extract the b0 images in the AP direction. For the PA_BVEC and PA_BVAL files, they should be in the follwing format (assuming you extract only one volume):
  69. # PA_BVEC: 0 0 0
  70. # PA_BVAL: 0
  71. dwiextract dwi_den.mif - -bzero | mrmath - mean mean_b0_AP.mif -axis 3
  72. # Average b0s for AP direction to calculate mean intensity
  73. mrconvert $REV_PHASE PA.mif # If the PA map contains only 1 image, add the option "-coord 3 0"
  74. # Average b0s for PA direction
  75. # mrconvert PA.mif -fslgrad $PA_BVEC $PA_BVAL - | mrmath - mean mean_b0_PA.mif -axis 3
  76. mrconvert PA.mif - | mrmath - mean mean_b0_PA.mif -axis 3
  77. # Concatenate the b0 images from AP and PA directions to create a paired b0 image
  78. mrcat mean_b0_AP.mif mean_b0_PA.mif -axis 3 b0_pair.mif
  79. # Run dwipreproc (wrapper for eddy and topup). *CHANGE READOUT TIME based on TotalReadoutTime in the AP/PA JSON file*
  80. # CogTE = 0.0476
  81. readout_time=0.0476
  82. dwifslpreproc dwi_den.mif dwi_den_preproc.mif -nocleanup -pe_dir AP -rpe_pair -se_epi b0_pair.mif -readout_time $readout_time -align_seepi # -align_seepi ensures that the first volume in the series provided to topup is also the first volume in the series provided to eddy to make sure the volumes are aligned.
  83. # Bias field correction. Needs ANTs to be installed in order to use the "ants" option
  84. dwibiascorrect ants dwi_den_preproc.mif dwi_den_preproc_unbiased.mif -bias bias.mif
  85. # Create whole brain mask from bias corrected data - mask should be tight
  86. dwi2mask dwi_den_preproc_unbiased.mif mask.mif
  87. # create new bvecs and bvals files for the preprocessed data
  88. mrinfo dwi_den_preproc_unbiased.mif -export_grad_fsl dwi/bvecs dwi/bvals
  89. mrconvert dwi_den_preproc_unbiased.mif dwi/dwi.nii.gz
  90. # ########################### STEP 2 ###################################
  91. # # Basis function for each tissue type #
  92. # ######################################################################
  93. # # Create a basis function from the subject's DWI data to estimate response function for spherical deconvolution. The "dhollander" function is best used for multi-shell acquisitions; it will estimate different basis functions for each tissue type. For single-shell acquisition, use the "tournier" function instead
  94. # # -------------- Type of shell-acquisition -------------- #
  95. # ### SINGLE-SHELL (single b-value)
  96. # # dwi2response tournier dwi_den_preproc_unbiased.mif response_wm.txt response_gm.txt response_csf.txt -voxels voxels.mif
  97. # dwi2response tournier dwi_den_preproc_unbiased.mif response_wm.txt -voxels voxels.mif
  98. # # move files into QC folder
  99. # # mv response_wm.txt response_gm.txt response_csf.txt $QC_DIR/response/$SBJ
  100. # # Estimate fibre orientation distributions (FOD) from diffusion data using spherical deconvolution, using the basis functions estimated above
  101. # # dwi2fod csd dwi.mif response_wm.txt wmfod.mif
  102. # dwi2fod csd dwi_den_preproc_unbiased.mif response_wm.txt wmfod.mif
  103. # # -------------------------------------------------------- #
  104. # ### MULTI-SHELL (more than 1 b-value)
  105. # # dwi2response dhollander dwi_den_preproc_unbiased.mif response_wm.txt response_gm.txt response_csf.txt-voxels voxels.mif
  106. # # Estimate fibre orientation distributions from diffusion data using spherical deconvolution, using the basis functions estimated above
  107. # # dwi2fod msmt_csd dwi_den_preproc_unbiased.mif -mask mask.mif response_wm.txt wmfod.mif response_gm.txt gmfod.mif response_csf.txt csffod.mif
  108. # # -------------------------------------------------------- #
  109. # # Create an image of the FODs overlaid onto the estimated tissues (Blue=WM; Green=GM; Red=CSF)
  110. # # mrconvert -coord 3 0 wmfod.mif - | mrcat csffod.mif gmfod.mif - vf.mif
  111. # # Normalize the FODs to enable comparison between subjects
  112. # # mtnormalise wmfod.mif wmfod_norm.mif gmfod.mif gmfod_norm.mif csffod.mif csffod_norm.mif -mask mask.mif
  113. # mtnormalise wmfod.mif wmfod_norm.mif -mask mask.mif
  114. # ########################### STEP 3 ###################################
  115. # # Create a GM/WM boundary for seed analysis #
  116. # ######################################################################
  117. # # Convert the anatomical image to .mif format, and then extract all five tissue catagories (1=GM; 2=Subcortical GM; 3=WM; 4=CSF; 5=Pathological tissue)
  118. # mrconvert $ANAT anat.mif
  119. # # Segment the anatomical image into the tissue types
  120. # 5ttgen fsl anat.mif 5tt_nocoreg.mif
  121. # # Average the b0 images
  122. # dwiextract dwi_den_preproc_unbiased.mif - -bzero | mrmath - mean mean_b0_processed.mif -axis 3
  123. # # Convert b0 average and 5tt image to NIFTI for co-registration
  124. # mrconvert mean_b0_processed.mif mean_b0_processed.nii.gz
  125. # mrconvert 5tt_nocoreg.mif 5tt_nocoreg.nii.gz
  126. # # Extract the first volume (gray matter) of the 5tt dataset
  127. # fslroi 5tt_nocoreg.nii.gz 5tt_vol0.nii.gz 0 1
  128. # # Use fsl to create a transformation matrix for co-registration between the tissue map and the b0 images
  129. # flirt -in mean_b0_processed.nii.gz -ref 5tt_vol0.nii.gz -interp nearestneighbour -dof 6 -omat diff2struct_fsl.mat
  130. # # Convert back to format that mrtrix can read
  131. # transformconvert diff2struct_fsl.mat mean_b0_processed.nii.gz 5tt_nocoreg.nii.gz flirt_import diff2struct_mrtrix.txt
  132. # # Co-register the anatomical image to the diffusion image
  133. # mrtransform 5tt_nocoreg.mif -linear diff2struct_mrtrix.txt -inverse 5tt_coreg.mif
  134. # # Create a seed region along the GM/WM boundary
  135. # 5tt2gmwmi 5tt_coreg.mif gmwmSeed_coreg.mif
  136. # ########################### STEP 4 ###################################
  137. # # QC #
  138. # ######################################################################
  139. # # --------------------------- QC for Step 1 --------------------------
  140. # # View original diffusion data overlaid on top of the eddy-corrected data and colored in red
  141. # # mrview dwi_den_preproc.mif -overlay.load raw_dwi.mif
  142. # ### Corrupt slices
  143. # # dwi_post_eddy.eddy_outlier_map indicates whether a slice is an outlier (1) or not (0), because of too much motion, eddy currents, or something else.
  144. # cd dwifslpreproc-tmp-*
  145. # # Calculate the total number of slices by multiplying the number of slices for a single volume by the total number of volumes in the dataset.
  146. # totalSlices=`mrinfo dwi.mif | grep Dimensions | awk '{print $6 * $8}'`
  147. # # The total number of 1?s in the outlier map is then calculated, and the percentage of outlier slices is generated by dividing the number of outlier slices by the total number of slices. If this number is greater than 10 - i.e., if more than 10 percent of the slices are flagged as outliers - you should consider removing the subject from further analyses.
  148. # totalOutliers=`awk '{ for(i=1;i<=NF;i++)sum+=$i } END { print sum }' dwi_post_eddy.eddy_outlier_map`
  149. # echo "scale=5; ($totalOutliers / $totalSlices * 100)/1" | bc | tee ${SBJ}_percentageOutliers.txt
  150. # mv ${SBJ}_percentageOutliers.txt $QC_DIR/corrupt
  151. # cd ..
  152. # # --------------------------- QC for Step 2 --------------------------
  153. # ### Brain mask
  154. # # mrview mask.mif -capture.folder $QC_DIR/mask -capture.prefix $SBJ -capture.grab
  155. # ### Basis function
  156. # # View the response functions for each tissue type (response_*.txt files in the QC folder). The WM function should flatten out at higher b-values, while the other tissues should remain spherical
  157. # ### FOD estimation
  158. # # View the voxels used for FOD estimation (Blue=WM; Green=GM; Red=CSF)"
  159. # # mrview dwi_den_preproc_unbiased.mif -overlay.load voxels.mif -capture.folder $QC_DIR/fod_voxels -capture.prefix $SBJ -capture.grab
  160. # # Views the FODs overlaid on the tissue types (Blue=WM; Green=GM; Red=CSF)
  161. # # mrview vf.mif -odf.load_sh wmfod.mif -capture.folder $QC_DIR/fod_overlay -capture.prefix $SBJ -capture.grab
  162. # # --------------------------- QC for Step 3 --------------------------
  163. # # Check alignment of the 5 tissue types before and after alignment (new alignment in red, old alignment in blue)
  164. # # mrview dwi_den_preproc_unbiased.mif -overlay.load 5tt_nocoreg.mif -overlay.colourmap 2 -overlay.load 5tt_coreg.mif -overlay.colourmap 1 -capture.folder $QC_DIR/tissue_align -capture.prefix $SBJ -capture.grab
  165. # # Check the seed region (should match up along the GM/WM boundary)
  166. # # mrview dwi_den_preproc_unbiased.mif -overlay.load gmwmSeed_coreg.mif -capture.folder $QC_DIR/gm_wm -capture.prefix $SBJ -capture.grab

dti_rpe_minimal_process.sh at commit 722e68c, no license · at the source

Overview

  1. Department of Biostatistics, Yale University,New Haven, CT USA
  2. Department of Statistics and Operations Research, University of North Carolina at Chapel Hill,Chapel Hill, NC USA
Institutions: Yale University (United States); University of North Carolina at Chapel Hill (United States)
Journal: Nature communications, volume 17, issue 1, article 4850
Dates: received 8 July 2025; accepted 23 March 2026; published online 9 April 2026
Type: Research article · Language: English
License: CC BY
Identifiers: DOI 10.1038/s41467-026-71719-y · PMID 41957369 · PMCID PMC13222887 · OpenAlex W4414321762
Open access: gold, a free copy (OpenAlex)
Status: code verified
Categories: human (organism), developmental (subfield)
Methods: Connectivity, Spectral & time-frequency, Statistics, Machine learning, Preprocessing, fMRI & imaging
Keywords: Neural patterning, Predictive markers
MeSH: Behavior*, Brain*, Adolescent, Adult, Child, Cognition, Connectome, Female, Humans, Magnetic Resonance Imaging, Male, Young Adult (* major topic)
Topic: Functional Brain Connectivity Studies (Cognitive Neuroscience, Neuroscience), according to OpenAlex
Funding: NIH (R01EB034720, RF1AG081413, R01AG068191)
Citations: cited by 1 paper (Europe PMC); 69 references in the paper

Abstract

Gradients provide low-dimensional representations of macroscale brain organization, yet how structural-functional gradient coupling develops and relates to behavioral and molecular features remains unclear. Here, we studied structural-functional gradient coupling across multiple metrics and spatial scales using high-resolution structural and functional connectivity from 5343 children in the Adolescent Brain Cognitive Development study and 875 adults from the Human Connectome Project. We find that gradient coupling shows developmental refinement from childhood to adulthood and distinct sex-specific patterns. Gradient coupling metrics are significantly associated with cognitive and mental health measures and enable robust out-of-sample prediction. Heritability analyses reveal that gradient coupling is strongly influenced by genetic factors. Transcriptomic analyses further demonstrate that highly heritable coupling patterns are enriched for genes expressed in deep-layer excitatory neurons. Together, our findings establish structural-functional gradient coupling as a biologically meaningful feature of brain organization that bridges macroscale connectivity, cognition, behavior, and molecular architecture.

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 20 matches between paragraphs and lines of code.

sbci-brain/SBCI_Pipeline

License: none: the authors keep all their rights
State: the link answers, verified on 29 September 2026
Evidence: files inventoried
Commit: 722e68c41001af38a02d6a867e8153cc9720b69f, 22 May 2026
Languages: Shell (133), Python (51), C/C++ (4), MATLAB (2)
Size: 583 files, 190 scripts
Software Heritage: not archived
Found in: “Code availability”
Holds: README
Not found: license file, CITATION.cff, environment file, tests, continuous integration, documentation
Tools: NumPy (48 files), MRtrix3 (35 files), FreeSurfer (30 files), FSL (27 files), SciPy (25 files), ANTs (15 files), NiBabel (14 files), DIPY (2 files), Dcm2Bids (1 file), dcm2niix (1 file), scikit-learn (1 file)
Availability: 1 check, the latest on 29 September 2026: the link answers
  • 29 September 2026: the link answers
191 files

mica-mni/brainspace

License: BSD-3-Clause
State: the link answers, verified on 26 September 2026
Evidence: files inventoried
Commit: 8730de88ae32c4f88eeaf16ef2a6e53c5c32dc34, 5 May 2026
Languages: Python (77), MATLAB (40), Jupyter (4)
Size: 373 files, 121 scripts
Software Heritage: not archived
Found in: “Code availability”
Holds: README, license file, CITATION.cff, environment (Dockerfile, requirements.txt, setup.cfg, setup.py, docs/requirements.txt), tests, continuous integration, documentation, 4 notebooks
Tools: NumPy (52 files), BrainSpace (32 files), SciPy (18 files), scikit-learn (13 files), Matplotlib (9 files), NiBabel (3 files), Nilearn (3 files), GIfTI library for MATLAB (1 file), Parallel Computing Toolbox (1 file), Statistics and Machine Learning Toolbox (1 file)
Availability: 1 check, the latest on 26 September 2026: the link answers
  • 26 September 2026: the link answers
123 files

Zhao-team/SF-Gradient-Coupling

License: MIT
State: the link answers, verified on 29 September 2026
Evidence: files inventoried
Commit: ff163a9e9aa01dc7600ac3cc3f0c4a0d0676574f, 15 December 2025
Languages: Python (6), R (1)
Size: 15 files, 7 scripts
Software Heritage: not archived
Found in: “Code availability”
Holds: README, license file
Not found: CITATION.cff, environment file, tests, continuous integration, documentation
Tools: NumPy (6 files), pandas (6 files), SciPy (5 files), scikit-learn (4 files), Matplotlib (3 files), BrainSpace (2 files), h5py (2 files), PyTorch (2 files), abagen (1 file), anndata (1 file), ggplot2 (1 file), limma (1 file), patchwork (1 file), Scanpy (1 file), tidyverse (1 file)
Availability: 1 check, the latest on 29 September 2026: the link answers
  • 29 September 2026: the link answers
9 files

Zenodo 18912522

License: CC-BY-4.0
State: the link answers, verified on 29 September 2026
Evidence: files inventoried
Size: 1 file
Software Heritage: not checked
Found in: “Code availability”
Not found: README, license file, CITATION.cff, environment file, tests, continuous integration, documentation
Availability: 1 check, the latest on 29 September 2026: the link answers (HTTP 200)
  • 29 September 2026: the link answers (HTTP 200)
At the source:

Code availability

The data preprocessing software FreeSurfer v6.0 is available at https://surfer.nmr.mgh.harvard.edu/. The Surface-Based Connectivity Integration pipeline can be accessed at https://github.com/sbci-brain/SBCI_Pipeline, and the BrainSpace toolbox is available at https://github.com/MICA-MNI/BrainSpace/tree/master. Python (2.7 and 3.12) and R v4.4 were used for data processing and analysis. Code used in this study is publicly available at https://github.com/Zhao-team/SF-Gradient-Coupling.git and stored at 10.5281/zenodo.1891252268.

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:

  • 4 repositories of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
  • 318 scripts, each with its path and the digest of its content;
  • 20 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

Datasets cited

Data availability

Neuroimaging and behavioral data from the ABCD Study can be obtained via the NIH Data Archive (https://nda.nih.gov/abcd) with approval from the ABCD consortium. Neuroimaging data and most behavioral measures from the HCP-YA are publicly available at https://db.humanconnectome.org; access to restricted data is subject to approval. Data from the HCP-D study are available through the NIH Data Archive (https://nda.nih.gov) and require approval for access. The raw data are protected and are not available due to data privacy laws and the terms of the original ethical approvals. Source data are provided with this paper. The data are supplied in multiple formats, including Python pickle (.pkl) files, which can be read using standard Python packages such as pickle or pandas. Source data are provided with this paper.

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, 29 September 2026: the first record

Recorded: type, language, journal, volume, issue, pages, dates, 7 authors, 2 keywords, 12 MeSH terms, 1 funder, 64 references.

Cite

This paper

Gao, S., Gu, Z., Ding, S., Wang, G., Zhang, Z., Zhao, H., & Zhao, Y. (2026). Brain functional-structural gradient coupling reflects development, behavior and genetic influences. Nature communications, 17(1), 4850. https://doi.org/10.1038/s41467-026-71719-y

BibTeX

@article{gao2026brain,
author = {Gao, Simiao and Gu, Zhiling and Ding, Shengxian and Wang, Gefei and Zhang, Zhengwu and Zhao, Hongyu and Zhao, Yize},
title = {{Brain functional-structural gradient coupling reflects development, behavior and genetic influences}},
journal = {Nature communications},
year = {2026},
month = apr,
volume = {17},
number = {1},
pages = {4850},
publisher = {Nature Publishing Group},
issn = {2041-1723},
doi = {10.1038/s41467-026-71719-y},
url = {https://doi.org/10.1038/s41467-026-71719-y},
pmid = {41957369},
pmcid = {PMC13222887}
}

RIS

TY - JOUR
AU - Gao, Simiao
AU - Gu, Zhiling
AU - Ding, Shengxian
AU - Wang, Gefei
AU - Zhang, Zhengwu
AU - Zhao, Hongyu
AU - Zhao, Yize
TI - Brain functional-structural gradient coupling reflects development, behavior and genetic influences
T2 - Nature communications
J2 - Nat Commun
PY - 2026
DA - 2026/04/09
VL - 17
IS - 1
SP - 4850
SN - 2041-1723
PB - Nature Publishing Group
DO - 10.1038/s41467-026-71719-y
UR - https://doi.org/10.1038/s41467-026-71719-y
LA - en
ER -

CSL-JSON

{
"id": "10.1038/s41467-026-71719-y",
"type": "article-journal",
"title": "Brain functional-structural gradient coupling reflects development, behavior and genetic influences",
"container-title": "Nature communications",
"author": [
{
"family": "Gao",
"given": "Simiao"
},
{
"family": "Gu",
"given": "Zhiling"
},
{
"family": "Ding",
"given": "Shengxian"
},
{
"family": "Wang",
"given": "Gefei"
},
{
"family": "Zhang",
"given": "Zhengwu"
},
{
"family": "Zhao",
"given": "Hongyu"
},
{
"family": "Zhao",
"given": "Yize"
}
],
"container-title-short": "Nat Commun",
"volume": "17",
"issue": "1",
"page": "4850",
"DOI": "10.1038/s41467-026-71719-y",
"PMID": "41957369",
"PMCID": "PMC13222887",
"ISSN": "2041-1723",
"publisher": "Nature Publishing Group",
"URL": "https://doi.org/10.1038/s41467-026-71719-y",
"language": "en",
"issued": {
"date-parts": [
[
2026,
4,
9
]
]
}
}

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.7554/elife.103097 [code]
Canonical neurodevelopmental trajectories of structural and functional manifolds.
Journal: eLife
In common: BrainSpace, Nilearn, h5py, 8 other tools, developmental, 14 references
[2] doi:10.1038/s41467-026-71270-w [code]
Spatiotemporal dynamics of the human cortical functional hierarchy across the lifespan.
Journal: Nature communications
In common: BrainSpace, GIfTI library for MATLAB, Parallel Computing Toolbox, 12 other tools, developmental, 8 references
[3] doi:10.1038/s41467-026-76011-7 [code]
Human cortex organizes dynamic co-fluctuations along the sensorimotor-association axis.
Journal: Nature communications
In common: BrainSpace, GIfTI library for MATLAB, ANTs, 10 other tools, 10 references
[4] doi:10.1038/s42003-025-09444-3 [code]
Decoupling of neurophysiological activity from structure mirrors global microarchitectural and neuromodulatory trends.
Journal: Communications biology
In common: abagen, Nilearn, NiBabel, 6 other tools, 14 references
[5] doi:10.1038/s41467-026-73668-y [code]
Convergent and divergent brain-cognition development in early adolescence.
Journal: Nature communications
In common: MRtrix3, GIfTI library for MATLAB, ANTs, 13 other tools, 6 references
[6] doi:10.1038/s41467-026-76812-w [code]
Assessing molecular, cellular and transcriptomic bases of laminar perfusion and cytoarchitecture coupling in the human cortex.
Journal: Nature communications
In common: BrainSpace, limma, ANTs, 12 other tools, 6 references
[7] doi:10.1038/s41586-026-10631-3 [code]
A prognostic human brain network for diffuse midline glioma.
Journal: Nature
In common: BrainSpace, GIfTI library for MATLAB, ANTs, 12 other tools, 5 references
[8] doi:10.1038/s41467-026-71568-9 [code]
Convergent and selective representations of pain, appetitive processes, aversive processes, and cognitive control in the insula.
Journal: Nature communications
In common: abagen, MRtrix3, ANTs, 11 other tools, 5 references
[9] doi:10.1016/j.isci.2026.116671 [code]
A high-resolution functional network-organized atlas of human superficial white matter from ultra-high-field diffusion MRI.
Journal: iScience
In common: DIPY, MRtrix3, ANTs, 11 other tools, 5 references
[10] doi:10.1038/s41398-026-04025-2 [code]
Brain energetic landscapes shape state dysregulation in major depressive disorder: a morphological network controllability perspective.
Journal: Translational psychiatry
In common: BrainSpace, GIfTI library for MATLAB, ANTs, 13 other tools, 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.

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.