Mapping tumor metabolic remodeling and their spatial crosstalk with surrounding non-tumor cells can fundamentally improve our understanding of tumor biology, facilitates the designing of advanced therapeutic strategies. Here, we present an integration of mass spectrometry imaging-based spatial metabolomics and lipidomics with microarray-based spatial transcriptomics to hierarchically visualize the intratumor metabolic heterogeneity and cell metabolic interactions in same gastric cancer sample. Tumor-associated metabolic reprogramming is imaged at metabolic-transcriptional levels, and maker metabolites, lipids, genes are connected in metabolic pathways and colocalized in the heterogeneous cancer tissues. Integrated data from spatial multi-omics approaches coherently identify cell types and distributions within the complex tumor microenvironment, and an immune cell-dominated “tumor-normal interface” region where tumor cells contact adjacent tissues are characterized with distinct transcriptional signatures and significant immunometabolic alterations. Our approach for mapping tissue molecular architecture provides highly integrated picture of intratumor heterogeneity, and transform the understanding of cancer metabolism at systemic level.

Background: Plant natural products have been co-opted for millennia by humans for various uses such as flavor, fragrances, and medicines. These compounds often are only produced in relatively low amounts and are difficult to chemically synthesize, limiting access. While elucidation of the underlying biosynthetic processes might help alleviate these issues (e.g., via metabolic engineering), investigation of this is hindered by the low levels of relevant gene expression and expansion of the corresponding enzymatic gene families. However, the often-inducible nature of such metabolic processes enables selection of those genes whose expression pattern indicates a role in production of the targeted natural product. Results: Here, we combine metabolomics and transcriptomics to investigate the inducible biosynthesis of the bioactive diterpenoid tanshinones from the Chinese medicinal herb, Salvia miltiorrhiza (Danshen). Untargeted metabolomics investigation of elicited hairy root cultures indicated that tanshinone production was a dominant component of the metabolic response, increasing at later time points. A transcriptomic approach was applied to not only define a comprehensive transcriptome (comprised of 20,972 non-redundant genes), but also its response to induction, revealing 6,358 genes that exhibited differential expression, with significant enrichment for up-regulation of genes involved in stress, stimulus and immune response processes. Consistent with our metabolomics analysis, there appears to be a slower but more sustained increased in transcript levels of known genes from diterpenoid and, more specifically, tanshinone biosynthesis. Among the co-regulated genes were 70 transcription factors and 8 cytochromes P450, providing targets for future investigation. Conclusions: Our results indicate a biphasic response of Danshen terpenoid metabolism to elicitation, with early induction of sesqui- and tri- terpenoid biosynthesis, followed by later and more sustained production of the diterpenoid tanshinones. Our data provides a firm foundation for further elucidation of tanshinone and other inducible natural product metabolism in Danshen.