{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,3,11]],"date-time":"2026-03-11T23:06:21Z","timestamp":1773270381119,"version":"3.50.1"},"reference-count":37,"publisher":"Oxford University Press (OUP)","issue":"17","license":[{"start":{"date-parts":[[2016,11,10]],"date-time":"2016-11-10T00:00:00Z","timestamp":1478736000000},"content-version":"vor","delay-in-days":73,"URL":"http:\/\/creativecommons.org\/licenses\/by-nc\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100001711","name":"Swiss National Science Foundation","doi-asserted-by":"publisher","award":["33CS30_148522"],"award-info":[{"award-number":["33CS30_148522"]}],"id":[{"id":"10.13039\/501100001711","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100001711","name":"Swiss National Science Foundation","doi-asserted-by":"publisher","award":["320030_159868"],"award-info":[{"award-number":["320030_159868"]}],"id":[{"id":"10.13039\/501100001711","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100001711","name":"Swiss National Science Foundation","doi-asserted-by":"publisher","award":["PZ00P3-142411"],"award-info":[{"award-number":["PZ00P3-142411"]}],"id":[{"id":"10.13039\/501100001711","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2016,9,1]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:p>The continuous time conjunctive Bayesian network (CT-CBN) is a graphical model for analyzing the waiting time process of the accumulation of genetic changes (mutations). CT-CBN models have been successfully used in several biological applications such as HIV drug resistance development and genetic progression of cancer. However, current approaches for parameter estimation and network structure learning of CBNs can only deal with a small number of mutations (&amp;lt;20). Here, we address this limitation by presenting an efficient and accurate approximate inference algorithm using a Monte Carlo expectation-maximization algorithm based on importance sampling. The new method can now be used for a large number of mutations, up to one thousand, an increase by two orders of magnitude. In simulation studies, we present the accuracy as well as the running time efficiency of the new inference method and compare it with a MLE method, expectation-maximization, and discrete time CBN model, i.e. a first-order approximation of the CT-CBN model. We also study the application of the new model on HIV drug resistance datasets for the combination therapy with zidovudine plus lamivudine (AZT\u2009+\u20093TC) as well as under no treatment, both extracted from the Swiss HIV Cohort Study database.<\/jats:p>\n               <jats:sec>\n                  <jats:title>Availability and implementation<\/jats:title>\n                  <jats:p>The proposed method is implemented as an R package available at https:\/\/github.com\/cbg-ethz\/MC-CBN.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Contact<\/jats:title>\n                  <jats:p>niko.beerenwinkel@bsse.ethz.ch<\/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\/btw459","type":"journal-article","created":{"date-parts":[[2016,9,1]],"date-time":"2016-09-01T07:53:39Z","timestamp":1472716419000},"page":"i727-i735","source":"Crossref","is-referenced-by-count":34,"title":["Large-scale inference of conjunctive Bayesian networks"],"prefix":"10.1093","volume":"32","author":[{"given":"Hesam","family":"Montazeri","sequence":"first","affiliation":[{"name":"Department of Biosystems Science and Engineering, ETH Zurich, Basel, Switzerland"},{"name":"SIB Swiss Institute of Bioinformatics, Basel, Switzerland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Jack","family":"Kuipers","sequence":"additional","affiliation":[{"name":"Department of Biosystems Science and Engineering, ETH Zurich, Basel, Switzerland"},{"name":"SIB Swiss Institute of Bioinformatics, Basel, Switzerland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Roger","family":"Kouyos","sequence":"additional","affiliation":[{"name":"Division of Infectious Diseases and Hospital Epidemiology, University Hospital Zurich, University of Zurich, Zurich, Switzerland"},{"name":"Institute of Medical Virology"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"J\u00fcrg","family":"B\u00f6ni","sequence":"additional","affiliation":[{"name":"Swiss National Center for Retroviruses, Institute of Medical Virology, University of Zurich, Zurich 8057, Switzerland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Sabine","family":"Yerly","sequence":"additional","affiliation":[{"name":"Laboratory of Virology, Division of Infectious Diseases, Geneva University Hospital, Geneva, Switzerland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Thomas","family":"Klimkait","sequence":"additional","affiliation":[{"name":"Department of Biomedicine-Petersplatz, University of Basel, Basel, Switzerland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Vincent","family":"Aubert","sequence":"additional","affiliation":[{"name":"Division of Immunology and Allergy, University Hospital Lausanne, Lausanne, Switzerland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Huldrych F","family":"G\u00fcnthard","sequence":"additional","affiliation":[{"name":"Division of Infectious Diseases and Hospital Epidemiology, University Hospital Zurich, University of Zurich, Zurich, Switzerland"},{"name":"Institute of Medical Virology"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Niko","family":"Beerenwinkel","sequence":"additional","affiliation":[{"name":"Department of Biosystems Science and Engineering, ETH Zurich, Basel, Switzerland"},{"name":"SIB Swiss Institute of Bioinformatics, Basel, Switzerland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"name":"The Swiss HIV Cohort Study","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2016,8,29]]},"reference":[{"key":"2023020113310869500_btw459-B1","doi-asserted-by":"crossref","first-page":"999","DOI":"10.1086\/597305","article-title":"Maurizio Zazzi, Rolf Kaiser, and Thomas Lengauer. Predicting the response to combination antiretroviral therapy: retrospective validation of geno2pheno-THEO on a large clinical database","volume":"199","author":"Altmann","year":"2009","journal-title":"J. Infect. Dis"},{"key":"2023020113310869500_btw459-B2","doi-asserted-by":"crossref","first-page":"645","DOI":"10.1093\/biomet\/asp023","article-title":"Markov models for accumulating mutations","volume":"96","author":"Beerenwinkel","year":"2009","journal-title":"Biometrika"},{"key":"2023020113310869500_btw459-B3","doi-asserted-by":"crossref","first-page":"3943","DOI":"10.1093\/bioinformatics\/bti654","article-title":"Computational methods for the design of effective therapies against drug resistant HIV strains","volume":"21","author":"Beerenwinkel","year":"2005","journal-title":"Bioinformatics"},{"key":"2023020113310869500_btw459-B4","doi-asserted-by":"crossref","first-page":"584","DOI":"10.1089\/cmb.2005.12.584","article-title":"Learning multiple evolutionary pathways from cross-sectional data","volume":"12","author":"Beerenwinkel","year":"2005","journal-title":"J. Comput. Biol"},{"key":"2023020113310869500_btw459-B5","doi-asserted-by":"crossref","first-page":"893","DOI":"10.3150\/07-BEJ6133","article-title":"Conjunctive Bayesian networks","volume":"13","author":"Beerenwinkel","year":"2007","journal-title":"Bernoulli"},{"key":"2023020113310869500_btw459-B6","doi-asserted-by":"crossref","first-page":"1","DOI":"10.2202\/1544-6115.1603","article-title":"Learning monotonic genotype\u2013phenotype maps","volume":"10","author":"Beerenwinkel","year":"2011","journal-title":"Stat. Appl. Genet. Mol. Biol"},{"key":"2023020113310869500_btw459-B7","doi-asserted-by":"crossref","first-page":"e1003203","DOI":"10.1371\/journal.pcbi.1003203","article-title":"The individualized genetic barrier predicts treatment response in a large cohort of HIV-1 infected patients","volume":"9","author":"Beerenwinkel","year":"2013","journal-title":"PLoS Comput. Biol"},{"key":"2023020113310869500_btw459-B8","volume-title":"Pattern Recognition and Machine Learning (Information Science and Statistics)","author":"Bishop","year":"2006"},{"key":"2023020113310869500_btw459-B9","doi-asserted-by":"crossref","DOI":"10.1017\/CBO9780511809088","volume-title":"Introduction to Lattices and Order","author":"Davey","year":"2002"},{"key":"2023020113310869500_btw459-B10","doi-asserted-by":"crossref","first-page":"399","DOI":"10.1177\/135965350801300316","article-title":"Modelled in vivo HIV fitness under drug selective pressure and estimated genetic barrier towards resistance are predictive for virological response","volume":"13","author":"Deforche","year":"2008","journal-title":"Antivir. Ther"},{"key":"2023020113310869500_btw459-B11","doi-asserted-by":"crossref","first-page":"37","DOI":"10.1089\/cmb.1999.6.37","article-title":"Inferring tree models for oncogenesis from comparative genome hybridization data","volume":"6","author":"Desper","year":"1999","journal-title":"J. Comput. Biol"},{"key":"2023020113310869500_btw459-B12","doi-asserted-by":"crossref","first-page":"1273","DOI":"10.1002\/nme.909","article-title":"Computation of the matrix exponential and its derivatives by scaling and squaring","volume":"59","author":"Fung","year":"2004","journal-title":"Int. J. Numer. Methods Eng"},{"key":"2023020113310869500_btw459-B13","doi-asserted-by":"crossref","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":"2023020113310869500_btw459-B14","doi-asserted-by":"crossref","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":"2023020113310869500_btw459-B15","doi-asserted-by":"crossref","first-page":"646","DOI":"10.1016\/j.cell.2011.02.013","article-title":"Hallmarks of cancer: the next generation","volume":"144","author":"Hanahan","year":"2011","journal-title":"Cell"},{"key":"2023020113310869500_btw459-B16","doi-asserted-by":"crossref","first-page":"545","DOI":"10.1093\/biostatistics\/kxh007","article-title":"Maximum likelihood estimation of oncogenetic tree models","volume":"5","author":"Heydebreck","year":"2004","journal-title":"Biostatistics (Oxford, England)"},{"key":"2023020113310869500_btw459-B17","doi-asserted-by":"crossref","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":"2023020113310869500_btw459-B18","doi-asserted-by":"crossref","first-page":"520.","DOI":"10.1186\/1471-2164-10-520","article-title":"The effect of sequencing errors on metagenomic gene prediction","volume":"10","author":"Hoff","year":"2009","journal-title":"BMC Genom"},{"key":"2023020113310869500_btw459-B19","first-page":"6503","article-title":"Construction of evolutionary tree models for renal cell carcinoma from comparative genomic hybridization data","volume":"60","author":"Jiang","year":"2000","journal-title":"Cancer Res"},{"key":"2023020113310869500_btw459-B20","doi-asserted-by":"crossref","first-page":"503","DOI":"10.1007\/BF01589116","article-title":"On the limited memory BFGS method for large scale optimization","volume":"45","author":"Liu","year":"1989","journal-title":"Math. Program"},{"key":"2023020113310869500_btw459-B21","doi-asserted-by":"crossref","first-page":"12025","DOI":"10.1073\/pnas.0905922106","article-title":"Stepwise acquisition of pyrimethamine resistance in the malaria parasite","volume":"106","author":"Lozovsky","year":"2009","journal-title":"Proc. Natl. Acad. Sci. U.S.A"},{"key":"2023020113310869500_btw459-B22","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1016\/j.cancergencyto.2003.12.019","article-title":"Dissecting karyotypic patterns in renal cell carcinoma: an analysis of the accumulated cytogenetic data","volume":"153","author":"Mattias","year":"2004","journal-title":"Cancer Genet. Cytogenet"},{"key":"2023020113310869500_btw459-B23","doi-asserted-by":"crossref","first-page":"924","DOI":"10.1038\/nrc2013","article-title":"Cancer as an evolutionary and ecological process","volume":"6","author":"Merlo","year":"2006","journal-title":"Nat. Rev. Cancer"},{"key":"2023020113310869500_btw459-B24","doi-asserted-by":"crossref","DOI":"10.1093\/biostatistics\/kxv019","article-title":"Estimating the dynamics and dependencies of accumulating mutations with applications to HIV drug resistance","author":"Montazeri","year":"2015","journal-title":"Biostatistics"},{"key":"2023020113310869500_btw459-B25","doi-asserted-by":"crossref","first-page":"457","DOI":"10.2307\/3318671","article-title":"The stochastic EM algorithm: estimation and asymptotic results","author":"Nielsen","year":"2000","journal-title":"Bernoulli"},{"key":"2023020113310869500_btw459-B26","doi-asserted-by":"crossref","first-page":"383","DOI":"10.1038\/nature05451","article-title":"Empirical fitness landscapes reveal accessible evolutionary paths","volume":"445","author":"Poelwijk","year":"2007","journal-title":"Nature"},{"key":"2023020113310869500_btw459-B27","doi-asserted-by":"crossref","first-page":"255","DOI":"10.1007\/0-387-30065-1_16","volume-title":"Large-Scale Nonlinear Optimization, Volume 83 of Nonconvex Optimization and Its Applications","author":"Powell","year":"2006"},{"key":"2023020113310869500_btw459-B28","doi-asserted-by":"crossref","first-page":"1040","DOI":"10.1093\/bioinformatics\/btn568","article-title":"Stochastic modelling of genotypic drug-resistance for human immunodeficiency virus towards long-term combination therapy optimization","volume":"25","author":"Prosperi","year":"2009","journal-title":"Bioinformatics"},{"key":"2023020113310869500_btw459-B29","doi-asserted-by":"crossref","first-page":"2438","DOI":"10.1093\/bioinformatics\/bti312","article-title":"Estimating cancer survival and clinical outcome based on genetic tumor progression scores","volume":"21","author":"Rahnenf\u00fchrer","year":"2005","journal-title":"Bioinformatics (Oxford, England)"},{"key":"2023020113310869500_btw459-B30","doi-asserted-by":"crossref","first-page":"2318","DOI":"10.1093\/bioinformatics\/bts433","article-title":"Efficient sampling for Bayesian inference of conjunctive Bayesian networks","volume":"28","author":"Sakoparnig","year":"2012","journal-title":"Bioinformatics"},{"key":"2023020113310869500_btw459-B31","doi-asserted-by":"crossref","first-page":"191","DOI":"10.1534\/genetics.114.172312","article-title":"A framework for inferring fitness landscapes of patient-derived viruses using quasispecies theory","volume":"199","author":"Seifert","year":"2015","journal-title":"Genetics"},{"key":"2023020113310869500_btw459-B32","volume-title":"Causation, Prediction, and Search","author":"Spirtes","year":"2000","edition":"2nd ed."},{"key":"2023020113310869500_btw459-B33","doi-asserted-by":"crossref","first-page":"1546","DOI":"10.1126\/science.1235122","article-title":"Cancer genome landscapes","volume":"339","author":"Vogelstein","year":"2013","journal-title":"Science"},{"key":"2023020113310869500_btw459-B34","doi-asserted-by":"crossref","first-page":"699","DOI":"10.1080\/01621459.1990.10474930","article-title":"A Monte Carlo implementation of the EM algorithm and the poor man\u2019s data augmentation algorithms","volume":"85","author":"Wei","year":"1990","journal-title":"J. Am. Stat. Assoc"},{"key":"2023020113310869500_btw459-B35","doi-asserted-by":"crossref","first-page":"111","DOI":"10.1126\/science.1123539","article-title":"Darwinian evolution can follow only very few mutational paths to fitter proteins","volume":"312","author":"Weinreich","year":"2006","journal-title":"Science"},{"key":"2023020113310869500_btw459-B36","first-page":"642","article-title":"Special contribution 2014 update of the drug resistance mutations in hiv-1","volume":"22","author":"Wensing","year":"2014","journal-title":"Topics Antivir. Med"},{"key":"2023020113310869500_btw459-B37","doi-asserted-by":"crossref","first-page":"4099","DOI":"10.1128\/JCM.37.12.4099-4106.1999","article-title":"Mutation patterns of the reverse transcriptase and protease genes in human immunodeficiency virus type 1-infected patients undergoing combination therapy: survey of 787 sequences","volume":"37","author":"Yahi","year":"1999","journal-title":"J. Clin. Microbiol"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/32\/17\/i727\/49023594\/bioinformatics_32_17_i727.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/32\/17\/i727\/49023594\/bioinformatics_32_17_i727.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,2,1]],"date-time":"2023-02-01T23:33:38Z","timestamp":1675294418000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/32\/17\/i727\/2450789"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2016,8,29]]},"references-count":37,"journal-issue":{"issue":"17","published-print":{"date-parts":[[2016,9,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btw459","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2016,9,1]]},"published":{"date-parts":[[2016,8,29]]}}}