Graph of graphs analysis for multiplexed data with application to imaging mass cytometry

Abstract

Hyper spectral imaging, sensor networks, spatial multiplexed proteomics, and spatial transcriptomics assays is a representative subset of distinct technologies from diverse domains of science and engineering that share common data structures. The data in all these modalities consist of high-dimensional multivariate observations (m-dimensional feature space) collected at different spatial positions and therefore can be analyzed using similar computational methodologies. Furthermore, in many studies practitioners collect datasets consisting of multiple spatial assays of this type, each capturing such data from a single biological sample, patient, or hyper spectral image, etc. Each of these spatial assays could be characterized by several regions of interest (ROIs). The focus of this paper is on a particular application, imaging mass cytometry (IMC), which falls into this problem setup. 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 mdimensional 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 92% accuracy. This serves as empirical evidence that the graph of graphs approach enables us to integrate multiple ROIs and the intrarelationships between the features at each ROI, giving rise to an informative representation that is strongly associated with the phenotypic state of the entire image. Importantly, this approach is applicable to other modalities such as spatial transcriptomics.

Managing Resistance to Immune Checkpoint Inhibitors in Lung Cancer: Treatment and Novel Strategies

“…The most novel analytic tools allow integrating genomic, transcriptomic, and spatial localization data, proposing a highly complex evaluation of multiparametric TME and biomarkers to select better and tailor immunotherapy approaches. 130-132…”

Procrustes Analysis on the Manifold of SPSD Matrices for Data Sets Alignment

“…We note that in contrast to the previous steps (centering and scaling), which are completely label-free, this step requires subsets of X and Y that are sufficiently large and allow for accurate estimations of M x c and M y c , respectively. Alternatively, an unsupervised rotation by matching the second-order moments of the data sets as in [5], [27], can be considered.…”

“…Extending this idea from shapes to high-dimensional point clouds facilitated an appealing data alignment approach, which is simple, efficient, and mathematically tractable, and it does not require any rigid a-priori model assumptions or estimates of the whole distribution of the data. Indeed, data alignment using PA has been successfully applied to various fields, including Brain-Computer Interface (BCI) [2], genetics and bioinformatics [3], [4], [5], indoor navigation [6], face recognition [7], and hierarchical representation [5], [8], to name but a few.…”