Graph of graphs analysis for multiplexed data with application to imaging mass cytometry.
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μm resolution. Each slide is a sp...
| Main Authors: | , , , , , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Public Library of Science (PLoS)
2021-03-01
|
| Series: | PLoS Computational Biology |
| Online Access: | https://doi.org/10.1371/journal.pcbi.1008741 |
| _version_ | 1818926846716149760 |
|---|---|
| author | Ya-Wei Eileen Lin Tal Shnitzer Ronen Talmon Franz Villarroel-Espindola Shruti Desai Kurt Schalper Yuval Kluger |
| author_facet | Ya-Wei Eileen Lin Tal Shnitzer Ronen Talmon Franz Villarroel-Espindola Shruti Desai Kurt Schalper Yuval Kluger |
| author_sort | Ya-Wei Eileen Lin |
| collection | DOAJ |
| description | 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μm resolution. Each slide is a spatial assay consisting of high-dimensional multivariate observations (m-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 m covariates and compute an m dimensional vector characterizing the steady state distribution among features. We then use all these m-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. |
| first_indexed | 2024-12-20T03:03:37Z |
| format | Article |
| id | doaj.art-1b79d21b07f74f38847553da6ed37362 |
| institution | Directory Open Access Journal |
| issn | 1553-734X 1553-7358 |
| language | English |
| last_indexed | 2024-12-20T03:03:37Z |
| publishDate | 2021-03-01 |
| publisher | Public Library of Science (PLoS) |
| record_format | Article |
| series | PLoS Computational Biology |
| spelling | doaj.art-1b79d21b07f74f38847553da6ed373622022-12-21T19:55:39ZengPublic Library of Science (PLoS)PLoS Computational Biology1553-734X1553-73582021-03-01173e100874110.1371/journal.pcbi.1008741Graph of graphs analysis for multiplexed data with application to imaging mass cytometry.Ya-Wei Eileen LinTal ShnitzerRonen TalmonFranz Villarroel-EspindolaShruti DesaiKurt SchalperYuval KlugerImaging 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μm resolution. Each slide is a spatial assay consisting of high-dimensional multivariate observations (m-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 m covariates and compute an m dimensional vector characterizing the steady state distribution among features. We then use all these m-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.https://doi.org/10.1371/journal.pcbi.1008741 |
| spellingShingle | Ya-Wei Eileen Lin Tal Shnitzer Ronen Talmon Franz Villarroel-Espindola Shruti Desai Kurt Schalper Yuval Kluger Graph of graphs analysis for multiplexed data with application to imaging mass cytometry. PLoS Computational Biology |
| title | Graph of graphs analysis for multiplexed data with application to imaging mass cytometry. |
| title_full | Graph of graphs analysis for multiplexed data with application to imaging mass cytometry. |
| title_fullStr | Graph of graphs analysis for multiplexed data with application to imaging mass cytometry. |
| title_full_unstemmed | Graph of graphs analysis for multiplexed data with application to imaging mass cytometry. |
| title_short | Graph of graphs analysis for multiplexed data with application to imaging mass cytometry. |
| title_sort | graph of graphs analysis for multiplexed data with application to imaging mass cytometry |
| url | https://doi.org/10.1371/journal.pcbi.1008741 |
| work_keys_str_mv | AT yaweieileenlin graphofgraphsanalysisformultiplexeddatawithapplicationtoimagingmasscytometry AT talshnitzer graphofgraphsanalysisformultiplexeddatawithapplicationtoimagingmasscytometry AT ronentalmon graphofgraphsanalysisformultiplexeddatawithapplicationtoimagingmasscytometry AT franzvillarroelespindola graphofgraphsanalysisformultiplexeddatawithapplicationtoimagingmasscytometry AT shrutidesai graphofgraphsanalysisformultiplexeddatawithapplicationtoimagingmasscytometry AT kurtschalper graphofgraphsanalysisformultiplexeddatawithapplicationtoimagingmasscytometry AT yuvalkluger graphofgraphsanalysisformultiplexeddatawithapplicationtoimagingmasscytometry |