Transcriptome-wide association study identifies new susceptibility genes and pathways for spondyloarthritis.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:37667381", "kind": "science-record", "revision": 3, "content_hash": "sha256:71bc04f39cf8bdaed4a3d019b1980535f9eefa291362ea9959eaec8cd142fda4", "created_at": "2026-11-01T19:02:27Z", "updated_at": "2026-11-01T21:31:28Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Transcriptome-wide association study identifies new susceptibility genes and pathways for spondyloarthritis.", "doi": "10.1186/s13018-023-04029-4", "pmid": "37667381", "authors": ["Xiaochen Su", "Anfa Chen", "Menghao Teng", "Wenchen Ji", "Yingang Zhang"], "journal": "Journal of orthopaedic surgery and research", "year": 2023, "abstract": "BACKGROUND: Spondyloarthritis (SpA) is a group of multifactorial bone diseases influenced by genetic factors, the environment and lifestyle. However, current studies have found a limited number of SpA-related genes, and the genetic and pathogenic mechanisms of SpA are still unclear.\nMETHODS: A tissue-specific transcriptome-wide association study (TWAS) of SpA was performed using GWAS (including 3966 SpA patients and 448,298 controls) summary data and gene expression weights of whole blood and skeletal muscle. The SpA-associated genes identified by TWAS were further compared with the differentially expressed genes (DEGs) identified in the SpA gene expression profile acquired from the Gene Expression Omnibus database (GEO, GSE58667). Finally, functional enrichment and annotation analyses of the identified genes were performed.\nRESULTS: The TWAS detected 499 suggestive genes associated with SpA in whole blood and skeletal muscle, such as CTNNAL1 (PSM = 3.04 × 10-2, PWB = 9.58 × 10-3). The gene expression profile of SpA identified 20 candidate genes that overlapped in the TWAS data, such as MCM4 (PTWAS = 1.32 × 10-2, PDEG = 2.75 × 10-2) and KIAA1109 (PTWAS = 3.71 × 10-2, PDEG = 4.67 × 10-2). Enrichment analysis of the genes identified by TWAS identified 93 significant GO terms and 33 KEGG pathways, such as mitochondrion organization (GO: 0007005) and axon guidance (hsa04360).\nCONCLUSION: We identified multiple candidate genes that were genetically related to SpA. Our study may provide novel clues regarding the genetic mechanism, diagnosis, and treatment of SpA.", "claims": [{"id": "claim-method", "statement": "The TWAS of SpA was carried out by functional summary-based imputation software (FUSION http://gusevlab.org/projects/fusion/).", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:37667381", "locator": "sections[4] \"TWAS of SpA\" paragraph 2", "quote": "The TWAS of SpA was carried out by functional summary-based imputation software (FUSION http://gusevlab.org/projects/fusion/)."}], "state": "complete"}, {"id": "claim-data-access", "statement": "The datasets analyzed during the current study are available from the Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/gds) accession number: GSE58667. The GWAS summary data of SpA obtained from UK biobank (https://nealelab.github.io/UKBB_ldsc/downloads.html).", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:37667381", "locator": "sections[20] \"notes\" paragraph 1", "quote": "The datasets analyzed during the current study are available from the Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/gds) accession number: GSE58667. The GWAS summary data of SpA obtained from UK biobank (https://nealelab.github.io/UKBB_ldsc/downloads.html)."}], "state": "complete"}, {"id": "claim-target-result", "statement": "The TWAS detected 499 suggestive genes associated with SpA in whole blood and skeletal muscle.", "kind": "qualitative", "target": null, "evidence": [{"asset": "source:pmid:37667381", "locator": "sections[0] \"Abstract\" paragraph 3", "quote": "The TWAS detected 499 suggestive genes associated with SpA in whole blood and skeletal muscle, such as CTNNAL1 (PSM = 3.04 × 10−2, PWB = 9.58 × 10−3)."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "complete", "objective": "The TWAS detected 499 suggestive genes associated with SpA in whole blood and skeletal muscle.", "steps": [{"id": "read_published_estimate", "order": 1, "name": "Read the published estimate from the source", "kind": "other", "state": "complete", "blocked_on": null, "description": "The claim, method and evidence span published with the source article.", "duration_hours": null, "parameters": {}, "inputs": ["source:pmid:37667381"], "outputs": ["the published estimate and its evidence"], "instruments": [], "reagents": []}, {"id": "obtain_source_data", "order": 2, "name": "Obtain the source data from the named repository", "kind": "data-preparation", "state": "complete", "blocked_on": null, "description": "The GEO series accession resolves and its sample accessions are listed. This lane did not download the dataset or read its licence.", "duration_hours": null, "parameters": {}, "inputs": ["data:geo:GSE58667"], "outputs": ["the source data with its repository provenance"], "instruments": [], "reagents": []}, {"id": "run_analysis", "order": 3, "name": "Run the published analysis", "kind": "analysis", "state": "complete", "blocked_on": null, "description": "Apply the paper's statistical method to the resolved series and compare the result with the published target.", "duration_hours": null, "parameters": {}, "inputs": ["data:geo:GSE58667"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:37667381", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1186/s13018-023-04029-4", "sha256": "sha256:e7fe43fc198702034bbbf854237a30d31702bb0f07f980c9e3cf0b0f8df732fd", "bytes": 29779, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1186/s13018-023-04029-4", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC10478464", "source": "europe_pmc_fulltext", "coverage": "whole_body_without_references", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE58667", "role": "raw", "format": "geo-series", "media_type": "application/vnd.ncbi.geo-series+json", "uri": "https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE58667", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE58667", "version": null, "note": "The datasets analyzed during the current study are available from the Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/gds) accession number: GSE58667. The GWAS summary data of SpA obtained from UK biobank (https://nealelab.github.io/UKBB_ldsc/downloads.html).", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The datasets analyzed during the current study are available from the Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/gds) accession number: GSE58667. The GWAS summary data of SpA obtained from UK biobank (https://nealelab.github.io/UKBB_ldsc/downloads.html).", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM1416673", "GSM1416674", "GSM1416675", "GSM1416676", "GSM1416677", "GSM1416678", "GSM1416679", "GSM1416680", "GSM1416681", "GSM1416682", "GSM1416683", "GSM1416684", "GSM1416685", "GSM1416686", "GSM1416687"], "sample_count": 15, "group_sizes": null, "group_sizes_basis": null, "platform": "570", "taxon": "Homo sapiens", "title": "Gene expression in juvenile spondyloarthritis", "resolution_checked_at": "2026-11-01T19:09:46Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}], "gaps": ["asset:data:geo:GSE58667 not downloaded", "asset:data:geo:GSE58667 licence not read"]}, "analysis": {"state": "partial", "methods": [{"id": "method-1", "family": "genetic_association", "specific": "gwas", "label": "Genetic association (GWAS / linkage / MR)", "dsl_operations": ["gwas_locus_scan"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["genotype_phenotype_table"], "evidence": {"score": 0.7, "band": "medium", "quote": "GWAS"}}, {"id": "method-2", "family": null, "specific": null, "label": "transcriptome-wide association study (TWAS) using functional summary-based imputation", "dsl_operations": [], "missing_operations": [], "expressible": null, "parameters": {}, "data_shapes": [], "evidence": {"score": 1.0, "band": "verbatim", "quote": "The TWAS of SpA was carried out by functional summary-based imputation software (FUSION http://gusevlab.org/projects/fusion/)."}}], "targets": [], "gaps": ["method.reader_named_unmapped"]}, "capability": {"state": "complete", "requirements": [{"capability_id": "capability:computational_analysis", "label": "Computational analysis platform", "assay_families": ["meta_analysis", "machine_learning", "survival_analysis", "statistical_regression", "pca", "molecular_docking", "network_analysis", "differential_expression", "gwas"], "availability_state": "owned", "tier_unlocked": "R1", "candidate": "Automated Assays Laravel application and simulated adapters"}], "gaps": []}, "readiness": {"state": "complete", "tier": "R1", "blocking_axis": "none", "blocked_on": "nothing:replicable_now", "axes": {"method": {"tier": "R1", "expressible": true, "required_operations": ["gwas_locus_scan"], "missing_operations": [], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: gwas."}, "capability": {"tier": "R1", "capability_id": "capability:computational_analysis", "availability_state": "owned", "basis": "The computational analysis is owned. The physical bench the paper used is not, so a physical validation cannot be run here."}, "data": {"tier": "R1", "reachable": true, "accessions": ["GSE58667"], "resolved_accessions": ["GSE58667"], "statements": ["The datasets analyzed during the current study are available from the Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/gds) accession number: GSE58667. The GWAS summary data of SpA obtained from UK biobank (https://nealelab.github.io/UKBB_ldsc/downloads.html)."], "basis": "The paper names GSE58667 in Gene Expression Omnibus (GEO) and the accession resolves in the repository. No dataset was downloaded and no licence was read."}}, "can_run": {"dsl_operations": ["gwas_locus_scan"], "assets": ["source:pmid:37667381", "data:geo:GSE58667"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE58667"]}, "gaps": []}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-11-01T18:52:29Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 9293, "completion_tokens": 1824, "cached_tokens": 256}}, "human_critique": null, "stages": [{"id": "extraction", "state": "complete", "reader": "MiniMax-M3.1-Flash", "provider": "minimax"}], "gaps": ["review.independent_reader", "review.machine_critique", "review.human_critique"]}, "execution": {"state": "absent", "status": "not_run", "mode": null, "runs": [], "gaps": ["execution.not_run"]}, "provenance": {"entities": [{"id": "pmid:37667381", "type": "paper", "label": "pmid:37667381"}, {"id": "source:pmid:37667381", "type": "manuscript", "label": "PMID 37667381 full text"}, {"id": "data:geo:GSE58667", "type": "dataset", "label": "GEO GSE58667"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:37667381", "type": "paper", "label": "pmid:37667381"}, {"id": "source:pmid:37667381", "type": "manuscript", "label": "PMID 37667381 full text"}, {"id": "data:geo:GSE58667", "type": "dataset", "label": "GEO GSE58667"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-37667381", "type": "extraction", "used": ["source:pmid:37667381"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE58667", "type": "accession_resolution", "used": ["data:geo:GSE58667"], "agent": "agent:corpus-assay-engine"}], "agents": [{"id": "agent:corpus-assay-engine", "type": "software", "label": "Corpus assay engine"}, {"id": "agent:minimax-reader", "type": "software", "label": "MiniMax literature reader"}], "source_run_id": "reading-sci-2026-11-01", "dsl_catalogue_fingerprint": "343fb30756ba6ec3", "capability_table_fingerprint": "2cb3e42f047d2ecf", "snapshot_id": "ba34b49490d9bfaf"}, "lineage": {"parent_id": "pmid:37667381", "parent_hash": "sha256:1c19662b4649feb3f5b93803d4956dea03ac45bf55bf3371427bdb0676ae8ef4", "fork_of": null, "generation": 2}, "gaps": ["asset:data:geo:GSE58667 licence not read", "asset:data:geo:GSE58667 not downloaded", "execution.not_run", "method.reader_named_unmapped", "review.human_critique", "review.independent_reader", "review.machine_critique"], "assays": [{"role": "analysis", "definition_id": "genetic_association", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": ["gwas_locus_scan"], "support_state": "supported", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 0.7, "evidence": ["GWAS"], "order_basis": "analysis_only", "id": "assay:1"}]}