Identification of moderate effect size genes in autism spectrum disorder through a novel gene pairing approach.
The 7 matches · 6 of them tie a paragraph to a whole file, not to given lines: weak matches, whose lines are not tinted
- [1] § Methods › Incorporation and assessment of non-neuro cohorts ↔ 2_MES_non_neuro/example_commands_to_annotate.sh, the whole file · a weak match · score 0.78 · SNPeff, dbNSFP, missense variants, SIFT score, Ensembl, neuro
- [2] § Methods › Simulation of discovery pipeline ↔ 2_MES_non_neuro/simulate_MES_genes.R, lines 66–143 · score 0.74 · variant assignments, twins, orphans, rbinom, MES gene, simulated
- [3] § Methods › Gene expression ↔ 4_Expression/query_gene_expression.R, the whole file · a weak match · score 0.72 · tissue consensus, nTPM, genes expressed, GTEx, command, MES
- [4] § Methods › Incorporation and assessment of non-neuro cohorts ↔ 3_MES_synonymous/non_neuro_example_commands_to_annotate.sh, the whole file · a weak match · score 0.71 · SNPeff, dbNSFP, SIFT score, Ensembl, deleterious variants, neuro
- [5] § Results › Enrichment of MES risk gene expression in the brain ↔ 4_Expression/query_gene_expression.R, the whole file · a weak match · score 0.58 · nTPM, adrenal, cerebral, consensus, MES gene, tissues
- [6] § Methods › Family-based search for MES risk genes in SPARK ↔ 2_MES_non_neuro/example_commands_to_annotate.sh, the whole file · a weak match · score 0.58 · missense variant, SIFT score, Ensembl, PTV, neuro, deleterious
- [7] § Results › Controlling false positives in candidate gene pair analysis ↔ 1_MES_screen/4_table.R, the whole file · a weak match · score 0.50 · SPARK parents, ultra rare deleterious, rank, sum, probability, inheritance
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 · 39 lines · 1.6 KB · GPL-3.0 · 2 matches
- # Moderate effect size genes project
- # Author: Madison Caballero
- # Desription: A set of steps to annotate and extract variants in the non-neuro cohorts.
- # These steps were used for 1kGP, BioMe, and All of Us
- # For speed, I just extract one sample. This keeps all loci even if 0/0
- bcftools view --threads 4 -s [random_sample_id] input.vcf.gz > singlesamp.vcf
- # SNPeff annotate for function
- java -Xmx8g -jar snpEff.jar -v GRCh38.86 -canon singlesamp.vcf > singlesamp.ann.vcf
- # SNPsift to get dbNSFP information
- java -jar SnpSift.jar dbnsfp -v -db dbNSFP4.1a.txt.gz \
- -f gnomAD_exomes_AF,SIFT_score,SIFT4G_score,Polyphen2_HDIV_score,Polyphen2_HVAR_score,genename,Ensembl_geneid \
- singlesamp.ann.vcf > singlesamp.ann.dbNSFP.vcf
- # Get gene IDs that are SNPeff compatible
- awk '{print "missense_variant|MODERATE|" $1}' gene_names_Snpeffcompat.txt > grab.txt
- awk '{print "HIGH|" $1}' gene_names_Snpeffcompat.txt >> grab.txt
- # Grab variants that are missense or highly damaging (PTVs) for those genes
- # There is a later filtering step in case grep grabs partial matches
- grep -f grab.txt singlesamp.ann.dbNSFP.vcf > singlesamp.ann.dbNSFP.functional.noheader.vcf
- perl filter.pl singlesamp.ann.dbNSFP.functional.noheader.vcf > singlesamp.ann.dbNSFP.deleterious.noheader.vcf
- # header check
- grep -P "^" singlesamp1.ann.dbNSFP.vcf > header.txt
- cat header.txt singlesamp.ann.dbNSFP.deleterious.noheader.vcf > singlesamp.ann.dbNSFP.deleterious.vcf
- # Grab regions
- bcftools view \
- -T singlesamp.ann.dbNSFP.deleterious.vcf \
- -O v --threads 8 \
- -o filtered.input.vcf \
- input.vcf.gz
- # Get counts per gene
- perl grab_counts.pl
example_commands_to_annotate.sh at commit 17c5b6a, under GPL-3.0 · at the source
Overview
- Seaver Autism Center for Research and Treatment, Icahn School of Medicine at Mount Sinai, New York, NY USA
- The Mindich Child Health and Development Institute, Icahn School of Medicine at Mount Sinai, New York, NY USA
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY USA
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA USA
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY USA
- Analytic and Translational Genetics Unit, Department of Medicine, Massachusetts General Hospital, Boston, MA USA
- Friedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, NY USA
- Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, MA USA
- Department of Neuroscience, Icahn School of Medicine at Mount Sinai, New York, NY USA
- Department of Artificial Intelligence and Human Health, Icahn School of Medicine at Mount Sinai, New York, NY USA
- Department of Molecular Medicine and Surgery, Karolinska Institutet, Stockholm, Sweden
- Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden
Abstract
Autism Spectrum Disorder (ASD) arises from complex genetic and environmental factors, with inherited genetic variation playing a substantial role. This study introduces a novel approach to uncover moderate effect size (MES) risk genes in ASD, which individually do not meet the ASD liability threshold but contribute to risk when paired with another MES risk gene. Analyzing 10,795 families from the SPARK dataset, we identified 97 MES risk genes forming 50 significant gene pairs, demonstrating a substantial association with ASD when considered jointly, but not individually. Our method leverages familial inheritance patterns and statistical analyses, refined by comparisons against control cohorts, to elucidate these gene pairs’ contribution to ASD liability. Furthermore, expression profile analyses of these genes in brain tissues underscore their relevance to ASD pathology. This study underscores the complexity of ASD’s genetic landscape, suggesting that gene combinations, beyond high impact single-gene mutations, significantly contribute to the disorder’s etiology and heterogeneity. Our findings pave the way for new avenues in understanding ASD’s genetic underpinnings and developing targeted therapeutic strategies.
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 7 matches between paragraphs and lines of code.
MahjaniLab/MES_Code
17c5b6ac0fbf57aac32245c7c7f8b6d3e04bbefa, 12 September 2024Availability: 1 check, the latest on 28 September 2026: the link answers
- 28 September 2026: the link answers
27 files
- 1_MES_screen/
1_identify_probands_batc , R, 114 linesh.R - 1_MES_screen/
2_grade.pl , Perl, 99 lines - 1_MES_screen/
3.5_collapse.R , R, 48 lines - 1_MES_screen/
3_elsewhere.R , R, 49 lines - 1_MES_screen/
4_table.R , R, 36 lines, 1 match - 1_MES_screen/
5_poll.pl , Perl, 40 lines - 1_MES_screen/
6_gather_all_SPARK.pl , Perl, 94 lines - 1_MES_screen/
7_SPARK_contingency_tabl , R, 63 lineses.R - 1_MES_screen/
Mendelian_inheritence_ch , R, 64 lineseck.R - 1_MES_screen/
parental_variants.pl , Perl, 71 lines - 2_MES_non_neuro/
equilibrium_test.R , R, 37 lines - 2_MES_non_neuro/
equilibrium_test_single_ , R, 31 linesgene.R - 2_MES_non_neuro/
example_commands_to_anno , Shell, 39 lines, 2 matchestate.sh - 2_MES_non_neuro/
filter.pl , Perl, 47 lines - 2_MES_non_neuro/
grab_counts.pl , Perl, 62 lines - 2_MES_non_neuro/
simulate_MES_genes.R , R, 143 lines, 1 match - 3_MES_synonymous/
2_grade.pl , Perl, 94 lines - 3_MES_synonymous/
6_gather_all_SPARK.pl , Perl, 101 lines - 3_MES_synonymous/
filter.pl , Perl, 27 lines - 3_MES_synonymous/
grab_counts.pl , Perl, 61 lines - 3_MES_synonymous/
non_neuro_example_comman , Shell, 43 lines, 1 matchds_to_annotate.sh - 4_Expression/
query_gene_expression.R , R, 47 lines, 2 matches - 5_Phenotypes/
Medical_query_ASD.R , R, 72 lines - 5_Phenotypes/
Medical_query_nonASD.R , R, 74 lines - 5_Phenotypes/
SCQ_query.R , R, 75 lines - LICENSE, License, 674 lines
- README, Text, 22 lines
Zenodo 19859234
Availability: 1 check, the latest on 28 September 2026: the link answers (HTTP 200)
- 28 September 2026: the link answers (HTTP 200)
27 files
- 1_MES_screen/
1_identify_probands_batc , R, 114 linesh.R - 1_MES_screen/
2_grade.pl , Perl, 99 lines - 1_MES_screen/
3.5_collapse.R , R, 48 lines - 1_MES_screen/
3_elsewhere.R , R, 49 lines - 1_MES_screen/
4_table.R , R, 36 lines - 1_MES_screen/
5_poll.pl , Perl, 40 lines - 1_MES_screen/
6_gather_all_SPARK.pl , Perl, 94 lines - 1_MES_screen/
7_SPARK_contingency_tabl , R, 63 lineses.R - 1_MES_screen/
Mendelian_inheritence_ch , R, 64 lineseck.R - 1_MES_screen/
parental_variants.pl , Perl, 71 lines - 2_MES_non_neuro/
equilibrium_test.R , R, 37 lines - 2_MES_non_neuro/
equilibrium_test_single_ , R, 31 linesgene.R - 2_MES_non_neuro/
example_commands_to_anno , Shell, 39 linestate.sh - 2_MES_non_neuro/
filter.pl , Perl, 47 lines - 2_MES_non_neuro/
grab_counts.pl , Perl, 62 lines - 2_MES_non_neuro/
simulate_MES_genes.R , R, 143 lines - 3_MES_synonymous/
2_grade.pl , Perl, 94 lines - 3_MES_synonymous/
6_gather_all_SPARK.pl , Perl, 101 lines - 3_MES_synonymous/
filter.pl , Perl, 27 lines - 3_MES_synonymous/
grab_counts.pl , Perl, 61 lines - 3_MES_synonymous/
non_neuro_example_comman , Shell, 43 linesds_to_annotate.sh - 4_Expression/
query_gene_expression.R , R, 47 lines - 5_Phenotypes/
Medical_query_ASD.R , R, 72 lines - 5_Phenotypes/
Medical_query_nonASD.R , R, 74 lines - 5_Phenotypes/
SCQ_query.R , R, 75 lines - LICENSE, License, 674 lines
- README, Text, 22 lines
Code availability
Code and resources used in this study is available both on GitHub at (https://
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:
- 2 repositories of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
- 50 scripts, each with its path and the digest of its content;
- 7 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
Whole genome sequencing from 1kGP is freely available from The International Genome Sample Resource (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, 28 September 2026: the first record
Recorded: type, language, journal, volume, issue, pages, dates, 5 authors, 2 keywords, 4 MeSH terms, 5 funders, 50 references.
Cite
This paper
Caballero, M., Satterstrom, F. K., de Rubeis, S., Buxbaum, J. D., & Mahjani, B. (2026). Identification of moderate effect size genes in autism spectrum disorder through a novel gene pairing approach. Communications biology, 9(1), 1146. https://
BibTeX
@article{caballero2026id
author = {Caballero, Madison and Satterstrom, F Kyle and de Rubeis, Silvia and Buxbaum, Joseph D and Mahjani, Behrang},
title = {{Identification of moderate effect size genes in autism spectrum disorder through a novel gene pairing approach}},
journal = {Communications biology},
year = {2026},
month = may,
volume = {9},
number = {1},
pages = {1146},
publisher = {Nature Publishing Group},
issn = {2399-3642},
doi = {10.1038/
url = {https://
pmid = {42218271},
pmcid = {PMC13518923}
}
RIS
TY - JOUR
AU - Caballero, Madison
AU - Satterstrom, F Kyle
AU - de Rubeis, Silvia
AU - Buxbaum, Joseph D
AU - Mahjani, Behrang
TI - Identification of moderate effect size genes in autism spectrum disorder through a novel gene pairing approach
T2 - Communications biology
J2 - Commun Biol
PY - 2026
DA - 2026/
VL - 9
IS - 1
SP - 1146
SN - 2399-3642
PB - Nature Publishing Group
DO - 10.1038/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1038/
"type": "article-journal",
"title": "Identification of moderate effect size genes in autism spectrum disorder through a novel gene pairing approach",
"container-title": "Communications biology",
"author": [
{
"family": "Caballero",
"given": "Madison"
},
{
"family": "Satterstrom",
"given": "F Kyle"
},
{
"family": "de Rubeis",
"given": "Silvia"
},
{
"family": "Buxbaum",
"given": "Joseph D"
},
{
"family": "Mahjani",
"given": "Behrang"
}
],
"container-title-short":
"volume": "9",
"issue": "1",
"page": "1146",
"DOI": "10.1038/
"PMID": "42218271",
"PMCID": "PMC13518923",
"ISSN": "2399-3642",
"publisher": "Nature Publishing Group",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
5,
30
]
]
}
}
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/s41467-026-72598-z [code]
- Functional impact of genetic background on variable expressivity in neurodevelopmental disorders.Journal: Nature communicationsIn common: BCFtools, data.table, ggplot2, 1 other tool, 6 references
- [2] doi:10.1016/j.xhgg.2026.100652 [code]
- CRISPR-engineered deletion of POGZ alters transcription factor binding at promoters of genes involved in synaptic signaling.Journal: HGG advancesIn common: data.table, ggplot2, tidyverse, autism, cellular / molecular, 6 references
- [3] doi:10.1038/s41588-026-02646-3 [code]
- Co-expression-based models improve eQTL predictions for transcriptome-wide association studies and highlight new schizophrenia-associated
genes. Journal: Nature geneticsIn common: data.table, tidyverse, genetics / omics, cellular / molecular, 2 references, author Joseph D Buxbaum - [4] doi:10.1038/s41467-026-76675-1 [code]
- Long-read proteogenomic atlas of human neuronal differentiation reveals isoform diversity informing neurodevelopmental risk mechanisms.Journal: Nature communicationsIn common: data.table, ggplot2, tidyverse, autism, genetics / omics, 4 references
- [5] doi:10.3390/genes17080874 [code]
- Multi-Omic Analysis of Cerebrospinal Fluid Metabolites in Autism Spectrum Disorder: Biomarker Identification, Metabolic Genetics Insights, and Network Toxicology.Journal: GenesIn common: autism, genetics / omics, 6 references
- [6] doi:10.1038/s41588-026-02570-6 [code]
- Somatic mosaicism in ALS and FTD identifies focal mutations associated with widespread degeneration.Journal: Nature geneticsIn common: ggplot2, tidyverse, genetics / omics, cellular / molecular, 5 references
- [7] doi:10.1038/s41467-026-76676-0 [code]
- Determinants of functional burden pleiotropy and gene dosage responses across human traits.Journal: Nature communicationsIn common: data.table, ggplot2, tidyverse, genetics / omics, cellular / molecular, 4 references
- [8] doi:10.1111/ejn.70638 [code]
- Autism Spectrum Disorder-Associated Genes Enrich in Discrete Cortical, Limbic and Cerebellar Brain Regions.Journal: The European journal of neuroscienceIn common: autism, cellular / molecular, 5 references
- [9] doi:10.1101/gr.280394.124 [code]
- De novo structural variants in autism spectrum disorder disrupt distal regulatory interactions of neuronal genes.Journal: Genome researchIn common: autism, cellular / molecular, 5 references
- [10] doi:10.1038/s42003-026-10059-5 [code]
- Transcriptomic analysis in autism spectrum disorder suggests three molecular subtypes with distinct phenotypic profiles and functional pathways.Journal: Communications biologyIn common: data.table, ggplot2, autism, genetics / omics, cellular / molecular, 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, 50 scripts, and 7 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:43fe13a7e7006118…
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.
