Genetic architectures of brain-related traits are shaped by strong selective constraints.
The 10 matches · 3 of them tie a paragraph to a whole file without given lines: weak matches, whose lines are not tinted
- [1] § Materials and Methods › Thresholding and Downsampling Quantitative Traits. ↔ scripts/gwas_deflation/05_thresholding_downsampling.R, lines 1 to 12 · score 0.89 · diastolic blood pressure, systolic blood pressure, HbA1c, heel bone mineral, BMI, LDL · Report a mismatch
- [2] § Materials and Methods › Thresholding and Downsampling Quantitative Traits. ↔ scripts/gwas_deflation/01_field_names.ipynb, lines 37 to 78 · score 0.85 · heel bone mineral, diastolic blood pressure, systolic blood pressure, HbA1c, medically relevant, fields · Report a mismatch
- [3] § Results › The Distinct Genetic Architectures of Brain-Related Traits. ↔ scripts/plotting/sldsc_plots.R, lines 1 to 48 · score 0.70 · Bonferroni correction, selected traits, DFP, FIS, IMR, AFS · Report a mismatch
- [4] § Materials and Methods › Distinct Subsets of (Approximately Uncorrelated) Brain-Related and Non-Brain-Related Traits. ↔ scripts/gwas_analyses/get_independent_traits.R, lines 1 to 69 · score 0.68 · Neale Lab, Genetic correlation, CNS traits, meta · Report a mismatch
- [5] § Results › Brain-Related Traits Retain Distinct Architectures After Accounting for Differences in Study Design. ↔ scripts/gwas_deflation/05_thresholding_downsampling.R, lines 22 to 67 · score 0.67 · reach equivalent power, original phenotypic, normal transformation, Deflation, downsampling, inverse · Report a mismatch
- [6] § Materials and Methods › Using Stratified LD Score Regression (S-LDSC) to Define Brain-Related Traits. ↔ scripts/gwas_analyses/07_sldsc_metaanalysis_ACAT.R, the whole file · a weak match · score 0.67 · Cauchy variables, Cauchy distribution, ACAT, meta, transformed, cell · Report a mismatch
- [7] § Materials and Methods › Thresholding and Downsampling Quantitative Traits. ↔ scripts/gwas_deflation/06_gwas_deflated.sh, the whole file · a weak match · score 0.64 · binary traits, fallback, firth, covar, glm, hide · Report a mismatch
- [8] § Materials and Methods › Identifying Approximately Independent Hits Using GCTA-COJO. ↔ scripts/gwas_analyses/02_cojo.sh, the whole file · a weak match · score 0.62 · cojo slct, GCTA COJO, GWAS · Report a mismatch
- [9] § Materials and Methods › Using Stratified LD Score Regression (S-LDSC) to Define Brain-Related Traits. ↔ ldsc.py, lines 547 to 617 · score 0.59 · munge_sumstats.py, LD scores, regression coefficient · Report a mismatch
- [10] § Materials and Methods › Thresholding and Downsampling Quantitative Traits. ↔ scripts/gwas_deflation/05_thresholding_downsampling.R, lines 22 to 67 · score 0.56 · normal transformed, nest, downsampling, inverse, rank, binary · Report a mismatch
Paper
Loaded from Europe PMC by your browser; OSCR does not store it: doi.org · Europe PMC
The paper is loaded when this pane is shown.
The authors' code
R · 67 lines · 2.6 KB · no license · 3 matches
- # Run this script on UKB-RAP
- set.seed(123)
- library(data.table)
- input<- fread('/mnt/project/Deflation/phenotypes_processed.tsv')
- data<- as.data.frame(input)
- trait_list<-setdiff(colnames(data), c("FID", "IID")) #List of traits to threhold on
- trait_names<- c("BMI","Creatinine in urine","HbA1c","LDL","Diastolic blood pressure","Systolic blood pressure","Heel bone mineral density T-score")
- k <- c(20,10,5,2,1) #Population prevalences
- calculate_N_prime <- function(K, omega, M) {
- # M: number of cases + number of controls included in case-control GWAS (M=N if omega=K)
- # K: population prevalence
- # omega: sample prevalence
- size <- M * omega * (1-omega) * (dnorm(qnorm(1-K)) **2) / ((K*(1-K))**2)
- return(round(size))
- }
- # x -- quantitative phenotypes
- # y -- binary phenotypes
- for (i in 1:length(trait_list)) {
- trait<-trait_list[i]
- trait_name<-trait_names[i]
- system(paste("dx mkdir -p /Deflation/",trait,"/plots", sep=""))
- filtered_data<-data[!is.na(data[[trait]]), c("FID", "IID", trait)] #Filter for non-NA entries
- my_data<-filtered_data[order(filtered_data[[trait]]),]
- N<-nrow(my_data) #Sample size
- setnames(my_data, old = trait, new = "x") #x: original phenotype, sorted in increasing order
- my_data$x_irnt<-qnorm((rank(my_data$x,na.last="keep")-0.5)/sum(!is.na(my_data$x))) #Inverse rank normal transformed phenotype
- previous_indices <- NULL
- for (prev in k){
- case_number <- round(N*prev/100,digits=0) # PLINK format: controls=1, controls=2
- # Left tail thresholding
- colname_left <- paste("y_left_k", prev, sep="")
- my_data[[colname_left]] <- rep(1,N)
- my_data[[colname_left]][1:case_number] <- rep(2, case_number)
- # Right tail thresholding
- colname_right <- paste("y_right_k", prev, sep="")
- my_data[[colname_right]] <- rep(1,N)
- my_data[[colname_right]][(N-(case_number-1)):N] <- rep(2, case_number)
- # Nested downsampling to reach equivalent power
- colname <- paste("x_irnt_k", prev, sep="")
- my_data[[colname]] <- rep(NaN,N)
- size<-calculate_N_prime(prev/100,prev/100,N) #Assuming omega=K and M=N
- if (is.null(previous_indices)) {
- indices <- sample(1:N, size)
- } else {
- indices <- sample(previous_indices, size)
- }
- new_x<- my_data$x[indices] #new sample
- my_data[[colname]][indices] <- qnorm((rank(new_x,na.last="keep")-0.5)/sum(!is.na(new_x))) #IRNT the new sample
- previous_indices<-indices
- }
- write.table(my_data, file = paste(trait,".txt", sep=""), sep = "\t", row.names = FALSE, quote = FALSE)
- system(paste("dx upload ", trait, ".txt --dest /Deflation/",trait, "/", sep=""))
- }
05_thresholding_downsampling.R at commit cff8309, no license · at the source
Overview
- Department of Biology, Stanford University, Stanford, CA 94305
- Department of Genetics, Stanford University, Stanford, CA 94305
- Section of Genetic Medicine, University of Chicago, Chicago, IL 60637
- Department of Human Genetics, University of Chicago, Chicago, IL 60637
- Institute for Human Genetics, University of California, San Francisco, CA 94158
- Department of Epidemiology and Biostatistics, University of California, San Francisco, CA 94158
- Department of Biological Sciences, Columbia University, New York City, NY 10027
- Program for Mathematical Genomics, Columbia University, New York City, NY 10032
Abstract
Genome-wide association studies (GWAS) have identified hundreds of significant loci for psychiatric disorders, yet the strength of these associations remains modest compared to other human complex traits with similar numbers of hits. Whether this pattern reflects statistical artifacts or real biological differences—and, if the latter, what underlies it—remains unclear. In addition to psychiatric disorders, we find that other traits with functional enrichment in the central nervous system (CNS), whether binary or quantitative, also share similar genetic architectures, characterized by GWAS hits of limited statistical significance and generally higher allele frequencies. In comparing the architecture of binary and quantitative traits, we adjust for statistical power in their respective studies. After this adjustment, we fit an evolutionary model of architecture and show that CNS-enriched traits have large mutational target sizes, with contributing variants and genes experiencing stronger selection than those for other traits. Our findings reveal heterogeneity among complex traits and provide insights into traits that more effectively capture fitness-relevant processes. More broadly, our results suggest that the genetic architectures of complex traits are shaped by the tissues through which these traits are mediated.
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 10 matches between paragraphs and lines of code.
bulik/ldsc
2fdeeb3b44379408794154993dbd6101b8946b7e, 16 January 2026Availability: 2 checks, the latest on 26 October 2026: the link answers
- 26 October 2026: the link answers
- 26 September 2026: the link answers
27 files
- ContinuousAnnotations/
quantile_M.pl , Perl, 241 lines - ContinuousAnnotations/
quantile_h2g.r , R, 76 lines - ldsc.py, Python, 660 lines, 1 match
- ldscore/
__init__.py , Python, 1 line - ldscore/
irwls.py , Python, 196 lines - ldscore/
jackknife.py , Python, 514 lines - ldscore/
ldscore.py , Python, 415 lines - ldscore/
parse.py , Python, 292 lines - ldscore/
regressions.py , Python, 743 lines - ldscore/
sumstats.py , Python, 581 lines - make_annot.py, Python, 56 lines
- munge_sumstats.py, Python, 745 lines
- setup.py, Python, 20 lines
- test/
parse_test/ , MATLAB, 1 linetest.l2.M - test/
parse_test/ , MATLAB, 1 linetest1.l2.M - test/
parse_test/ , MATLAB, 1 linetest2.l2.M - test/
parse_test/ , MATLAB, 1 linetest_bad.l2.M - test/
simulate.py , Python, 81 lines - test/
test_irwls.py , Python, 69 lines - test/
test_jackknife.py , Python, 267 lines - test/
test_ldscore.py , Python, 111 lines - test/
test_munge_sumstats.py , Python, 358 lines - test/
test_parse.py , Python, 129 lines - test/
test_regressions.py , Python, 342 lines - test/
test_sumstats.py , Python, 487 lines - LICENSE, License, 675 lines
- README.md, Text, 122 lines
huishengz/cns-selection
cff83096f0ca242980b30769973e0b3b222bd7eb, 20 March 2026Availability: 1 check, the latest on 26 September 2026: the link answers
- 26 September 2026: the link answers
45 files
- scripts/
gwas_analyses/ , R, 43 lines01_prep_gcta.R - scripts/
gwas_analyses/ , Shell, 47 lines01_prep_submit_cojo.sh - scripts/
gwas_analyses/ , Shell, 65 lines, 1 match02_cojo.sh - scripts/
gwas_analyses/ , R, 57 lines03_collect_cojo.R - scripts/
gwas_analyses/ , Shell, 46 lines03_collect_hits.sh - scripts/
gwas_analyses/ , Shell, 39 lines04_munge_sumstats.sh - scripts/
gwas_analyses/ , Shell, 24 lines05_sldsc.sh - scripts/
gwas_analyses/ , Shell, 15 lines05_submit_sldsc.sh - scripts/
gwas_analyses/ , Shell, 38 lines06_collect_sldsc_z.sh - scripts/
gwas_analyses/ , R, 44 lines, 1 match07_sldsc_metaanalysis_AC AT.R - scripts/
gwas_analyses/ , R, 24 linesget_Nprime.R - scripts/
gwas_analyses/ , R, 169 lines, 1 matchget_independent_traits.R - scripts/
gwas_deflation/ , Jupyter, 97 lines, 1 match01_field_names.ipynb - scripts/
gwas_deflation/ , Jupyter, 82 lines02_extract_phenotypes.ip ynb - scripts/
gwas_deflation/ , Jupyter, 110 lines03_process_phenotypes.ip ynb - scripts/
gwas_deflation/ , Shell, 24 lines04_gwas_original.sh - scripts/
gwas_deflation/ , R, 67 lines, 3 matches05_thresholding_downsamp ling.R - scripts/
gwas_deflation/ , Shell, 53 lines, 1 match06_gwas_deflated.sh - scripts/
inference_and_simulation , Shell, 61 liness/ 01_inference.sh - scripts/
inference_and_simulation , Shell, 22 liness/ 02_job.sh - scripts/
inference_and_simulation , Shell, 24 liness/ 02_submit_jobs.sh - scripts/
inference_and_simulation , Shell, 21 liness/ 11_simulations.sh - scripts/
inference_and_simulation , MATLAB, 87 liness/ collate_boot.m - scripts/
inference_and_simulation , MATLAB, 132 liness/ collate_jack.m - scripts/
inference_and_simulation , MATLAB, 191 liness/ do_inference.m - scripts/
inference_and_simulation , MATLAB, 175 liness/ infer_this_data.m - scripts/
inference_and_simulation , MATLAB, 117 liness/ make_datafile_from_cojo. m - scripts/
inference_and_simulation , MATLAB, 53 liness/ preprocess.m - scripts/
inference_and_simulation , MATLAB, 29 liness/ run_simulation_changeL.m - scripts/
inference_and_simulation , MATLAB, 40 liness/ run_simulation_different fs.m - scripts/
inference_and_simulation , MATLAB, 45 liness/ run_simulation_gwas_pred iction.m - scripts/
inference_and_simulation , MATLAB, 99 liness/ simulate_a_no_pleiotropy _trait.m - scripts/
inference_and_simulation , MATLAB, 97 liness/ simulate_a_trait.m - scripts/
plotting/ , R, 184 linesburden_bootstrap.R - scripts/
plotting/ , R, 316 linescdf_plots.R - scripts/
plotting/ , R, 253 linesgenetic_architecture_2D_ plots.R - scripts/
plotting/ , R, 139 linesgwas_deflation_plots.R - scripts/
plotting/ , R, 157 linesgwas_prediction.R - scripts/
plotting/ , R, 289 linesinference_plot.R - scripts/
plotting/ , R, 95 linesmanhattan_plots.R - scripts/
plotting/ , R, 454 linessimulation_plots.R - scripts/
plotting/ , R, 145 lines, 1 matchsldsc_plots.R - supplement/
trait_specific_fs.R , R, 122 lines - supplement/
wgs_h2_comparison.R , R, 108 lines - README.md, Text, 1 line
The paper's code and data availability statement is in the Data section.
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;
- 69 scripts, each with its path and the digest of its content;
- 10 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
- ukbiobank.ac.uk/
use-our-data/ , at UK Biobank; found in “Data, Materials, and Software Availability”apply-for-access - zenodo:19154957, at Zenodo; found in “Data, Materials, and Software Availability”
Data, Materials, and Software Availability
Codes used for this article are available 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, 5 authors, 5 keywords, 10 MeSH terms, 6 funders, 93 references.
Cite
This paper
Zhu, H., Simons, Y. B., Spence, J. P., Sella, G., & Pritchard, J. K. (2026). Genetic architectures of brain-related traits are shaped by strong selective constraints. Proceedings of the National Academy of Sciences of the United States of America, 123(36), e2609814123. https://
BibTeX
@article{zhu2026genetic,
author = {Zhu, Huisheng and Simons, Yuval B and Spence, Jeffrey P and Sella, Guy and Pritchard, Jonathan K},
title = {{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 America},
year = {2026},
month = sep,
volume = {123},
number = {36},
pages = {e2609814123},
publisher = {National Academy of Sciences},
issn = {0027-8424},
doi = {10.1073/
url = {https://
pmid = {42679042},
pmcid = {PMC13552918}
}
RIS
TY - JOUR
AU - Zhu, Huisheng
AU - Simons, Yuval B
AU - Spence, Jeffrey P
AU - Sella, Guy
AU - Pritchard, Jonathan K
TI - Genetic architectures of brain-related traits are shaped by strong selective constraints
T2 - Proceedings of the National Academy of Sciences of the United States of America
J2 - Proc Natl Acad Sci U S A
PY - 2026
DA - 2026/
VL - 123
IS - 36
SP - e2609814123
SN - 0027-8424
PB - National Academy of Sciences
DO - 10.1073/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1073/
"type": "article-journal",
"title": "Genetic architectures of brain-related traits are shaped by strong selective constraints",
"container-title": "Proceedings of the National Academy of Sciences of the United States of America",
"author": [
{
"family": "Zhu",
"given": "Huisheng"
},
{
"family": "Simons",
"given": "Yuval B"
},
{
"family": "Spence",
"given": "Jeffrey P"
},
{
"family": "Sella",
"given": "Guy"
},
{
"family": "Pritchard",
"given": "Jonathan K"
}
],
"container-title-short":
"volume": "123",
"issue": "36",
"page": "e2609814123",
"DOI": "10.1073/
"PMID": "42679042",
"PMCID": "PMC13552918",
"ISSN": "0027-8424",
"publisher": "National Academy of Sciences",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
9,
1
]
]
}
}
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/s43587-025-00962-7 [code]
- Genome-wide analysis of brain age identifies 59 associated loci and unveils relationships with mental and physical health.Journal: Nature agingIn common: patchwork, cowplot, data.table, 6 other tools, genetics / omics, cellular / molecular, 9 references
- [2] doi:10.1016/j.cell.2024.12.002 [code]
- Trans-ancestry genome-wide study of depression identifies 697 associations implicating cell types and pharmacotherapies.Journal: CellIn common: patchwork, cowplot, data.table, 3 other tools, genetics / omics, cellular / molecular, 10 references
- [3] doi:10.1038/s41562-025-02145-1 [code]
- Genome-wide association meta-analysis of age at onset of walking in over 70,000 infants of European ancestry.Journal: Nature human behaviourIn common: data.table, ggplot2, tidyverse, 4 other tools, genetics / omics, 11 references
- [4] doi:10.1038/s41562-024-02051-y [code]
- The shared genetic architecture and evolution of human language and musical rhythm.Journal: Nature human behaviourIn common: BEDTools, cowplot, data.table, 5 other tools, genetics / omics, 8 references
- [5] doi:10.1038/s41588-025-02189-z [code]
- Genome-wide analyses identify 30 loci associated with obsessive-compulsive disorder.Journal: Nature geneticsIn common: BEDTools, data.table, ggplot2, 4 other tools, genetics / omics, cellular / molecular, 9 references
- [6] doi:10.1038/s41467-026-76676-0 [code]
- Determinants of functional burden pleiotropy and gene dosage responses across human traits.Journal: Nature communicationsIn common: patchwork, data.table, ggplot2, 5 other tools, genetics / omics, cellular / molecular, 7 references
- [7] doi:10.1016/j.biopsych.2023.06.022 [code]
- The Genetic Architecture of Amygdala Nuclei.Journal: Biological psychiatryIn common: BEDTools, patchwork, data.table, 7 other tools, genetics / omics, 7 references
- [8] doi:10.1186/s13073-023-01212-4 [code]
- Shared genetic architecture between irritable bowel syndrome and psychiatric disorders reveals molecular pathways of the gut-brain axis.Journal: Genome medicineIn common: data.table, Statistics and Machine Learning Toolbox, pandas, 3 other tools, genetics / omics, cellular / molecular, 10 references
- [9] doi:10.1038/s41586-025-09820-3 [code]
- Mapping the genetic landscape across 14 psychiatric disorders.Journal: NatureIn common: data.table, ggplot2, tidyverse, 4 other tools, genetics / omics, 9 references
- [10] doi:10.1038/s43856-026-01510-z [code]
- Mapping genetic convergence across brain structure, mental health, and cardiometabolic disease.Journal: Communications medicineIn common: BEDTools, data.table, ggplot2, 6 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: 2 repositories of the authors' code, each at its verified commit and with its license, 69 scripts, and 10 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:2117d7a980933bfc…
Add the badge to the code's README
If you are an author or a maintainer of this paper's code, you can add this badge. It says what OSCR knows of the paper (its code, its validated tracing map) and links the code to this page. Paste one of these snippets under the README's title. Keep the two marks around it, so that the badge can be replaced in place and is never added twice. OSCR changes nothing in your repository.
Markdown
<!-- oscr-badge:start -->
[](https://openscicode.org/badge/doi/10.1073/pnas.2609814123)
<!-- oscr-badge:end -->reStructuredText
.. oscr-badge:start
.. image:: https://openscicode.org/badge/doi/10.1073/pnas.2609814123.svg
:target: https://openscicode.org/badge/doi/10.1073/pnas.2609814123
:alt: OSCR
.. oscr-badge:endHTML
<!-- oscr-badge:start -->
<a href="https://openscicode.org/badge/doi/10.1073/pnas.2609814123"><img src="https://openscicode.org/badge/doi/10.1073/pnas.2609814123.svg" alt="OSCR"></a>
<!-- oscr-badge:end -->On GitHub, open the README, choose Edit (the pencil), paste the snippet under the title, then "Commit changes".
If you choose to, you can have the badge kept up to date for you, with oscr badge --write from the command line or with the GitHub Action in your repository's workflows. Every other form and option (a card, AsciiDoc, LaTeX, shields.io) is on the badge 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: code errors and mismatches between the code and the paper about this paper and its code, from readers and its authors, beside the paper and its tracing map.
Report a mismatch between the code and the paper · Report a code error (signed in)
Reproductions: reports from readers who ran the authors' code and did not get the paper's result, what they ran, with which environment, commit and data.
Report a reproduction failure (signed in)
Versions of its code: the releases of its code tied to a version of this paper (the preprint, the accepted manuscript, the version of record…), each with the tracing map versioned with it.
Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.
