{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,10,21]],"date-time":"2026-10-21T14:43:44Z","timestamp":1792593824092,"version":"4.5.0"},"reference-count":50,"publisher":"Oxford University Press (OUP)","issue":"1","license":[{"start":{"date-parts":[[2026,1,1]],"date-time":"2026-01-01T00:00:00Z","timestamp":1767225600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/pages\/standard-publication-reuse-rights"}],"funder":[{"DOI":"10.13039\/100000072","name":"National Institute of Dental and Craniofacial Research","doi-asserted-by":"publisher","award":["R21DE026872"],"award-info":[{"award-number":["R21DE026872"]}],"id":[{"id":"10.13039\/100000072","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000057","name":"National Institute of General Medical Sciences","doi-asserted-by":"publisher","award":["R01GM126257"],"award-info":[{"award-number":["R01GM126257"]}],"id":[{"id":"10.13039\/100000057","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2026,1,20]]},"abstract":"<jats:title>Summary<\/jats:title>\n                  <jats:p>Advances in cellular imaging technologies, especially those based on fluorescence in situ hybridization (FISH), now allow detailed visualization of the spatial organization of human or bacterial cells. Quantifying this spatial organization is crucial for understanding the function of multicellular tissues or biofilms, with implications for human health and disease. To address the need for better methods to achieve such quantification, we propose a flexible multivariate point process model that characterizes and estimates complex spatial interactions among multiple cell types. The proposed Bayesian framework is appealing due to its unified estimation process and the ability to directly quantify uncertainty in key estimates of interest, such as those of inter-type correlation and the proportion of variance due to inter-type relationships. To ensure stable and interpretable estimation, we consider shrinkage priors for coefficients associated with latent processes that induce dependencies among point patterns. Model selection and comparison are conducted using a deviance information criterion designed for models with latent variables, providing a practical criterion for balancing model fit and complexity. Furthermore, we use a Bayesian hierarchical pooling model to synthesize image-specific posterior summaries, allowing inference at both the global- (across subjects) and subject-specific levels. An R package, mspatPPM, implements an efficient computational scheme based on Hamiltonian Monte Carlo and adaptive Metropolis-Hastings algorithms. Numerical studies evaluate the practical performance of the proposed framework for model selection and for estimating quantities that characterize the multivariate spatial distribution of cell types. We apply the proposed method to microbial biofilm image data from the human tongue dorsum and find that specific taxon pairs, such as Streptococcus mitis-Streptococcus salivarius and S. mitis-Veillonella, exhibit strong positive spatial correlations, while others, such as Actinomyces-Rothia, show slight negative correlations. For most of the taxa, a substantial portion of spatial variance can be attributed to inter-taxon relationships.<\/jats:p>","DOI":"10.1093\/biostatistics\/kxag028","type":"journal-article","created":{"date-parts":[[2026,7,15]],"date-time":"2026-07-15T11:47:41Z","timestamp":1784116061000},"source":"Crossref","is-referenced-by-count":0,"title":["A Bayesian multivariate spatial point pattern model: application to oral microbiome FISH image data"],"prefix":"10.1093","volume":"27","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-5313-8375","authenticated-orcid":false,"given":"Kyu Ha","family":"Lee","sequence":"first","affiliation":[{"name":"Harvard T.H. Chan School of Public Health Department of Nutrition, , 677 Huntington Avenue , Boston, MA 02115,","place":["United States"]},{"name":"Harvard T.H. Chan School of Public Health Department of Epidemiology, , 677 Huntington Avenue , Boston, MA 02115,","place":["United States"]},{"name":"Harvard T.H. Chan School of Public Health Department of Biostatistics, , 677 Huntington Avenue , Boston, MA 02115,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-1808-4156","authenticated-orcid":false,"given":"Brent A","family":"Coull","sequence":"additional","affiliation":[{"name":"Harvard T.H. Chan School of Public Health Department of Biostatistics, , 677 Huntington Avenue , Boston, MA 02115,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-3678-5642","authenticated-orcid":false,"given":"Suman","family":"Majumder","sequence":"additional","affiliation":[{"name":"Indian Statistical Institute Interdisciplinary Statistical Research Unit, , 203 Barrackpore Trunk Road , Kolkata-700108,","place":["India"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Patrick J","family":"La Riviere","sequence":"additional","affiliation":[{"name":"University of Chicago Department of Radiology, , 5841 South Maryland Avenue, MC2026 , Chicago, IL 60637,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Jessica L","family":"Mark Welch","sequence":"additional","affiliation":[{"name":"American Dental Association Forsyth Institute Department of Microbiology, , 100 Chestnut Street , Somerville, MA 02143,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Jacqueline R","family":"Starr","sequence":"additional","affiliation":[{"name":"Brigham and Women\u2019s Hospital Channing Division of Network Medicine, , 399 Revolution Drive , Somerville, MA 02145,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2026,8,18]]},"reference":[{"key":"2026102110221891500_kxag028-B1","doi-asserted-by":"publisher","first-page":"99","DOI":"10.1111\/j.2517-6161.1974.tb00989.x","article-title":"Scale mixtures of normal distributions","volume":"36","author":"Andrews","year":"1974","journal-title":"J R Stat Soc Ser B (Methodol)"},{"key":"2026102110221891500_kxag028-B2","doi-asserted-by":"publisher","first-page":"5283","DOI":"10.1038\/s41467-021-25557-9","article-title":"Computational challenges and opportunities in spatially resolved transcriptomic data analysis","volume":"12","author":"Atta","year":"2021","journal-title":"Nat Commun"},{"key":"2026102110221891500_kxag028-B3","doi-asserted-by":"publisher","first-page":"825","DOI":"10.1111\/j.1467-9868.2008.00663.x","article-title":"Gaussian predictive process models for large spatial data sets","volume":"70","author":"Banerjee","year":"2008","journal-title":"J R Stat Soc Ser B Stat Methodol"},{"key":"2026102110221891500_kxag028-B4","doi-asserted-by":"publisher","first-page":"471","DOI":"10.1111\/1467-9469.00249","article-title":"Space-time multi type log Gaussian Cox processes with a view to modelling weeds","volume":"28","author":"Brix","year":"2001","journal-title":"Scand J Stat"},{"key":"2026102110221891500_kxag028-B5","doi-asserted-by":"publisher","first-page":"651","DOI":"10.1214\/06-BA122","article-title":"Deviation information criteria for missing data models","volume":"1","author":"Celeux","year":"2006","journal-title":"Bayesian Anal"},{"key":"2026102110221891500_kxag028-B6","doi-asserted-by":"publisher","first-page":"183","DOI":"10.1016\/0304-4076(94)90063-9","article-title":"Bayes inference in regression models with ARMA (p, q) errors","volume":"64","author":"Chib","year":"1994","journal-title":"J Econometrics"},{"key":"2026102110221891500_kxag028-B7","doi-asserted-by":"publisher","first-page":"649","DOI":"10.1007\/s11222-019-09911-y","article-title":"Regularized estimation for highly multivariate log Gaussian Cox processes","volume":"30","author":"Choiruddin","year":"2020","journal-title":"Stat Comput"},{"key":"2026102110221891500_kxag028-B8","doi-asserted-by":"publisher","first-page":"404","DOI":"10.1111\/insr.12205","article-title":"A tutorial on Palm distributions for spatial point processes","volume":"85","author":"Coeurjolly","year":"2017","journal-title":"Int Stat Rev"},{"key":"2026102110221891500_kxag028-B9","doi-asserted-by":"publisher","first-page":"253","DOI":"10.1198\/1085711032291","article-title":"A Bayesian hierarchical model for risk assessment of methylmercury","volume":"8","author":"Coull","year":"2003","journal-title":"J Agric Biol Environ Stat"},{"key":"2026102110221891500_kxag028-B10","doi-asserted-by":"publisher","first-page":"200","DOI":"10.1111\/j.1462-2920.2005.00880.x","article-title":"Daime, a novel image analysis program for microbial ecology and biofilm research","volume":"8","author":"Daims","year":"2006","journal-title":"Environ Microbiol"},{"key":"2026102110221891500_kxag028-B11","doi-asserted-by":"publisher","first-page":"5002","DOI":"10.1128\/JB.00542-10","article-title":"The human oral microbiome","volume":"192","author":"Dewhirst","year":"2010","journal-title":"J Bacteriol"},{"key":"2026102110221891500_kxag028-B12","first-page":"423","article-title":"Point process methodology for on-line spatio-temporal disease surveillance","volume":"16","author":"Diggle","year":"2005","journal-title":"Environ Off J Int Environ Soc"},{"key":"2026102110221891500_kxag028-B13","doi-asserted-by":"publisher","first-page":"542","DOI":"10.1214\/13-STS441","article-title":"Spatial and spatio-temporal log-Gaussian cox processes: extending the geostatistical paradigm","author":"Diggle","year":"2013","journal-title":"Stat Sci"},{"key":"2026102110221891500_kxag028-B14","doi-asserted-by":"publisher","first-page":"299","DOI":"10.1111\/1467-9876.00113","article-title":"Model-based geostatistics","volume":"47","author":"Diggle","year":"1998","journal-title":"J R Stat Soc Ser C Appl Stat"},{"key":"2026102110221891500_kxag028-B15","doi-asserted-by":"publisher","first-page":"263","DOI":"10.1111\/1467-985X.00170","article-title":"Combining evidence on air pollution and daily mortality from the 20 largest us cities: a hierarchical modelling strategy","volume":"163","author":"Dominici","year":"2000","journal-title":"J R Stat Soc Ser A Stat Soc"},{"key":"2026102110221891500_kxag028-B16","doi-asserted-by":"publisher","first-page":"1128","DOI":"10.1080\/02664763.2021.2023116","article-title":"The integrated nested Laplace approximation applied to spatial log-Gaussian Cox process models","volume":"50","author":"Flagg","year":"2023","journal-title":"J Appl Stat"},{"key":"2026102110221891500_kxag028-B17","doi-asserted-by":"publisher","first-page":"216","DOI":"10.1109\/JPROC.2004.840301","article-title":"The design and implementation of FFTW3","volume":"93","author":"Frigo","year":"2005","journal-title":"Proc IEEE"},{"key":"2026102110221891500_kxag028-B18","doi-asserted-by":"publisher","first-page":"775","DOI":"10.1081\/STA-120018828","article-title":"Note on the sampling distribution for the Metropolis-Hastings algorithm","volume":"32","author":"Geweke","year":"2003","journal-title":"Commun Stat Theory Methods"},{"key":"2026102110221891500_kxag028-B19","doi-asserted-by":"publisher","first-page":"1871","DOI":"10.1137\/17M1149730","article-title":"Analysis of circulant embedding methods for sampling stationary random fields","volume":"56","author":"Graham","year":"2018","journal-title":"SIAM J Numer Anal"},{"key":"2026102110221891500_kxag028-B20","doi-asserted-by":"publisher","first-page":"12375","DOI":"10.1073\/pnas.1919099117","article-title":"Spatial mapping of polymicrobial communities reveals a precise biogeography associated with human dental caries","volume":"117","author":"Kim","year":"2020","journal-title":"Proc Natl Acad Sci USA"},{"key":"2026102110221891500_kxag028-B21","doi-asserted-by":"publisher","first-page":"0000102202155746791301","DOI":"10.2202\/1557-4679.1301","article-title":"Bayesian variable selection in semiparametric proportional hazards model for high dimensional survival data","volume":"7","author":"Lee","year":"2011","journal-title":"Int J Biostat"},{"key":"2026102110221891500_kxag028-B22","doi-asserted-by":"publisher","first-page":"114","DOI":"10.1002\/sam.11266","article-title":"Survival prediction and variable selection with simultaneous shrinkage and grouping priors","volume":"8","author":"Lee","year":"2015","journal-title":"Stat Anal Data Min ASA Data Sci J"},{"key":"2026102110221891500_kxag028-B23","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1016\/j.csda.2017.02.014","article-title":"Variable selection for high-dimensional genomic data with censored outcomes using group lasso prior","volume":"112","author":"Lee","year":"2017","journal-title":"Comput Stat Data Anal"},{"key":"2026102110221891500_kxag028-B24","doi-asserted-by":"publisher","first-page":"5446","DOI":"10.1002\/sim.10261","article-title":"Multivariate cluster point process to quantify and explore multi-entity configurations: Application to biofilm image data","volume":"43","author":"Majumder","year":"2024","journal-title":"Stat Med"},{"key":"2026102110221891500_kxag028-B25","doi-asserted-by":"publisher","first-page":"2104","DOI":"10.1128\/JB.02512-14","article-title":"The interaction between Streptococcus spp. and Veillonella tobetsuensis in the early stages of oral biofilm formation.","volume":"197","author":"Mashima","year":"2015","journal-title":"J Bacteriol"},{"key":"2026102110221891500_kxag028-B26","doi-asserted-by":"publisher","first-page":"e32219","DOI":"10.1371\/journal.pone.0032219","article-title":"Identifying low pH active and lactate-utilizing taxa within oral microbiome communities from healthy children using stable isotope probing techniques","volume":"7","author":"McLean","year":"2012","journal-title":"PLoS One"},{"key":"2026102110221891500_kxag028-B27","doi-asserted-by":"publisher","first-page":"10839","DOI":"10.1073\/pnas.1100292108","article-title":"Social evolution in multispecies biofilms","volume":"108","author":"Mitri","year":"2011","journal-title":"Proc Natl Acad Sci USA"},{"key":"2026102110221891500_kxag028-B28","doi-asserted-by":"crossref","DOI":"10.1201\/9780203496930","volume-title":"Statistical inference and simulation for spatial point processes","author":"Moller","year":"2003"},{"key":"2026102110221891500_kxag028-B29","doi-asserted-by":"publisher","first-page":"451","DOI":"10.1111\/1467-9469.00115","article-title":"Log Gaussian Cox processes","volume":"25","author":"M\u00f8ller","year":"1998","journal-title":"Scand J Stat"},{"key":"2026102110221891500_kxag028-B30","volume-title":"Bayesian learning for neural networks","author":"Neal","year":"2012"},{"key":"2026102110221891500_kxag028-B31","doi-asserted-by":"publisher","first-page":"681","DOI":"10.1198\/016214508000000337","article-title":"The Bayesian lasso","volume":"103","author":"Park","year":"2008","journal-title":"J Am Stat Assoc"},{"key":"2026102110221891500_kxag028-B32","doi-asserted-by":"publisher","first-page":"e1002012","DOI":"10.1371\/journal.ppat.1002012","article-title":"Metabolite cross-feeding enhances virulence in a model polymicrobial infection","volume":"7","author":"Ramsey","year":"2011","journal-title":"PLoS Pathogens"},{"key":"2026102110221891500_kxag028-B33","doi-asserted-by":"publisher","first-page":"1412","DOI":"10.1177\/09622802241262523","article-title":"Group lasso priors for Bayesian accelerated failure time models with left-truncated and interval-censored data","volume":"33","author":"Reeder","year":"2024","journal-title":"Stat Methods Med Res"},{"key":"2026102110221891500_kxag028-B34","doi-asserted-by":"crossref","DOI":"10.1201\/9780203492024","volume-title":"Gaussian Markov random fields: theory and applications","author":"Rue","year":"2005"},{"key":"2026102110221891500_kxag028-B35","doi-asserted-by":"publisher","first-page":"217","DOI":"10.1111\/rssc.12295","article-title":"Bayesian log-Gaussian Cox process regression: applications to meta-analysis of neuroimaging working memory studies","volume":"68","author":"Samartsidis","year":"2019","journal-title":"J R Stat Soc Ser C (Appl Stat)"},{"key":"2026102110221891500_kxag028-B36","doi-asserted-by":"publisher","first-page":"e37583","DOI":"10.1371\/journal.pone.0037583","article-title":"Co-localized or randomly distributed? pair cross correlation of in vivo grown subgingival biofilm bacteria quantified by digital image analysis","volume":"7","author":"Schillinger","year":"2012","journal-title":"PLoS One"},{"key":"2026102110221891500_kxag028-B37","doi-asserted-by":"publisher","first-page":"676","DOI":"10.1038\/nmeth.2019","article-title":"Fiji: an open-source platform for biological-image analysis","volume":"9","author":"Schindelin","year":"2012","journal-title":"Nat Methods"},{"key":"2026102110221891500_kxag028-B38","doi-asserted-by":"publisher","first-page":"e00536\u201323","DOI":"10.1128\/iai.00536-23","article-title":"Oral streptococci S. anginosus and S. mitis induce distinct morphological, inflammatory, and metabolic signatures in macrophages","volume":"92","author":"Senthil Kumar","year":"2024","journal-title":"Infect Immunity"},{"key":"2026102110221891500_kxag028-B39","doi-asserted-by":"crossref","first-page":"543","DOI":"10.1111\/rssc.12321","article-title":"Careful prior specification avoids incautious inference for log-Gaussian cox point processes","volume":"68","author":"Sorbye","year":"2019","journal-title":"J R Stat Soc Ser C (Appl Stat)"},{"key":"2026102110221891500_kxag028-B40","doi-asserted-by":"publisher","first-page":"583","DOI":"10.1111\/1467-9868.00353","article-title":"Bayesian measures of model complexity and fit","volume":"64","author":"Spiegelhalter","year":"2002","journal-title":"J R Stat Soc Ser B (Stat Methodol)"},{"key":"2026102110221891500_kxag028-B41","doi-asserted-by":"publisher","first-page":"1","DOI":"10.18637\/jss.v063.i07","article-title":"Bayesian inference and data augmentation schemes for spatial, spatiotemporal and multivariate log-Gaussian cox processes in R","volume":"63","author":"Taylor","year":"2015","journal-title":"J Stat Softw"},{"key":"2026102110221891500_kxag028-B42","doi-asserted-by":"publisher","first-page":"2227","DOI":"10.1080\/00949655.2017.1326117","article-title":"Bayesian computation for Log-Gaussian Cox processes: a comparative analysis of methods","volume":"87","author":"Teng","year":"2017","journal-title":"J Stat Comput Simulation"},{"key":"2026102110221891500_kxag028-B43","doi-asserted-by":"publisher","first-page":"4152","DOI":"10.1073\/pnas.1101134108","article-title":"Systems-level analysis of microbial community organization through combinatorial labeling and spectral imaging","volume":"108","author":"Valm","year":"2011","journal-title":"Proc Natl Acad Sci USA"},{"key":"2026102110221891500_kxag028-B44","doi-asserted-by":"publisher","first-page":"77","DOI":"10.1111\/rssc.12108","article-title":"Analysis of multispecies point patterns by using multivariate log-Gaussian Cox processes","volume":"65","author":"Waagepetersen","year":"2016","journal-title":"J R Stat Soc Ser C (Appl Stat)"},{"key":"2026102110221891500_kxag028-B45","doi-asserted-by":"publisher","first-page":"4003","DOI":"10.1016\/j.celrep.2020.02.097","article-title":"Spatial ecology of the human tongue dorsum microbiome","volume":"30","author":"Wilbert","year":"2020","journal-title":"Cell Rep"},{"key":"2026102110221891500_kxag028-B46","doi-asserted-by":"publisher","first-page":"409","DOI":"10.2307\/1390903","article-title":"Simulation of stationary Gaussian processes in [0, 1] d","volume":"3","author":"Wood","year":"1994","journal-title":"J Comput Graph Stat"},{"key":"2026102110221891500_kxag028-B47","first-page":"10","article-title":"Method of the year 2020: spatially resolved transcriptomics","volume":"18","author":"Xiaowei","year":"2021","journal-title":"Nat Methods"},{"key":"2026102110221891500_kxag028-B48","doi-asserted-by":"publisher","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":"2026102110221891500_kxag028-B49","doi-asserted-by":"publisher","first-page":"921","DOI":"10.1093\/biomet\/92.4.921","article-title":"Towards reconciling two asymptotic frameworks in spatial statistics","volume":"92","author":"Zhang","year":"2005","journal-title":"Biometrika"},{"key":"2026102110221891500_kxag028-B50","doi-asserted-by":"publisher","first-page":"774115","DOI":"10.3389\/froh.2021.774115","article-title":"Veillonellae: beyond bridging species in oral biofilm ecology","volume":"2","author":"Zhou","year":"2021","journal-title":"Front Oral Health"}],"container-title":["Biostatistics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/biostatistics\/article-pdf\/27\/1\/kxag028\/70687319\/kxag028.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/biostatistics\/article-pdf\/27\/1\/kxag028\/70687319\/kxag028.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,10,21]],"date-time":"2026-10-21T14:22:22Z","timestamp":1792592542000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/biostatistics\/article\/doi\/10.1093\/biostatistics\/kxag028\/8764052"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,1]]},"references-count":50,"journal-issue":{"issue":"1","published-print":{"date-parts":[[2026,1,20]]}},"URL":"https:\/\/doi.org\/10.1093\/biostatistics\/kxag028","relation":{},"ISSN":["1465-4644","1468-4357"],"issn-type":[{"value":"1465-4644","type":"print"},{"value":"1468-4357","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2026,1]]},"published":{"date-parts":[[2026,1]]},"article-number":"kxag028"}}