aipoch/go-kegg-enrichment
scientific-skills/Data Analysis/go-kegg-enrichment/SKILL.md
Performs GO (Gene Ontology) and KEGG pathway enrichment analysis on gene lists. Trigger when: - User provides a list of genes (symbols or IDs) and asks for enrichment analysis - User mentions "GO enrichment", "KEGG enrichment", "pathway analysis" - User wants to understand biological functions of gene sets - User provides differentially expressed genes (DEGs) and asks for interpretation - Input: gene list (file or inline), organism (human/mouse/rat), background gene set (optional) - Output: enriched terms, statistics, visualizations (barplot, dotplot, enrichment map)
$ install --global
skillsauthnpx skillsauth add aipoch/medical-research-skills go-kegg-enrichmentInstall this skill globally with one command. Works with Claude Code, Cursor, and Windsurf.
Security Scan Results
4 of 9 scanners reported clean
Some scanners were skipped, did not run, or reported a non-clean status. Review each row below.
4
Scanners Passed
9
Scanners in report
Clean
TrivyContainer and dependency vulnerability scanner
95%
Clean
SemgrepStatic code analysis for vulnerabilities
95%
Clean
mcp-scan (Snyk)Model Context Protocol security validation
95%
Skipped
Snyk (dep)Open source security scanning
50%
Skipped
Socket.devSupply chain security analysis
50%
Skipped
CrowdStrikeAdvanced threat intelligence
50%
Skipped
OSV-ScannerOpen Source Vulnerability database check
50%
Skipped
OWASP Dep-Check
50%
Last scanned: Apr 5, 2026, 12:52 AM41.1s1 file scanned
SKILL.md
- name:
- go-kegg-enrichment
- description:
- Performs GO (Gene Ontology) and KEGG pathway enrichment analysis on\
- \ gene lists.\nTrigger when:
- \n- User provides a list of genes (symbols or IDs)\
- - Input:
- gene list (file or inline), organism (human/mouse/rat), background gene\
- \ set (optional)\n- Output:
- enriched terms, statistics, visualizations (barplot,\
- version:
- 1.0.0
- category:
- Bioinfo
- tags:
- []
- author:
- AIPOCH
- license:
- MIT
- status:
- Draft
- risk_level:
- Medium
- skill_type:
- Tool/Script
- owner:
- AIPOCH
- last_updated:
- 2026-02-06
GO/KEGG Enrichment Analysis
Automated pipeline for Gene Ontology and KEGG pathway enrichment analysis with result interpretation and visualization.
Features
- GO Enrichment: Biological Process (BP), Molecular Function (MF), Cellular Component (CC)
- KEGG Pathway: Pathway enrichment with organism-specific mapping
- Multiple ID Support: Gene symbols, Entrez IDs, Ensembl IDs, RefSeq
- Statistical Methods: Hypergeometric test, Fisher's exact test, GSEA support
- Visualizations: Bar plots, dot plots, enrichment maps, cnet plots
- Result Interpretation: Automatic biological significance summary
Supported Organisms
| Common Name | Scientific Name | KEGG Code | OrgDB Package | |-------------|-----------------|-----------|---------------| | Human | Homo sapiens | hsa | org.Hs.eg.db | | Mouse | Mus musculus | mmu | org.Mm.eg.db | | Rat | Rattus norvegicus | rno | org.Rn.eg.db | | Zebrafish | Danio rerio | dre | org.Dr.eg.db | | Fly | Drosophila melanogaster | dme | org.Dm.eg.db | | Yeast | Saccharomyces cerevisiae | sce | org.Sc.sgd.db |
Usage
Basic Usage
# Run enrichment analysis with gene list
python scripts/main.py --genes gene_list.txt --organism human --output results/
Parameters
| Parameter | Description | Default | Required |
|-----------|-------------|---------|----------|
| --genes | Path to gene list file (one gene per line) | - | Yes |
| --organism | Organism code (human/mouse/rat/zebrafish/fly/yeast) | human | No |
| --id-type | Gene ID type (symbol/entrez/ensembl/refseq) | symbol | No |
| --background | Background gene list file | all genes | No |
| --pvalue-cutoff | P-value cutoff for significance | 0.05 | No |
| --qvalue-cutoff | Adjusted p-value (q-value) cutoff | 0.2 | No |
| --analysis | Analysis type (go/kegg/all) | all | No |
| --output | Output directory | ./enrichment_results | No |
| --format | Output format (csv/tsv/excel/all) | all | No |
Advanced Usage
# GO enrichment only with specific ontology
python scripts/main.py \
--genes deg_upregulated.txt \
--organism mouse \
--analysis go \
--go-ontologies BP,MF \
--pvalue-cutoff 0.01 \
--output go_results/
# KEGG enrichment with custom background
python scripts/main.py \
--genes treatment_genes.txt \
--background all_expressed_genes.txt \
--organism human \
--analysis kegg \
--qvalue-cutoff 0.05 \
--output kegg_results/
Input Format
Gene List File
TP53
BRCA1
EGFR
MYC
KRAS
PTEN
With Expression Values (for GSEA)
gene,log2FoldChange
TP53,2.5
BRCA1,-1.8
EGFR,3.2
Output Files
output/
├── go_enrichment/
│ ├── GO_BP_results.csv # Biological Process results
│ ├── GO_MF_results.csv # Molecular Function results
│ ├── GO_CC_results.csv # Cellular Component results
│ ├── GO_BP_barplot.pdf # Visualization
│ ├── GO_MF_dotplot.pdf
│ └── GO_summary.txt # Interpretation summary
├── kegg_enrichment/
│ ├── KEGG_results.csv # Pathway results
│ ├── KEGG_barplot.pdf
│ ├── KEGG_dotplot.pdf
│ └── KEGG_pathview/ # Pathway diagrams
└── combined_report.html # Interactive report
Result Interpretation
The tool automatically generates biological interpretation including:
- Top Enriched Terms: Significant GO terms/pathways ranked by enrichment ratio
- Functional Themes: Clustered biological themes from enriched terms
- Key Genes: Core genes driving enrichment in significant terms
- Network Relationships: Gene-term relationship visualization
- Clinical Relevance: Disease associations (for human genes)
Technical Difficulty: HIGH
⚠️ AI自主验收状态: 需人工检查
This skill requires:
- R/Bioconductor environment with clusterProfiler
- Multiple annotation databases (org.*.eg.db)
- KEGG REST API access
- Complex visualization dependencies
Dependencies
Required R Packages
install.packages(c("BiocManager", "ggplot2", "dplyr", "readr"))
BiocManager::install(c(
"clusterProfiler",
"org.Hs.eg.db", "org.Mm.eg.db", "org.Rn.eg.db",
"enrichplot", "pathview", "DOSE"
))
Python Dependencies
pip install pandas numpy matplotlib seaborn rpy2
Example Workflow
- Prepare Input: Create gene list from DEG analysis
- Run Analysis: Execute main.py with appropriate parameters
- Review Results: Check generated CSV files and visualizations
- Interpret: Read auto-generated summary for biological insights
References
See references/ for:
- clusterProfiler documentation
- KEGG API guide
- Statistical methods explanation
- Visualization examples
Limitations
- Requires internet connection for KEGG database queries
- Large gene lists (>5000) may require increased memory
- Some pathways may not be available for all organisms
- KEGG API has rate limits (max 3 requests/second)
Risk Assessment
| Risk Indicator | Assessment | Level | |----------------|------------|-------| | Code Execution | Python/R scripts executed locally | Medium | | Network Access | No external API calls | Low | | File System Access | Read input files, write output files | Medium | | Instruction Tampering | Standard prompt guidelines | Low | | Data Exposure | Output files saved to workspace | Low |
Security Checklist
- [ ] No hardcoded credentials or API keys
- [ ] No unauthorized file system access (../)
- [ ] Output does not expose sensitive information
- [ ] Prompt injection protections in place
- [ ] Input file paths validated (no ../ traversal)
- [ ] Output directory restricted to workspace
- [ ] Script execution in sandboxed environment
- [ ] Error messages sanitized (no stack traces exposed)
- [ ] Dependencies audited
Prerequisites
# Python dependencies
pip install -r requirements.txt
Evaluation Criteria
Success Metrics
- [ ] Successfully executes main functionality
- [ ] Output meets quality standards
- [ ] Handles edge cases gracefully
- [ ] Performance is acceptable
Test Cases
- Basic Functionality: Standard input → Expected output
- Edge Case: Invalid input → Graceful error handling
- Performance: Large dataset → Acceptable processing time
Lifecycle Status
- Current Stage: Draft
- Next Review Date: 2026-03-06
- Known Issues: None
- Planned Improvements:
- Performance optimization
- Additional feature support
Related Skills
aipoch/conventional-oncology-hub-gene
tools
VerifiedTrustedCommunity
Generates complete conventional oncology bulk-transcriptome biomarker and hub-gene research designs from a user-provided cancer type and study direction. Always use this skill whenever a user wants to design, plan, or build a tumor bioinformatics study centered on differential expression, prognostic filtering or risk modeling, PPI-based hub-gene prioritization, diagnostic/prognostic evaluation, clinical association, immune infiltration context, methylation context, and optional tissue or cell validation. Covers five study patterns (signature-first prognostic workflow, hub-gene-first biomarker workflow, hybrid signature-to-hub workflow, immune-context biomarker workflow, translational validation workflow) and always outputs four workload configs (Lite / Standard / Advanced / Publication+) with recommended primary plan, step-by-step workflow, figure plan, validation strategy, minimal executable version, publication upgrade path...
348SKILL.mdUpdated Apr 28, 2026
aipoch/conventional-non-oncology-hub-gene
development
VerifiedTrustedCommunity
Generates complete conventional non-oncology bioinformatics research designs from a user-provided disease context, process-related gene family or biological theme, and validation direction. Use when a study centers on multi-dataset bulk transcriptome integration, DEG analysis, process-gene intersection, enrichment analysis, GSEA, PPI hub-gene prioritization, TF/miRNA regulatory networks, ROC-based biomarker evaluation, and immune infiltration analysis. Covers five study patterns (process-DEG discovery, enrichment/GSEA interpretation, hub-gene prioritization, regulatory-network and immune interpretation, multi-layer public validation) and always outputs Lite / Standard / Advanced / Publication+ with a recommended primary plan, stepwise workflow, figure plan, validation hierarchy, minimal executable version, publication upgrade path, and strictly verified literature retrieval.
348SKILL.mdUpdated Apr 28, 2026
aipoch/confounder-and-bias-control-planner
tools
VerifiedTrustedCommunity
Plans confounder control, variable adjustment logic, and bias mitigation strategies at the protocol stage for clinical, epidemiologic, translational, observational, and biomarker studies. Always use this skill when a user needs to identify major confounders, decide which variables should or should not be adjusted for, compare matching/stratification/weighting approaches, anticipate selection or measurement bias, or pressure-test a study design before execution. Focus on bias sensing, causal structure awareness, variable-role classification, and critical design review rather than generic statistical advice.
348SKILL.mdUpdated Apr 28, 2026
aipoch/comparative-network-toxicology-shared-mechanism-reference-grounded
testing
VerifiedTrustedCommunity
Generates complete comparative network-toxicology research designs from a user-provided exposure pair, shared toxic phenotype, and validation direction. Use when a study centers on two related exposures under one outcome and needs target collection, shared-vs-specific target decomposition, enrichment, PPI hub prioritization, docking, optional transcriptomic cross-checks, and conservative mechanistic synthesis. Covers five study patterns and always outputs Lite / Standard / Advanced / Publication+ with a recommended primary plan, stepwise workflow, figure plan, validation hierarchy, minimal executable version, publication upgrade path, and strictly verified literature retrieval.
348SKILL.mdUpdated Apr 28, 2026