{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,3,26]],"date-time":"2026-03-26T12:11:10Z","timestamp":1774527070071,"version":"3.50.1"},"reference-count":49,"publisher":"Oxford University Press (OUP)","issue":"1","license":[{"start":{"date-parts":[[2025,11,19]],"date-time":"2025-11-19T00:00:00Z","timestamp":1763510400000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/pages\/standard-publication-reuse-rights"}],"funder":[{"DOI":"10.13039\/100000002","name":"U.S. National Institutes of Health","doi-asserted-by":"crossref","award":["R01HG003747"],"award-info":[{"award-number":["R01HG003747"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"crossref"}]},{"DOI":"10.13039\/100000002","name":"U.S. National Institutes of Health","doi-asserted-by":"crossref","award":["R21HG012881"],"award-info":[{"award-number":["R21HG012881"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"crossref"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2026,2,13]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:p>Motivated by the challenge of estimating effects of DNA methylation on 3D genomic contacts captured by multimodal single-cell Hi-C data, we consider the tensor-response partial least squares model with $ \\mathcal{Y}=\\mathcal{B}\\times_{1}X+\\mathcal{F} $, where the correlated high-dimensional predictors $ X\\in\\mathbb{R}^{n\\times d_{1}} $ and the sparse and noisy high-dimensional responses $ \\mathcal{Y}\\in\\mathbb{R}^{n\\times\\prod_{m}d_{m}} $ are observed, the low-rank and sparse partial least squares coefficient tensor $ \\mathcal{B}\\in\\mathbb{R}^{\\prod_{m}d_{m}} $ is unknown, and $ \\mathcal{F}\\in\\mathbb{R}^{n\\times\\prod_{m}d_{m}} $ is the noise tensor. In this work, we study the problem of estimating the partial least squares coefficient $ \\mathcal{B} $ and identifying its active entries in the tensor partial least squares framework. We show that the consistency of the existing tensor partial least squares estimator (Zhao et al., 2012) cannot be guaranteed under a high-dimensional regime in which both the number of predictors and the response tensor dimensions grow faster than the sample size. To address this high dimensionality, we propose the sparse higher-order partial least squares estimator and an accompanying algorithm that can simultaneously perform variable selection, dimension reduction and tensor-response denoising. We establish asymptotic guarantees for the proposed estimator in the high-dimensional regime and validate these results through comprehensive simulation studies, demonstrating our method\u2019s advantages over baseline approaches. Finally, application of the proposed estimator to the motivating multimodal single-cell Hi-C data provides novel biological insights into gene regulation by multiple regulatory elements.<\/jats:p>","DOI":"10.1093\/biomet\/asaf082","type":"journal-article","created":{"date-parts":[[2025,11,7]],"date-time":"2025-11-07T13:02:13Z","timestamp":1762520533000},"source":"Crossref","is-referenced-by-count":0,"title":["Sparse higher-order partial least squares for simultaneous variable selection, dimension reduction and tensor denoising"],"prefix":"10.1093","volume":"113","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-0987-923X","authenticated-orcid":false,"given":"Kwangmoon","family":"Park","sequence":"first","affiliation":[{"name":"Department of Statistics, University of Wisconsin-Madison , 1205 University Avenue , Madison, Wisconsin 53706, \u00a0 kpark243@wisc.edu","place":["U.S.A"]}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-9048-0922","authenticated-orcid":false,"given":"S\u00fcnd\u00fcz","family":"Kele\u015f","sequence":"additional","affiliation":[{"name":"Department of Biostatistics and Medical Informatics, University of Wisconsin , 425 Henry Mall , Madison, Wisconsin 53706,","place":["U.S.A"]}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2025,11,19]]},"reference":[{"key":"2026032607084828000_asaf082-B1","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1186\/s13059-020-02015-1","article-title":"MOFA+: A statistical framework for comprehensive integration of multi-modal single-cell data","volume":"21","author":"Argelaguet","year":"2020","journal-title":"Genome Biol"},{"key":"2026032607084828000_asaf082-B2","first-page":"1044","article-title":"Perturbation bounds for (nearly) orthogonally decomposable tensors with statistical applications","volume":"12","author":"Auddy","year":"2023","journal-title":"Info. Infer"},{"key":"2026032607084828000_asaf082-B3","doi-asserted-by":"crossref","first-page":"1382","DOI":"10.1016\/j.jmva.2005.08.003","article-title":"Eigenvalues of large sample covariance matrices of spiked population models","volume":"97","author":"Baik","year":"2006","journal-title":"J. Mult. Anal"},{"key":"2026032607084828000_asaf082-B4","doi-asserted-by":"crossref","first-page":"3243","DOI":"10.1038\/s41467-021-23142-8","article-title":"Impact of DNA methylation on 3D genome structure","volume":"12","author":"Buitrago","year":"2021","journal-title":"Nature Commun"},{"key":"2026032607084828000_asaf082-B5","first-page":"3074","article-title":"Sparse PCA: Optimal rates and adaptive estimation","volume":"41","author":"Cai","year":"2012","journal-title":"Ann. Statist."},{"key":"2026032607084828000_asaf082-B6","doi-asserted-by":"crossref","first-page":"3","DOI":"10.1111\/j.1467-9868.2009.00723.x","article-title":"Sparse partial least squares regression for simultaneous dimension reduction and variable selection","volume":"72","author":"Chun","year":"2010","journal-title":"J. R. Statist. Soc. B"},{"key":"2026032607084828000_asaf082-B7","doi-asserted-by":"crossref","first-page":"17","DOI":"10.2202\/1544-6115.1492","article-title":"Sparse partial least squares classification for high dimensional data","volume":"9","author":"Chung","year":"2010","journal-title":"Statist. Appl. Genet. Molec. Biol"},{"key":"2026032607084828000_asaf082-B8","doi-asserted-by":"crossref","first-page":"851","DOI":"10.1111\/rssb.12018","article-title":"Envelopes and partial least squares regression","volume":"75","author":"Cook","year":"2013","journal-title":"J. R. Statist. Soc. B"},{"key":"2026032607084828000_asaf082-B9","doi-asserted-by":"crossref","first-page":"251","DOI":"10.1016\/0169-7439(93)85002-X","article-title":"SIMPLS: An alternative approach to partial least squares regression","volume":"18","author":"De Jong","year":"1993","journal-title":"Chemometr. Intel. Lab. Syst"},{"key":"2026032607084828000_asaf082-B10","doi-asserted-by":"crossref","first-page":"1253","DOI":"10.1137\/S0895479896305696","article-title":"A multilinear singular value decomposition","volume":"21","author":"De Lathauwer","year":"2000","journal-title":"SIAM J. Matrix Anal. Applic"},{"key":"2026032607084828000_asaf082-B11","doi-asserted-by":"crossref","first-page":"1084","DOI":"10.1137\/06066518X","article-title":"Tensor rank and the ill-posedness of the best low-rank approximation problem","volume":"30","author":"De Silva","year":"2008","journal-title":"SIAM J. Matrix Anal. Applic"},{"key":"2026032607084828000_asaf082-B12","first-page":"4913","article-title":"Sparse PCA via covariance thresholding","volume":"17","author":"Deshpande","year":"2014","journal-title":"J. Mach. Learn. Res"},{"key":"2026032607084828000_asaf082-B13","doi-asserted-by":"crossref","first-page":"977","DOI":"10.1137\/18M1211350","article-title":"Sparse principal component analysis via variable projection","volume":"80","author":"Erichson","year":"2020","journal-title":"SIAM J. Appl. Math."},{"key":"2026032607084828000_asaf082-B14","author":"Friedman","year":"2004"},{"key":"2026032607084828000_asaf082-B15","doi-asserted-by":"crossref","first-page":"2074","DOI":"10.1214\/16-AOS1519","article-title":"Sparse CCA: Adaptive estimation and computational barriers","volume":"45","author":"Gao","year":"2017","journal-title":"Ann. Statist."},{"key":"2026032607084828000_asaf082-B16","doi-asserted-by":"crossref","first-page":"3573","DOI":"10.1016\/j.cell.2021.04.048","article-title":"Integrated analysis of multimodal single-cell data","volume":"184","author":"Hao","year":"2021","journal-title":"Cell"},{"key":"2026032607084828000_asaf082-B17","doi-asserted-by":"crossref","DOI":"10.1201\/b18401","volume-title":"Statistical Learning with Sparsity: The Lasso and Generalizations","author":"Hastie","year":"2015"},{"key":"2026032607084828000_asaf082-B18","first-page":"97","article-title":"Partial least squares regression and statistical models","volume":"17","author":"Helland","year":"1990","journal-title":"Scand. J. Statist"},{"key":"2026032607084828000_asaf082-B19","doi-asserted-by":"crossref","first-page":"583","DOI":"10.1080\/01621459.1994.10476783","article-title":"Comparison of prediction methods when only a few components are relevant","volume":"89","author":"Helland","year":"1994","journal-title":"J. Am. Statist. Assoc"},{"key":"2026032607084828000_asaf082-B20","doi-asserted-by":"crossref","first-page":"295","DOI":"10.1214\/aos\/1009210544","article-title":"On the distribution of the largest eigenvalue in principal components analysis","volume":"29","author":"Johnstone","year":"2001","journal-title":"Ann. Statist."},{"key":"2026032607084828000_asaf082-B21","doi-asserted-by":"crossref","first-page":"682","DOI":"10.1198\/jasa.2009.0121","article-title":"On consistency and sparsity for principal components analysis in high dimensions","volume":"104","author":"Johnstone","year":"2009","journal-title":"J. Am. Statist. Assoc"},{"key":"2026032607084828000_asaf082-B22","doi-asserted-by":"crossref","first-page":"455","DOI":"10.1137\/07070111X","article-title":"Tensor decompositions and applications","volume":"51","author":"Kolda","year":"2009","journal-title":"SIAM Rev."},{"key":"2026032607084828000_asaf082-B23","doi-asserted-by":"crossref","first-page":"1455","DOI":"10.1214\/21-AOAS1475","article-title":"Tensor quantile regression with application to association between neuroimages and human intelligence","volume":"15","author":"Li","year":"2021","journal-title":"Ann. Appl. Statist"},{"key":"2026032607084828000_asaf082-B24","first-page":"1","article-title":"Feature selection: A data perspective","volume":"50","author":"Li","year":"2017","journal-title":"ACM Comp. Surv"},{"key":"2026032607084828000_asaf082-B25","doi-asserted-by":"crossref","first-page":"1131","DOI":"10.1080\/01621459.2016.1193022","article-title":"Parsimonious tensor response regression","volume":"112","author":"Li","year":"2017","journal-title":"J. Am. Statist. Assoc"},{"key":"2026032607084828000_asaf082-B26","doi-asserted-by":"crossref","first-page":"129","DOI":"10.1038\/s41586-021-03604-1","article-title":"An atlas of gene regulatory elements in adult mouse cerebrum","volume":"598","author":"Li","year":"2021","journal-title":"Nature"},{"key":"2026032607084828000_asaf082-B27","doi-asserted-by":"crossref","first-page":"e2303647120","DOI":"10.1073\/pnas.2303647120","article-title":"Quantifying common and distinct information in single-cell multimodal data with tilted canonical correlation analysis","volume":"120","author":"Lin","year":"2023","journal-title":"Proc. Nat. Acad. Sci."},{"key":"2026032607084828000_asaf082-B28","doi-asserted-by":"crossref","first-page":"120","DOI":"10.1038\/s41586-020-03182-8","article-title":"DNA methylation atlas of the mouse brain at single-cell resolution","volume":"598","author":"Liu","year":"2021","journal-title":"Nature"},{"key":"2026032607084828000_asaf082-B29","doi-asserted-by":"crossref","first-page":"1436","DOI":"10.1214\/009053606000000281","article-title":"High-dimensional graphs and variable selection with the lasso","volume":"34","author":"Meinshausen","year":"2006","journal-title":"Ann. Statist."},{"key":"2026032607084828000_asaf082-B30","doi-asserted-by":"crossref","first-page":"404","DOI":"10.1038\/s41594-024-01241-6","article-title":"The impact of DNA methylation on CTCF-mediated 3D genome organization","volume":"31","author":"Monteagudo-S\u00e1nchez","year":"2024","journal-title":"Nature Struct. Molec. Biol"},{"key":"2026032607084828000_asaf082-B31","doi-asserted-by":"crossref","first-page":"763","DOI":"10.1111\/1467-9868.00262","article-title":"Partial least squares estimator for single-index models","volume":"62","author":"Naik","year":"2000","journal-title":"J. R. Statist. Soc. B"},{"key":"2026032607084828000_asaf082-B32","doi-asserted-by":"crossref","first-page":"2464","DOI":"10.1080\/01621459.2024.2358557","article-title":"Joint tensor modeling of single cell 3D genome and epigenetic data with Muscle","volume":"119","author":"Park","year":"2024","journal-title":"J. Am. Statist. Assoc"},{"key":"2026032607084828000_asaf082-B33","doi-asserted-by":"crossref","first-page":"1194","DOI":"10.1016\/j.cell.2009.06.001","article-title":"CTCF: Master weaver of the genome","volume":"137","author":"Phillips","year":"2009","journal-title":"Cell"},{"key":"2026032607084828000_asaf082-B34","doi-asserted-by":"crossref","first-page":"858","DOI":"10.1016\/j.molcel.2018.06.044","article-title":"Cicero predicts cis-regulatory DNA interactions from single-cell chromatin accessibility data","volume":"71","author":"Pliner","year":"2018","journal-title":"Molec. Cell"},{"key":"2026032607084828000_asaf082-B35","doi-asserted-by":"crossref","first-page":"6976","DOI":"10.1109\/TIT.2011.2165799","article-title":"Minimax rates of estimation for high-dimensional linear regression over $ l_{q} $-balls","volume":"57","author":"Raskutti","year":"2011","journal-title":"IEEE Trans. Info. Theory"},{"key":"2026032607084828000_asaf082-B36","doi-asserted-by":"crossref","first-page":"1554","DOI":"10.1214\/18-AOS1725","article-title":"Convex regularization for high-dimensional multiresponse tensor regression","volume":"47","author":"Raskutti","year":"2019","journal-title":"Ann. Statist."},{"key":"2026032607084828000_asaf082-B37","doi-asserted-by":"crossref","first-page":"437","DOI":"10.1038\/s41576-019-0128-0","article-title":"Long-range enhancer\u2013promoter contacts in gene expression control","volume":"20","author":"Schoenfelder","year":"2019","journal-title":"Nature Rev. Genet"},{"key":"2026032607084828000_asaf082-B38","first-page":"4908","article-title":"STORE: Sparse tensor response regression and neuroimaging analysis","volume":"18","author":"Sun","year":"2017","journal-title":"J. Mach. Learn. Res"},{"key":"2026032607084828000_asaf082-B39","doi-asserted-by":"crossref","first-page":"899","DOI":"10.1111\/rssb.12190","article-title":"Provable sparse tensor decomposition","volume":"79","author":"Sun","year":"2017","journal-title":"J. R. Statist. Soc. B"},{"key":"2026032607084828000_asaf082-B40","doi-asserted-by":"crossref","first-page":"267","DOI":"10.1111\/j.2517-6161.1996.tb02080.x","article-title":"Regression shrinkage and selection via the lasso","volume":"58","author":"Tibshirani","year":"1996","journal-title":"J. R. Statist. Soc. B"},{"key":"2026032607084828000_asaf082-B41","doi-asserted-by":"crossref","first-page":"279","DOI":"10.1007\/BF02289464","article-title":"Some mathematical notes on three-mode factor analysis","volume":"31","author":"Tucker","year":"1966","journal-title":"Psychometrika"},{"key":"2026032607084828000_asaf082-B42","volume-title":"Multivariate Analysis","author":"Wold","year":"1966"},{"key":"2026032607084828000_asaf082-B43","doi-asserted-by":"crossref","first-page":"49","DOI":"10.1111\/j.1467-9868.2005.00532.x","article-title":"Model selection and estimation in regression with grouped variables","volume":"68","author":"Yuan","year":"2006","journal-title":"J. R. Statist. Soc. B"},{"key":"2026032607084828000_asaf082-B44","doi-asserted-by":"crossref","first-page":"1708","DOI":"10.1080\/01621459.2018.1527227","article-title":"Optimal sparse singular value decomposition for high-dimensional high-order data","volume":"114","author":"Zhang","year":"2019","journal-title":"J. Am. Statist. Assoc"},{"key":"2026032607084828000_asaf082-B45","doi-asserted-by":"crossref","first-page":"7311","DOI":"10.1109\/TIT.2018.2841377","article-title":"Tensor SVD: Statistical and computational limits","volume":"64","author":"Zhang","year":"2018","journal-title":"IEEE Trans. Info. Theory"},{"key":"2026032607084828000_asaf082-B46","doi-asserted-by":"crossref","first-page":"426","DOI":"10.1080\/00401706.2016.1272495","article-title":"Tensor envelope partial least-squares regression","volume":"59","author":"Zhang","year":"2017","journal-title":"Technometrics"},{"key":"2026032607084828000_asaf082-B47","doi-asserted-by":"crossref","first-page":"1660","DOI":"10.1109\/TPAMI.2012.254","article-title":"Higher order partial least squares (HOPLS): A generalized multilinear regression method","volume":"35","author":"Zhao","year":"2012","journal-title":"IEEE Trans. Pat. Anal. Mach. Intel"},{"key":"2026032607084828000_asaf082-B48","doi-asserted-by":"crossref","first-page":"161","DOI":"10.1214\/18-AOS1796","article-title":"Envelope-based sparse partial least squares","volume":"48","author":"Zhu","year":"2020","journal-title":"Ann. Statist."},{"key":"2026032607084828000_asaf082-B49","doi-asserted-by":"crossref","first-page":"301","DOI":"10.1111\/j.1467-9868.2005.00503.x","article-title":"Regularization and variable selection via the elastic net","volume":"67","author":"Zou","year":"2005","journal-title":"J. R. Statist. Soc. B"}],"container-title":["Biometrika"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/biomet\/advance-article-pdf\/doi\/10.1093\/biomet\/asaf082\/65394805\/asaf082.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/biomet\/article-pdf\/113\/1\/asaf082\/65394805\/asaf082.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/biomet\/article-pdf\/113\/1\/asaf082\/65394805\/asaf082.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,3,26]],"date-time":"2026-03-26T11:09:00Z","timestamp":1774523340000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/biomet\/article\/doi\/10.1093\/biomet\/asaf082\/8328024"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2025,11,19]]},"references-count":49,"journal-issue":{"issue":"1","published-print":{"date-parts":[[2026,2,13]]}},"URL":"https:\/\/doi.org\/10.1093\/biomet\/asaf082","relation":{},"ISSN":["0006-3444","1464-3510"],"issn-type":[{"value":"0006-3444","type":"print"},{"value":"1464-3510","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2026]]},"published":{"date-parts":[[2025,11,19]]},"article-number":"asaf082"}}