A Dirichlet-multinomial mixed model for determining differential abundance of mutational signatures Journal Article
| Authors: | Morrill Gavarró, L.; Couturier, D. L.; Markowetz, F. |
| Article Title: | A Dirichlet-multinomial mixed model for determining differential abundance of mutational signatures |
| Abstract: | Background: Mutational processes of diverse origin leave their imprints in the genome during tumour evolution. These imprints are called mutational signatures and they have been characterised for point mutations, structural variants and copy number changes. Each signature has an exposure, 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 cell vulnerabilities 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 conditions or time-points when comparing groups of samples. In general, the data consists of multivariate count mutational data (e.g. signature exposures) with two observations per patient, each reflecting a group. Results: We propose a mixed-effects Dirichlet-multinomial model: within-patient correlations are taken into account with random effects, possible correlations between signatures by making such random effects multivariate, and a group-specific dispersion parameter can deal with particularities of the groups. Moreover, the model is flexible in its fixed-effects structure, so that the two-group comparison can be generalised to several groups, or to a regression setting. 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. Conclusions: Mutational signature analysis is an expanding field and we envision our framework to be used widely to detect global changes in mutational process activity. Our methodology is available in the R package CompSign and offers an ample toolkit for the analysis and visualisation of differential abundance of compositional data such as, but not restricted to, mutational signatures. © The Author(s) 2025. |
| Keywords: | genes; cancer genomics; cancer evolution; clonal evolution; mutational processes; mutational signatures; mutational signature; mixed effects models; compositional data; cancer aetiology; compositional regression; differential abundance; mixed-effects model; cancer etiology; mutational process |
| Journal Title: | BMC Bioinformatics |
| Volume: | 26 |
| ISSN: | 1471-2105 |
| Publisher: | Biomed Central Ltd |
| Date Published: | 2025-02-18 |
| Start Page: | 59 |
| Language: | English |
| DOI: | 10.1186/s12859-025-06055-x |
| PROVIDER: | scopus |
| PMCID: | PMC11837616 |
| PUBMED: | 39966709 |
| DOI/URL: | |
| Notes: | Article -- Source: Scopus |
Altmetric
Citation Impact
BMJ Impact Analytics
MSK Authors
Related MSK Work