Abstract:Heterogeneity in characteristics from one region (sub-population) to another, commonly observed in complex systems, such as glasses and a collection of cells, is hard to describe theoretically. In the context of cancer, intra-tumor heterogeneity (ITH), characterized by cells with genetic and phenotypic variability that co-exist within a single tumor, is often the cause of ineffective therapy and recurrence of cancer. Multi-region sequencing (M-Seq), obtained by sampling multiple regions of a single tumor, has vividly demonstrated the pervasive nature of ITH, raising the need for a theory that accounts for evolution of tumor heterogeneity. Here, we develop a statistical mechanical theory to quantify ITH, using the Hamming distance, between genetic mutations in distinct regions within a single tumor. An analytic expression for ITH, expressed in terms of cell division probability ($\alpha$) and mutation probability ($p$), is validated using cellular-automaton type simulations. Application of the theory successfully captures ITH extracted from M-seq data in patients with exogenous cancers (melanoma and lung). The theory, based on punctuated evolution at the early stages of the tumor followed by neutral evolution, is accurate provided the spatial variation in the tumor mutation burden is not large. We show that there are substantial variations in ITH in distinct regions of a single solid tumor, which supports the notion that distinct subclones could co-exist. The simulations show that there are substantial variations in the sub-populations, with the ITH increasing as the distance between the regions increases. The analytical and simulation framework developed here could be used in the quantitative analyses of the experimental (M-Seq) data. More broadly, our theory is likely to be useful in analyzing dynamic heterogeneity in complex systems such as super-cooled liquids.

The Context Sensitivity Problem in Biological Sequence Segmentation

Extending the Recursive Jensen-Shannon Segmentation of Biological Sequences

New Stopping Criteria for Segmenting DNA Sequences

Bayesian Approach to DNA Segmentation into Regions with Different Average Nucleotide Composition

A (possible) mathematical model to describe biological "context-dependence" : case study with protein structure

Disentangling Cobionts and Contamination in Long-Read Genomic Data using Sequence Composition

Metagenome2Vec: Building Contextualized Representations for Scalable Metagenome Analysis

Genotype Sequence Segmentation: Handling Constraints and Noise

A Statistical Perspective on the Challenges in Molecular Microbial Biology

Sequence-based multiscale modeling for high-throughput chromosome conformation capture (Hi-C) data analysis

Computational Genome Analysis

Exploring Spatial Segregation Induced by Competition Avoidance as Driving Mechanism for Emergent Coexistence in Microbial Communities

Sequence Clustering in Bioinformatics: an Empirical Study.

Categorizing Sequences of Concern by Function To Better Assess Mechanisms of Microbial Pathogenesis

iSeg: an algorithm for segmentation of genomic data

Towards a core genome: pairwise similarity searches on interspecific genomic data

Microbial species exist and are maintained by ecological cohesiveness coupled to high homologous recombination.

Identification of Transcriptional Regulatory Elements in Chemosensory Receptor Genes by Probabilistic Segmentation

Statistical Mechanical theory for spatio-temporal evolution of Intra-tumor heterogeneity in cancers: Analysis of Multiregion sequencing data

Efficient Analysis of Annotation Colocalization Accounting for Genomic Contexts

Conserved intergenic sequences revealed by CTAG-profiling in Salmonella : thermodynamic modeling for function prediction