name: 'tumor-clonal-evolution-agent'
description: 'AI-powered analysis of tumor clonal architecture, subclonal dynamics, and evolutionary trajectories from multi-region sequencing and longitudinal liquid biopsy data.'
measurable_outcome: Execute skill workflow successfully with valid output within 15 minutes.
allowed-tools:
- read_file
- run_shell_command
Tumor Clonal Evolution Agent
The Tumor Clonal Evolution Agent analyzes intratumoral heterogeneity (ITH), reconstructs tumor phylogenies, and tracks clonal dynamics over time. It integrates multi-region sequencing data, longitudinal liquid biopsies, and mathematical modeling to predict treatment response and resistance emergence.
When to Use This Skill
- When analyzing multi-region tumor sequencing to map spatial heterogeneity.
- To reconstruct tumor phylogenetic trees and identify ancestral mutations.
- For tracking clonal evolution through serial liquid biopsy samples.
- To predict time to treatment failure using evolutionary modeling.
- When identifying resistance-conferring subclones before clinical progression.
Core Capabilities
Clonal Deconvolution: Identifies tumor subpopulations and estimates their cellular fractions using variant allele frequencies (VAF) from bulk sequencing.
Phylogenetic Reconstruction: Builds tumor evolutionary trees showing relationships between subclones and their mutational acquisition order.
Longitudinal Tracking: Monitors subclone dynamics over time using ctDNA variant frequencies from serial blood draws.
Resistance Prediction: Applies Bayesian evolutionary frameworks to forecast emergence of resistant clones and time to progression.
Spatial ITH Mapping: Integrates multi-region data to visualize spatial distribution of subclones across tumor sites.
Fitness Estimation: Calculates subclone fitness parameters to identify aggressive populations driving tumor progression.
Workflow
Input: Multi-region or longitudinal mutation data (VCF/MAF), tumor purity estimates, copy number profiles.
Clustering: Cluster mutations into subclones using PyClone, SciClone, or MOBSTER.
Phylogeny: Reconstruct evolutionary trees using CITUP, PhyloWGS, or CALDER.
Modeling: Apply mathematical models (Lotka-Volterra, birth-death) to estimate dynamics.
Prediction: Forecast treatment response and resistance timeline.
Output: Phylogenetic trees, subclone trajectories, resistance predictions, actionable insights.
Example Usage
User: "Analyze the clonal evolution from these 6 longitudinal ctDNA samples and predict time to progression."
Agent Action:
python3 Skills/Oncology/Tumor_Clonal_Evolution_Agent/clonal_evolution.py \
--input longitudinal_ctdna_variants.maf \
--timepoints 0,4,8,12,16,20 \
--tumor_burden cea_values.csv \
--method bayesian_evolution \
--predict_ttp true \
--output evolution_analysis/
Key Methods and Algorithms
| Tool/Method |
Application |
Reference |
| PyClone-VI |
Bayesian clustering of mutations |
Nature Methods 2014 |
| MOBSTER |
Subclonal deconvolution with selection |
Nature Genetics 2020 |
| PhyloWGS |
Phylogenetic tree reconstruction |
Genome Biology 2015 |
| CALDER |
Copy-number aware phylogeny |
Nature Methods 2019 |
| CHESS |
Cancer heterogeneity from single samples |
Cell Systems 2019 |
Mathematical Framework
The agent applies evolutionary dynamics models:
Lotka-Volterra Competition:
dNi/dt = ri * Ni * (1 - sum(aij * Nj) / Ki)
Where:
- Ni = population of subclone i
- ri = growth rate (fitness)
- aij = competition coefficient
- Ki = carrying capacity
VAF Dynamics Modeling:
- Serial ctDNA VAF measurements enable real-time fitness estimation
- Bayesian inference updates subclone parameters with each sample
- Monte Carlo simulations generate prediction intervals
Prerequisites
- Python 3.10+
- PyClone-VI, PhyloWGS, or MOBSTER
- Copy number calling tools (ASCAT, Sequenza)
- Statistical modeling (PyMC, Stan)
Related Skills
- ctDNA_Analysis - For cfDNA variant calling
- Liquid_Biopsy_Analysis - For blood-based biomarker detection
- Variant_Interpretation - For mutation annotation
Clinical Applications
- Treatment Selection: Identify dominant subclones to target
- Resistance Monitoring: Detect emerging resistant populations early
- Prognosis: Predict time to treatment failure
- Combination Therapy: Design strategies targeting multiple subclones
Author
AI Group - Biomedical AI Platform
1---2name: tumor-clonal-evolution-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-clonal-evolution-agent'20description: 'AI-powered analysis of tumor clonal architecture, subclonal dynamics, and evolutionary trajectories from multi-region sequencing and longitudinal liquid biopsy data.'21measurable_outcome: Execute skill workflow successfully with valid output within 15 minutes.22allowed-tools:23 - read_file24 - run_shell_command25---262728# Tumor Clonal Evolution Agent2930The **Tumor Clonal Evolution Agent** analyzes intratumoral heterogeneity (ITH), reconstructs tumor phylogenies, and tracks clonal dynamics over time. It integrates multi-region sequencing data, longitudinal liquid biopsies, and mathematical modeling to predict treatment response and resistance emergence.3132## When to Use This Skill3334* When analyzing multi-region tumor sequencing to map spatial heterogeneity.35* To reconstruct tumor phylogenetic trees and identify ancestral mutations.36* For tracking clonal evolution through serial liquid biopsy samples.37* To predict time to treatment failure using evolutionary modeling.38* When identifying resistance-conferring subclones before clinical progression.3940## Core Capabilities41421. **Clonal Deconvolution**: Identifies tumor subpopulations and estimates their cellular fractions using variant allele frequencies (VAF) from bulk sequencing.43442. **Phylogenetic Reconstruction**: Builds tumor evolutionary trees showing relationships between subclones and their mutational acquisition order.45463. **Longitudinal Tracking**: Monitors subclone dynamics over time using ctDNA variant frequencies from serial blood draws.47484. **Resistance Prediction**: Applies Bayesian evolutionary frameworks to forecast emergence of resistant clones and time to progression.49505. **Spatial ITH Mapping**: Integrates multi-region data to visualize spatial distribution of subclones across tumor sites.51526. **Fitness Estimation**: Calculates subclone fitness parameters to identify aggressive populations driving tumor progression.5354## Workflow55561. **Input**: Multi-region or longitudinal mutation data (VCF/MAF), tumor purity estimates, copy number profiles.57582. **Clustering**: Cluster mutations into subclones using PyClone, SciClone, or MOBSTER.59603. **Phylogeny**: Reconstruct evolutionary trees using CITUP, PhyloWGS, or CALDER.61624. **Modeling**: Apply mathematical models (Lotka-Volterra, birth-death) to estimate dynamics.63645. **Prediction**: Forecast treatment response and resistance timeline.65666. **Output**: Phylogenetic trees, subclone trajectories, resistance predictions, actionable insights.6768## Example Usage6970**User**: "Analyze the clonal evolution from these 6 longitudinal ctDNA samples and predict time to progression."7172**Agent Action**:73```bash74python3 Skills/Oncology/Tumor_Clonal_Evolution_Agent/clonal_evolution.py \75 --input longitudinal_ctdna_variants.maf \76 --timepoints 0,4,8,12,16,20 \77 --tumor_burden cea_values.csv \78 --method bayesian_evolution \79 --predict_ttp true \80 --output evolution_analysis/81```8283## Key Methods and Algorithms8485| Tool/Method | Application | Reference |86|-------------|-------------|-----------|87| PyClone-VI | Bayesian clustering of mutations | Nature Methods 2014 |88| MOBSTER | Subclonal deconvolution with selection | Nature Genetics 2020 |89| PhyloWGS | Phylogenetic tree reconstruction | Genome Biology 2015 |90| CALDER | Copy-number aware phylogeny | Nature Methods 2019 |91| CHESS | Cancer heterogeneity from single samples | Cell Systems 2019 |9293## Mathematical Framework9495The agent applies evolutionary dynamics models:9697**Lotka-Volterra Competition**:98```99dNi/dt = ri * Ni * (1 - sum(aij * Nj) / Ki)100```101102Where:103- Ni = population of subclone i104- ri = growth rate (fitness)105- aij = competition coefficient106- Ki = carrying capacity107108**VAF Dynamics Modeling**:109- Serial ctDNA VAF measurements enable real-time fitness estimation110- Bayesian inference updates subclone parameters with each sample111- Monte Carlo simulations generate prediction intervals112113## Prerequisites114115* Python 3.10+116* PyClone-VI, PhyloWGS, or MOBSTER117* Copy number calling tools (ASCAT, Sequenza)118* Statistical modeling (PyMC, Stan)119120## Related Skills121122* ctDNA_Analysis - For cfDNA variant calling123* Liquid_Biopsy_Analysis - For blood-based biomarker detection124* Variant_Interpretation - For mutation annotation125126## Clinical Applications1271281. **Treatment Selection**: Identify dominant subclones to target1292. **Resistance Monitoring**: Detect emerging resistant populations early1303. **Prognosis**: Predict time to treatment failure1314. **Combination Therapy**: Design strategies targeting multiple subclones132133## Author134135AI Group - Biomedical AI Platform136137138<!-- AUTHOR_SIGNATURE: 9a7f3c2e-MD-BABU-MIA-2026-MSSM-SECURE -->