tooluniverse-polygenic-risk-score

Build and interpret polygenic risk scores from GWAS summary statistics and genotype data.

1.7k|254|Updated Mar 3, 2025
One-click install
npx skills add https://github.com/mims-harvard/ToolUniverse --skill tooluniverse-polygenic-risk-score
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: tooluniverse-polygenic-risk-score
Source: https://github.com/mims-harvard/ToolUniverse/tree/main/plugins/tooluniverse/skills/tooluniverse-polygenic-risk-score
Command: npx skills add https://github.com/mims-harvard/ToolUniverse --skill tooluniverse-polygenic-risk-score

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill requires tooluniverse.

What problem does it solve?

Estimating genetic predisposition to complex diseases requires aggregating thousands of small-effect variants from GWAS data, a process that is error-prone when done manually and easily misinterpreted without population context.

Core Features & Use Cases

  • PRS Model Construction: Extract genome-wide significant variants (p < 5e-8) from GWAS Catalog and build weighted scores using beta coefficients or log odds ratios.
  • Individual Score Calculation: Compute personal PRS from genotype data (23andMe, AncestryDNA, VCF) and standardize to z-scores and population percentiles.
  • Risk Interpretation: Classify results into risk categories (low, average, elevated, high) with ancestry-aware caveats and LDSC-based heritability and genetic correlation analysis.
  • Use Case: A researcher queries GWAS Catalog for coronary artery disease associations, builds a weighted PRS model, then calculates a study participant's percentile relative to the population distribution.

Quick Start

Ask the agent to build a polygenic risk score for type 2 diabetes using genome-wide significant variants from GWAS Catalog and interpret the resulting percentile.

Frequently Asked Questions about tooluniverse-polygenic-risk-score

High-intent search queries and answers about installing and using this skill.

FAQPage Schema
How do I build a polygenic risk score from GWAS data?

Query GWAS Catalog for genome-wide significant associations (p < 5e-8) for your trait, extract beta coefficients or convert odds ratios to log-odds, then compute a weighted sum of effect allele dosages. Production PRS should also apply LD clumping to avoid double-counting correlated variants.

How to calculate PRS from 23andMe raw data?

Export your raw genotype file, map each rsID to its allele pair, then count effect allele dosages (0, 1, or 2) for each SNP in the PRS model. Multiply dosages by effect sizes and sum across SNPs, then standardize against population mean and standard deviation.

Does PRS accuracy depend on ancestry?

Yes, PRS accuracy drops substantially when GWAS discovery ancestry differs from the target population. A European-derived PRS can lose 50-70% of predictive power in West African populations due to differences in allele frequencies and linkage disequilibrium patterns.

What is the difference between PRS percentile and absolute disease risk?

PRS percentile indicates relative position in the population genetic risk distribution, not absolute disease probability. A high percentile means elevated risk compared to average, but most high-PRS individuals never develop the disease because environmental factors contribute substantially.

When should I use PGS Catalog models instead of building a PRS?

Use validated PGS Catalog models for clinical-grade applications because they include LD clumping, ancestry-specific weights, and external validation. Building a PRS from raw GWAS associations is appropriate for research exploration and educational purposes.

Why does my PRS explain so little of disease heritability?

PRS captures only common variant effects, typically explaining 10-20% of trait variance while total heritability may be 50% or higher. Rare variants, epistasis, and gene-environment interactions are not captured by standard PRS models.