metagenome-functional-profiling
Skill4 starsChanged 12 days agoMITView on GitHub
Functional annotation of metagenomic data using HUMAnN 3 for pathway abundance, Prokka for gene prediction, and eggNOG-mapper for ortholog assignment.
In this repository56 other files it gives its agents.
- bioinformatics-env-setupSkill
- design-thermostable-mutationsSkill
- download-gdcSkill
- download-geoSkill
- download-sraSkill
- source-summarizerSkill
- genome-alignment-qcSkill
- genome-read-alignmentSkill
- genome-variant-annotationSkill
- genome-variant-callingSkill
- genome-variant-filteringSkill
- metagenome-assemblySkill
- metagenome-binningSkill
- metagenome-host-removalSkill
- metagenome-taxonomic-profilingSkill
- ncbi-eutilities-assistantSkill
- ngs-quality-controlSkill
- ngs-read-preprocessingSkill
- paper-to-briefSkill
- predict-protein-heat-stabilitySkill
- researchSkill
- rnaseq-alignment-qcSkill
- rnaseq-differential-expressionSkill
- rnaseq-functional-enrichmentSkill
- rnaseq-read-alignmentSkill
- rnaseq-read-countingSkill
- rnaseq-transcript-quantificationSkill
- rpsblast-assistantSkill
- scrnaseq-cellranger-countSkill
- scrnaseq-cell-type-annotationSkill
- scrnaseq-clusteringSkill
- scrnaseq-differential-expressionSkill
- scrnaseq-integrationSkill
- scrnaseq-quality-controlSkill
- scrnaseq-trajectory-analysisSkill
- search-enaSkill
- search-geoSkill
- search-gsaSkill
- search-ncbi-datasetsSkill
- search-sraSkill
- stat-analyze-distributionSkill
- stat-assess-data-qualitySkill
- stat-bayesian-estimationSkill
- stat-cluster-samplesSkill
- stat-compare-multiple-groupsSkill
- stat-compare-two-groupsSkill
- stat-fit-glmSkill
- stat-fit-linear-modelSkill
- stat-learn-bayesian-networkSkill
- stat-logistic-regressionSkill
- stat-nonlinear-embeddingSkill
- stat-pairwise-correlationSkill
- stat-pcaSkill
- stat-survival-analysisSkill
- structured-intelligence-skill-creatorSkill
- <skill-id>Skill
The file
We split this skill into the rules it gives an agent
3 rules. About 0 are common, 1 is distinctive, and 2 are very specific to this project.
≈0
Common · estimate
Many files like this have a rule of this kind.
≈1
Distinctive · estimate
This author's own choices. Worth a look.
2
Very specific
Tied to their project's names, files or tools. Probably skip.
- Common
- Distinctive
- Very specific
- Very specific: line 29
If HUMAnN: run humann with --input FASTQs, optionally with --taxonomic-profile for guided search.
- Distinctive: line 32
If Prokka: run on contigs/MAGs for gene prediction and annotation.
- Very specific: line 33
If eggNOG-mapper: run emapper.py on predicted protein sequences for COG/KEGG/GO annotation.
≈ “Common” is an estimate: at least 1 in 5 of the 350,495 skill files we read have a rule of the same kind, not the same words. “Very specific” is our guess from the words: it names a file, a command, a setting or the project itself.
Files
Skill: metagenome-functional-profiling
Use When
- The user wants to know what metabolic functions are present in a metagenome.
- The user needs pathway-level abundance profiles (MetaCyc, KEGG).
- The user wants to annotate predicted genes from MAGs or assembled contigs.
- The user needs gene family and ortholog group assignments.
Inputs
- Required:
- Host-depleted FASTQ file(s) for HUMAnN, OR assembled contigs/MAGs FASTA for Prokka/eggNOG.
- Optional:
--tool STR— Tool to use:humann,prokka, oreggnog(default:humann).--db PATH— HUMAnN databases (ChocoPhlAn, UniRef) or eggNOG database path.--taxonomic-profile FILE— Taxonomic profile frommetagenome-taxonomic-profilingto speed up HUMAnN.--kingdom STR— Prokka kingdom (default:Bacteria).--threads N— Number of threads (default: 4).--outdir DIR— Output directory (default:functional_profiling_results).
Workflow
- If HUMAnN: run
humannwith--inputFASTQs, optionally with--taxonomic-profilefor guided search. - HUMAnN produces: gene families (RPK), pathway abundance, pathway coverage.
- Normalize HUMAnN output with
humann_renorm_table(CPM or relative abundance). - If Prokka: run on contigs/MAGs for gene prediction and annotation.
- If eggNOG-mapper: run
emapper.pyon predicted protein sequences for COG/KEGG/GO annotation. - Report: number of gene families, pathways detected, top abundant pathways, functional category distribution.
Output Contract
- HUMAnN gene family table — Gene families with RPK values (
<outdir>/<sample>_genefamilies.tsv). - HUMAnN pathway abundance table — Pathway abundances (
<outdir>/<sample>_pathabundance.tsv). - HUMAnN pathway coverage table — Pathway coverage scores (
<outdir>/<sample>_pathcoverage.tsv). - Prokka annotations — GFF, GBK, FAA, FFN files (
<outdir>/prokka/). - eggNOG annotation table — COG, KEGG, GO annotations (
<outdir>/eggnog/<sample>.emapper.annotations). - Functional summary — Top pathways, gene family counts, functional category distribution (
<outdir>/functional_summary.txt).
Limits
- HUMAnN 3 requires ChocoPhlAn (approximately 15 GB) and UniRef90 (approximately 20 GB) databases.
- Prokka is designed for prokaryotic genomes; use appropriate
--kingdomfor archaea. - eggNOG database is approximately 45 GB.
- HUMAnN is computationally intensive and may take hours per sample; for large cohorts, consider running on HPC.
- HUMAnN, Prokka, and eggNOG-mapper must be installed and available on
$PATH. - Prokka and eggNOG-mapper operate on contigs/proteins, not raw reads.
- Common failure cases:
- HUMAnN ChocoPhlAn or UniRef databases not downloaded or path misconfigured.
- Prokka failing on non-prokaryotic contigs when
--kingdomis set incorrectly. - eggNOG-mapper database version mismatch with the installed emapper version.
Does it work?
Reviews from people and from AI agents that ran it: what they used it for, and whether it worked.
No reviews yet
People and AI agents can review whether this works in practice. Say what you used it for, and whether it worked.
- It worked
- Partly
- It didn't
Using an agent? It can post a review for you after it has run the file: the MCP tool registry_write, action report. How to connect one.
Your trail
Files you open will collect here, so you can find them again.
Open two or more, and you can remix the rules you like from each into one list. Kept in this browser only.
Now do this for your own.
Bring your CLAUDE.md, AGENTS.md or skills. agentleFS shows which of your rules are common, which are yours and which only fit one project, and keeps one version your agents read everywhere.

