{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,12,31]],"date-time":"2025-12-31T10:26:37Z","timestamp":1767176797519,"version":"build-2238731810"},"update-to":[{"DOI":"10.1371\/journal.pcbi.1013241","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2025,7,16]],"date-time":"2025-07-16T00:00:00Z","timestamp":1752624000000}}],"reference-count":65,"publisher":"Public Library of Science (PLoS)","issue":"7","license":[{"start":{"date-parts":[[2025,7,7]],"date-time":"2025-07-07T00:00:00Z","timestamp":1751846400000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/100010661","name":"Horizon 2020 Framework Programme","doi-asserted-by":"publisher","award":["955708"],"award-info":[{"award-number":["955708"]}],"id":[{"id":"10.13039\/100010661","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":["www.ploscompbiol.org"],"crossmark-restriction":false},"short-container-title":["PLoS Comput Biol"],"abstract":"<jats:p>\n                    Genetic mutations are footprints of cancer evolution and reveal critical dynamic parameters of tumour growth, which otherwise are hard to measure\n                    <jats:italic>in vivo<\/jats:italic>\n                    . The mutation accumulation in tumour cell populations has been described by various statistics, such as site frequency spectra (SFS), single-cell division distributions (DD) and mutational burden distributions (MBD). While DD and SFS have been intensively studied in phylogenetics especially after the development of whole genome sequencing technology of bulk samples, MBD has drawn attention more recently with the single-cell sequencing data. Although those statistics all arise from the same somatic evolutionary process, an integrated understanding of these distributions is missing and requires novel mathematical tools to better inform the ecological and evolutionary dynamics of tumours. Here we introduce dynamical matrices to analyse and unite the SFS, DD and MBD and derive recurrence relations for the expectations of these three distributions. While we successfully recover classic exact results in pure-birth cases for the SFS and the DD through our new framework, we derive a new expression for the MBD and approximate all three distributions when death is introduced. We demonstrate a natural link between the SFS and the single-cell MBD, and show that the MBD can be regenerated through the DD. Counter-intuitively, the single-cell MBD is mainly driven by the stochasticity arising in the DD, rather than the extra stochasticity in the number of mutations at each cell division.\n                  <\/jats:p>","DOI":"10.1371\/journal.pcbi.1013241","type":"journal-article","created":{"date-parts":[[2025,7,7]],"date-time":"2025-07-07T14:06:41Z","timestamp":1751897201000},"page":"e1013241","update-policy":"https:\/\/doi.org\/10.1371\/journal.pcbi.corrections_policy","source":"Crossref","is-referenced-by-count":0,"title":["Single-cell mutational burden distributions in birth\u2013death processes"],"prefix":"10.1371","volume":"21","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-9350-7833","authenticated-orcid":true,"given":"Christo","family":"Morison","sequence":"first","affiliation":[]},{"ORCID":"https:\/\/orcid.org\/0000-0001-8527-3996","authenticated-orcid":true,"given":"Dudley","family":"Stark","sequence":"additional","affiliation":[]},{"ORCID":"https:\/\/orcid.org\/0000-0002-9016-2665","authenticated-orcid":true,"given":"Weini","family":"Huang","sequence":"additional","affiliation":[]}],"member":"340","published-online":{"date-parts":[[2025,7,7]]},"reference":[{"key":"pcbi.1013241.ref001","doi-asserted-by":"crossref","DOI":"10.1201\/9780429258794","volume-title":"The biology of cancer","author":"RA Weinberg","year":"2013"},{"issue":"12","key":"pcbi.1013241.ref002","doi-asserted-by":"crossref","first-page":"1083","DOI":"10.1016\/j.tree.2021.08.011","article-title":"Evolution via somatic genetic variation in modular species","volume":"36","author":"TBH Reusch","year":"2021","journal-title":"Trends Ecol Evol"},{"issue":"6375","key":"pcbi.1013241.ref003","doi-asserted-by":"crossref","first-page":"550","DOI":"10.1126\/science.aan8690","article-title":"Different mutational rates and mechanisms in human cells at pregastrulation and neurogenesis","volume":"359","author":"T Bae","year":"2018","journal-title":"Science"},{"issue":"7724","key":"pcbi.1013241.ref004","doi-asserted-by":"crossref","first-page":"473","DOI":"10.1038\/s41586-018-0497-0","article-title":"Population dynamics of normal human blood inferred from somatic mutations","volume":"561","author":"H Lee-Six","year":"2018","journal-title":"Nature"},{"issue":"1","key":"pcbi.1013241.ref005","doi-asserted-by":"crossref","first-page":"1035","DOI":"10.1038\/s41467-020-14844-6","article-title":"Measuring single cell divisions in human tissues from multi-region sequencing data","volume":"11","author":"B Werner","year":"2020","journal-title":"Nat Commun"},{"issue":"3","key":"pcbi.1013241.ref006","doi-asserted-by":"crossref","first-page":"291","DOI":"10.1002\/bies.20000","article-title":"Problems of somatic mutation and cancer","volume":"26","author":"SA Frank","year":"2004","journal-title":"Bioessays"},{"issue":"6","key":"pcbi.1013241.ref007","doi-asserted-by":"crossref","first-page":"494","DOI":"10.1016\/j.semcancer.2005.07.003","article-title":"Cancer, aging and the optimal tissue design","volume":"15","author":"NL Komarova","year":"2005","journal-title":"Semin Cancer Biol"},{"issue":"7467","key":"pcbi.1013241.ref008","doi-asserted-by":"crossref","first-page":"338","DOI":"10.1038\/nature12625","article-title":"The causes and consequences of genetic heterogeneity in cancer evolution","volume":"501","author":"RA Burrell","year":"2013","journal-title":"Nature"},{"issue":"43","key":"pcbi.1013241.ref009","doi-asserted-by":"crossref","first-page":"18545","DOI":"10.1073\/pnas.1010978107","article-title":"Accumulation of driver and passenger mutations during tumor progression","volume":"107","author":"I Bozic","year":"2010","journal-title":"Proc Natl Acad Sci U S A"},{"issue":"47","key":"pcbi.1013241.ref010","doi-asserted-by":"crossref","DOI":"10.1073\/pnas.1519556112","article-title":"Extremely high genetic diversity in a single tumor points to prevalence of non-Darwinian cell evolution","volume":"112","author":"S Ling","year":"2015","journal-title":"Proc Natl Acad Sci U S A"},{"issue":"3","key":"pcbi.1013241.ref011","doi-asserted-by":"crossref","first-page":"209","DOI":"10.1038\/ng.3214","article-title":"A Big Bang model of human colorectal tumor growth","volume":"47","author":"A Sottoriva","year":"2015","journal-title":"Nat Genet"},{"issue":"3","key":"pcbi.1013241.ref012","doi-asserted-by":"crossref","first-page":"238","DOI":"10.1038\/ng.3489","article-title":"Identification of neutral tumor evolution across cancer types","volume":"48","author":"MJ Williams","year":"2016","journal-title":"Nat Genet"},{"issue":"46","key":"pcbi.1013241.ref013","doi-asserted-by":"crossref","first-page":"18135","DOI":"10.1073\/pnas.0705778104","article-title":"Clonal interference in large populations","volume":"104","author":"S-C Park","year":"2007","journal-title":"Proc Natl Acad Sci U S A"},{"issue":"7964","key":"pcbi.1013241.ref014","doi-asserted-by":"crossref","first-page":"383","DOI":"10.1038\/s41586-023-06102-8","article-title":"Deterministic evolution and stringent selection during preneoplasia","volume":"618","author":"K Karlsson","year":"2023","journal-title":"Nature"},{"issue":"1","key":"pcbi.1013241.ref015","doi-asserted-by":"crossref","first-page":"141","DOI":"10.2307\/3212065","article-title":"The pseudo-transient distribution and its uses in genetics","volume":"1","author":"WJ Ewens","year":"1964","journal-title":"J Appl Prob"},{"issue":"3","key":"pcbi.1013241.ref016","doi-asserted-by":"crossref","first-page":"247","DOI":"10.1017\/S0016672300011459","article-title":"Genetic variability maintained in a finite population due to mutational production of neutral and nearly neutral isoalleles","volume":"11","author":"M Kimura","year":"1968","journal-title":"Genet Res"},{"issue":"2","key":"pcbi.1013241.ref017","doi-asserted-by":"crossref","first-page":"172","DOI":"10.1006\/tpbi.1995.1025","article-title":"Statistical properties of segregating sites","volume":"48","author":"YX Fu","year":"1995","journal-title":"Theor Popul Biol"},{"issue":"2","key":"pcbi.1013241.ref018","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pcbi.1004731","article-title":"Quantifying clonal and subclonal passenger mutations in cancer evolution","volume":"12","author":"I Bozic","year":"2016","journal-title":"PLoS Comput Biol"},{"key":"pcbi.1013241.ref019","doi-asserted-by":"crossref","DOI":"10.7554\/eLife.89780","article-title":"Measures of genetic diversification in somatic tissues at bulk and single-cell resolution","volume":"12","author":"ME Moeller","year":"2024","journal-title":"Elife"},{"issue":"1","key":"pcbi.1013241.ref020","doi-asserted-by":"crossref","first-page":"230","DOI":"10.1214\/11-AAP824","article-title":"Population genetics of neutral mutations in exponentially growing cancer cell populations","volume":"23","author":"R Durrett","year":"2013","journal-title":"Ann Appl Probab"},{"key":"pcbi.1013241.ref021","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1007\/978-3-319-16065-8_1","article-title":"Branching process models of cancer.","volume-title":"Branching process models of cancer","author":"R. Durrett","year":"2015"},{"key":"pcbi.1013241.ref022","doi-asserted-by":"crossref","first-page":"67","DOI":"10.1016\/j.tpb.2021.09.004","article-title":"Exact site frequency spectra of neutrally evolving tumors: a transition between power laws reveals a signature of cell viability","volume":"142","author":"EB Gunnarsson","year":"2021","journal-title":"Theor Popul Biol"},{"issue":"9","key":"pcbi.1013241.ref023","doi-asserted-by":"crossref","first-page":"618","DOI":"10.1038\/nrg3542","article-title":"Single-cell sequencing-based technologies will revolutionize whole-organism science","volume":"14","author":"E Shapiro","year":"2013","journal-title":"Nat Rev Genet"},{"issue":"4","key":"pcbi.1013241.ref024","doi-asserted-by":"crossref","first-page":"598","DOI":"10.1016\/j.molcel.2015.05.005","article-title":"Advances and applications of single-cell sequencing technologies","volume":"58","author":"Y Wang","year":"2015","journal-title":"Mol Cell"},{"issue":"7859","key":"pcbi.1013241.ref025","doi-asserted-by":"crossref","first-page":"405","DOI":"10.1038\/s41586-021-03477-4","article-title":"Somatic mutation landscapes at single-molecule resolution","volume":"593","author":"F Abascal","year":"2021","journal-title":"Nature"},{"issue":"1","key":"pcbi.1013241.ref026","doi-asserted-by":"crossref","first-page":"34","DOI":"10.1186\/s13073-017-0424-2","article-title":"Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden","volume":"9","author":"ZR Chalmers","year":"2017","journal-title":"Genome Med"},{"key":"pcbi.1013241.ref027","article-title":"Cancer-specific thresholds adjust for whole exome sequencing-based tumor mutational burden distribution","volume":"3","author":"EM Fernandez","year":"2019","journal-title":"JCO Precis Oncol"},{"issue":"1","key":"pcbi.1013241.ref028","doi-asserted-by":"crossref","first-page":"31","DOI":"10.1038\/s41698-021-00169-0","article-title":"Panels and models for accurate prediction of tumor mutation burden in tumor samples","volume":"5","author":"E Mart\u00ednez-P\u00e9rez","year":"2021","journal-title":"NPJ Precis Oncol"},{"key":"pcbi.1013241.ref029","doi-asserted-by":"crossref","first-page":"044407","DOI":"10.1103\/PhysRevE.109.044407","article-title":"How mutation accumulation depends on the structure of the cell lineage tree","volume":"109","author":"I Der\u00e9nyi","year":"2024","journal-title":"Phys Rev E"},{"key":"pcbi.1013241.ref030","volume-title":"Molecular evolution: a phylogenetic approach","author":"RD Page","year":"2009"},{"issue":"1","key":"pcbi.1013241.ref031","doi-asserted-by":"crossref","first-page":"91","DOI":"10.1016\/S0025-5564(00)00061-4","article-title":"Properties of phylogenetic trees generated by Yule-type speciation models","volume":"170","author":"M Steel","year":"2001","journal-title":"Math Biosci"},{"issue":"4","key":"pcbi.1013241.ref032","doi-asserted-by":"crossref","first-page":"299","DOI":"10.1093\/comjnl\/7.4.299","article-title":"More combinatorial properties of certain trees","volume":"7","author":"WC Lynch","year":"1965","journal-title":"Comput J"},{"issue":"4","key":"pcbi.1013241.ref033","doi-asserted-by":"crossref","first-page":"893","DOI":"10.1093\/genetics\/61.4.893","article-title":"The number of heterozygous nucleotide sites maintained in a finite population due to steady flux of mutations","volume":"61","author":"M Kimura","year":"1969","journal-title":"Genetics"},{"issue":"11","key":"pcbi.1013241.ref034","doi-asserted-by":"crossref","first-page":"1885","DOI":"10.1101\/gr.220707.117","article-title":"Single-cell sequencing data reveal widespread recurrence and loss of mutational hits in the life histories of tumors","volume":"27","author":"J Kuipers","year":"2017","journal-title":"Genome Res"},{"issue":"6","key":"pcbi.1013241.ref035","doi-asserted-by":"crossref","DOI":"10.1214\/18-AAP1413","article-title":"Mutation frequencies in a birth\u2013death branching process","volume":"28","author":"D Cheek","year":"2018","journal-title":"Ann Appl Probab"},{"issue":"11","key":"pcbi.1013241.ref036","doi-asserted-by":"crossref","first-page":"6580","DOI":"10.1016\/j.spa.2020.06.003","article-title":"Genetic composition of an exponentially growing cell population","volume":"130","author":"D Cheek","year":"2020","journal-title":"Stochast Process Appl"},{"issue":"1","key":"pcbi.1013241.ref037","doi-asserted-by":"crossref","first-page":"181","DOI":"10.1534\/genetics.113.152587","article-title":"Learning natural selection from the site frequency spectrum","volume":"195","author":"R Ronen","year":"2013","journal-title":"Genetics"},{"issue":"19","key":"pcbi.1013241.ref038","doi-asserted-by":"crossref","first-page":"4297","DOI":"10.1016\/j.physa.2009.06.032","article-title":"Growing networks with preferential deletion and addition of edges","volume":"388","author":"M Deijfen","year":"2009","journal-title":"Phys A: Statist Mech Appl"},{"key":"pcbi.1013241.ref039","volume-title":"Probability models for DNA sequence evolution","author":"R Durrett"},{"issue":"1","key":"pcbi.1013241.ref040","doi-asserted-by":"crossref","first-page":"128","DOI":"10.1073\/pnas.1516123113","article-title":"Deep sequencing as a probe of normal stem cell fate and preneoplasia in human epidermis","volume":"113","author":"BD Simons","year":"2016","journal-title":"Proc Natl Acad Sci U S A"},{"issue":"52","key":"pcbi.1013241.ref041","doi-asserted-by":"crossref","first-page":"26863","DOI":"10.1073\/pnas.1910301116","article-title":"Extensive subclonal mutational diversity in human colorectal cancer and its significance","volume":"116","author":"LA Loeb","year":"2019","journal-title":"Proc Natl Acad Sci U S A"},{"issue":"6485","key":"pcbi.1013241.ref042","doi-asserted-by":"crossref","first-page":"1449","DOI":"10.1126\/science.aay9333","article-title":"The evolutionary dynamics and fitness landscape of clonal hematopoiesis","volume":"367","author":"CJ Watson","year":"2020","journal-title":"Science"},{"issue":"11","key":"pcbi.1013241.ref043","doi-asserted-by":"crossref","first-page":"1597","DOI":"10.1038\/s41588-021-00957-1","article-title":"Synonymous mutations reveal genome-wide levels of positive selection in healthy tissues","volume":"53","author":"GYP Poon","year":"2021","journal-title":"Nat Genet"},{"issue":"2","key":"pcbi.1013241.ref044","article-title":"Signatures of neutral evolution in exponentially growing tumors: a theoretical perspective","volume":"17","author":"H-R Tung","year":"2021","journal-title":"PLoS Comput Biol"},{"key":"pcbi.1013241.ref045","doi-asserted-by":"crossref","first-page":"889438","DOI":"10.3389\/fevo.2022.889438","article-title":"Modes of selection in tumors as reflected by two mathematical models and site frequency spectra","volume":"10","author":"MK Kurpas","year":"2022","journal-title":"Front Ecol Evol"},{"issue":"8","key":"pcbi.1013241.ref046","doi-asserted-by":"crossref","first-page":"8505","DOI":"10.3934\/mbe.2022395","article-title":"Comparison of cell state models derived from single-cell RNA sequencing data: graph versus multi-dimensional space","volume":"19","author":"H Cho","year":"2022","journal-title":"Math Biosci Eng"},{"issue":"1","key":"pcbi.1013241.ref047","doi-asserted-by":"crossref","DOI":"10.1214\/24-AAP2110","article-title":"Asymptotics for the site frequency spectrum associated with the genealogy of a birth and death process","volume":"35","author":"J Schweinsberg","year":"2025","journal-title":"Ann Appl Probab"},{"issue":"1","key":"pcbi.1013241.ref048","article-title":"Statistical Inference for the Evolutionary History of Cancer Genomes","volume":"35","author":"KN Dinh","year":"2020","journal-title":"Statist Sci"},{"issue":"4","key":"pcbi.1013241.ref049","doi-asserted-by":"crossref","DOI":"10.1214\/105051604000000486","article-title":"Asymptotic genealogy of a critical branching process","volume":"14","author":"L Popovic","year":"2004","journal-title":"Ann Appl Probab"},{"issue":"3","key":"pcbi.1013241.ref050","first-page":"359","article-title":"The allelic partition for coalescent point processes","volume":"15","author":"A Lambert","year":"2009","journal-title":"Markov Process Relat Fields"},{"key":"pcbi.1013241.ref051","doi-asserted-by":"crossref","first-page":"113","DOI":"10.1016\/j.tpb.2013.10.002","article-title":"Birth-death models and coalescent point processes: the shape and probability of reconstructed phylogenies","volume":"90","author":"A Lambert","year":"2013","journal-title":"Theor Popul Biol"},{"key":"pcbi.1013241.ref052","doi-asserted-by":"crossref","first-page":"30","DOI":"10.1016\/j.tpb.2018.04.005","article-title":"The coalescent of a sample from a binary branching process","volume":"122","author":"A Lambert","year":"2018","journal-title":"Theor Popul Biol"},{"issue":"9","key":"pcbi.1013241.ref053","doi-asserted-by":"crossref","DOI":"10.1093\/bioinformatics\/btad561","article-title":"cloneRate: fast estimation of single-cell clonal dynamics using coalescent theory","volume":"39","author":"B Johnson","year":"2023","journal-title":"Bioinformatics"},{"issue":"6","key":"pcbi.1013241.ref054","doi-asserted-by":"crossref","first-page":"895","DOI":"10.1038\/s41588-018-0128-6","article-title":"Quantification of subclonal selection in cancer from bulk sequencing data","volume":"50","author":"MJ Williams","year":"2018","journal-title":"Nat Genet"},{"issue":"7702","key":"pcbi.1013241.ref055","doi-asserted-by":"crossref","first-page":"457","DOI":"10.1038\/s41586-018-0024-3","article-title":"Intra-tumour diversification in colorectal cancer at the single-cell level","volume":"556","author":"SF Roerink","year":"2018","journal-title":"Nature"},{"issue":"1","key":"pcbi.1013241.ref056","doi-asserted-by":"crossref","first-page":"31","DOI":"10.1186\/s13059-020-1926-6","article-title":"Eleven grand challenges in single-cell data science","volume":"21","author":"D L\u00e4hnemann","year":"2020","journal-title":"Genome Biol"},{"key":"pcbi.1013241.ref057","doi-asserted-by":"crossref","first-page":"2978","DOI":"10.1016\/j.csbj.2022.06.013","article-title":"Somatic variant calling from single-cell DNA sequencing data","volume":"20","author":"M Valecha","year":"2022","journal-title":"Comput Struct Biotechnol J"},{"issue":"7913","key":"pcbi.1013241.ref058","doi-asserted-by":"crossref","first-page":"343","DOI":"10.1038\/s41586-022-04786-y","article-title":"Clonal dynamics of haematopoiesis across the human lifespan","volume":"606","author":"E Mitchell","year":"2022","journal-title":"Nature"},{"key":"pcbi.1013241.ref059","article-title":"On the patterns of genetic intra-tumour heterogeneity before and after treatment","author":"A Stein","year":"2025","journal-title":"Genetics"},{"issue":"07","key":"pcbi.1013241.ref060","doi-asserted-by":"crossref","DOI":"10.1088\/1742-5468\/2010\/07\/P07028","article-title":"Exact solution of a two-type branching process: clone size distribution in cell division kinetics","volume":"2010","author":"T Antal","year":"2010","journal-title":"J Stat Mech"},{"issue":"08","key":"pcbi.1013241.ref061","doi-asserted-by":"crossref","DOI":"10.1088\/1742-5468\/2011\/08\/P08018","article-title":"Exact solution of a two-type branching process: models of tumor progression","volume":"2011","author":"T Antal","year":"2011","journal-title":"J Stat Mech"},{"issue":"10","key":"pcbi.1013241.ref062","doi-asserted-by":"crossref","first-page":"640","DOI":"10.1038\/nchembio.2007.38","article-title":"Telomere length, stem cells and aging","volume":"3","author":"MA Blasco","year":"2007","journal-title":"Nat Chem Biol"},{"key":"pcbi.1013241.ref063","doi-asserted-by":"crossref","DOI":"10.3389\/fgene.2019.00008","article-title":"The impact of heterogeneity on single-cell sequencing","volume":"10","author":"SL Goldman","year":"2019","journal-title":"Front Genet"},{"issue":"4","key":"pcbi.1013241.ref064","doi-asserted-by":"crossref","first-page":"456","DOI":"10.1016\/j.ccell.2020.03.008","article-title":"Advancing cancer research and medicine with single-cell genomics","volume":"37","author":"B Lim","year":"2020","journal-title":"Cancer Cell"},{"issue":"4","key":"pcbi.1013241.ref065","doi-asserted-by":"crossref","first-page":"403","DOI":"10.1016\/0021-9991(76)90041-3","article-title":"A general method for numerically simulating the stochastic time evolution of coupled chemical reactions","volume":"22","author":"DT Gillespie","year":"1976","journal-title":"J Comput Phys"}],"updated-by":[{"DOI":"10.1371\/journal.pcbi.1013241","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2025,7,16]],"date-time":"2025-07-16T00:00:00Z","timestamp":1752624000000}}],"container-title":["PLOS Computational Biology"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1013241","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,7,16]],"date-time":"2025-07-16T17:46:14Z","timestamp":1752687974000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1013241"}},"subtitle":[],"editor":[{"given":"Ivana","family":"Bozic","sequence":"first","affiliation":[]}],"short-title":[],"issued":{"date-parts":[[2025,7,7]]},"references-count":65,"journal-issue":{"issue":"7","published-online":{"date-parts":[[2025,7,7]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pcbi.1013241","relation":{},"ISSN":["1553-7358"],"issn-type":[{"value":"1553-7358","type":"electronic"}],"subject":[],"published":{"date-parts":[[2025,7,7]]}}}