{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,15]],"date-time":"2026-09-15T19:35:17Z","timestamp":1789500917725,"version":"build-2803163510"},"reference-count":33,"publisher":"Oxford University Press (OUP)","issue":"9","license":[{"start":{"date-parts":[[2026,5,11]],"date-time":"2026-05-11T00:00:00Z","timestamp":1778457600000},"content-version":"vor","delay-in-days":6,"URL":"https:\/\/creativecommons.org\/licenses\/by-nc\/4.0\/"}],"funder":[{"DOI":"10.13039\/100000002","name":"National Institutes of Health","doi-asserted-by":"publisher","award":["R35GM146586"],"award-info":[{"award-number":["R35GM146586"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000002","name":"National Institutes of Health","doi-asserted-by":"publisher","award":["P20GM130454"],"award-info":[{"award-number":["P20GM130454"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000002","name":"National Institutes of Health","doi-asserted-by":"publisher","award":["P30CA023108"],"award-info":[{"award-number":["P30CA023108"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2026,5,5]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:p>High-throughput spatial transcriptomics (ST) now profiles hundreds of thousands of cells or locations per section, creating computational bottlenecks for routine analysis. Sketching, or intelligent sub-sampling, addresses scale by selecting small, representative subsets. While effective for single-cell RNA sequencing data, existing sketching methods, which optimize coverage in expression space but ignore physical location, can introduce spatial bias when applied to ST data. To explore the impact of sketching on ST analysis, we systematically benchmarked uniform sampling, leverage-score sampling, Geosketch (minimax\/Hausdorff), and scSampler (maximin) across multiple real ST datasets (mouse ovary, MERFISH brain, human breast cancer, lung) and simulations, using three input representations: Principle Component Analysis (PCA) embeddings, spatial coordinates, and spatially smoothed embeddings. We show that expression-only designs capture global transcriptomic heterogeneity but distort tissue architecture by over-sampling high-variability regions and under-sampling homogeneous areas. Coordinate-only sampling restores tissue coverage but misses transcriptional extremes. A simple spatially aware extension, computing leverage scores from a randomized singular value decomposition (SVD) basis smoothed by a spatial weights matrix, strikes a favorable balance, recovering rare cell states while maintaining uniform tissue coverage and avoiding edge effects. Across robust Hausdorff distances, clustering stability (Adjusted Rand Index), PCA loading drift, and local cell-type mean squared error (MSE), spatially smoothed leverage scores match or outperform alternatives. These results motivate joint spatial\u2013transcriptomic sketching objectives to enable fast, unbiased analyses of increasingly large ST datasets.<\/jats:p>","DOI":"10.1093\/nar\/gkag434","type":"journal-article","created":{"date-parts":[[2026,4,24]],"date-time":"2026-04-24T11:42:32Z","timestamp":1777030952000},"source":"Crossref","is-referenced-by-count":1,"title":["Benchmarking sketching methods on spatial transcriptomics data"],"prefix":"10.1093","volume":"54","author":[{"ORCID":"https:\/\/orcid.org\/0000-0003-0881-5258","authenticated-orcid":false,"given":"Ian K","family":"Gingerich","sequence":"first","affiliation":[{"name":"Department of Biomedical Data Science, Geisel School of Medicine, Dartmouth College , Hanover, NH 03755 ,","place":["United States"]},{"name":"Thayer School of Engineering, Dartmouth College , Hanover, NH 03755 ,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Brittany A","family":"Goods","sequence":"additional","affiliation":[{"name":"Thayer School of Engineering, Dartmouth College , Hanover, NH 03755 ,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Hildreth Robert","family":"Frost","sequence":"additional","affiliation":[{"name":"Department of Biomedical Data Science, Geisel School of Medicine, Dartmouth College , Hanover, NH 03755 ,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2026,5,11]]},"reference":[{"key":"2026051105204999900_B1","doi-asserted-by":"publisher","first-page":"28","DOI":"10.1080\/08982112.2015.1100447","article-title":"Space-filling designs for computer experiments: a review","volume":"28","author":"Joseph","year":"2016","journal-title":"Qual Eng"},{"key":"2026051105204999900_B2","doi-asserted-by":"publisher","DOI":"10.48550\/arXiv.1203.3472","article-title":"Super-samples from Kernel Herding. Proceedings of the Twent-Sixth Conference on Uncertainty in Artificial Intelligence","author":"Chen","year":"2012"},{"key":"2026051105204999900_B3","doi-asserted-by":"publisher","first-page":"68","DOI":"10.1016\/j.spasta.2017.08.004","article-title":"Spatial subsemble estimator for large geostatistical data","volume":"22","author":"Barbian","year":"2017","journal-title":"Spat Stat"},{"key":"2026051105204999900_B4","doi-asserted-by":"crossref","first-page":"131","DOI":"10.1016\/0378-3758(90)90122-B","article-title":"Minimax and maximin distance designs","volume":"26","author":"Johnson","year":"1990","journal-title":"J Stat Plan Inference"},{"key":"2026051105204999900_B5","doi-asserted-by":"publisher","first-page":"54:1","DOI":"10.1145\/3019134","article-title":"Low-rank approximation and regression in input sparsity time","volume":"63","author":"Clarkson","year":"2017","journal-title":"J ACM"},{"key":"2026051105204999900_B6","doi-asserted-by":"publisher","first-page":"621","DOI":"10.1038\/s41592-020-0825-9","article-title":"Putative cell type discovery from single-cell gene expression data","volume":"17","author":"Miao","year":"2020","journal-title":"Nat Methods"},{"key":"2026051105204999900_B7","doi-asserted-by":"publisher","first-page":"e8746","DOI":"10.15252\/msb.20188746","article-title":"Current best practices in single-cell RNA-seq analysis: a tutorial","volume":"15","author":"Luecken","year":"2019","journal-title":"Mol Syst Biol"},{"key":"2026051105204999900_B8","doi-asserted-by":"publisher","first-page":"40","DOI":"10.1186\/s12859-019-2599-6","article-title":"Comparative analysis of differential gene expression analysis tools for single-cell RNA sequencing data","volume":"20","author":"Wang","year":"2019","journal-title":"BMC Bioinf"},{"key":"2026051105204999900_B9","doi-asserted-by":"publisher","first-page":"3126","DOI":"10.1093\/bioinformatics\/btac271","article-title":"scSampler: fast diversity-preserving subsampling of large-scale single-cell transcriptomic data","volume":"38","author":"Song","year":"2022","journal-title":"Bioinformatics"},{"key":"2026051105204999900_B10","doi-asserted-by":"publisher","first-page":"483","DOI":"10.1016\/j.cels.2019.05.003","article-title":"Geometric sketching compactly summarizes the single-cell transcriptomic landscape","volume":"8","author":"Hie","year":"2019","journal-title":"Cell Syst"},{"key":"2026051105204999900_B11","doi-asserted-by":"publisher","first-page":"293","DOI":"10.1038\/s41587-023-01767-y","article-title":"Dictionary learning for integrative, multimodal and scalable single-cell analysis","volume":"42","author":"Hao","year":"2024","journal-title":"Nat Biotechnol"},{"key":"2026051105204999900_B12","doi-asserted-by":"publisher","first-page":"i236","DOI":"10.1093\/bioinformatics\/btaa408","article-title":"Hopper: a mathematically optimal algorithm for sketching biological data","volume":"36","author":"DeMeo","year":"2020","journal-title":"Bioinformatics"},{"key":"2026051105204999900_B13","doi-asserted-by":"publisher","first-page":"3","DOI":"10.1093\/bib\/bbaf279","article-title":"scValue: value-based subsampling of large-scale single-cell transcriptomic data for machine and deep learning tasks","volume-title":"Briefings in Bioinformatics","author":"Huang","year":"2025"},{"key":"2026051105204999900_B14","doi-asserted-by":"publisher","first-page":"101126","DOI":"10.1016\/j.isci.2020.101126","article-title":"Spherical thresholding improves sketching of single-cell transcriptomic heterogeneity","volume":"23","author":"Do","year":"2020","journal-title":"iScience"},{"key":"2026051105204999900_B15","doi-asserted-by":"publisher","first-page":"104093","DOI":"10.1016\/j.jbi.2022.104093","article-title":"A downsampling method enables robust clustering and integration of single-cell transcriptome data","volume":"130","author":"Ren","year":"2022","journal-title":"J Biomed Inform"},{"key":"2026051105204999900_B16","doi-asserted-by":"publisher","first-page":"btae546","DOI":"10.1093\/bioinformatics\/btae546","article-title":"Bin2cell reconstructs cells from high resolution Visium HD data","volume":"40","author":"Polan\u00b4ski","year":"2024","journal-title":"Bioinformatics"},{"key":"2026051105204999900_B17","doi-asserted-by":"publisher","first-page":"313","DOI":"10.1038\/s41587-020-0739-1","article-title":"Highly sensitive spatial transcriptomics at near-cellular resolution with Slide-seqV2","volume":"39","author":"Stickels","year":"2021","journal-title":"Nat Biotechnol"},{"key":"2026051105204999900_B18","doi-asserted-by":"publisher","DOI":"10.1101\/2022.10.06.510405v1","article-title":"High resolution mapping of the breast cancer tumor microenvironment using integrated single cell, spatial and in situ analysis of FFPE tissue","author":"Janesick","year":"2022"},{"key":"2026051105204999900_B19","doi-asserted-by":"publisher","DOI":"10.1101\/2023.03.06.531121v1","article-title":"A high-resolution transcriptomic and spatial atlas of cell types in the whole mouse brain","author":"Yao","year":"2023"},{"key":"2026051105204999900_B20","doi-asserted-by":"publisher","first-page":"1777","DOI":"10.1016\/j.cell.2022.04.003","article-title":"Spatiotemporal transcriptomic atlas of mouse organogenesis using DNA nanoball-patterned arrays. en","volume":"185","author":"Chen","year":"2022","journal-title":"Cell"},{"key":"2026051105204999900_B21","doi-asserted-by":"crossref","first-page":"D1053","DOI":"10.1093\/nar\/gkad933","article-title":"STOmicsDB: a comprehensive database for spatial transcriptomics data sharing, analysis and visualization","volume":"52","author":"Xu","year":"2024","journal-title":"Nucleic Acids Res"},{"key":"2026051105204999900_B22","doi-asserted-by":"publisher","first-page":"5892","DOI":"10.1016\/j.cell.2023.11.016","article-title":"Spatiotemporal transcriptome atlas reveals the regional specification of the developing human brain","volume":"186","author":"Li","year":"2023","journal-title":"Cell"},{"key":"2026051105204999900_B23","doi-asserted-by":"publisher","DOI":"10.1101\/2023.04.22.537909v1","article-title":"Spatiotemporal transcriptome atlas of human embryos after gastrulation","author":"Pan","year":"2023"},{"key":"2026051105204999900_B24","doi-asserted-by":"publisher","first-page":"552","DOI":"10.1038\/s41586-023-06569-5","article-title":"Spatial atlas of the mouse central nervous system at molecular resolution","volume":"622","author":"Shi","year":"2023","journal-title":"Nature"},{"key":"2026051105204999900_B25","doi-asserted-by":"crossref","first-page":"2141","DOI":"10.1016\/j.neuron.2025.02.015","article-title":"Single-cell spatial transcriptomic atlas of the whole mouse brain","volume":"113","author":"Han","year":"2025","journal-title":"Neuron"},{"key":"2026051105204999900_B26","doi-asserted-by":"publisher","first-page":"343","DOI":"10.1038\/s41586-023-06808-9","article-title":"Molecularly defined and spatially resolved cell atlas of the whole mouse brain","volume":"624","author":"Zhang","year":"2023","journal-title":"Nature"},{"key":"2026051105204999900_B27","doi-asserted-by":"publisher","first-page":"163","DOI":"10.1038\/s41587-021-01206-w","article-title":"A Python library for probabilistic analysis of single-cell omics data","volume":"40","author":"Gayoso","year":"2022","journal-title":"Nat Biotechnol"},{"key":"2026051105204999900_B28","doi-asserted-by":"publisher","first-page":"2525","DOI":"10.1038\/s41592-025-02814-z","article-title":"Nicheformer: a foundation model for single-cell and spatial omics","volume":"22","author":"Tejada-Lapuerta","year":"2025","journal-title":"Nat Methods"},{"key":"2026051105204999900_B29","doi-asserted-by":"publisher","DOI":"10.1101\/2025.02.05.636714v1","article-title":"scGPT-spatial: continual pretraining of single-cell foundation model for spatial transcriptomics","author":"Wang","year":"2025"},{"key":"2026051105204999900_B30","doi-asserted-by":"publisher","first-page":"e3003193","DOI":"10.1371\/journal.pbio.3003193","article-title":"Single-cell and spatiotemporal profile of ovulation in the mouse ovary","volume":"23","author":"Huang","year":"2025","journal-title":"PLoS Biol"},{"key":"2026051105204999900_B31","doi-asserted-by":"publisher","first-page":"99","DOI":"10.1186\/s13059-024-03241-7","article-title":"Library size confounds biology in spatial transcriptomics data","volume":"25","author":"Bhuva","year":"2024","journal-title":"Genome Biol"},{"key":"2026051105204999900_B32","doi-asserted-by":"publisher","first-page":"39","DOI":"10.1186\/s13059-023-02879-z","article-title":"SRTsim: spatial pattern preserving simulations for spatially resolved transcrip-tomics","volume":"24","author":"Zhu","year":"2023","journal-title":"Genome Biol"},{"key":"2026051105204999900_B33","doi-asserted-by":"publisher","first-page":"e1013759","DOI":"10.1371\/journal.pcbi.1013759","article-title":"Randomized Spatial PCA (RASP): a computationally efficient method for dimensionality reduction of high-resolution spatial transcriptomics data","volume":"21","author":"Gingerich","year":"2025","journal-title":"PLoS Comput Biol"}],"container-title":["Nucleic Acids Research"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/nar\/article-pdf\/54\/9\/gkag434\/68264987\/gkag434.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/nar\/article-pdf\/54\/9\/gkag434\/68264987\/gkag434.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,5,11]],"date-time":"2026-05-11T09:21:00Z","timestamp":1778491260000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/nar\/article\/doi\/10.1093\/nar\/gkag434\/8675578"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,5,5]]},"references-count":33,"journal-issue":{"issue":"9","published-print":{"date-parts":[[2026,5,5]]}},"URL":"https:\/\/doi.org\/10.1093\/nar\/gkag434","relation":{},"ISSN":["0305-1048","1362-4962"],"issn-type":[{"value":"0305-1048","type":"print"},{"value":"1362-4962","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2026,5,22]]},"published":{"date-parts":[[2026,5,5]]},"article-number":"gkag434"}}