{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,23]],"date-time":"2026-09-23T21:54:40Z","timestamp":1790200480286,"version":"4.1.0"},"reference-count":55,"publisher":"Oxford University Press (OUP)","issue":"2","license":[{"start":{"date-parts":[[2026,2,2]],"date-time":"2026-02-02T00:00:00Z","timestamp":1769990400000},"content-version":"vor","delay-in-days":1,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100001711","name":"Swiss National Science Foundation","doi-asserted-by":"publisher","award":["179170"],"award-info":[{"award-number":["179170"]}],"id":[{"id":"10.13039\/501100001711","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2026,2,2]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:p>Phylodynamic methods are widely used to infer the population dynamics of viruses between and within hosts. For HIV-1, these methods have been used to estimate migration rates between different anatomical compartments within a host. These methods typically assume that the genomic regions used for reconstruction are evolving without selective pressure, even though other parts of the viral genome are known to experience strong selection. In this study, we investigate how selection affects phylodynamic migration rate estimates. To this end, we developed a novel agent-based simulation tool, virolution, to simulate the evolution of virus within two anatomical compartments of a host. Using this tool, we generated viral sequences and genealogies assuming both, neutral evolution and selection governed by an empirically-supported distribution of fitness effects that is concordant in both compartments. We found that, under the selection regime, migration rates are significantly overestimated with a stochastic mixture model and a structured coalescent model in the Bayesian inference framework BEAST2. Our results reveal that commonly used phylogeographic methods, which assume neutral evolution, can significantly bias migration rate estimates in selective regimes. This study underscores the need for assessing the robustness of phylodynamic analysis with respect to more realistic selection regimes.<\/jats:p>","DOI":"10.1093\/molbev\/msag014","type":"journal-article","created":{"date-parts":[[2026,1,12]],"date-time":"2026-01-12T12:21:35Z","timestamp":1768220495000},"source":"Crossref","is-referenced-by-count":2,"title":["Viral Simulation Reveals Overestimation Bias in Within-Host Phylodynamic Migration Rate Estimates Under Selection"],"prefix":"10.1093","volume":"43","author":[{"ORCID":"https:\/\/orcid.org\/0009-0009-0796-1502","authenticated-orcid":false,"given":"Nicolas","family":"Ochsner","sequence":"first","affiliation":[{"name":"Institute of Integrative Biology, ETH Zurich , Zurich ,","place":["Switzerland"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-3781-0675","authenticated-orcid":false,"given":"Judith","family":"Bouman","sequence":"additional","affiliation":[{"name":"Institute of Social and Preventive Medicine, University of Bern , Bern ,","place":["Switzerland"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-6220-2239","authenticated-orcid":false,"given":"Timothy","family":"Vaughan","sequence":"additional","affiliation":[{"name":"Department of Biosystems Science and Engineering, ETH Zurich , Basel ,","place":["Switzerland"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-6431-535X","authenticated-orcid":false,"given":"Tanja","family":"Stadler","sequence":"additional","affiliation":[{"name":"Department of Biosystems Science and Engineering, ETH Zurich , Basel ,","place":["Switzerland"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-8052-3925","authenticated-orcid":false,"given":"Sebastian","family":"Bonhoeffer","sequence":"additional","affiliation":[{"name":"Institute of Integrative Biology, ETH Zurich , Zurich ,","place":["Switzerland"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-8319-5293","authenticated-orcid":false,"given":"Roland","family":"Regoes","sequence":"additional","affiliation":[{"name":"Institute of Integrative Biology, ETH Zurich , Zurich ,","place":["Switzerland"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2026,2,2]]},"reference":[{"key":"2026021712080309400_msag014-B1","doi-asserted-by":"publisher","first-page":"e90","DOI":"10.1371\/journal.pbio.0040090","article-title":"Inefficient cytotoxic t lymphocyte\u2013mediated killing of HIV-1\u2013infected cells in vivo","volume":"4","author":"Asquith","year":"2006","journal-title":"PLoS Biol"},{"key":"2026021712080309400_msag014-B2","doi-asserted-by":"publisher","first-page":"220","DOI":"10.1186\/1471-2148-11-220","article-title":"Strength and tempo of selection revealed in viral gene genealogies","volume":"11","author":"Bedford","year":"2011","journal-title":"BMC Evol Biol"},{"key":"2026021712080309400_msag014-B3","doi-asserted-by":"publisher","first-page":"2400","DOI":"10.1093\/molbev\/msz155","article-title":"Parallel evolution of HIV-1 in a long-term experiment","volume":"36","author":"Bertels","year":"2019","journal-title":"Mol Biol Evol"},{"key":"2026021712080309400_msag014-B4","doi-asserted-by":"publisher","DOI":"10.24072\/pcjournal.6","article-title":"Convergent evolution as an indicator for selection during acute HIV-1 infection","volume":"1","author":"Bertels","year":"2021","journal-title":"Peer Community J"},{"key":"2026021712080309400_msag014-B6","doi-asserted-by":"publisher","first-page":"vey029","DOI":"10.1093\/ve\/vey029","article-title":"Estimating the mutational fitness effects distribution during early HIV infection","volume":"4","author":"Bons","year":"2018","journal-title":"Virus Evol"},{"key":"2026021712080309400_msag014-B7","doi-asserted-by":"publisher","first-page":"134","DOI":"10.1111\/imr.2018.285.issue-1","article-title":"Virus dynamics and phyloanatomy: merging population dynamic and phylogenetic approaches","volume":"285","author":"Bons","year":"2018","journal-title":"Immunol Rev"},{"key":"2026021712080309400_msag014-B8","doi-asserted-by":"publisher","first-page":"e1006650","DOI":"10.1371\/journal.pcbi.1006650","article-title":"BEAST 2.5: an advanced software platform for Bayesian evolutionary analysis","volume":"15","author":"Bouckaert","year":"2019","journal-title":"PLoS Comput Biol"},{"key":"2026021712080309400_msag014-B9","doi-asserted-by":"publisher","first-page":"8986","DOI":"10.1073\/pnas.0900233106","article-title":"Evidence for an ancient adaptive episode of convergent molecular evolution","volume":"106","author":"Castoe","year":"2009","journal-title":"Proc Natl Acad Sci U S A"},{"key":"2026021712080309400_msag014-B10","doi-asserted-by":"publisher","first-page":"1642","DOI":"10.1093\/infdis\/jit678","article-title":"Hiv migration between blood and cerebrospinal fluid or semen over time","volume":"209","author":"Chaillon","year":"2014","journal-title":"J Infect Dis"},{"key":"2026021712080309400_msag014-B11","doi-asserted-by":"publisher","first-page":"20120206","DOI":"10.1098\/rstb.2012.0206","article-title":"Graph hierarchies for phylogeography","volume":"368","author":"Cybis","year":"2013","journal-title":"Philos Trans R Soc Lond B Biol Sci"},{"key":"2026021712080309400_msag014-B12","doi-asserted-by":"publisher","first-page":"e88","DOI":"10.1371\/journal.pbio.0040088","article-title":"Relaxed phylogenetics and dating with confidence","volume":"4","author":"Drummond","year":"2006","journal-title":"PLoS Biol"},{"key":"2026021712080309400_msag014-B13","doi-asserted-by":"publisher","first-page":"1441","DOI":"10.1534\/genetics.105.052019","article-title":"Evolution of the human immunodeficiency virus envelope gene is dominated by purifying selection","volume":"174","author":"Edwards","year":"2006","journal-title":"Genetics"},{"key":"2026021712080309400_msag014-B14","doi-asserted-by":"publisher","first-page":"2407","DOI":"10.1534\/genetics.104.030411","article-title":"Using temporally spaced sequences to simultaneously estimate migration rates, mutation rate and population sizes in measurably evolving populations","volume":"168","author":"Ewing","year":"2004","journal-title":"Genetics"},{"key":"2026021712080309400_msag014-B15","doi-asserted-by":"publisher","first-page":"56","DOI":"10.1126\/science.1256739","article-title":"The early spread and epidemic ignition of HIV-1 in human populations","volume":"346","author":"Faria","year":"2014","journal-title":"Science"},{"key":"2026021712080309400_msag014-B16","doi-asserted-by":"publisher","first-page":"5721","DOI":"10.1128\/JVI.79.9.5721-5731.2005","article-title":"Rapid viral escape at an immunodominant simian-human immunodeficiency virus cytotoxic t-lymphocyte epitope exacts a dramatic fitness cost","volume":"79","author":"Fernandez","year":"2005","journal-title":"J Virol"},{"key":"2026021712080309400_msag014-B17","doi-asserted-by":"publisher","first-page":"1557","DOI":"10.1126\/science.1176062","article-title":"Pandemic potential of a strain of influenza a (h1n1): early findings","volume":"324","author":"Fraser","year":"2009","journal-title":"Science"},{"key":"2026021712080309400_msag014-B18","doi-asserted-by":"publisher","first-page":"327","DOI":"10.1126\/science.1090727","article-title":"Unifying the epidemiological and evolutionary dynamics of pathogens","volume":"303","author":"Grenfell","year":"2004","journal-title":"Science"},{"key":"2026021712080309400_msag014-B19","doi-asserted-by":"publisher","first-page":"905","DOI":"10.1093\/molbev\/msp296","article-title":"Isolation with migration models for more than two populations","volume":"27","author":"Hey","year":"2010","journal-title":"Mol Biol Evol"},{"key":"2026021712080309400_msag014-B20","doi-asserted-by":"publisher","first-page":"7552","DOI":"10.1073\/pnas.0802203105","article-title":"Identification and characterization of transmitted and early founder virus envelopes in primary HIV-1 infection","volume":"105","author":"Keele","year":"2008","journal-title":"Proc Natl Acad Sci U S A"},{"key":"2026021712080309400_msag014-B21","doi-asserted-by":"publisher","first-page":"27","DOI":"10.2307\/3213548","article-title":"On the genealogy of large populations","volume":"19","author":"Kingman","year":"1982","journal-title":"J Appl Probab"},{"key":"2026021712080309400_msag014-B22","doi-asserted-by":"publisher","first-page":"20230314","DOI":"10.1098\/rstb.2023.0314","article-title":"Phylodynamics beyond neutrality: the impact of incomplete purifying selection on viral phylogenies and inference","volume":"380","author":"Koelle","year":"2025","journal-title":"Philos Trans R Soc Lond B Biol Sci"},{"key":"2026021712080309400_msag014-B23","doi-asserted-by":"publisher","first-page":"2102","DOI":"10.1093\/molbev\/msw064","article-title":"Phylodynamics with migration: a computational framework to quantify population structure from genomic data","volume":"33","author":"K\u00fchnert","year":"2016","journal-title":"Mol Biol Evol"},{"key":"2026021712080309400_msag014-B24","doi-asserted-by":"publisher","first-page":"341","DOI":"10.1016\/j.jtbi.2009.07.038","article-title":"Modeling sequence evolution in acute HIV-1 infection","volume":"261","author":"Lee","year":"2009","journal-title":"J Theor Biol"},{"key":"2026021712080309400_msag014-B25","doi-asserted-by":"publisher","first-page":"e1000520","DOI":"10.1371\/journal.pcbi.1000520","article-title":"Bayesian phylogeography finds its roots","volume":"5","author":"Lemey","year":"2009","journal-title":"PLoS Comput Biol"},{"key":"2026021712080309400_msag014-B26","doi-asserted-by":"publisher","first-page":"1877","DOI":"10.1093\/molbev\/msq067","article-title":"Phylogeography takes a relaxed random walk in continuous space and time","volume":"27","author":"Lemey","year":"2010","journal-title":"Mol Biol Evol"},{"key":"2026021712080309400_msag014-B27","doi-asserted-by":"publisher","first-page":"51","DOI":"10.1038\/nature16933","article-title":"Persistent HIV-1 replication maintains the tissue reservoir during therapy","volume":"530","author":"Lorenzo-Redondo","year":"2016","journal-title":"Nature"},{"key":"2026021712080309400_msag014-B28","doi-asserted-by":"publisher","first-page":"1053","DOI":"10.1093\/bioinformatics\/btx701","article-title":"Efficient comparative phylogenetics on large trees","volume":"34","author":"Louca","year":"2018","journal-title":"Bioinformatics"},{"key":"2026021712080309400_msag014-B29","doi-asserted-by":"publisher","first-page":"11982","DOI":"10.1128\/JVI.00946-07","article-title":"Estimating the effectiveness of simian immunodeficiency virus-specific cd8+ t cells from the dynamics of viral immune escape","volume":"81","author":"Mandl","year":"2007","journal-title":"J Virol"},{"key":"2026021712080309400_msag014-B30","doi-asserted-by":"publisher","first-page":"e1013421","DOI":"10.1371\/journal.pcbi.1013421","article-title":"MASCOT-skyline integrates population and migration dynamics to enhance phylogeographic reconstructions","volume":"21","author":"M\u00fcller","year":"2025","journal-title":"PLoS Comput Biol"},{"key":"2026021712080309400_msag014-B31","doi-asserted-by":"publisher","first-page":"3843","DOI":"10.1093\/bioinformatics\/bty406","article-title":"MASCOT: parameter and state inference under the marginal structured coalescent approximation","volume":"34","author":"M\u00fcller","year":"2018","journal-title":"Bioinformatics"},{"key":"2026021712080309400_msag014-B32","doi-asserted-by":"publisher","first-page":"2970","DOI":"10.1093\/molbev\/msx186","article-title":"The structured coalescent and its approximations","volume":"34","author":"M\u00fcller","year":"2017","journal-title":"Mol Biol Evol"},{"key":"2026021712080309400_msag014-B33","doi-asserted-by":"publisher","first-page":"5540","DOI":"10.1128\/JVI.77.9.5540-5546.2003","article-title":"Evolutionary indicators of human immunodeficiency virus type 1 reservoirs and compartments","volume":"77","author":"Nickle","year":"2003","journal-title":"J Virol"},{"key":"2026021712080309400_msag014-B34","doi-asserted-by":"publisher","first-page":"595","DOI":"10.1534\/genetics.114.172791","article-title":"Inferring epidemiological dynamics with Bayesian coalescent inference: the merits of deterministic and stochastic models","volume":"199","author":"Popinga","year":"2015","journal-title":"Genetics"},{"key":"2026021712080309400_msag014-B35","doi-asserted-by":"publisher","first-page":"8240","DOI":"10.1128\/JVI.72.10.8240-8251.1998","article-title":"Evolution of envelope sequences from the genital tract and peripheral blood of women infected with clade a human immunodeficiency virus type 1","volume":"72","author":"Poss","year":"1998","journal-title":"J Virol"},{"key":"2026021712080309400_msag014-B36","doi-asserted-by":"publisher","first-page":"1232","DOI":"10.1038\/s41592-023-01917-9","article-title":"CherryML: scalable maximum likelihood estimation of phylogenetic models","volume":"20","author":"Prillo","year":"2023","journal-title":"Nat Methods"},{"key":"2026021712080309400_msag014-B37","doi-asserted-by":"publisher","first-page":"e1002136","DOI":"10.1371\/journal.pcbi.1002136","article-title":"Inference for nonlinear epidemiological models using genealogies and time series","volume":"7","author":"Rasmussen","year":"2011","journal-title":"PLoS Comput Biol"},{"key":"2026021712080309400_msag014-B38","doi-asserted-by":"publisher","first-page":"e45562","DOI":"10.7554\/eLife.45562","article-title":"Coupling adaptive molecular evolution to phylodynamics using fitness-dependent birth-death models","volume":"8","author":"Rasmussen","year":"2019","journal-title":"Elife"},{"key":"2026021712080309400_msag014-B39","doi-asserted-by":"publisher","first-page":"185","DOI":"10.1093\/molbev\/msw217","article-title":"Phylogenetic tools for generalized HIV-1 epidemics: findings from the PANGEA-HIV methods comparison","volume":"34","author":"Ratmann","year":"2017","journal-title":"Mol Biol Evol"},{"key":"2026021712080309400_msag014-B40","doi-asserted-by":"publisher","first-page":"e87428","DOI":"10.1371\/journal.pone.0087428","article-title":"Evaluating the effects of non-neutral molecular markers on phylogeny inference","volume":"9","author":"Roje","year":"2014","journal-title":"PLoS One"},{"key":"2026021712080309400_msag014-B41","doi-asserted-by":"publisher","first-page":"vex042","DOI":"10.1093\/ve\/vex042","article-title":"TreeTime: maximum-likelihood phylodynamic analysis","volume":"4","author":"Sagulenko","year":"2018","journal-title":"Virus Evol"},{"key":"2026021712080309400_msag014-B42","doi-asserted-by":"publisher","first-page":"110","DOI":"10.2174\/1570162X13666151029102413","article-title":"Phylogenetics and phyloanatomy of HIV\/SIV intra-host compartments and reservoirs: the key role of the central nervous system","volume":"14","author":"Salemi","year":"2016","journal-title":"Curr HIV Res"},{"key":"2026021712080309400_msag014-B43","author":"Shirreff"},{"key":"2026021712080309400_msag014-B44","doi-asserted-by":"publisher","first-page":"676","DOI":"10.1093\/sysbio\/syr029","article-title":"Simulating trees with a fixed number of extant species","volume":"60","author":"Stadler","year":"2011","journal-title":"Syst Biol"},{"key":"2026021712080309400_msag014-B45","doi-asserted-by":"publisher","first-page":"228","DOI":"10.1073\/pnas.1207965110","article-title":"Birth\u2013death skyline plot reveals temporal changes of epidemic spread in HIV and hepatitis c virus (HCV)","volume":"110","author":"Stadler","year":"2013","journal-title":"Proc Natl Acad Sci U S A"},{"key":"2026021712080309400_msag014-B46","doi-asserted-by":"publisher","DOI":"10.1371\/currents.outbreaks.02bc6d927ecee7bbd33532ec8ba6a25f","article-title":"Insights into the early epidemic spread of ebola in sierra leone provided by viral sequence data","volume":"6","author":"Stadler","year":"2014","journal-title":"PLoS Curr"},{"key":"2026021712080309400_msag014-B47","doi-asserted-by":"publisher","first-page":"2298","DOI":"10.1093\/oxfordjournals.molbev.a003776","article-title":"Exploring the demographic history of DNA sequences using the generalized skyline plot","volume":"18","author":"Strimmer","year":"2001","journal-title":"Mol Biol Evol"},{"key":"2026021712080309400_msag014-B48","doi-asserted-by":"publisher","first-page":"1370","DOI":"10.1093\/bioinformatics\/btp244","article-title":"Many-core algorithms for statistical phylogenetics","volume":"25","author":"Suchard","year":"2009","journal-title":"Bioinformatics"},{"key":"2026021712080309400_msag014-B49","doi-asserted-by":"publisher","first-page":"1480","DOI":"10.1093\/molbev\/mst057","article-title":"A stochastic simulator of birth\u2013death master equations with application to phylodynamics","volume":"30","author":"Vaughan","year":"2013","journal-title":"Mol Biol Evol"},{"key":"2026021712080309400_msag014-B50","doi-asserted-by":"publisher","first-page":"2272","DOI":"10.1093\/bioinformatics\/btu201","article-title":"Efficient Bayesian inference under the structured coalescent","volume":"30","author":"Vaughan","year":"2014","journal-title":"Bioinformatics"},{"key":"2026021712080309400_msag014-B51","doi-asserted-by":"publisher","first-page":"e2308125121","DOI":"10.1073\/pnas.2308125121","article-title":"Estimates of early outbreak-specific SARS-CoV-2 epidemiological parameters from genomic data","volume":"121","author":"Vaughan","year":"2024","journal-title":"Proc Natl Acad Sci U S A"},{"key":"2026021712080309400_msag014-B52","doi-asserted-by":"publisher","first-page":"e1002947","DOI":"10.1371\/journal.pcbi.1002947","article-title":"Viral phylodynamics","volume":"9","author":"Volz","year":"2013","journal-title":"PLoS Comput Biol"},{"key":"2026021712080309400_msag014-B53","doi-asserted-by":"publisher","first-page":"1421","DOI":"10.1534\/genetics.109.106021","article-title":"Phylodynamics of infectious disease epidemics","volume":"183","author":"Volz","year":"2009","journal-title":"Genetics"},{"key":"2026021712080309400_msag014-B54","doi-asserted-by":"publisher","first-page":"1376","DOI":"10.1093\/oxfordjournals.molbev.a004199","article-title":"The genealogy of a sequence subject to purifying selection at multiple sites","volume":"19","author":"Williamson","year":"2002","journal-title":"Mol Biol Evol"},{"key":"2026021712080309400_msag014-B55","doi-asserted-by":"publisher","first-page":"vex003","DOI":"10.1093\/ve\/vex003","article-title":"In vivo mutation rates and the landscape of fitness costs of HIV-1","volume":"3","author":"Zanini","year":"2017","journal-title":"Virus Evol"},{"key":"2026021712080309400_msag014-B56","doi-asserted-by":"publisher","first-page":"6643","DOI":"10.1128\/JVI.02268-06","article-title":"Comparative study of methods for detecting sequence compartmentalization in human immunodeficiency virus type 1","volume":"81","author":"Z\u00e1rate","year":"2007","journal-title":"J Virol"}],"container-title":["Molecular Biology and Evolution"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/mbe\/advance-article-pdf\/doi\/10.1093\/molbev\/msag014\/66683962\/msag014.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/mbe\/article-pdf\/43\/2\/msag014\/66683962\/msag014.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/mbe\/article-pdf\/43\/2\/msag014\/66683962\/msag014.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,2,17]],"date-time":"2026-02-17T17:08:14Z","timestamp":1771348094000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/mbe\/article\/doi\/10.1093\/molbev\/msag014\/8451695"}},"subtitle":[],"editor":[{"given":"Yuseob","family":"Kim","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2026,2]]},"references-count":55,"journal-issue":{"issue":"2","published-print":{"date-parts":[[2026,2,2]]}},"URL":"https:\/\/doi.org\/10.1093\/molbev\/msag014","relation":{},"ISSN":["0737-4038","1537-1719"],"issn-type":[{"value":"0737-4038","type":"print"},{"value":"1537-1719","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2026,2]]},"published":{"date-parts":[[2026,2]]},"article-number":"msag014"}}