{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,5,16]],"date-time":"2026-05-16T13:33:40Z","timestamp":1778938420870,"version":"3.51.4"},"update-to":[{"DOI":"10.1371\/journal.pcbi.1008741","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2021,4,8]],"date-time":"2021-04-08T00:00:00Z","timestamp":1617840000000}}],"reference-count":43,"publisher":"Public Library of Science (PLoS)","issue":"3","license":[{"start":{"date-parts":[[2021,3,29]],"date-time":"2021-03-29T00:00:00Z","timestamp":1616976000000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["www.ploscompbiol.org"],"crossmark-restriction":false},"short-container-title":["PLoS Comput Biol"],"abstract":"<jats:p>\n                    Imaging Mass Cytometry (IMC) combines laser ablation and mass spectrometry to quantitate metal-conjugated primary antibodies incubated in intact tumor tissue slides. This strategy allows spatially-resolved multiplexing of dozens of simultaneous protein targets with 1\n                    <jats:italic>\u03bcm<\/jats:italic>\n                    resolution. Each slide is a spatial assay consisting of high-dimensional multivariate observations (\n                    <jats:italic>m<\/jats:italic>\n                    -dimensional feature space) collected at different spatial positions and capturing data from a single biological sample or even representative spots from multiple samples when using tissue microarrays. Often, each of these spatial assays could be characterized by several regions of interest (ROIs). To extract meaningful information from the multi-dimensional observations recorded at different ROIs across different assays, we propose to analyze such datasets using a two-step graph-based approach. We first construct for each ROI a graph representing the interactions between the\n                    <jats:italic>m<\/jats:italic>\n                    covariates and compute an\n                    <jats:italic>m<\/jats:italic>\n                    dimensional vector characterizing the steady state distribution among features. We then use all these\n                    <jats:italic>m<\/jats:italic>\n                    -dimensional vectors to construct a graph between the ROIs from all assays. This second graph is subjected to a nonlinear dimension reduction analysis, retrieving the intrinsic geometric representation of the ROIs. Such a representation provides the foundation for efficient and accurate organization of the different ROIs that correlates with their phenotypes. Theoretically, we show that when the ROIs have a particular bi-modal distribution, the new representation gives rise to a better distinction between the two modalities compared to the maximum a posteriori (MAP) estimator. We applied our method to predict the sensitivity to PD-1 axis blockers treatment of lung cancer subjects based on IMC data, achieving 97.3% average accuracy on two IMC datasets. This serves as empirical evidence that the graph of graphs approach enables us to integrate multiple ROIs and the intra-relationships between the features at each ROI, giving rise to an informative representation that is strongly associated with the phenotypic state of the entire image.\n                  <\/jats:p>","DOI":"10.1371\/journal.pcbi.1008741","type":"journal-article","created":{"date-parts":[[2021,3,29]],"date-time":"2021-03-29T13:41:48Z","timestamp":1617025308000},"page":"e1008741","update-policy":"https:\/\/doi.org\/10.1371\/journal.pcbi.corrections_policy","source":"Crossref","is-referenced-by-count":10,"title":["Graph of graphs analysis for multiplexed data with application to imaging mass cytometry"],"prefix":"10.1371","volume":"17","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-4592-0244","authenticated-orcid":true,"given":"Ya-Wei Eileen","family":"Lin","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-9058-1870","authenticated-orcid":true,"given":"Tal","family":"Shnitzer","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-6838-1423","authenticated-orcid":true,"given":"Ronen","family":"Talmon","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Franz","family":"Villarroel-Espindola","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-5234-5890","authenticated-orcid":true,"given":"Shruti","family":"Desai","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Kurt","family":"Schalper","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Yuval","family":"Kluger","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"340","published-online":{"date-parts":[[2021,3,29]]},"reference":[{"issue":"5500","key":"pcbi.1008741.ref001","doi-asserted-by":"crossref","first-page":"2319","DOI":"10.1126\/science.290.5500.2319","article-title":"A global geometric framework for nonlinear dimensionality reduction","volume":"290","author":"JB Tenenbaum","year":"2000","journal-title":"science"},{"issue":"5500","key":"pcbi.1008741.ref002","doi-asserted-by":"crossref","first-page":"2323","DOI":"10.1126\/science.290.5500.2323","article-title":"Nonlinear dimensionality reduction by locally linear embedding","volume":"290","author":"ST Roweis","year":"2000","journal-title":"science"},{"issue":"6","key":"pcbi.1008741.ref003","doi-asserted-by":"crossref","first-page":"1373","DOI":"10.1162\/089976603321780317","article-title":"Laplacian eigenmaps for dimensionality reduction and data representation","volume":"15","author":"M Belkin","year":"2003","journal-title":"Neural computation"},{"key":"pcbi.1008741.ref004","doi-asserted-by":"crossref","first-page":"5","DOI":"10.1016\/j.acha.2006.04.006","article-title":"Diffusion Maps","volume":"21","author":"R Coifman","year":"2006","journal-title":"Appl Comput Harmon Anal"},{"key":"pcbi.1008741.ref005","volume-title":"Machine learning: a probabilistic perspective","author":"KP Murphy","year":"2012"},{"issue":"4","key":"pcbi.1008741.ref006","doi-asserted-by":"crossref","first-page":"545","DOI":"10.1093\/bioinformatics\/btu674","article-title":"Statistical significance of variables driving systematic variation in high-dimensional data","volume":"31","author":"NC Chung","year":"2014","journal-title":"Bioinformatics"},{"key":"pcbi.1008741.ref007","article-title":"Learning the geometry of common latent variables using alternating-diffusion","author":"RR Lederman","year":"2015","journal-title":"Applied and Computational Harmonic Analysis"},{"issue":"3","key":"pcbi.1008741.ref008","doi-asserted-by":"crossref","first-page":"848","DOI":"10.1016\/j.acha.2017.12.006","article-title":"Latent common manifold learning with alternating diffusion: analysis and applications","volume":"47","author":"R Talmon","year":"2019","journal-title":"Applied and Computational Harmonic Analysis"},{"issue":"3","key":"pcbi.1008741.ref009","doi-asserted-by":"crossref","first-page":"588","DOI":"10.1137\/18M1218157","article-title":"Recovering hidden components in multimodal data with composite diffusion operators","volume":"1","author":"T Shnitzer","year":"2019","journal-title":"SIAM Journal on Mathematics of Data Science"},{"key":"pcbi.1008741.ref010","doi-asserted-by":"crossref","first-page":"346","DOI":"10.1016\/j.inffus.2018.01.007","article-title":"Alternating diffusion maps for multimodal data fusion","volume":"45","author":"O Katz","year":"2019","journal-title":"Information Fusion"},{"key":"pcbi.1008741.ref011","doi-asserted-by":"crossref","first-page":"127","DOI":"10.1016\/j.inffus.2019.08.005","article-title":"Multi-view diffusion maps","volume":"55","author":"O Lindenbaum","year":"2020","journal-title":"Information Fusion"},{"issue":"12","key":"pcbi.1008741.ref012","doi-asserted-by":"crossref","first-page":"2505","DOI":"10.1109\/TPAMI.2015.2408348","article-title":"Multimodal manifold analysis by simultaneous diagonalization of laplacians","volume":"37","author":"D Eynard","year":"2015","journal-title":"IEEE transactions on pattern analysis and machine intelligence"},{"issue":"9","key":"pcbi.1008741.ref013","doi-asserted-by":"crossref","first-page":"858","DOI":"10.1038\/nbt.2317","article-title":"Multiplexed mass cytometry profiling of cellular states perturbed by small-molecule regulators","volume":"30","author":"B Bodenmiller","year":"2012","journal-title":"Nature biotechnology"},{"issue":"4","key":"pcbi.1008741.ref014","doi-asserted-by":"crossref","first-page":"417","DOI":"10.1038\/nmeth.2869","article-title":"Highly multiplexed imaging of tumor tissues with subcellular resolution by mass cytometry","volume":"11","author":"C Giesen","year":"2014","journal-title":"Nature methods"},{"issue":"2","key":"pcbi.1008741.ref015","doi-asserted-by":"crossref","first-page":"160","DOI":"10.1002\/cyto.a.23053","article-title":"Imaging mass cytometry","volume":"91","author":"Q Chang","year":"2017","journal-title":"Cytometry part A"},{"key":"pcbi.1008741.ref016","doi-asserted-by":"crossref","unstructured":"Sun J, Ovsjanikov M, Guibas L. A concise and provably informative multi-scale signature based on heat diffusion. In: Computer graphics forum. vol. 28. Wiley Online Library; 2009. p. 1383\u20131392.","DOI":"10.1111\/j.1467-8659.2009.01515.x"},{"key":"pcbi.1008741.ref017","doi-asserted-by":"crossref","unstructured":"Aubry M, Schlickewei U, Cremers D. The wave kernel signature: A quantum mechanical approach to shape analysis. In: 2011 IEEE international conference on computer vision workshops (ICCV workshops). IEEE; 2011. p. 1626\u20131633.","DOI":"10.1109\/ICCVW.2011.6130444"},{"key":"pcbi.1008741.ref018","unstructured":"Krause H. Localization theory for triangulated categories. arXiv preprint arXiv:08061324. 2008;."},{"issue":"1","key":"pcbi.1008741.ref019","first-page":"1","article-title":"Random walks on graphs: A survey","volume":"2","author":"L Lov\u00e1sz","year":"1993","journal-title":"Combinatorics, Paul erdos is eighty"},{"issue":"1","key":"pcbi.1008741.ref020","doi-asserted-by":"crossref","first-page":"132","DOI":"10.1109\/TASL.2012.2215593","article-title":"Single-channel transient interference suppression with diffusion maps","volume":"21","author":"R Talmon","year":"2012","journal-title":"IEEE transactions on audio, speech, and language processing"},{"issue":"4","key":"pcbi.1008741.ref021","doi-asserted-by":"crossref","first-page":"732","DOI":"10.1109\/TASLP.2015.2405481","article-title":"Audio-visual voice activity detection using diffusion maps","volume":"23","author":"D Dov","year":"2015","journal-title":"IEEE\/ACM Transactions on Audio, Speech, and Language Processing"},{"issue":"1","key":"pcbi.1008741.ref022","doi-asserted-by":"crossref","first-page":"111","DOI":"10.1109\/JSTSP.2012.2232279","article-title":"Multiscale anomaly detection using diffusion maps","volume":"7","author":"G Mishne","year":"2012","journal-title":"IEEE Journal of selected topics in signal processing"},{"issue":"5","key":"pcbi.1008741.ref023","doi-asserted-by":"crossref","first-page":"1065","DOI":"10.1109\/TPAMI.2010.210","article-title":"Shape recognition with spectral distances","volume":"33","author":"MM Bronstein","year":"2010","journal-title":"IEEE Transactions on Pattern Analysis and Machine Intelligence"},{"key":"pcbi.1008741.ref024","unstructured":"Cheng X, Mishne G. Spectral Embedding Norm: Looking Deep into the Spectrum of the Graph Laplacian. arXiv preprint arXiv:181010695. 2018;."},{"key":"pcbi.1008741.ref025","unstructured":"Devroye L, Mehrabian A, Reddad T. The total variation distance between high-dimensional Gaussians. arXiv preprint arXiv:181008693. 2018;."},{"key":"pcbi.1008741.ref026","doi-asserted-by":"crossref","DOI":"10.3389\/fimmu.2019.02657","article-title":"Mass Cytometry Imaging for the Study of Human Diseases\u2014Applications and Data Analysis Strategies","volume":"10","author":"H Baharlou","year":"2019","journal-title":"Frontiers in Immunology"},{"issue":"2","key":"pcbi.1008741.ref027","doi-asserted-by":"crossref","first-page":"163","DOI":"10.1038\/s43018-020-0026-6","article-title":"Imaging mass cytometry and multiplatform genomics define the phenogenomic landscape of breast cancer","volume":"1","author":"HR Ali","year":"2020","journal-title":"Nature Cancer"},{"issue":"Nov","key":"pcbi.1008741.ref028","first-page":"2579","article-title":"Visualizing data using t-SNE","volume":"9","author":"Lvd Maaten","year":"2008","journal-title":"Journal of machine learning research"},{"issue":"9","key":"pcbi.1008741.ref029","doi-asserted-by":"crossref","first-page":"1475","DOI":"10.1016\/S0031-3203(99)00138-7","article-title":"Combining multiple classifiers by averaging or by multiplying?","volume":"33","author":"DM Tax","year":"2000","journal-title":"Pattern recognition"},{"key":"pcbi.1008741.ref030","volume-title":"The PageRank citation ranking: Bringing order to the web","author":"L Page","year":"1999"},{"issue":"6434","key":"pcbi.1008741.ref031","doi-asserted-by":"crossref","first-page":"1463","DOI":"10.1126\/science.aaw1219","article-title":"Slide-seq: A scalable technology for measuring genome-wide expression at high spatial resolution","volume":"363","author":"SG Rodriques","year":"2019","journal-title":"Science"},{"issue":"10","key":"pcbi.1008741.ref032","doi-asserted-by":"crossref","first-page":"987","DOI":"10.1038\/s41592-019-0548-y","article-title":"High-definition spatial transcriptomics for in situ tissue profiling","volume":"16","author":"S Vickovic","year":"2019","journal-title":"Nature methods"},{"key":"pcbi.1008741.ref033","first-page":"563338","volume-title":"bioRxiv","author":"S Vickovic","year":"2019"},{"issue":"10","key":"pcbi.1008741.ref034","doi-asserted-by":"crossref","first-page":"eaax5851","DOI":"10.1126\/sciadv.aax5851","article-title":"MIBI-TOF: A multiplexed imaging platform relates cellular phenotypes and tissue structure","volume":"5","author":"L Keren","year":"2019","journal-title":"Science Advances"},{"key":"pcbi.1008741.ref035","doi-asserted-by":"crossref","unstructured":"Liu Y, Yang M, Deng Y, Su G, Guo C, Zhang D, et al. High-Spatial-Resolution Multi-Omics Atlas Sequencing of Mouse Embryos via Deterministic Barcoding in Tissue. Available at SSRN 3466428. 2019;.","DOI":"10.2139\/ssrn.3466428"},{"issue":"4","key":"pcbi.1008741.ref036","doi-asserted-by":"crossref","first-page":"373","DOI":"10.1007\/BF01896401","article-title":"Embedding Riemannian manifolds by their heat kernel","volume":"4","author":"P B\u00e9rard","year":"1994","journal-title":"Geometric & Functional Analysis GAFA"},{"key":"pcbi.1008741.ref037","doi-asserted-by":"crossref","DOI":"10.4135\/9781412985475","volume-title":"Principal components analysis","author":"GH Dunteman","year":"1989"},{"key":"pcbi.1008741.ref038","doi-asserted-by":"crossref","unstructured":"Raviv D, Bronstein MM, Bronstein AM, Kimmel R. Volumetric heat kernel signatures. In: Proceedings of the ACM workshop on 3D object retrieval. ACM; 2010. p. 39\u201344.","DOI":"10.1145\/1877808.1877817"},{"key":"pcbi.1008741.ref039","unstructured":"Rustamov RM. Laplace-Beltrami eigenfunctions for deformation invariant shape representation. In: Proceedings of the fifth Eurographics symposium on Geometry processing. Eurographics Association; 2007. p. 225\u2013233."},{"issue":"3","key":"pcbi.1008741.ref040","doi-asserted-by":"crossref","first-page":"242","DOI":"10.4153\/CJM-1949-021-5","article-title":"Some properties of the eigenfunctions of the Laplace-operator on Riemannian manifolds","volume":"1","author":"S Minakshisundaram","year":"1949","journal-title":"Canadian Journal of Mathematics"},{"issue":"5","key":"pcbi.1008741.ref041","doi-asserted-by":"crossref","first-page":"604","DOI":"10.1145\/324133.324140","article-title":"Authoritative sources in a hyperlinked environment","volume":"46","author":"JM Kleinberg","year":"1999","journal-title":"Journal of the ACM (JACM)"},{"key":"pcbi.1008741.ref042","unstructured":"Wu Z, Song S, Khosla A, Yu F, Zhang L, Tang X, et al. 3d shapenets: A deep representation for volumetric shapes. In: Proceedings of the IEEE conference on computer vision and pattern recognition; 2015. p. 1912\u20131920."},{"key":"pcbi.1008741.ref043","doi-asserted-by":"crossref","unstructured":"Ma G, Lu CT, He L, Philip SY, Ragin AB. Multi-view graph embedding with hub detection for brain network analysis. In: 2017 IEEE International Conference on Data Mining (ICDM). IEEE; 2017. p. 967\u2013972.","DOI":"10.1109\/ICDM.2017.123"}],"updated-by":[{"DOI":"10.1371\/journal.pcbi.1008741","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2021,4,8]],"date-time":"2021-04-08T00:00:00Z","timestamp":1617840000000}}],"container-title":["PLOS Computational Biology"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1008741","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2021,4,8]],"date-time":"2021-04-08T17:10:07Z","timestamp":1617901807000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1008741"}},"subtitle":[],"editor":[{"given":"Min","family":"Xu","sequence":"first","affiliation":[],"role":[{"role":"editor","vocabulary":"crossref"}]}],"short-title":[],"issued":{"date-parts":[[2021,3,29]]},"references-count":43,"journal-issue":{"issue":"3","published-online":{"date-parts":[[2021,3,29]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pcbi.1008741","relation":{"has-preprint":[{"id-type":"doi","id":"10.1101\/2020.08.23.263467","asserted-by":"object"}]},"ISSN":["1553-7358"],"issn-type":[{"value":"1553-7358","type":"electronic"}],"subject":[],"published":{"date-parts":[[2021,3,29]]}}}