Skip to content
cccskills
FreePublished on GitHub

tooluniverse-metabolomics-analysis

Analyze metabolomics data end-to-end — metabolite identification, quantification (TIC normalization, batch correction), differential analysis, and pathway interpretation. Use for processing mass-spec metabolomics output, normalization choice, untargeted metabolomics workflows, and integrating with other omics layers.

The English summary is being prepared. Showing the original description.

See how to install

Files included (3)

  • SKILL.md8.1 KB
  • code_examples.md8.0 KB
  • report_template.md3.9 KB

SKILL.md (original)

Review the instructions the agent will receive before you install.

Metabolomics Analysis

Comprehensive analysis of metabolomics data from metabolite identification through quantification, statistical analysis, pathway interpretation, and integration with other omics layers.

Domain Reasoning

Metabolomics quantification depends critically on normalization. Total ion current (TIC) normalization corrects for sample-loading variation and works well for global abundance changes; internal standard normalization is more accurate for targeted analysis where specific metabolite concentrations matter. Missing values in a peak table may reflect signal below the detection limit — not true absence — and should be imputed or handled explicitly rather than treated as zero. Failing to account for batch effects across instrument runs is a frequent source of spurious differential metabolites.

LOOK UP DON'T GUESS

  • Metabolite identities: use Metabolite_search and Metabolite_get_info to confirm names, CIDs, and HMDB IDs; never assume identity from m/z alone.
  • Pathway memberships: query KEGG, MetaCyc, or Reactome tools; do not list pathways from memory.
  • Disease associations: retrieve from CTD via Metabolite_get_diseases; do not infer clinical relevance without database evidence.
  • CV thresholds and QC criteria: apply the values defined in this workflow (CV < 30%, blank ratio > 3x); do not override with guesses.

When to Use This Skill

Triggers:

  • User has metabolomics data (LC-MS, GC-MS, NMR)
  • Questions about metabolite abundance or concentrations
  • Differential metabolite analysis requests
  • Metabolic pathway analysis
  • Multi-omics integration with metabolomics
  • Metabolic biomarker discovery
  • Flux balance analysis or metabolic modeling
  • Metabolite-enzyme correlation

Example Questions:

  1. "Analyze this LC-MS metabolomics data for differential metabolites"
  2. "Which metabolic pathways are dysregulated between conditions?"
  3. "Identify metabolite biomarkers for disease classification"
  4. "Correlate metabolite levels with enzyme expression"
  5. "Perform pathway enrichment for differential metabolites"
  6. "Integrate metabolomics with transcriptomics data"

Core Capabilities

CapabilityDescription
Data ImportLC-MS, GC-MS, NMR, targeted/untargeted platforms
Metabolite IdentificationMatch to HMDB, KEGG, PubChem, spectral libraries
Quality ControlPeak quality, blank subtraction, internal standard normalization
NormalizationProbabilistic quotient, total ion current, internal standards
Statistical AnalysisUnivariate and multivariate (PCA, PLS-DA, OPLS-DA)
Differential AnalysisIdentify significant metabolite changes
Pathway EnrichmentKEGG, Reactome, BioCyc metabolic pathway analysis
Metabolite-Enzyme IntegrationCorrelate with expression data
Flux AnalysisMetabolic flux balance analysis (FBA)
Biomarker DiscoveryMulti-metabolite signatures

Workflow Overview

Input: Metabolomics Data (Peak Table or Spectra)
    |
    v
Phase 1: Data Import & Metabolite Identification
    |-- Load peak table or process raw spectra
    |-- Match features to HMDB, KEGG (accurate mass +/- 5 ppm)
    |-- Confidence scoring (Level 1-4)
    |
    v
Phase 2: Quality Control & Filtering
    |-- CV in QC samples (<30%)
    |-- Blank subtraction (sample/blank > 3)
    |-- Remove features with >50% missing
    |
    v
Phase 3: Normalization
    |-- Sample-wise: TIC, PQN, or internal standards
    |-- Transformation: log2, Pareto, or auto-scaling
    |-- Batch effect correction (if multi-batch)
    |
    v
Phase 4: Exploratory Analysis
    |-- PCA for sample clustering
    |-- PLS-DA for supervised separation
    |-- Outlier detection
    |
    v
Phase 5: Differential Analysis
    |-- t-test / ANOVA / Wilcoxon
    |-- Fold change + FDR correction
    |-- Volcano plots, heatmaps
    |
    v
Phase 6: Pathway Analysis
    |-- Metabolite set enrichment (MSEA)
    |-- KEGG/Reactome pathway mapping
    |-- Pathway topology (hub/bottleneck metabolites)
    |
    v
Phase 7: Multi-Omics Integration
    |-- Metabolite-enzyme Spearman correlation
    |-- Pathway-level concordance scoring
    |-- Metabolic flux inference
    |
    v
Phase 8: Generate Report
    |-- Summary statistics, differential metabolites
    |-- Pathway diagrams, biomarker panel

Phase Summaries

Phase 1: Data Import & Identification

Load peak tables (CSV/TSV) or process raw spectra (mzML). Match features to HMDB by accurate mass (+/- 5 ppm). Assign confidence levels: L1 (standard match), L2 (MS/MS), L3 (mass only), L4 (unknown).

Phase 2: Quality Control

Assess CV in QC samples (reject >30%), compute blank ratios (keep >3x blank), filter features with >50% missing values. Check internal standard recovery (95-105% acceptable).

Phase 3: Normalization

Three methods available: TIC (simple, assumes similar total abundance), PQN (robust to large changes, recommended), Internal Standard (most accurate with spiked standards). Follow with log2 transform or Pareto scaling.

Phase 4: Exploratory Analysis

PCA reveals sample grouping and batch effects. PLS-DA provides supervised separation (report R2 and Q2 for model quality). Flag and investigate outliers.

Phase 5: Differential Analysis

Welch's t-test (two groups) or ANOVA (multiple groups) with Benjamini-Hochberg FDR correction. Significance thresholds: adj. p < 0.05 and |log2FC| > 1.0.

Phase 6: Pathway Analysis

Map differential metabolites to KEGG compound IDs. Perform MSEA for pathway enrichment. Consider topology: metabolites at pathway hubs (high degree/betweenness centrality) have greater impact.

Phase 7: Multi-Omics Integration

Correlate metabolite levels with enzyme expression (Spearman). Expected: substrate-enzyme negative correlation (consumption), product-enzyme positive correlation (production). Score pathway dysregulation using combined metabolite + gene evidence.

Phase 8: Report

See report_template.md for full example output.


Integration with ToolUniverse

SkillUsed ForPhase
tooluniverse-gene-enrichmentPathway enrichmentPhase 6
tooluniverse-rnaseq-deseq2Enzyme expression for integrationPhase 7
tooluniverse-proteomics-analysisProtein levels for integrationPhase 7
tooluniverse-multi-omics-integrationComprehensive integrationPhase 7

Quantified Minimums

ComponentRequirement
MetabolitesAt least 50 identified metabolites
ReplicatesAt least 3 per condition
QCCV < 30% in QC samples, blank subtraction
Statistical testt-test or Wilcoxon with FDR correction
Pathway analysisMSEA with KEGG or Reactome
ReportQC, differential metabolites, pathways, visualizations

Limitations

  • Identification: Many features remain unidentified (Level 4)
  • Coverage: Cannot detect all metabolites (depends on method)
  • Quantification: Relative abundance (not absolute without standards)
  • Isomers: Difficult to distinguish structural isomers
  • Ion suppression: Matrix effects can affect quantification
  • Dynamic range: Limited compared to targeted methods

References

Methods:

Databases:


Reference Files

  • code_examples.md - Python code for all phases (data loading, QC, normalization, statistics, pathway analysis)
  • report_template.md - Full example report (LC-MS disease vs control)

Reviews

No reviews yet. Tell others what you thought.

Comments

No login needed. Share how it worked for you, ask a question, or note where it did not run.

No comments yet. Be the first to write one.

Comments are published right away. Do not include personal information such as email addresses or phone numbers. Up to 2 URLs.

0 / 1,000 characters

This form is protected by reCAPTCHA, and the Google Privacy Policy and Terms of Service apply.

More skills from this repository

Simplest distributed training API. 4 lines to add distributed support to any PyTorch script. Unified API for DeepSpeed/FSDP/Megatron/DDP. Automatic device placement, mixed precision (FP16/BF16/FP8). Interactive config, single launch command. HuggingFace ecosystem standard.

The English summary is being prepared. Showing the original description.

huang-sh/DeepScience4Updated Jul 15, 2026

ADMET property prediction for drug candidates. Full pharmacokinetic panel (Caco-2, PPB, clearance, CYP), toxicity (hERG, AMES, DILI), drug-likeness (Lipinski, QED), using RDKit descriptors and TDC models.

The English summary is being prepared. Showing the original description.

huang-sh/DeepScience4Updated Jul 15, 2026

Interpretable ADMET analysis with mechanistic reasoning. Maps liabilities to structural causes and biological pathways. Based on CoTox (Park 2025) and DrugR (Liu 2026).

The English summary is being prepared. Showing the original description.

huang-sh/DeepScience4Updated Jul 15, 2026

aeon

Free

This skill should be used for time series machine learning tasks including classification, regression, clustering, forecasting, anomaly detection, segmentation, and similarity search. Use when working with temporal data, sequential patterns, or time-indexed observations requiring specialized algorithms beyond standard ML approaches. Particularly suited for univariate and multivariate time series analysis with scikit-learn compatible APIs.

The English summary is being prepared. Showing the original description.

huang-sh/DeepScience4Updated Jul 15, 2026

Infer gene regulatory networks (GRNs) from gene expression data using scalable algorithms (GRNBoost2, GENIE3). Use when analyzing transcriptomics data (bulk RNA-seq, single-cell RNA-seq) to identify transcription factor-target gene relationships and regulatory interactions. Supports distributed computation for large-scale datasets.

The English summary is being prepared. Showing the original description.

huang-sh/DeepScience4Updated Jul 15, 2026

Comprehensive Python library for astronomy and astrophysics. This skill should be used when working with astronomical data including celestial coordinates, physical units, FITS files, cosmological calculations, time systems, tables, world coordinate systems (WCS), and astronomical data analysis. Use when tasks involve coordinate transformations, unit conversions, FITS file manipulation, cosmological distance calculations, time scale conversions, or astronomical data processing.

The English summary is being prepared. Showing the original description.

huang-sh/DeepScience4Updated Jul 15, 2026

See all skills by huang-sh

Report a problem with this skill