Agent skill
tumor-heterogeneity-agent
Install this agent skill to your Project
npx add-skill https://github.com/FreedomIntelligence/OpenClaw-Medical-Skills/tree/main/skills/tumor-heterogeneity-agent
SKILL.md
name: 'tumor-heterogeneity-agent' description: 'AI-powered intratumor heterogeneity analysis for clonal architecture reconstruction, subclonal evolution tracking, and therapy resistance prediction using multi-region and longitudinal sequencing.' measurable_outcome: Execute skill workflow successfully with valid output within 15 minutes. allowed-tools:
- read_file
- run_shell_command
Tumor Heterogeneity Agent
The Tumor Heterogeneity Agent provides comprehensive analysis of intratumor heterogeneity (ITH) for understanding clonal architecture, tracking subclonal evolution, and predicting therapy resistance. It integrates multi-region sequencing, single-cell data, and longitudinal samples to reconstruct tumor phylogenies and identify actionable subclones.
When to Use This Skill
- When analyzing multi-region tumor sequencing for clonal architecture.
- For tracking clonal evolution under treatment pressure.
- To predict resistance emergence from subclonal populations.
- When assessing tumor heterogeneity impact on treatment response.
- For integrating single-cell and bulk sequencing for ITH analysis.
Core Capabilities
-
Clonal Deconvolution: Infer clonal populations and their frequencies.
-
Phylogeny Reconstruction: Build tumor evolutionary trees from variants.
-
Subclonal Tracking: Monitor subclone dynamics over time.
-
Resistance Prediction: Identify pre-existing resistant subclones.
-
Multi-Region Integration: Combine spatial heterogeneity data.
-
Single-Cell ITH: Integrate scDNA-seq for ground-truth clones.
Heterogeneity Metrics
| Metric | Definition | Clinical Relevance |
|---|---|---|
| MATH Score | Mutant-allele tumor heterogeneity | ITH quantification |
| Shannon Index | Clonal diversity | Evolutionary potential |
| Clone Count | Number of distinct clones | Complexity |
| Truncal Fraction | % truncal mutations | Targetability |
| ITH Score | Composite heterogeneity | Prognosis |
Workflow
-
Input: Multi-region/longitudinal WES/WGS, copy number, tumor purity.
-
Preprocessing: Variant calling, CNV calling, purity estimation.
-
CCF Estimation: Calculate cancer cell fraction for each mutation.
-
Clustering: Group mutations into clonal populations.
-
Phylogeny: Reconstruct evolutionary tree.
-
Temporal Analysis: Track clone dynamics over time.
-
Output: Clone structures, phylogenies, heterogeneity metrics.
Example Usage
User: "Analyze the clonal architecture of this multi-region lung tumor sequencing to understand heterogeneity and identify resistant subclones."
Agent Action:
python3 Skills/Oncology/Tumor_Heterogeneity_Agent/ith_analysis.py \
--multi_region_vcfs region1.vcf,region2.vcf,region3.vcf \
--cnv_segments cnv_calls.seg \
--purity 0.7,0.65,0.72 \
--sample_names Primary,Met1,Met2 \
--method pyclone-vi \
--phylogeny_method citup \
--output ith_analysis/
Deconvolution Methods
| Method | Approach | Best For |
|---|---|---|
| PyClone-VI | Variational inference | Large datasets |
| SciClone | Kernel density | High purity |
| EXPANDS | Probabilistic | Multi-region |
| Canopy | EM algorithm | CNV integration |
| Clonevol | Phylogeny-aware | Longitudinal |
| CITUP | Integer programming | Tree optimization |
Input Requirements
| Input | Format | Required |
|---|---|---|
| Somatic Variants | VCF with depth | Yes |
| Copy Number | SEG file | Yes |
| Tumor Purity | Float (0-1) | Yes |
| Sample Metadata | TSV | Yes |
| Normal BAM | BAM | Recommended |
Output Components
| Output | Description | Format |
|---|---|---|
| Clone Assignments | Mutation-to-clone mapping | .csv |
| Clone Frequencies | Per-sample clone fractions | .csv |
| Phylogenetic Tree | Newick and visualization | .nwk, .pdf |
| ITH Metrics | Heterogeneity scores | .json |
| Subclone Variants | Clone-specific mutations | .vcf |
| Evolution Plot | Clone dynamics over time | .png |
| Actionable Subclones | Druggable clone mutations | .csv |
Clonal Classification
| Clone Type | Definition | Implications |
|---|---|---|
| Truncal | Present in all samples | Ideal targets |
| Branch | Present in subset | Regional targets |
| Private | Single sample only | Local significance |
| Resistant | Expand under therapy | Resistance mechanism |
AI/ML Components
Clone Inference:
- Variational autoencoders for CCF estimation
- Dirichlet process mixture models
- Graph neural networks for phylogeny
Resistance Prediction:
- Time-series models for clone trajectories
- Classification of resistant signatures
- Drug-clone interaction prediction
Multi-Region Integration:
- Multi-task learning across regions
- Spatial models for regional patterns
- Transfer learning across cancers
Clinical Applications
| Application | ITH Insight | Clinical Action |
|---|---|---|
| Treatment Selection | Truncal vs branch targets | Prioritize truncal targets |
| Resistance Monitoring | Pre-existing resistant clones | Early combination therapy |
| Prognosis | ITH score | Risk stratification |
| Biomarker Development | Clonal biomarkers | Robust biomarker selection |
Cancer-Specific Patterns
| Cancer Type | Typical ITH | Key Drivers |
|---|---|---|
| Lung (NSCLC) | High | EGFR, KRAS subclonal |
| Breast | Moderate-High | PIK3CA, ESR1 evolution |
| Colorectal | Moderate | KRAS, BRAF clonal |
| Renal | Very High | VHL truncal, diverse branches |
| Melanoma | High | BRAF/NRAS truncal |
Prerequisites
- Python 3.10+
- PyClone-VI, SciClone
- CITUP, Clonevol
- CNVkit/FACETS for CNV
- R with clonal evolution packages
Related Skills
- ctDNA_Dynamics_MRD_Agent - Liquid biopsy tracking
- Single_Cell_CNV_Agent - scDNA-seq analysis
- HRD_Analysis_Agent - Genomic instability
- Pan_Cancer_MultiOmics_Agent - Multi-omic integration
Phylogeny Visualization
| View Type | Shows | Best For |
|---|---|---|
| Fish Plot | Clone dynamics over time | Longitudinal |
| Tree Diagram | Branching evolution | Multi-region |
| Muller Plot | Population dynamics | Treatment response |
| Clone Map | Spatial distribution | Multi-region spatial |
Special Considerations
- Sampling Bias: Multi-region captures more heterogeneity
- Purity Effects: Low purity reduces clone resolution
- CNV Complexity: High CNV burden complicates CCF
- Single-Cell Validation: Ground truth from scDNA-seq
- Temporal Resolution: Frequent sampling improves tracking
Resistance Mechanisms
| Mechanism | Detection | Intervention |
|---|---|---|
| Pre-existing resistant clone | Subclonal at baseline | Combination therapy |
| Acquired resistance | New clone emerges | Switch therapy |
| Phenotypic plasticity | Expression change | Monitor phenotype |
| Microenvironment | TME evolution | Immunotherapy |
Author
AI Group - Biomedical AI Platform
Recommended Agent Skills
Expand your agent's capabilities with these related and highly-rated skills.
vcf-annotator
Annotate VCF variants with VEP, ClinVar, gnomAD frequencies, and ancestry-aware context. Generates prioritised variant reports.
chemist-analyst
Analyzes events through chemistry lens using molecular structure, reaction mechanisms, thermodynamics, kinetics, and analytical techniques (spectroscopy, chromatography, mass spectrometry). Provides insights on chemical processes, material properties, reaction pathways, synthesis, and analytical methods. Use when: Chemical reactions, material analysis, synthesis planning, process optimization, environmental chemistry. Evaluates: Molecular structure, reaction mechanisms, yield, selectivity, safety, environmental impact.
bio-alignment-io
Read, write, and convert multiple sequence alignment files using Biopython Bio.AlignIO. Supports Clustal, PHYLIP, Stockholm, FASTA, Nexus, and other alignment formats for phylogenetics and conservation analysis. Use when reading, writing, or converting alignment file formats.
sleep-analyzer
分析睡眠数据、识别睡眠模式、评估睡眠质量,并提供个性化睡眠改善建议。支持与其他健康数据的关联分析。
metabolomics-workbench-database
Access NIH Metabolomics Workbench via REST API (4,200+ studies). Query metabolites, RefMet nomenclature, MS/NMR data, m/z searches, study metadata, for metabolomics and biomarker discovery.
bio-hi-c-analysis-matrix-operations
Balance, normalize, and transform Hi-C contact matrices using cooler and cooltools. Apply iterative correction (ICE), compute expected values, and generate observed/expected matrices. Use when normalizing or transforming Hi-C matrices.
Didn't find tool you were looking for?