| name | gcell-pathway |
| description | Pathway enrichment analysis using gcell. Use this skill when users ask about:
- Gene set enrichment analysis
- GO (Gene Ontology) enrichment
- KEGG pathway analysis
- Reactome pathway enrichment
- Custom pathway/gene set analysis
Triggers: pathway enrichment, GO enrichment, KEGG, Reactome, gene set analysis, functional enrichment, ontology
|
Pathway Enrichment Analysis
Quick Enrichment with gprofiler
from gcell.ontology.pathway import gprofiler_enrichment
gene_list = ['TP53', 'BRCA1', 'MYC', 'EGFR', 'KRAS']
results = gprofiler_enrichment(gene_list, organism='hsapiens')
results = gprofiler_enrichment(
gene_list,
organism='hsapiens',
sources=['GO:BP', 'GO:MF', 'GO:CC', 'KEGG', 'REAC']
)
Working with Results
print(results.columns)
significant = results[results['p_value'] < 0.05]
top_terms = results.sort_values('p_value').head(20)
for _, row in top_terms.iterrows():
print(f"{row['term_name']}: {row['intersections']}")
Mouse and Other Organisms
results = gprofiler_enrichment(gene_list, organism='mmusculus')
results = gprofiler_enrichment(gene_list, organism='rnorvegicus')
Custom Pathways from GMT Files
from gcell.ontology.pathway import Pathways
pathways = Pathways.from_gmt('custom_pathways.gmt')
background_genes = [...]
enriched = pathways.enrichment(gene_list, background_genes)
Key Functions and Classes
| Name | Purpose |
|---|
gprofiler_enrichment() | Quick enrichment via g:Profiler |
Pathways | Custom pathway collections |
Pathways.from_gmt() | Load GMT format gene sets |
Pathways.enrichment() | Run enrichment analysis |
Tips
- Always use appropriate background genes when possible
- Multiple testing correction is applied automatically
- Use specific sources (e.g., just 'GO:BP') to reduce multiple testing burden
- Gene symbols should match the organism (human: HUGO symbols)