Abstract:Mutational processes of diverse origin leave their imprints in the genome during tumour evolution. These imprints are called and they have been characterised for point mutations, structural variants and copy number changes. Each signature has an , or abundance, per sample, which indicates how much a process has contributed to the overall genomic change. Mutational processes are not static, and a better understanding of their dynamics is key to characterise tumour evolution and identify cancer weaknesses that can be exploited during treatment. However, the structure of the data typically collected in this context makes it difficult to test whether signature exposures differ between samples or time-points. In general, the data consist of (1) patient-dependent vectors of counts for each sample and clonality group (2) generated from a covariate-dependent and compositional vector of probabilities with (3) a possibly group-dependent over-dispersion level. To model these data, we build on the Dirichlet-multinomial model to be able to model multivariate overdispersed vectors of counts as well as within-sample dependence and positive correlations between signatures. To estimate the model parameters, we implement a maximum likelihood estimator with a Laplace approximation of the random effect high-dimensional integrals and assess its bias and coverage by means of Monte Carlo simulations. We apply our approach to characterise differences of mutational processes between clonal and subclonal mutations across 23 cancer types of the PCAWG cohort. We find ubiquitous differential abundance of clonal and subclonal signatures across cancer types, and higher dispersion of signatures in the subclonal group, indicating higher variability between patients at subclonal level, possibly due to the presence of different clones with distinct active mutational processes. Mutational signature analysis is an expanding field and we envision our framework to be used widely to detect global changes in mutational process activity.

Single-cell mutational burden distributions in birth-death processes

Phylodynamic inference on single-cell data

Measures of genetic diversification in somatic tissues at bulk and single-cell resolution

Extreme value statistics of mutation accumulation in renewing cell populations

A Dirichlet-multinomial mixed model for determining differential abundance of mutational signatures

Sequential mutations in exponentially growing populations

How mutation accumulation depends on the structure of the cell lineage tree

Closed-form stochastic solutions for non-equilibrium dynamics and inheritance of cellular components over many cell divisions

An Exact, Time-Independent Approach to Clone Size Distributions in Normal and Mutated Cells

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

Inference of mutability landscapes of tumors from single cell sequencing data

Evolutionary dynamics of tumor progression with random fitness values

Single-cell genomic variation induced by mutational processes in cancer

Use of high throughput sequencing to observe genome dynamics at a single cell level

Modelling the order of driver mutations and metabolic mutations as structures in cancer dynamics

Stochastic modelling, Bayesian inference, and new in vivo measurements elucidate the debated mtDNA bottleneck mechanism

Age-dependent Branching Processes and Applications to the Luria-Delbrück Experiment

Hierarchical Bayesian analysis of somatic mutation data in cancer

Single-Cell Technologies for Studying the Evolution and Function of Mitochondrial DNA Heteroplasmy in Cancer

Mathematical modeling and quantitative analysis of phenotypic plasticity during tumor evolution based on single-cell data

Genetic Composition of Supercritical Branching Populations under Power Law Mutation Rates