{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,4,3]],"date-time":"2026-04-03T20:14:35Z","timestamp":1775247275190,"version":"3.50.1"},"reference-count":53,"publisher":"Oxford University Press (OUP)","issue":"10","license":[{"start":{"date-parts":[[2018,10,13]],"date-time":"2018-10-13T00:00:00Z","timestamp":1539388800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/journals\/pages\/open_access\/funder_policies\/chorus\/standard_publication_model"}],"funder":[{"DOI":"10.13039\/100007225","name":"Ministry of Science and Technology","doi-asserted-by":"publisher","award":["MOST 106-2221-E-126-008"],"award-info":[{"award-number":["MOST 106-2221-E-126-008"]}],"id":[{"id":"10.13039\/100007225","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100007225","name":"Ministry of Science and Technology","doi-asserted-by":"publisher","award":["MOST 106-2118-M-007-004"],"award-info":[{"award-number":["MOST 106-2118-M-007-004"]}],"id":[{"id":"10.13039\/100007225","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2019,5,15]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:sec>\n                    <jats:title>Motivation<\/jats:title>\n                    <jats:p>High-throughput sequencing technology has revolutionized the study of metagenomics and cancer evolution. In a relatively simple environment, a metagenomics sequencing data is dominated by a few species. By analyzing the alignment of reads from microbial species, single nucleotide polymorphisms can be discovered and the evolutionary history of the populations can be reconstructed. The ever-increasing read length will allow more detailed analysis about the evolutionary history of microbial or tumor cell population. A simulator of shotgun sequences from such populations will be helpful in the development or evaluation of analysis algorithms.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>Here, we described an efficient algorithm, MetaSMC, which simulates reads from evolving microbial populations. Based on the coalescent theory, our simulator supports all evolutionary scenarios supported by other coalescent simulators. In addition, the simulator supports various substitution models, including Jukes-Cantor, HKY85 and generalized time-reversible models. The simulator also supports mutator phenotypes by allowing different mutation rates and substitution models in different subpopulations. Our algorithm ignores unnecessary chromosomal segments and thus is more efficient than standard coalescent when recombination is frequent. We showed that the process behind our algorithm is equivalent to Sequentially Markov Coalescent with an incomplete sample. The accuracy of our algorithm was evaluated by summary statistics and likelihood curves derived from Monte Carlo integration over large number of random genealogies.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Availability and implementation<\/jats:title>\n                    <jats:p>MetaSMC is written in C. The source code is available at https:\/\/github.com\/tarjxvf\/metasmc.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Supplementary information<\/jats:title>\n                    <jats:p>Supplementary data are available at Bioinformatics online.<\/jats:p>\n                  <\/jats:sec>","DOI":"10.1093\/bioinformatics\/bty840","type":"journal-article","created":{"date-parts":[[2018,10,11]],"date-time":"2018-10-11T07:45:08Z","timestamp":1539243908000},"page":"1677-1685","source":"Crossref","is-referenced-by-count":1,"title":["MetaSMC: a coalescent-based shotgun sequence simulator for evolving microbial populations"],"prefix":"10.1093","volume":"35","author":[{"given":"Ki-Hok","family":"Liao","sequence":"first","affiliation":[{"name":"Department of Computer Science, National Tsing-Hua University, Hsinchu, Taiwan"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Wing-Kai","family":"Hon","sequence":"additional","affiliation":[{"name":"Department of Computer Science, National Tsing-Hua University, Hsinchu, Taiwan"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Chuan-Yi","family":"Tang","sequence":"additional","affiliation":[{"name":"Department of Computer Science, National Tsing-Hua University, Hsinchu, Taiwan"},{"name":"Department of Computer Science and Information Engineering, Providence University, Taichung, Taiwan"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Wen-Ping","family":"Hsieh","sequence":"additional","affiliation":[{"name":"Institute of Statistics, National Tsing-Hua University, Hsinchu, Taiwan"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2018,10,13]]},"reference":[{"key":"2023013107481477000_bty840-B1","doi-asserted-by":"crossref","first-page":"i420","DOI":"10.1093\/bioinformatics\/btq365","article-title":"Characteristics of 454 pyrosequencing data-enabling realistic simulation with flowsim","volume":"26","author":"Balzer","year":"2010","journal-title":"Bioinformatics"},{"key":"2023013107481477000_bty840-B2","doi-asserted-by":"crossref","first-page":"e1002745","DOI":"10.1371\/journal.ppat.1002745","article-title":"A high-resolution view of genome-wide pneumococcal transformation","volume":"8","author":"Croucher","year":"2012","journal-title":"PLoS Pathogens"},{"key":"2023013107481477000_bty840-B3","doi-asserted-by":"crossref","first-page":"820","DOI":"10.1111\/j.1365-2958.2006.05150.x","article-title":"Evolution of mutation rates in bacteria","volume":"60","author":"Denamur","year":"2006","journal-title":"Mol. Microbiol"},{"key":"2023013107481477000_bty840-B4","doi-asserted-by":"crossref","first-page":"997","DOI":"10.1534\/genetics.111.128116","article-title":"The balance between mutators and nonmutators in asexual populations","volume":"188","author":"Desai","year":"2011","journal-title":"Genetics"},{"key":"2023013107481477000_bty840-B5","doi-asserted-by":"crossref","first-page":"150","DOI":"10.1038\/nrmicro.2015.13","article-title":"Within-host evolution of bacterial pathogens","volume":"14","author":"Didelot","year":"2016","journal-title":"Nat. Rev. Microbiol"},{"key":"2023013107481477000_bty840-B6","doi-asserted-by":"crossref","first-page":"1053","DOI":"10.1093\/jac\/dkn044","article-title":"Increased mutability of Pseudomonas aeruginosa in biofilms","volume":"61","author":"Driffield","year":"2008","journal-title":"J. Antimicrob. Chemother"},{"key":"2023013107481477000_bty840-B7","doi-asserted-by":"crossref","first-page":"398","DOI":"10.1186\/1471-2105-8-398","article-title":"Strainer: software for analysis of population variation in community genomic datasets","volume":"8","author":"Eppley","year":"2007","journal-title":"BMC Bioinformatics"},{"key":"2023013107481477000_bty840-B8","doi-asserted-by":"crossref","first-page":"e1004651.","DOI":"10.1371\/journal.pgen.1004651","article-title":"Coexistence and within-host evolution of diversified lineages of hypermutable pseudomonas aeruginosa in long-term cystic fibrosis infections","volume":"10","author":"Feliziani","year":"2014","journal-title":"PLoS Genet"},{"key":"2023013107481477000_bty840-B9","doi-asserted-by":"crossref","first-page":"368","DOI":"10.1007\/BF01734359","article-title":"Evolutionary trees from DNA sequences: a maximum likelihood approach","volume":"17","author":"Felsenstein","year":"1981","journal-title":"J. Mol. Evol"},{"key":"2023013107481477000_bty840-B10","doi-asserted-by":"crossref","first-page":"e1000208","DOI":"10.1371\/journal.pgen.1000208","article-title":"Mutability and importance of a hypermutable cell subpopulation that produces stress-induced mutants in Escherichia coli","volume":"4","author":"Gonzalez","year":"2008","journal-title":"PLoS Genet"},{"key":"2023013107481477000_bty840-B11","doi-asserted-by":"crossref","first-page":"1249","DOI":"10.1534\/genetics.116.193565","article-title":"Evolution of mutation rates in rapidly adapting asexual populations","volume":"204","author":"Good","year":"2016","journal-title":"Genetics"},{"key":"2023013107481477000_bty840-B12","doi-asserted-by":"crossref","first-page":"479","DOI":"10.1089\/cmb.1996.3.479","article-title":"Ancestral inference from samples of DNA sequences with recombination","volume":"3","author":"Griffiths","year":"1996","journal-title":"J. Comput. Biol"},{"key":"2023013107481477000_bty840-B13","doi-asserted-by":"crossref","DOI":"10.1098\/rstb.1994.0079","article-title":"Sampling theory for neutral alleles in a varying environment","volume":"344","author":"Griffiths","year":"1994","journal-title":"Philos. Trans. Biol. Sci"},{"key":"2023013107481477000_bty840-B14","doi-asserted-by":"crossref","first-page":"277","DOI":"10.1111\/j.1365-294X.2009.04482.x","article-title":"mlRho a program for estimating the population mutation and recombination rates from shotgun-sequenced diploid genomes","volume":"19","author":"Haubold","year":"2010","journal-title":"Mol. Ecol"},{"key":"2023013107481477000_bty840-B15","doi-asserted-by":"crossref","first-page":"1020","DOI":"10.1101\/gr.074187.107","article-title":"Population genetic analysis of shotgun assemblies of genomic sequences from multiple individuals","volume":"18","author":"Hellmann","year":"2008","journal-title":"Genome Res"},{"key":"2023013107481477000_bty840-B16","doi-asserted-by":"crossref","first-page":"593","DOI":"10.1093\/bioinformatics\/btr708","article-title":"ART: a next-generation sequencing read simulator","volume":"28","author":"Huang","year":"2012","journal-title":"Bioinformatics"},{"key":"2023013107481477000_bty840-B17","doi-asserted-by":"crossref","first-page":"183","DOI":"10.1016\/0040-5809(83)90013-8","article-title":"Properties of a neutral allele model with intragenic recombination","volume":"23","author":"Hudson","year":"1983","journal-title":"Theor. Popul. Biol"},{"key":"2023013107481477000_bty840-B18","doi-asserted-by":"crossref","first-page":"337","DOI":"10.1093\/bioinformatics\/18.2.337","article-title":"Generating samples under a Wright-Fisher neutral model of genetic variation","volume":"18","author":"Hudson","year":"2002","journal-title":"Bioinformatics"},{"key":"2023013107481477000_bty840-B19","doi-asserted-by":"crossref","first-page":"e1000674.","DOI":"10.1371\/journal.pgen.1000674","article-title":"Inference of microbial recombination rates from metagenomic data","volume":"5","author":"Johnson","year":"2009","journal-title":"PLoS Genet"},{"key":"2023013107481477000_bty840-B20","doi-asserted-by":"crossref","first-page":"1320","DOI":"10.1101\/gr.5431206","article-title":"Inference of population genetic parameters in metagenomics: a clean look at messy data","volume":"16","author":"Johnson","year":"2006","journal-title":"Genome Res"},{"key":"2023013107481477000_bty840-B21","doi-asserted-by":"crossref","first-page":"887","DOI":"10.1093\/genetics\/123.4.887","article-title":"The \u201cHitchhiking Effect\u201d revisited","volume":"123","author":"Kaplan","year":"1989","journal-title":"Genetics"},{"key":"2023013107481477000_bty840-B22","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1371\/journal.pcbi.1004842","article-title":"Efficient coalescent simulation and genealogical analysis for large sample sizes","volume":"12","author":"Kelleher","year":"2016","journal-title":"PLoS Comput. Biol"},{"key":"2023013107481477000_bty840-B23","doi-asserted-by":"crossref","first-page":"247.","DOI":"10.1186\/1471-2105-10-247","article-title":"Accurate and fast methods to estimate the population mutation rate from error prone sequences","volume":"10","author":"Knudsen","year":"2009","journal-title":"BMC Bioinformatics"},{"key":"2023013107481477000_bty840-B24","doi-asserted-by":"crossref","first-page":"86","DOI":"10.1016\/j.tree.2008.09.007","article-title":"Coalescent genealogy samplers: windows into population history","volume":"24","author":"Kuhner","year":"2009","journal-title":"Trends Ecol. Evol"},{"key":"2023013107481477000_bty840-B25","doi-asserted-by":"crossref","first-page":"557","DOI":"10.1128\/MMBR.00009-08","article-title":"A bioinformatician\u2019s guide to metagenomics","volume":"72","author":"Kunin","year":"2008","journal-title":"Microbiol. Mol. Biol. Rev"},{"key":"2023013107481477000_bty840-B26","doi-asserted-by":"crossref","first-page":"375","DOI":"10.1038\/nrc3039","article-title":"An analogy between the evolution of drug resistance in bacterial communities and malignant tissues","volume":"11","author":"Lambert","year":"2011","journal-title":"Nat. Rev.. Cancer"},{"key":"2023013107481477000_bty840-B27","doi-asserted-by":"crossref","first-page":"493","DOI":"10.1038\/nature10231","article-title":"Inference of human population history from individual whole-genome sequences","volume":"475","author":"Li","year":"2011","journal-title":"Nature"},{"key":"2023013107481477000_bty840-B28","doi-asserted-by":"crossref","first-page":"e37919","DOI":"10.1371\/journal.pone.0037919","article-title":"Deep sequencing of the oral microbiome reveals signatures of periodontal disease","volume":"7","author":"Liu","year":"2012","journal-title":"PloS One"},{"key":"2023013107481477000_bty840-B29","doi-asserted-by":"crossref","first-page":"1479","DOI":"10.1093\/molbev\/msp059","article-title":"Inferring population mutation rate and sequencing error rate using the SNP frequency spectrum in a sample of DNA sequences","volume":"26","author":"Liu","year":"2009","journal-title":"Mol. Biol. Evol"},{"key":"2023013107481477000_bty840-B30","doi-asserted-by":"crossref","first-page":"101","DOI":"10.1101\/gr.097543.109","article-title":"Estimating population genetic parameters and comparing model goodness-of-fit using DNA sequences with error","volume":"20","author":"Liu","year":"2010","journal-title":"Genome Res"},{"key":"2023013107481477000_bty840-B31","doi-asserted-by":"crossref","first-page":"16","DOI":"10.1186\/1471-2156-7-16","article-title":"Fast \u201ccoalescent\u201d simulation","volume":"7","author":"Marjoram","year":"2006","journal-title":"BMC Genet"},{"key":"2023013107481477000_bty840-B32","doi-asserted-by":"crossref","first-page":"1.","DOI":"10.1186\/2042-5783-4-1","article-title":"Deep sequencing of evolving pathogen populations: applications, errors, and bioinformatic solutions","volume":"4","author":"McElroy","year":"2014","journal-title":"Microb. Inform. Exp"},{"key":"2023013107481477000_bty840-B33","doi-asserted-by":"crossref","first-page":"74","DOI":"10.1186\/1471-2164-13-74","article-title":"GemSIM: general, error-model based simulator of next-generation sequencing data","volume":"13","author":"McElroy","year":"2012","journal-title":"BMC Genomics"},{"key":"2023013107481477000_bty840-B34","doi-asserted-by":"crossref","first-page":"1387","DOI":"10.1098\/rstb.2005.1673","article-title":"Approximating the coalescent with recombination","volume":"360","author":"McVean","year":"2005","journal-title":"Philos. Trans. R. Soc. Lond. B Biol. Sci"},{"key":"2023013107481477000_bty840-B35","doi-asserted-by":"crossref","first-page":"717","DOI":"10.1534\/g3.113.009597","article-title":"Extensive cotransformation of natural variation into chromosomes of naturally competent haemophilus influenzae","volume":"4","author":"Mell","year":"2014","journal-title":"G3"},{"key":"2023013107481477000_bty840-B36","doi-asserted-by":"crossref","first-page":"59","DOI":"10.1007\/BF00173909","article-title":"The coalescent and the genealogical process in geographically structured population","volume":"29","author":"Notohara","year":"1990","journal-title":"J. Math. Biol"},{"key":"2023013107481477000_bty840-B37","doi-asserted-by":"crossref","first-page":"321","DOI":"10.1534\/genetics.110.117986","article-title":"A principled approach to deriving approximate conditional sampling distributions in population genetics models with recombination","volume":"186","author":"Paul","year":"2010","journal-title":"Genetics"},{"key":"2023013107481477000_bty840-B38","doi-asserted-by":"crossref","first-page":"1115","DOI":"10.1534\/genetics.110.125534","article-title":"An accurate sequentially Markov conditional sampling distribution for the coalescent with recombination","volume":"187","author":"Paul","year":"2011","journal-title":"Genetics"},{"key":"2023013107481477000_bty840-B39","doi-asserted-by":"crossref","first-page":"752","DOI":"10.1111\/mmi.13661","article-title":"The processive kinetics of gene conversion in bacteria","volume":"104","author":"Paulsson","year":"2017","journal-title":"Mol. Microbiol"},{"key":"2023013107481477000_bty840-B40","doi-asserted-by":"crossref","first-page":"55","DOI":"10.1038\/nature11450","article-title":"A metagenome-wide association study of gut microbiota in type 2 diabetes","volume":"490","author":"Qin","year":"2012","journal-title":"Nature"},{"key":"2023013107481477000_bty840-B41","first-page":"235","article-title":"Seq-Gen: an application for the Monte Carlo simulation of DNA sequence evolution along phylogenetic trees","volume":"13","author":"Rambaut","year":"1997","journal-title":"Comput. Appl. Biosci"},{"key":"2023013107481477000_bty840-B42","doi-asserted-by":"crossref","first-page":"e3373","DOI":"10.1371\/journal.pone.0003373","article-title":"MetaSim-A sequencing simulator for genomics and metagenomics","volume":"3","author":"Richter","year":"2008","journal-title":"Plos One"},{"key":"2023013107481477000_bty840-B43","doi-asserted-by":"crossref","first-page":"919","DOI":"10.1038\/ng.3015","article-title":"Inferring human population size and separation history from multiple genome sequences","volume":"46","author":"Schiffels","year":"2014","journal-title":"Nat. Genet"},{"key":"2023013107481477000_bty840-B44","doi-asserted-by":"crossref","first-page":"647","DOI":"10.1534\/genetics.112.149096","article-title":"Estimating variable effective population sizes from multiple genomes: a sequentially markov conditional sampling distribution approach","volume":"194","author":"Sheehan","year":"2013","journal-title":"Genetics"},{"key":"2023013107481477000_bty840-B45","doi-asserted-by":"crossref","first-page":"e177","DOI":"10.1371\/journal.pbio.0060177","article-title":"Population genomic analysis of strain variation in Leptospirillum group II bacteria involved in acid mine drainage formation","volume":"6","author":"Simmons","year":"2008","journal-title":"PLoS Biol"},{"key":"2023013107481477000_bty840-B46","doi-asserted-by":"crossref","first-page":"703","DOI":"10.1038\/42701","article-title":"Evolution of high mutation rates in experimental populations of E. coli","volume":"387","author":"Sniegowski","year":"1997","journal-title":"Nature"},{"key":"2023013107481477000_bty840-B47","doi-asserted-by":"crossref","first-page":"51","DOI":"10.1016\/j.tpb.2013.11.004","article-title":"Within a sample from a population, the distribution of the number of descendants of a subsample\u2019s most recent common ancestor","volume":"92","author":"Spouge","year":"2014","journal-title":"Theor. Popul. Biol"},{"key":"2023013107481477000_bty840-B48","doi-asserted-by":"crossref","first-page":"51","DOI":"10.1016\/j.tpb.2012.08.004","article-title":"A sequentially Markov conditional sampling distribution for structured populations with migration and recombination","volume":"87","author":"Steinr\u00fccken","year":"2012","journal-title":"Theor. Popul. Biol"},{"key":"2023013107481477000_bty840-B49","doi-asserted-by":"crossref","first-page":"213","DOI":"10.1017\/S0016672300027683","article-title":"The coalescent in two partially isolated diffusion populations","volume":"52","author":"Takahata","year":"1988","journal-title":"Genet. Res"},{"key":"2023013107481477000_bty840-B50","doi-asserted-by":"crossref","first-page":"37","DOI":"10.1038\/nature02340","article-title":"Community structure and metabolism through reconstruction of microbial genomes from the environment","volume":"428","author":"Tyson","year":"2004","journal-title":"Nature"},{"key":"2023013107481477000_bty840-B51","doi-asserted-by":"crossref","first-page":"248","DOI":"10.1006\/tpbi.1998.1403","article-title":"Recombination as a point process along sequences","volume":"55","author":"Wiuf","year":"1999","journal-title":"Theor. Popul. Biol"},{"key":"2023013107481477000_bty840-B52","doi-asserted-by":"crossref","first-page":"451","DOI":"10.1093\/genetics\/155.1.451","article-title":"The coalescent with gene conversion","volume":"155","author":"Wiuf","year":"2000","journal-title":"Genetics"},{"key":"2023013107481477000_bty840-B53","doi-asserted-by":"crossref","first-page":"86","DOI":"10.1093\/jhered\/ess074","article-title":"Coalescent entanglement and the conditional dependence of the times to common ancestry of mutually exclusive pairs of individuals","volume":"104","author":"Wu","year":"2013","journal-title":"J. Hered"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/35\/10\/1677\/48969892\/bioinformatics_35_10_1677.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/35\/10\/1677\/48969892\/bioinformatics_35_10_1677.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,4,3]],"date-time":"2026-04-03T19:31:50Z","timestamp":1775244710000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/35\/10\/1677\/5128939"}},"subtitle":[],"editor":[{"given":"Russell","family":"Schwartz","sequence":"additional","affiliation":[],"role":[{"role":"editor","vocabulary":"crossref"}]}],"short-title":[],"issued":{"date-parts":[[2018,10,13]]},"references-count":53,"journal-issue":{"issue":"10","published-print":{"date-parts":[[2019,5,15]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/bty840","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2019,5,15]]},"published":{"date-parts":[[2018,10,13]]}}}