{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,8,4]],"date-time":"2026-08-04T23:50:20Z","timestamp":1785887420044,"version":"3.56.0"},"reference-count":40,"publisher":"Oxford University Press (OUP)","issue":"8","license":[{"start":{"date-parts":[[2026,7,24]],"date-time":"2026-07-24T00:00:00Z","timestamp":1784851200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100001659","name":"German Research Foundation","doi-asserted-by":"publisher","id":[{"id":"10.13039\/501100001659","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100001711","name":"Swiss National Science Foundation","doi-asserted-by":"publisher","award":["179518"],"award-info":[{"award-number":["179518"]}],"id":[{"id":"10.13039\/501100001711","id-type":"DOI","asserted-by":"publisher"}]},{"name":"Swiss Cancer League","award":["KFS-2977-08-2012"],"award-info":[{"award-number":["KFS-2977-08-2012"]}]},{"name":"Free State of Bavaria [Marianne-Plehn-Program"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2026,8,3]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:sec>\n                    <jats:title>Motivation<\/jats:title>\n                    <jats:p>Cancer progresses through the accumulation of genomic events. Cancer progression models such as Mutual Hazard Networks (MHNs) describe this dynamic, enabling prediction of temporal event positions and patient-specific risks of acquiring mutations. However, current MHN analyses rely on single most likely models and do not quantify the uncertainty inherent to parameter estimation. Assessing forecast stability is essential before using them to anticipate treatment-relevant mutations, adapt targeted therapies, or prioritize monitoring of patients at elevated progression risk.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>We address a key prerequisite for the responsible clinical use of cancer progression models by making MHN-derived predictions uncertainty-aware. We present a Bayesian framework for MHN that uses Markov Chain Monte Carlo to sample from the posterior distributions of model parameters and derived predictions. For practical use we implemented the Random-Walk Metropolis, Metropolis-Adjusted Langevin Algorithm (MALA), and simplified manifold MALA samplers as part of the existing mhn Python package. Only MALA and smMALA were successful in sampling from MHN posteriors, with MALA performing best. While most MHN parameters and predictions showed low posterior variance, a small subset displayed greater variability across the posterior distribution. This differentiation cannot be obtained from a single most likely model, emphasizing the need for uncertainty quantification, especially in clinical contexts. As an illustrative example, posterior sampling identified a subgroup of STK11$-$, KRAS$+$ lung adenocarcinoma patients with a high predicted short-term risk\u2014with low variance across posterior samples\u2014to develop an STK11 mutation. This subgroup exhibited poorer survival under immunotherapy, resembling patterns observed in STK11+ patients.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Availability and implementation<\/jats:title>\n                    <jats:p>Our implementation is part of version 1.2.0 of the mhn package (https:\/\/github.com\/spang-lab\/LearnMHN). All analyses including the code to produce all figures in this article can be found under https:\/\/github.com\/huy29433\/MCMC-sampling-for-MHN (https:\/\/doi.org\/10.5281\/zenodo.21160219).<\/jats:p>\n                  <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btag526","type":"journal-article","created":{"date-parts":[[2026,7,24]],"date-time":"2026-07-24T11:41:55Z","timestamp":1784893315000},"source":"Crossref","is-referenced-by-count":0,"title":["Quantifying uncertainty of predictions from cancer progression models"],"prefix":"10.1093","volume":"42","author":[{"ORCID":"https:\/\/orcid.org\/0009-0002-7263-3635","authenticated-orcid":false,"given":"Yanren Linda","family":"Hu","sequence":"first","affiliation":[{"name":"Department for Statistical Bioinformatics, University of Regensburg , Regensburg 93053,","place":["Germany"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0009-0001-7364-4005","authenticated-orcid":false,"given":"Simon","family":"Pfahler","sequence":"additional","affiliation":[{"name":"Department of Physics, University of Regensburg , Regensburg 93040,","place":["Germany"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0009-0000-9100-2367","authenticated-orcid":false,"given":"Andreas","family":"L\u00f6sch","sequence":"additional","affiliation":[{"name":"Department for Statistical Bioinformatics, University of Regensburg , Regensburg 93053,","place":["Germany"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0009-0007-9077-1376","authenticated-orcid":false,"given":"Stefan","family":"Vocht","sequence":"additional","affiliation":[{"name":"Department for Statistical Bioinformatics, University of Regensburg , Regensburg 93053,","place":["Germany"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-8458-4829","authenticated-orcid":false,"given":"Stefan","family":"Hansch","sequence":"additional","affiliation":[{"name":"Department for Statistical Bioinformatics, University of Regensburg , Regensburg 93053,","place":["Germany"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-4639-4717","authenticated-orcid":false,"given":"Kevin","family":"Rupp","sequence":"additional","affiliation":[{"name":"Department of Biosystems Science and Engineering, ETH Z\u00fcrich , Basel 4056,","place":["Switzerland"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-0573-6119","authenticated-orcid":false,"given":"Niko","family":"Beerenwinkel","sequence":"additional","affiliation":[{"name":"Department of Biosystems Science and Engineering, ETH Z\u00fcrich , Basel 4056,","place":["Switzerland"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-6732-9204","authenticated-orcid":false,"given":"Tilo","family":"Wettig","sequence":"additional","affiliation":[{"name":"Department of Physics, University of Regensburg , Regensburg 93040,","place":["Germany"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-4649-7190","authenticated-orcid":false,"given":"Rudolf","family":"Schill","sequence":"additional","affiliation":[{"name":"Department of Biosystems Science and Engineering, ETH Z\u00fcrich , Basel 4056,","place":["Switzerland"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-1326-4297","authenticated-orcid":false,"given":"Rainer","family":"Spang","sequence":"additional","affiliation":[{"name":"Department for Statistical Bioinformatics, University of Regensburg , Regensburg 93053,","place":["Germany"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2026,7,24]]},"reference":[{"key":"2026080419211954200_btag526-B1","author":"Abadi","year":"2015"},{"key":"2026080419211954200_btag526-B2","doi-asserted-by":"publisher","first-page":"e1012393","DOI":"10.1371\/journal.pcbi.1012393","article-title":"HyperTraPS-CT: inference and prediction for accumulation pathways with flexible data and model structures","volume":"20","author":"Aga","year":"2024","journal-title":"PLoS Comput Biol"},{"key":"2026080419211954200_btag526-B3","doi-asserted-by":"publisher","first-page":"754","DOI":"10.1093\/bioinformatics\/btab717","article-title":"PMCE: efficient inference of expressive models of cancer evolution with high prognostic power","volume":"38","author":"Angaroni","year":"2022","journal-title":"Bioinformatics"},{"key":"2026080419211954200_btag526-B4","doi-asserted-by":"publisher","first-page":"e0283004","DOI":"10.1371\/journal.pone.0283004","article-title":"Timed hazard networks: incorporating temporal difference for oncogenetic analysis","volume":"18","author":"Chen","year":"2023","journal-title":"PLoS One"},{"key":"2026080419211954200_btag526-B5","doi-asserted-by":"publisher","first-page":"23","DOI":"10.5334\/jors.252","article-title":"Probabilistic inference on noisy time series (PINTS)","volume":"7","author":"Clerx","year":"2019","journal-title":"JORS"},{"key":"2026080419211954200_btag526-B6","doi-asserted-by":"publisher","first-page":"e1009055","DOI":"10.1371\/journal.pcbi.1009055","article-title":"Conditional prediction of consecutive tumor evolution using cancer progression models: what genotype comes next?","volume":"17","author":"Diaz-Colunga","year":"2021","journal-title":"PLoS Comput Biol"},{"key":"2026080419211954200_btag526-B7","doi-asserted-by":"publisher","first-page":"823","DOI":"10.1016\/j.trecan.2021.04.009","article-title":"Co-occurrence and mutual exclusivity: what cross-cancer mutation patterns can tell us","volume":"7","author":"El Tekle","year":"2021","journal-title":"Trends Cancer"},{"key":"2026080419211954200_btag526-B8","doi-asserted-by":"publisher","first-page":"457","DOI":"10.1214\/ss\/1177011136","article-title":"Inference from iterative simulation using multiple sequences","volume":"7","author":"Gelman","year":"1992","journal-title":"Statist Sci"},{"key":"2026080419211954200_btag526-B9","doi-asserted-by":"publisher","first-page":"2809","DOI":"10.1093\/bioinformatics\/btp505","article-title":"Quantifying cancer progression with conjunctive Bayesian networks","volume":"25","author":"Gerstung","year":"2009","journal-title":"Bioinformatics"},{"key":"2026080419211954200_btag526-B10","doi-asserted-by":"publisher","first-page":"e27136","DOI":"10.1371\/journal.pone.0027136","article-title":"The temporal order of genetic and pathway alterations in tumorigenesis","volume":"6","author":"Gerstung","year":"2011","journal-title":"PLoS One"},{"key":"2026080419211954200_btag526-B11","doi-asserted-by":"publisher","first-page":"2340","DOI":"10.1021\/j100540a008","article-title":"Exact stochastic simulation of coupled chemical reactions","volume":"81","author":"Gillespie","year":"1977","journal-title":"J Phys Chem"},{"key":"2026080419211954200_btag526-B12","doi-asserted-by":"publisher","first-page":"123","DOI":"10.1111\/j.1467-9868.2010.00765.x","article-title":"Riemann manifold Langevin and Hamiltonian Monte Carlo methods","volume":"73","author":"Girolami","year":"2011","journal-title":"J R Stat Soc Ser B Stat Methodol"},{"key":"2026080419211954200_btag526-B13","doi-asserted-by":"publisher","first-page":"39","DOI":"10.1016\/j.cels.2019.10.009","article-title":"Hypertraps: inferring probabilistic patterns of trait acquisition in evolutionary and disease progression pathways","volume":"10","author":"Greenbury","year":"2020","journal-title":"Cell Syst"},{"key":"2026080419211954200_btag526-B14","doi-asserted-by":"publisher","DOI":"10.1007\/978-0-387-84858-7","volume-title":"The Elements of Statistical Learning. Springer Series in Statistics","author":"Hastie","year":"2009"},{"key":"2026080419211954200_btag526-B15","doi-asserted-by":"publisher","first-page":"853","DOI":"10.1089\/cmb.2006.13.853","article-title":"New probabilistic network models and algorithms for oncogenesis","volume":"13","author":"Hjelm","year":"2006","journal-title":"J Comput Biol"},{"key":"2026080419211954200_btag526-B16","doi-asserted-by":"publisher","first-page":"111272","DOI":"10.1016\/j.celrep.2022.111272","article-title":"Pervasive conditional selection of driver mutations and modular epistasis networks in cancer","volume":"40","author":"Iranzo","year":"2022","journal-title":"Cell Rep"},{"key":"2026080419211954200_btag526-B17","doi-asserted-by":"publisher","first-page":"776","DOI":"10.7150\/jca.41296","article-title":"Pan-cancer analysis of ARID1A alterations as biomarkers for immunotherapy outcomes","volume":"11","author":"Jiang","year":"2020","journal-title":"J Cancer"},{"key":"2026080419211954200_btag526-B18","doi-asserted-by":"crossref","first-page":"63","DOI":"10.1038\/s41746-019-0140-y","article-title":"Precision identification of high-risk phenotypes and progression pathways in severe malaria without requiring longitudinal data","volume":"2","author":"Johnston","year":"2019","journal-title":"NPJ Digit Med"},{"key":"2026080419211954200_btag526-B19","doi-asserted-by":"publisher","first-page":"101","DOI":"10.1016\/j.cels.2016.01.013","article-title":"Evolutionary inference across eukaryotes identifies specific pressures favoring mitochondrial gene retention","volume":"2","author":"Johnston","year":"2016","journal-title":"Cell Syst"},{"key":"2026080419211954200_btag526-B20","doi-asserted-by":"publisher","first-page":"1143","DOI":"10.21105\/joss.01143","article-title":"Arviz a unified library for exploratory analysis of Bayesian models in python","volume":"4","author":"Kumar","year":"2019","journal-title":"JOSS"},{"key":"2026080419211954200_btag526-B21","doi-asserted-by":"publisher","first-page":"3676","DOI":"10.1038\/s41467-023-39400-w","article-title":"Joint inference of exclusivity patterns and recurrent trajectories from tumor","volume":"14","author":"Luo","year":"2023","journal-title":"Nat Commun"},{"key":"2026080419211954200_btag526-B22","doi-asserted-by":"publisher","first-page":"551","DOI":"10.1007\/s11673-019-09944-y","article-title":"Tissue vs liquid biopsies for cancer detection: ethical issues","volume":"16","author":"Mannelli","year":"2019","journal-title":"J Bioeth Inq"},{"key":"2026080419211954200_btag526-B23","doi-asserted-by":"publisher","first-page":"btac803","DOI":"10.1093\/bioinformatics\/btac803","article-title":"Hyperhmm: efficient inference of evolutionary and progressive dynamics on hypercubic transition graphs","volume":"39","author":"Moen","year":"2023","journal-title":"Bioinformatics"},{"key":"2026080419211954200_btag526-B24","doi-asserted-by":"publisher","first-page":"i727","DOI":"10.1093\/bioinformatics\/btw459","article-title":"Large-scale inference of conjunctive Bayesian networks","volume":"32","author":"Montazeri","year":"2016","journal-title":"Bioinformatics"},{"key":"2026080419211954200_btag526-B25","doi-asserted-by":"publisher","first-page":"563","DOI":"10.1016\/j.cell.2022.01.003","article-title":"Genomic characterization of metastatic patterns from prospective clinical sequencing of 25,000 patients","volume":"185","author":"Nguyen","year":"2022","journal-title":"Cell"},{"key":"2026080419211954200_btag526-B26","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1136\/jitc-2019-000438","article-title":"ARID1A alterations function as a biomarker for longer progression-free survival after anti-pd-1\/pd-l1 immunotherapy","volume":"8","author":"Okamura","year":"2020","journal-title":"J Immunother Cancer"},{"key":"2026080419211954200_btag526-B27","doi-asserted-by":"publisher","author":"Pfahler","year":"2026","DOI":"10.64898\/2026.01.20.700556"},{"key":"2026080419211954200_btag526-B28","doi-asserted-by":"publisher","first-page":"399","DOI":"10.1016\/j.jtho.2021.10.013","article-title":"Diminished efficacy of programmed death-(ligand)1 inhibition in STK11- and KEAP1-mutant lung adenocarcinoma is affected by KRAS mutation status","volume":"17","author":"Ricciuti","year":"2022","journal-title":"J Thorac Oncol"},{"key":"2026080419211954200_btag526-B29","doi-asserted-by":"publisher","first-page":"11","DOI":"10.1214\/ss\/1015346320","article-title":"Optimal scaling for various metropolis-hastings algorithms","volume":"16","author":"Roberts","year":"2001","journal-title":"Statist Sci"},{"key":"2026080419211954200_btag526-B30","doi-asserted-by":"publisher","first-page":"255","DOI":"10.1111\/1467-9868.00123","article-title":"Optimal scaling of discrete approximations to Langevin diffusions","volume":"60","author":"Roberts","year":"1998","journal-title":"J R Stat Soc Ser B Stat Methodol"},{"key":"2026080419211954200_btag526-B31","doi-asserted-by":"publisher","first-page":"i140","DOI":"10.1093\/bioinformatics\/btae250","article-title":"Modeling metastatic progression from cross-sectional cancer genomics data","volume":"40","author":"Rupp","year":"2024","journal-title":"Bioinformatics"},{"key":"2026080419211954200_btag526-B32","doi-asserted-by":"publisher","first-page":"927","DOI":"10.1089\/cmb.2024.0666","article-title":"Correcting for observation bias in cancer progression modeling","volume":"31","author":"Schill","year":"2024","journal-title":"J Comput Biol"},{"key":"2026080419211954200_btag526-B33","doi-asserted-by":"publisher","first-page":"241","DOI":"10.1093\/bioinformatics\/btz513","article-title":"Modelling cancer progression using mutual hazard networks","volume":"36","author":"Schill","year":"2019","journal-title":"Bioinformatics"},{"key":"2026080419211954200_btag526-B34","doi-asserted-by":"publisher","first-page":"822","DOI":"10.1158\/2159-8290.CD-18-0099","article-title":"STK11\/LKB1 mutations and PD-1 inhibitor resistance in KRAS-mutant lung adenocarcinoma","volume":"8","author":"Skoulidis","year":"2018","journal-title":"Cancer Discov"},{"key":"2026080419211954200_btag526-B35","doi-asserted-by":"publisher","first-page":"818","DOI":"10.1158\/2159-8290.CD-17-0151","article-title":"Aacr project genie: powering precision medicine through an international consortium","volume":"7","author":"Sweeney","year":"2017","journal-title":"Cancer Discov"},{"key":"2026080419211954200_btag526-B36","doi-asserted-by":"publisher","first-page":"1701","DOI":"10.1214\/aos\/1176325750","article-title":"Markov chains for exploring posterior distributions","volume":"22","author":"Tierney","year":"1994","journal-title":"Ann Statist"},{"key":"2026080419211954200_btag526-B37","doi-asserted-by":"publisher","first-page":"119","DOI":"10.1016\/j.ejca.2020.09.006","article-title":"The impact of arid1a mutation on molecular characteristics in colorectal cancer","volume":"140","author":"Tokunaga","year":"2020","journal-title":"Eur J Cancer"},{"key":"2026080419211954200_btag526-B38","doi-asserted-by":"publisher","first-page":"667","DOI":"10.1214\/20-BA1221","article-title":"Rank-normalization, folding, and localization: an improved R^ for assessing convergence of MCMC (with discussion)","volume":"16","author":"Vehtari","year":"2021","journal-title":"Bayesian Anal"},{"key":"2026080419211954200_btag526-B39","doi-asserted-by":"publisher","first-page":"vbaf283","DOI":"10.1093\/bioadv\/vbaf283","article-title":"mhn: a python package for analyzing cancer progression with mutual hazard networks","volume":"6","author":"Vocht","year":"2026","journal-title":"Bioinform Adv"},{"key":"2026080419211954200_btag526-B40","doi-asserted-by":"publisher","first-page":"1546","DOI":"10.1126\/science.1235122","article-title":"Cancer genome landscapes","volume":"339","author":"Vogelstein","year":"2013","journal-title":"Science"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btag526\/69950952\/btag526.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/42\/8\/btag526\/69950952\/btag526.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/42\/8\/btag526\/69950952\/btag526.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,8,4]],"date-time":"2026-08-04T23:21:30Z","timestamp":1785885690000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btag526\/8741438"}},"subtitle":[],"editor":[{"given":"Anthony","family":"Mathelier","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2026,7,24]]},"references-count":40,"journal-issue":{"issue":"8","published-print":{"date-parts":[[2026,8,3]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btag526","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2026,8]]},"published":{"date-parts":[[2026,7,24]]},"article-number":"btag526"}}