{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,6,28]],"date-time":"2026-06-28T01:06:21Z","timestamp":1782608781029,"version":"3.54.5"},"reference-count":52,"publisher":"Oxford University Press (OUP)","issue":"19","license":[{"start":{"date-parts":[[2019,3,9]],"date-time":"2019-03-09T00:00:00Z","timestamp":1552089600000},"content-version":"vor","delay-in-days":1,"URL":"http:\/\/creativecommons.org\/licenses\/by-nc\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100002341","name":"Academy of Finland","doi-asserted-by":"publisher","award":["317255"],"award-info":[{"award-number":["317255"]}],"id":[{"id":"10.13039\/501100002341","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2019,10,1]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:sec><jats:title>Motivation<\/jats:title><jats:p>Recent advances in high dimensional phenotyping bring time as an extra dimension into the phenotypes. This promotes the quantitative trait locus (QTL) studies of function-valued traits such as those related to growth and development. Existing approaches for analyzing functional traits utilize either parametric methods or semi-parametric approaches based on splines and wavelets. However, very limited choices of software tools are currently available for practical implementation of functional QTL mapping and variable selection.<\/jats:p><\/jats:sec><jats:sec><jats:title>Results<\/jats:title><jats:p>We propose a Bayesian Gaussian process (GP) approach for functional QTL mapping. We use GPs to model the continuously varying coefficients which describe how the effects of molecular markers on the quantitative trait are changing over time. We use an efficient gradient based algorithm to estimate the tuning parameters of GPs. Notably, the GP approach is directly applicable to the incomplete datasets having even larger than 50% missing data rate (among phenotypes). We further develop a stepwise algorithm to search through the model space in terms of genetic variants, and use a minimal increase of Bayesian posterior probability as a stopping rule to focus on only a small set of putative QTL. We also discuss the connection between GP and penalized B-splines and wavelets. On two simulated and three real datasets, our GP approach demonstrates great flexibility for modeling different types of phenotypic trajectories with low computational cost. The proposed model selection approach finds the most likely QTL reliably in tested datasets.<\/jats:p><\/jats:sec><jats:sec><jats:title>Availability and implementation<\/jats:title><jats:p>Software and simulated data are available as a MATLAB package \u2018GPQTLmapping\u2019, and they can be downloaded from GitHub (https:\/\/github.com\/jpvanhat\/GPQTLmapping). Real datasets used in case studies are publicly available at QTL Archive.<\/jats:p><\/jats:sec><jats:sec><jats:title>Supplementary information<\/jats:title><jats:p>Supplementary data are available at Bioinformatics online.<\/jats:p><\/jats:sec>","DOI":"10.1093\/bioinformatics\/btz164","type":"journal-article","created":{"date-parts":[[2019,3,7]],"date-time":"2019-03-07T04:22:53Z","timestamp":1551932573000},"page":"3684-3692","source":"Crossref","is-referenced-by-count":15,"title":["A Gaussian process model and Bayesian variable selection for mapping function-valued quantitative traits with incomplete phenotypic data"],"prefix":"10.1093","volume":"35","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-6831-0211","authenticated-orcid":false,"given":"Jarno","family":"Vanhatalo","sequence":"first","affiliation":[{"name":"Department of Mathematics and Statistics and Organismal and Evolutionary Biology Research Programme, University of Helsinki , Helsinki, Finland"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-8469-7295","authenticated-orcid":false,"given":"Zitong","family":"Li","sequence":"additional","affiliation":[{"name":"CSIRO Agriculture & Food , GPO Box 1600, Canberra, ACT 2601, Australia"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-2808-2768","authenticated-orcid":false,"given":"Mikko J","family":"Sillanp\u00e4\u00e4","sequence":"additional","affiliation":[{"name":"Department of Mathematical Sciences, Biocenter Oulu and Infotech Oulu University of Oulu , Oulu FI-90014, Finland"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2019,3,8]]},"reference":[{"key":"2023013108143166500_btz164-B1","doi-asserted-by":"crossref","first-page":"722","DOI":"10.1002\/gepi.20623","article-title":"Stability selection for genome-wide asssociation","volume":"35","author":"Alexander","year":"2011","journal-title":"Genet. Epidemiol"},{"key":"2023013108143166500_btz164-B2","doi-asserted-by":"crossref","first-page":"113","DOI":"10.1093\/bioinformatics\/btu274","article-title":"Methods for time series analysis of RNA-seq data with application to human Th17 cell differentiation","volume":"30","author":"\u00c4ij\u00f6","year":"2014","journal-title":"Bioinformatics"},{"key":"2023013108143166500_btz164-B3","first-page":"1745","author":"Andersen","year":"2014"},{"key":"2023013108143166500_btz164-B4","first-page":"1","article-title":"Bayesian inference for spatio-temporal spike-and-slab priors","volume":"18","author":"Andersen","year":"2017","journal-title":"J. Mach. Learn. Res"},{"key":"2023013108143166500_btz164-B5","doi-asserted-by":"crossref","first-page":"1493","DOI":"10.1093\/bioinformatics\/btw018","article-title":"FINEMAP: efficient variable selection using summary data from genome-wide association studies","volume":"32","author":"Benner","year":"2016","journal-title":"Bioinformatics"},{"key":"2023013108143166500_btz164-B6","doi-asserted-by":"crossref","first-page":"859","DOI":"10.1080\/01621459.2017.1285773","article-title":"Varitional inference: a review for statisticians","volume":"112","author":"Blei","year":"2017","journal-title":"J. Am. Stat. Assoc"},{"key":"2023013108143166500_btz164-B7","doi-asserted-by":"crossref","first-page":"S1.","DOI":"10.1186\/1753-6561-4-S1-S3","article-title":"QTLMAS 2009: simulated dataset","volume":"4","author":"Coster","year":"2010","journal-title":"BMC Proc"},{"key":"2023013108143166500_btz164-B8","doi-asserted-by":"crossref","DOI":"10.1137\/1.9780898718881","volume-title":"Direct Methods for Sparse Linear Systems","author":"Davis","year":"2006"},{"key":"2023013108143166500_btz164-B9","doi-asserted-by":"crossref","first-page":"856","DOI":"10.1002\/gepi.21673","article-title":"Longitudinal association analysis of quantitative traits","volume":"36","author":"Fan","year":"2012","journal-title":"Genet. Epidemiol"},{"key":"2023013108143166500_btz164-B10","doi-asserted-by":"crossref","DOI":"10.1093\/acprof:oso\/9780199533022.001.0001","volume-title":"Bayesian Smoothing and Regression for Longitudinal, Spatial and Event History Data","author":"Fahrmeir","year":"2011"},{"key":"2023013108143166500_btz164-B11","doi-asserted-by":"crossref","first-page":"213","DOI":"10.1534\/genetics.115.177790","article-title":"Genetics of rapid and extreme size evolution in island mice","volume":"201","author":"Gray","year":"2015","journal-title":"Genetics"},{"key":"2023013108143166500_btz164-B12","doi-asserted-by":"crossref","first-page":"82.","DOI":"10.1186\/1471-2156-9-82","article-title":"Multiple-trait quantitative trait locus mapping with incomplete phenotypic data","volume":"9","author":"Guo","year":"2008","journal-title":"BMC Genetics"},{"key":"2023013108143166500_btz164-B13","doi-asserted-by":"crossref","first-page":"315","DOI":"10.1038\/hdy.1992.131","article-title":"A simple regression method for mapping quantitative trait loci in line crosses using flanking markers","volume":"69","author":"Haley","year":"1992","journal-title":"Heredity"},{"key":"2023013108143166500_btz164-B14","doi-asserted-by":"crossref","first-page":"1375","DOI":"10.1214\/17-AOAS1029","article-title":"Gaussian process framework for temporal dependence and discrepancy functions in Ricker-type population growth models","volume":"11","author":"Hartmann","year":"2017","journal-title":"Ann. Appl. Stat"},{"key":"2023013108143166500_btz164-B15","doi-asserted-by":"crossref","DOI":"10.1007\/b98818","volume-title":"Matrix Algebra From a Statistician\u2019s Perspective","author":"Harville","year":"1997"},{"key":"2023013108143166500_btz164-B16","doi-asserted-by":"crossref","first-page":"1026","DOI":"10.1093\/bioinformatics\/btr057","article-title":"tigre: transcription factor inference through Gaussian process reconstruction of expression for bioconductor","volume":"27","author":"Honkela","year":"2011","journal-title":"Bioinformatics"},{"key":"2023013108143166500_btz164-B17","doi-asserted-by":"crossref","first-page":"13115","DOI":"10.1073\/pnas.1420404112","article-title":"Genome-wide modelling of transcription kinetics reveals patterns of RNA production delays","volume":"112","author":"Honkela","year":"2015","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023013108143166500_btz164-B18","doi-asserted-by":"crossref","first-page":"913","DOI":"10.1093\/genetics\/156.2.913","article-title":"Statistical models for estimating the genetic basis of repeated measures and other function-valued traits","volume":"156","author":"Jaffr\u00e9zic","year":"2000","journal-title":"Genetics"},{"key":"2023013108143166500_btz164-B19","doi-asserted-by":"crossref","first-page":"1409","DOI":"10.1534\/genetics.114.166306","article-title":"A simple regression-based method to map quantitative trait loci underlying function-valued phenotypes","volume":"197","author":"Kwak","year":"2014","journal-title":"Genetics"},{"key":"2023013108143166500_btz164-B20","doi-asserted-by":"crossref","first-page":"79","DOI":"10.1534\/g3.115.024133","article-title":"Mapping quantitative trait loci underlying function-valued traits using functional principal component analysis and multi-trait mapping","volume":"6","author":"Kwak","year":"2016","journal-title":"G3 (Bethesda)"},{"key":"2023013108143166500_btz164-B21","doi-asserted-by":"crossref","first-page":"997","DOI":"10.1534\/genetics.113.152736","article-title":"A Bayesian nonparametric approach for mapping dynamic quantitative traits","volume":"194","author":"Li","year":"2013","journal-title":"Genetics"},{"key":"2023013108143166500_btz164-B22","doi-asserted-by":"crossref","first-page":"2365","DOI":"10.1534\/g3.114.014068","article-title":"Functional multi-locus QTL mapping of temporal trends in scots pine wood traits","volume":"4","author":"Li","year":"2014","journal-title":"G3"},{"key":"2023013108143166500_btz164-B23","doi-asserted-by":"crossref","first-page":"822","DOI":"10.1016\/j.tplants.2015.08.012","article-title":"Dynamic quantitative trait locus analysis of plant phenomic data","volume":"20","author":"Li","year":"2015","journal-title":"Trends Plant Sci"},{"key":"2023013108143166500_btz164-B24","doi-asserted-by":"crossref","first-page":"640","DOI":"10.1214\/15-AOAS808","article-title":"Bayesian group LASSO for nonparametric varying-coefficient models with application to functional genome-wide studies","volume":"9","author":"Li","year":"2015","journal-title":"Ann. Appl. Stat"},{"key":"2023013108143166500_btz164-B25","doi-asserted-by":"crossref","first-page":"286","DOI":"10.1111\/tpj.13777","article-title":"Two-stage identification of SNP effects on dynamic poplar growth","volume":"93","author":"Liu","year":"2018","journal-title":"Plant J"},{"key":"2023013108143166500_btz164-B26","doi-asserted-by":"crossref","first-page":"266","DOI":"10.1080\/01621459.2013.850086","article-title":"Feature selection for varying coefficient models with ultrahigh dimensional covariates","volume":"109","author":"Liu","year":"2014","journal-title":"J. Am. Stat. Assoc"},{"key":"2023013108143166500_btz164-B27","doi-asserted-by":"crossref","first-page":"1751","DOI":"10.1093\/genetics\/161.4.1751","article-title":"Functional mapping of quantitative trait loci underlying the character process: a theoretical framework","volume":"161","author":"Ma","year":"2002","journal-title":"Genetics"},{"key":"2023013108143166500_btz164-B28","doi-asserted-by":"crossref","first-page":"417","DOI":"10.1111\/j.1467-9868.2010.00740.x","article-title":"Stability selection","volume":"72","author":"Meinshausen","year":"2010","journal-title":"J. R. Stat. Soc. Series B"},{"key":"2023013108143166500_btz164-B29","doi-asserted-by":"crossref","DOI":"10.1017\/CBO9780511977176","volume-title":"Numerical Methods of Statistics","author":"Monahan","year":"2011"},{"key":"2023013108143166500_btz164-B30","doi-asserted-by":"crossref","first-page":"1077","DOI":"10.1534\/genetics.113.153346","article-title":"High-throughput computer vision introduces the time axis to a quantitative trait map of a plant growth response","volume":"195","author":"Moore","year":"2013","journal-title":"Genetics"},{"key":"2023013108143166500_btz164-B31","doi-asserted-by":"crossref","first-page":"199","DOI":"10.1002\/wics.199","article-title":"The Bayesian information criterion: background, derivation, and applications","volume":"4","author":"Neath","year":"2012","journal-title":"WIREs Comput. Stat"},{"key":"2023013108143166500_btz164-B32","doi-asserted-by":"crossref","first-page":"e0164766.","DOI":"10.1371\/journal.pone.0164766","article-title":"RNA-Seq count data modelling by grey relational analysis and nonparametric Gaussian process","volume":"11","author":"Nguyen","year":"2016","journal-title":"PLoS One"},{"key":"2023013108143166500_btz164-B33","doi-asserted-by":"crossref","first-page":"590.","DOI":"10.1038\/s41598-017-00638-2","article-title":"Perfomance gains in genome-wide association studies for longitudinal traits via modeling time-varied effects","volume":"7","author":"Ning","year":"2017","journal-title":"Sci. Rep"},{"key":"2023013108143166500_btz164-B34","doi-asserted-by":"crossref","first-page":"797","DOI":"10.1080\/10618600.2012.679897","article-title":"Regression density estimation with variational methods and stochastic approximation","volume":"21","author":"Nott","year":"2013","journal-title":"J. Comput. Graph. Stat"},{"key":"2023013108143166500_btz164-B35","volume-title":"Kendals Advanced Theory of Statistics, Volume 2B: Bayesian Inference","author":"O\u2019Hagan","year":"2004","edition":"2nd edn"},{"key":"2023013108143166500_btz164-B36","first-page":"85","article-title":"A review of Bayesian variable selection methods: what, how and which","volume":"4","author":"O\u2019Hara","year":"2009","journal-title":"Bayesian Anal"},{"key":"2023013108143166500_btz164-B37","doi-asserted-by":"crossref","first-page":"919","DOI":"10.1038\/ng.3609","article-title":"Genome-wide association study of behavioral, physiological and gene expression traits in commercially available outbred CFW mice","volume":"48","author":"Parker","year":"2016","journal-title":"Nat. Genet"},{"key":"2023013108143166500_btz164-B38","doi-asserted-by":"crossref","first-page":"711","DOI":"10.1007\/s11222-016-9649-y","article-title":"Comparison of Bayesian predictive methods for model selection","volume":"27","author":"Piironen","year":"2017","journal-title":"Stat. Comput"},{"key":"2023013108143166500_btz164-B39","doi-asserted-by":"crossref","first-page":"825","DOI":"10.1093\/genetics\/153.2.825","article-title":"The genetic analysis of age-dependent traits: modelling the character process","volume":"153","author":"Pletcher","year":"1999","journal-title":"Genetics"},{"key":"2023013108143166500_btz164-B40","volume-title":"Gaussian Processes for Machine Learning","author":"Rasmussen","year":"2006"},{"key":"2023013108143166500_btz164-B41","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1214\/16-STS576","article-title":"Penalising model component complexity: a principled, practical approach to constructing priors","volume":"32","author":"Simpson","year":"2017","journal-title":"Stat. Sci"},{"key":"2023013108143166500_btz164-B42","doi-asserted-by":"crossref","first-page":"1324","DOI":"10.1111\/j.1365-294X.2011.05021.x","article-title":"On statistical methods for estimating heritability in wild populations","volume":"20","author":"Sillanp\u00e4\u00e4","year":"2011","journal-title":"Mol. Ecol"},{"key":"2023013108143166500_btz164-B43","doi-asserted-by":"crossref","first-page":"134","DOI":"10.1038\/hdy.2011.56","article-title":"Simultaneous estimation of multiple quantitative trait loci and growth curve parameters through hierarchical Bayesian modeling","volume":"108","author":"Sillanp\u00e4\u00e4","year":"2012","journal-title":"Heredity"},{"key":"2023013108143166500_btz164-B44","doi-asserted-by":"crossref","first-page":"1580","DOI":"10.1002\/sim.3895","article-title":"Approximate inference for disease mapping with sparse Gaussian processes","volume":"2010","author":"Vanhatalo","year":"2010","journal-title":"Stat. Med"},{"key":"2023013108143166500_btz164-B45","first-page":"1175","article-title":"GPstuff: Bayesian modeling with Gaussian processes","volume":"14","author":"Vanhatalo","year":"2013","journal-title":"J. Mach. Learn. Res"},{"key":"2023013108143166500_btz164-B46","doi-asserted-by":"crossref","first-page":"229","DOI":"10.1038\/nrg1804","article-title":"Functional mapping\u2014how to map and study the genetic architecture of dynamical complex traits","volume":"7","author":"Wu","year":"2006","journal-title":"Nat. Rev. Genet"},{"key":"2023013108143166500_btz164-B47","doi-asserted-by":"crossref","first-page":"305","DOI":"10.1534\/genetics.111.129221","article-title":"A flexible estimating equations approach for mapping function valued traits","volume":"189","author":"Xiong","year":"2011","journal-title":"Genetics"},{"key":"2023013108143166500_btz164-B48","doi-asserted-by":"crossref","first-page":"1103","DOI":"10.1534\/genetics.113.155309","article-title":"Genetic mapping and genomic selection using recombination breakpoint data","volume":"195","author":"Xu","year":"2013","journal-title":"Genetics"},{"key":"2023013108143166500_btz164-B49","doi-asserted-by":"crossref","first-page":"30","DOI":"10.1111\/j.1541-0420.2008.01063.x","article-title":"Nonparametric functional mapping of quantitative trait loci","volume":"65","author":"Yang","year":"2009","journal-title":"Biometrics"},{"key":"2023013108143166500_btz164-B50","doi-asserted-by":"crossref","first-page":"1169","DOI":"10.1534\/genetics.106.064279","article-title":"Bayesian shrinkage analysis of quantitative trait loci for dynamic traits","volume":"176","author":"Yang","year":"2007","journal-title":"Genetics"},{"key":"2023013108143166500_btz164-B51","doi-asserted-by":"crossref","first-page":"250","DOI":"10.1198\/016214504000000241","article-title":"Inconsistent estimation and asymptotically equal interpolations in model-Based geostatistics","volume":"99","author":"Zhang","year":"2004","journal-title":"J. Am. Stat. Assoc"},{"key":"2023013108143166500_btz164-B52","doi-asserted-by":"crossref","first-page":"385","DOI":"10.1534\/genetics.109.113688","article-title":"Nonparametric Bayesian variable selection with applications to multiple quantitative trait loci mapping with epistasis and gene-environment interaction","volume":"186","author":"Zou","year":"2010","journal-title":"Genetics"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"http:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btz164\/28696096\/btz164.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/35\/19\/3684\/48976302\/bioinformatics_35_19_3684.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/35\/19\/3684\/48976302\/bioinformatics_35_19_3684.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,7,15]],"date-time":"2024-07-15T12:53:34Z","timestamp":1721048014000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/35\/19\/3684\/5372341"}},"subtitle":[],"editor":[{"given":"Russell","family":"Schwartz","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2019,3,8]]},"references-count":52,"journal-issue":{"issue":"19","published-print":{"date-parts":[[2019,10,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btz164","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2019,10,1]]},"published":{"date-parts":[[2019,3,8]]}}}