Integrating multiplex immunofluorescent and mass spectrometry imaging to map myeloid heterogeneity in its metabolic and cellular context.

Cells often adopt different phenotypes, dictated by tissue-specific or local signals such as cell-cell and cell-matrix contacts or molecular micro-environment. This holds in extremis for macrophages with their high phenotypic plasticity. Their broad range of functions, some even opposing, reflects their heterogeneity, and a multitude of subsets has been described in different tissues and diseases. Such micro-environmental imprint cannot be adequately studied by single-cell applications, as cells are detached from their context, while histology-based assessment lacks the phenotypic depth due to limitations in marker combination. Here, we present a novel, integrative approach in which 15-color multispectral imaging allows comprehensive cell classification based on multi-marker expression patterns, followed by downstream analysis pipelines to link their phenotypes to contextual, micro-environmental cues, such as their cellular ("community") and metabolic ("local lipidome") niches in complex tissue. The power of this approach is illustrated for myeloid subsets and associated lipid signatures in murine atherosclerotic plaque.

Integrative multiomics analysis of human atherosclerosis reveals a serum response factor-driven network associated with intraplaque hemorrhage.

BACKGROUND: While single-omics analyses on human atherosclerotic plaque have been very useful to map stage- or disease-related differences in expression, they only partly capture the array of changes in this tissue and suffer from scale-intrinsic limitations. In order to better identify processes associated with intraplaque hemorrhage and plaque instability, we therefore combined multiple omics into an integrated model. METHODS: In this study, we compared protein and gene makeup of low- versus high-risk atherosclerotic lesion segments from carotid endarterectomy patients, as judged from the absence or presence of intraplaque hemorrhage, respectively. Transcriptomic, proteomic, and peptidomic data of this plaque cohort were aggregated and analyzed by DIABLO, an integrative multivariate classification and feature selection method. RESULTS: We identified a protein-gene associated multiomics model able to segregate stable, nonhemorrhaged from vulnerable, hemorrhaged lesions at high predictive performance (AUC >0.95). The dominant component of this model correlated with αSMA- PDGFRα+ fibroblast-like cell content (p = 2.4E-05) and Arg1+ macrophage content (p = 2.2E-04) and was driven by serum response factor (SRF), possibly in a megakaryoblastic leukemia-1/2 (MKL1/2) dependent manner. Gene set overrepresentation analysis on the selected key features of this model pointed to a clear cardiovascular disease signature, with overrepresentation of extracellular matrix synthesis and organization, focal adhesion, and cholesterol metabolism terms, suggestive of the model's relevance for the plaque vulnerability. Finally, we were able to corroborate the predictive power of the selected features in several independent mRNA and proteomic plaque cohorts. CONCLUSIONS: In conclusion, our integrative omics study has identified an intraplaque hemorrhage-associated cardiovascular signature that provides excellent stratification of low- from high-risk carotid artery plaques in several independent cohorts. Further study revealed suppression of an SRF-regulated disease network, controlling lesion stability, in vulnerable plaque, which can serve as a scaffold for the design of targeted intervention in plaque destabilization.