A site specific model and analysis of the neutral somatic mutation rate in whole-genome cancer data

Research output: Contribution to journal/Conference contribution in journal/Contribution to newspaperJournal articleResearchpeer-review

Documents

DOI

Background: Detailed modelling of the neutral mutational process in cancer cells is crucial for identifying driver mutations and understanding the mutational mechanisms that act during cancer development. The neutral mutational process is very complex: whole-genome analyses have revealed that the mutation rate differs between cancer types, between patients and along the genome depending on the genetic and epigenetic context. Therefore, methods that predict the number of different types of mutations in regions or specific genomic elements must consider local genomic explanatory variables. A major drawback of most methods is the need to average the explanatory variables across the entire region or genomic element. This procedure is particularly problematic if the explanatory variable varies dramatically in the element under consideration.
Results: Totakeintoaccountthefinescaleoftheexplanatoryvariables,wemodeltheprobabilitiesofdifferenttypes of mutations for each position in the genome by multinomial logistic regression. We analyse 505 cancer genomes from 14 different cancer types and compare the performance in predicting mutation rate for both regional based models and site-specific models. We show that for 1000 randomly selected genomic positions, the site-specific model predicts the mutation rate much better than regional based models.
We use a forward selection procedure to identify the most important explanatory variables. The procedure identifies site-specific conservation (phyloP), replication timing, and expression level as the best predictors for the mutation rate. Finally, our model confirms and quantifies certain well-known mutational signatures.
Conclusion: Wefindthatoursite-specificmultinomialregressionmodeloutperformstheregionalbasedmodels. The possibility of including genomic variables on different scales and patient specific variables makes it a versatile framework for studying different mutational mechanisms. Our model can serve as the neutral null model for the mutational process; regions that deviate from the null model are candidates for elements that drive cancer development.
Keywords: Multinomiallogisticregression,Site-specificmodel,Somaticcancermutations
Original languageEnglish
JournalBMC Bioinformatics
Volume19
Issue number147
Number of pages15
ISSN1471-2105
DOIs
StatePublished - 19 Apr 2018

See relations at Aarhus University Citationformats

ID: 128633465