{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,5,21]],"date-time":"2026-05-21T04:25:38Z","timestamp":1779337538529,"version":"3.51.4"},"reference-count":19,"publisher":"Oxford University Press (OUP)","issue":"21","content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2016,11,1]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:p>Motivation: Cell differentiation is steered by extracellular signals that activate a cell type specific transcriptional program. Molecular mechanisms that drive the differentiation can be analyzed by combining mathematical modeling with population average data. For standard mathematical models, the population average data is informative only if the measurements come from a homogeneous cell culture. In practice, however, the differentiation efficiencies are always imperfect. Consequently, cell cultures are inherently mixtures of several cell types, which have different molecular mechanisms and exhibit quantitatively different dynamics. There is an urgent need for data-driven mathematical modeling approaches that can detect possible heterogeneity and, further, recover the molecular mechanisms from heterogeneous data.<\/jats:p>\n               <jats:p>Results: We develop a novel method that models a heterogeneous population using homogeneous subpopulations that evolve in parallel. Different subpopulations can represent different cell types and each subpopulation can have cell type specific molecular mechanisms. We present statistical methodology that can be used to quantify the effect of heterogeneity and to infer the subpopulation specific molecular interactions. After a proof of principle study with simulated data, we apply our methodology to analyze the differentiation of human Th17 cells using time-course RNA sequencing data. We construct putative molecular networks driving the T cell activation and Th17 differentiation and allow the cell populations to be split into two subpopulations in the case of heterogeneous samples. Our analysis shows that the heterogeneity indeed has a statistically significant effect on observed dynamics and, furthermore, our statistical methodology can infer both the subpopulation specific molecular mechanisms and the effect of heterogeneity.<\/jats:p>\n               <jats:p>Availability and Implementation: An implementation of the method is available at http:\/\/research.ics.aalto.fi\/csb\/software\/subpop\/ .<\/jats:p>\n               <jats:p>Contact: \u00a0jukka.intosalmi@aalto.fi or harri.lahdesmaki@aalto.fi<\/jats:p>\n               <jats:p>Supplementary information: \u00a0Supplementary data are available at Bioinformatics online.<\/jats:p>","DOI":"10.1093\/bioinformatics\/btw395","type":"journal-article","created":{"date-parts":[[2016,7,12]],"date-time":"2016-07-12T07:21:31Z","timestamp":1468308091000},"page":"3306-3313","source":"Crossref","is-referenced-by-count":6,"title":["A subpopulation model to analyze heterogeneous cell differentiation dynamics"],"prefix":"10.1093","volume":"32","author":[{"given":"Yat Hin","family":"Chan","sequence":"first","affiliation":[{"name":"1 Department of Computer Science, Aalto University, 00076 Aalto, Finland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Jukka","family":"Intosalmi","sequence":"additional","affiliation":[{"name":"1 Department of Computer Science, Aalto University, 00076 Aalto, Finland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Sini","family":"Rautio","sequence":"additional","affiliation":[{"name":"1 Department of Computer Science, Aalto University, 00076 Aalto, Finland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Harri","family":"L\u00e4hdesm\u00e4ki","sequence":"additional","affiliation":[{"name":"1 Department of Computer Science, Aalto University, 00076 Aalto, Finland"},{"name":"2 Turku Centre for Biotechnology, University of Turku and \u00c5bo Akademi, 20521 Turku, Finland"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2016,7,10]]},"reference":[{"key":"2023020113513670900_btw395-B1","doi-asserted-by":"crossref","first-page":"166","DOI":"10.1093\/bioinformatics\/btu638","article-title":"HTSeq \u2013 a Python framework to work with high-throughput sequencing data","volume":"31","author":"Anders","year":"2015","journal-title":"Bioinformatics"},{"key":"2023020113513670900_btw395-B2","doi-asserted-by":"crossref","first-page":"235","DOI":"10.1038\/nature04753","article-title":"Reciprocal developmental pathways for the generation of pathogenic effector TH17 and regulatory T cells","volume":"441","author":"Bettelli","year":"2006","journal-title":"Nature"},{"key":"2023020113513670900_btw395-B3","doi-asserted-by":"crossref","first-page":"4028","DOI":"10.1016\/j.csda.2009.07.025","article-title":"Estimating Bayes factors via thermodynamic integration and population MCMC","volume":"53","author":"Calderhead","year":"2009","journal-title":"Comput. Stat. Data Anal"},{"key":"2023020113513670900_btw395-B4","doi-asserted-by":"crossref","first-page":"289","DOI":"10.1016\/j.cell.2012.09.016","article-title":"A validated regulatory network for Th17 cell specification","volume":"151","author":"Ciofani","year":"2012","journal-title":"Cell"},{"key":"2023020113513670900_btw395-B5","doi-asserted-by":"crossref","first-page":"709","DOI":"10.1007\/s11222-013-9397-1","article-title":"Improving power posterior estimation of statistical evidence","volume":"24","author":"Friel","year":"2014","journal-title":"Stat. Comput"},{"key":"2023020113513670900_btw395-B6","doi-asserted-by":"crossref","first-page":"589","DOI":"10.1111\/j.1467-9868.2007.00650.x","article-title":"Marginal likelihood estimation via power posteriors","volume":"70","author":"Friel","year":"2008","journal-title":"J. R. Stat. Soc. Ser. B Stat. Methodol"},{"key":"2023020113513670900_btw395-B7","doi-asserted-by":"crossref","DOI":"10.1201\/b16018","volume-title":"Bayesian Data Analysis","author":"Gelman","year":"2013","edition":"3rd edn."},{"key":"2023020113513670900_btw395-B8","doi-asserted-by":"crossref","first-page":"125","DOI":"10.1186\/1471-2105-12-125","article-title":"Identification of models of heterogeneous cell populations from population snapshot data","volume":"12","author":"Hasenauer","year":"2011","journal-title":"BMC Bioinformatics"},{"key":"2023020113513670900_btw395-B9","doi-asserted-by":"crossref","first-page":"e1003686","DOI":"10.1371\/journal.pcbi.1003686","article-title":"ODE constrained mixture modelling: a method for unraveling subpopulation structures and dynamics","volume":"10","author":"Hasenauer","year":"2014","journal-title":"PLoS Comput. Biol"},{"key":"2023020113513670900_btw395-B10","doi-asserted-by":"crossref","first-page":"81","DOI":"10.1186\/s12918-015-0223-6","article-title":"Analyzing Th17 cell differentiation dynamics using a novel integrative modeling framework for time-course RNA sequencing data","volume":"9","author":"Intosalmi","year":"2015","journal-title":"BMC Syst. Biol"},{"key":"2023020113513670900_btw395-B11","doi-asserted-by":"crossref","first-page":"263","DOI":"10.1007\/s11222-007-9028-9","article-title":"On population-based simulation for static inference","volume":"17","author":"Jasra","year":"2007","journal-title":"Stat. Comput"},{"key":"2023020113513670900_btw395-B12","doi-asserted-by":"crossref","first-page":"485","DOI":"10.1146\/annurev.immunol.021908.132710","article-title":"IL-17 and Th17 Cells","volume":"27","author":"Korn","year":"2009","journal-title":"Annu. Rev. Immunol"},{"key":"2023020113513670900_btw395-B13","doi-asserted-by":"crossref","first-page":"674","DOI":"10.1038\/ni.1899","article-title":"Effector T cell plasticity: flexibility in the face of changing circumstances","volume":"11","author":"Murphy","year":"2010","journal-title":"Nat. Immunol"},{"key":"2023020113513670900_btw395-B14","doi-asserted-by":"crossref","first-page":"139","DOI":"10.1093\/bioinformatics\/btp616","article-title":"edgeR: a Bioconductor package for differential expression analysis of digital gene expression data","volume":"26","author":"Robinson","year":"2010","journal-title":"Bioinformatics"},{"key":"2023020113513670900_btw395-B15","doi-asserted-by":"crossref","first-page":"673","DOI":"10.1016\/j.immuni.2009.03.013","article-title":"Sequential polarization and imprinting of type 1 T helper lymphocytes by interferon-gamma and interleukin-12","volume":"30","author":"Schulz","year":"2009","journal-title":"Immunity"},{"key":"2023020113513670900_btw395-B16","doi-asserted-by":"crossref","first-page":"13416","DOI":"10.18632\/oncotarget.7963","article-title":"Comparative analysis of human and mouse transcriptomes of Th17 cell priming","volume":"7","author":"Tuomela","year":"2016","journal-title":"Oncotarget"},{"key":"2023020113513670900_btw395-B17","doi-asserted-by":"crossref","first-page":"833","DOI":"10.1093\/bioinformatics\/btm607","article-title":"Bayesian ranking of biochemical system models","volume":"24","author":"Vyshemirsky","year":"2008","journal-title":"Bioinformatics"},{"key":"2023020113513670900_btw395-B18","first-page":"134","article-title":"Inferring signaling pathway topologies from multiple perturbation measurements of specific biochemical species","volume":"3","author":"Xu","year":"2010","journal-title":"Sci. Sig"},{"key":"2023020113513670900_btw395-B19","doi-asserted-by":"crossref","first-page":"461","DOI":"10.1038\/nature11981","article-title":"Dynamic regulatory network controlling TH17 cell differentiation","volume":"496","author":"Yosef","year":"2013","journal-title":"Nature"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/32\/21\/3306\/49021529\/bioinformatics_32_21_3306.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/32\/21\/3306\/49021529\/bioinformatics_32_21_3306.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,2,1]],"date-time":"2023-02-01T23:54:23Z","timestamp":1675295663000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/32\/21\/3306\/2415099"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2016,7,10]]},"references-count":19,"journal-issue":{"issue":"21","published-print":{"date-parts":[[2016,11,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btw395","relation":{},"ISSN":["1367-4811","1367-4803"],"issn-type":[{"value":"1367-4811","type":"electronic"},{"value":"1367-4803","type":"print"}],"subject":[],"published-other":{"date-parts":[[2016,11,1]]},"published":{"date-parts":[[2016,7,10]]}}}