qzzqzzb/variant-gwas-associations

EBI GWAS Catalog — Trait Associations

Usage

Tool Description

Query EBI GWAS Catalog REST API to get GWAS association information for a variant.
Database: EBI GWAS Catalog (https://www.ebi.ac.uk/gwas/)
API: GET https://www.ebi.ac.uk/gwas/rest/api/associations/search/findByRsId?rsId={rsid}
Args:
    rs_id (str): dbSNP rsID (e.g. "rs7412")
Return:
    associations (list): Each contains pvalue, risk allele, trait/disease, source study.

Query Example

import requests

rs_id = "rs7412"
url = f"https://www.ebi.ac.uk/gwas/rest/api/associations/search/findByRsId?rsId={rs_id}"
resp = requests.get(url, headers={"Accept": "application/json"}, timeout=30).json()
associations = resp.get("_embedded", {}).get("associations", [])
print(f"[GWAS Catalog] {rs_id} 关联数: {len(associations)}")

for i, assoc in enumerate(associations[:10]):  # 展示前10条
    pval = assoc.get("pvalue", "N/A")
    beta = assoc.get("betaNum", "N/A")
    risk_freq = assoc.get("riskFrequency", "N/A")
    # 获取 trait 名称需要额外请求 _links
    trait_link = assoc.get("_links", {}).get("efoTraits", {}).get("href", "")
    if trait_link:
        trait_resp = requests.get(trait_link, headers={"Accept": "application/json"}, timeout=15).json()
        traits = [t.get("trait", "") for t in trait_resp.get("_embedded", {}).get("efoTraits", [])]
    else:
        traits = ["unknown"]
    print(f"  [{i+1}] p={pval}, beta={beta}, traits={traits}")