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
1---2name: tumor-heterogeneity-agent3description: <!--4---5<!--6# COPYRIGHT NOTICE7# This file is part of the "Universal Biomedical Skills" project.8# Copyright (c) 2026 MD BABU MIA, PhD <md.babu.mia@mssm.edu>9# All Rights Reserved.10#11# This code is proprietary and confidential.12# Unauthorized copying of this file, via any medium is strictly prohibited.13#14# Provenance: Authenticated by MD BABU MIA1516-->1718---19name: 'tumor-heterogeneity-agent'20description: 'AI-powered intratumor heterogeneity analysis for clonal architecture reconstruction, subclonal evolution tracking, and therapy resistance prediction using multi-region and longitudinal sequencing.'21measurable_outcome: Execute skill workflow successfully with valid output within 15 minutes.22allowed-tools:23 - read_file24 - run_shell_command25---262728# Tumor Heterogeneity Agent2930The **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.3132## When to Use This Skill3334* When analyzing multi-region tumor sequencing for clonal architecture.35* For tracking clonal evolution under treatment pressure.36* To predict resistance emergence from subclonal populations.37* When assessing tumor heterogeneity impact on treatment response.38* For integrating single-cell and bulk sequencing for ITH analysis.3940## Core Capabilities41421. **Clonal Deconvolution**: Infer clonal populations and their frequencies.43442. **Phylogeny Reconstruction**: Build tumor evolutionary trees from variants.45463. **Subclonal Tracking**: Monitor subclone dynamics over time.47484. **Resistance Prediction**: Identify pre-existing resistant subclones.49505. **Multi-Region Integration**: Combine spatial heterogeneity data.51526. **Single-Cell ITH**: Integrate scDNA-seq for ground-truth clones.5354## Heterogeneity Metrics5556| Metric | Definition | Clinical Relevance |57|--------|------------|-------------------|58| MATH Score | Mutant-allele tumor heterogeneity | ITH quantification |59| Shannon Index | Clonal diversity | Evolutionary potential |60| Clone Count | Number of distinct clones | Complexity |61| Truncal Fraction | % truncal mutations | Targetability |62| ITH Score | Composite heterogeneity | Prognosis |6364## Workflow65661. **Input**: Multi-region/longitudinal WES/WGS, copy number, tumor purity.67682. **Preprocessing**: Variant calling, CNV calling, purity estimation.69703. **CCF Estimation**: Calculate cancer cell fraction for each mutation.71724. **Clustering**: Group mutations into clonal populations.73745. **Phylogeny**: Reconstruct evolutionary tree.75766. **Temporal Analysis**: Track clone dynamics over time.77787. **Output**: Clone structures, phylogenies, heterogeneity metrics.7980## Example Usage8182**User**: "Analyze the clonal architecture of this multi-region lung tumor sequencing to understand heterogeneity and identify resistant subclones."8384**Agent Action**:85```bash86python3 Skills/Oncology/Tumor_Heterogeneity_Agent/ith_analysis.py \87 --multi_region_vcfs region1.vcf,region2.vcf,region3.vcf \88 --cnv_segments cnv_calls.seg \89 --purity 0.7,0.65,0.72 \90 --sample_names Primary,Met1,Met2 \91 --method pyclone-vi \92 --phylogeny_method citup \93 --output ith_analysis/94```9596## Deconvolution Methods9798| Method | Approach | Best For |99|--------|----------|----------|100| PyClone-VI | Variational inference | Large datasets |101| SciClone | Kernel density | High purity |102| EXPANDS | Probabilistic | Multi-region |103| Canopy | EM algorithm | CNV integration |104| Clonevol | Phylogeny-aware | Longitudinal |105| CITUP | Integer programming | Tree optimization |106107## Input Requirements108109| Input | Format | Required |110|-------|--------|----------|111| Somatic Variants | VCF with depth | Yes |112| Copy Number | SEG file | Yes |113| Tumor Purity | Float (0-1) | Yes |114| Sample Metadata | TSV | Yes |115| Normal BAM | BAM | Recommended |116117## Output Components118119| Output | Description | Format |120|--------|-------------|--------|121| Clone Assignments | Mutation-to-clone mapping | .csv |122| Clone Frequencies | Per-sample clone fractions | .csv |123| Phylogenetic Tree | Newick and visualization | .nwk, .pdf |124| ITH Metrics | Heterogeneity scores | .json |125| Subclone Variants | Clone-specific mutations | .vcf |126| Evolution Plot | Clone dynamics over time | .png |127| Actionable Subclones | Druggable clone mutations | .csv |128129## Clonal Classification130131| Clone Type | Definition | Implications |132|------------|------------|--------------|133| Truncal | Present in all samples | Ideal targets |134| Branch | Present in subset | Regional targets |135| Private | Single sample only | Local significance |136| Resistant | Expand under therapy | Resistance mechanism |137138## AI/ML Components139140**Clone Inference**:141- Variational autoencoders for CCF estimation142- Dirichlet process mixture models143- Graph neural networks for phylogeny144145**Resistance Prediction**:146- Time-series models for clone trajectories147- Classification of resistant signatures148- Drug-clone interaction prediction149150**Multi-Region Integration**:151- Multi-task learning across regions152- Spatial models for regional patterns153- Transfer learning across cancers154155## Clinical Applications156157| Application | ITH Insight | Clinical Action |158|-------------|-------------|-----------------|159| Treatment Selection | Truncal vs branch targets | Prioritize truncal targets |160| Resistance Monitoring | Pre-existing resistant clones | Early combination therapy |161| Prognosis | ITH score | Risk stratification |162| Biomarker Development | Clonal biomarkers | Robust biomarker selection |163164## Cancer-Specific Patterns165166| Cancer Type | Typical ITH | Key Drivers |167|-------------|-------------|-------------|168| Lung (NSCLC) | High | EGFR, KRAS subclonal |169| Breast | Moderate-High | PIK3CA, ESR1 evolution |170| Colorectal | Moderate | KRAS, BRAF clonal |171| Renal | Very High | VHL truncal, diverse branches |172| Melanoma | High | BRAF/NRAS truncal |173174## Prerequisites175176* Python 3.10+177* PyClone-VI, SciClone178* CITUP, Clonevol179* CNVkit/FACETS for CNV180* R with clonal evolution packages181182## Related Skills183184* ctDNA_Dynamics_MRD_Agent - Liquid biopsy tracking185* Single_Cell_CNV_Agent - scDNA-seq analysis186* HRD_Analysis_Agent - Genomic instability187* Pan_Cancer_MultiOmics_Agent - Multi-omic integration188189## Phylogeny Visualization190191| View Type | Shows | Best For |192|-----------|-------|----------|193| Fish Plot | Clone dynamics over time | Longitudinal |194| Tree Diagram | Branching evolution | Multi-region |195| Muller Plot | Population dynamics | Treatment response |196| Clone Map | Spatial distribution | Multi-region spatial |197198## Special Considerations1992001. **Sampling Bias**: Multi-region captures more heterogeneity2012. **Purity Effects**: Low purity reduces clone resolution2023. **CNV Complexity**: High CNV burden complicates CCF2034. **Single-Cell Validation**: Ground truth from scDNA-seq2045. **Temporal Resolution**: Frequent sampling improves tracking205206## Resistance Mechanisms207208| Mechanism | Detection | Intervention |209|-----------|-----------|--------------|210| Pre-existing resistant clone | Subclonal at baseline | Combination therapy |211| Acquired resistance | New clone emerges | Switch therapy |212| Phenotypic plasticity | Expression change | Monitor phenotype |213| Microenvironment | TME evolution | Immunotherapy |214215## Author216217AI Group - Biomedical AI Platform218219220<!-- AUTHOR_SIGNATURE: 9a7f3c2e-MD-BABU-MIA-2026-MSSM-SECURE -->