by Gergely Tibély, Dominik Schrempf, Imre Derényi, Gergely J. Szöllősi
Tumors often harbor orders of magnitude more mutations than healthy tissues. The increased number of mutations may be due to an elevated mutation rate or frequent cell death and correspondingly rapid cell turnover, or a combination of the two. It is difficult to disentangle these two mechanisms based on widely available bulk sequencing data, where sequences from individual cells are intermixed and, thus, the cell lineage tree of the tumor cannot be resolved. Here we present a method that can simultaneously estimate the cell turnover rate and the rate of mutations from bulk sequencing data. Our method works by simulating tumor growth and finding the parameters with which the observed data can be reproduced with maximum likelihood. Applying this method to a real tumor sample, we find that both the mutation rate and the frequency of death may be high.
Unless otherwise indicated, articles and accompanying materials published by PLOS on the PLOS Sites, including peer reviews, are licensed by the respective authors for use and distribution by you subject to citation of the original source in accordance with the Creative Commons Attribution (CC BY) license.