Single-Cell Metabolomics: A Comprehensive Overview
Metabolomics aims to systematically identify and quantify small-molecule metabolites to
elucidate metabolic states and regulatory mechanisms
within biological systems. While
nuclear magnetic resonance (NMR)
spectroscopy and
mass spectrometry (MS) are the primary
analytical
platforms, MS has become the method of choice due to its superior sensitivity, broad coverage,
and structural elucidation capability.
Conventional
LC-MS-based metabolomics typically analyzes bulk
cell populations, yielding averaged signals that obscure intercellular differences.
In reality, each cell operates as an independent metabolic unit; this metabolic variability,
known as cellular heterogeneity, underpins critical
processes such as drug resistance, metastasis, and cell fate decisions. Furthermore, rare cell
populations like circulating tumor cells, cancer
stem cells, and antigen-specific T cells, play pivotal roles in disease progression despite
their low abundance, making them challenging to capture
via bulk analysis. Consequently, performing metabolomic analysis at single-cell resolution has
become increasingly critical
[1].
Single-Cell Sampling and Ionization Techniques
MS is the cornerstone of single-cell metabolomics. A typical single-cell MS experiment
involves several key steps, including sampling,
ionization, detection, and data analysis, among which the precise extraction and ionization of
metabolites from an individual cell represent
the primary technical hurdles. As illustrated in Figure 1, a variety of complementary strategies
have been developed to address these
challenges. Vacuum-based techniques, such as secondary ion MS (SIMS) and matrix-assisted laser
desorption/ionization (MALDI), enable
simultaneous sampling and ionization with sub-micrometer spatial resolution. In contrast,
ambient-pressure approaches, including laser
ablation electrospray ionization (LAESI) and the Single-probe technique, allow for the in situ
analysis of living cells. Additionally,
methods such as capillary electrophoresis–electrospray ionization–mass spectrometry (CE–ESI–MS)
enhance metabolite identification by
incorporating an electrophoretic separation dimension. Collectively, these techniques constitute
a versatile and complementary toolbox
for single-cell metabolic analysis.
Figure 1. Examples of single cell MS metabolomics techniques[1].
(A) SIMS; (B) MALDI MS; (C) LAESI MS; (D) MALDI/C₆₀-SIMS hybrid Q-TOF; (E) Single-probe MS; (F)
T-probe MS; (G) Dean flow assisted cell ordering; (H) CE-ESI-MS interface.
Single-Cell Metabolomics Data Analysis
Single-cell MS data are characterized by low signal intensity, high noise levels, and
substantial variability, necessitating rigorous
computational pipelines (Figure 2). A standard workflow typically includes raw data
preprocessing (denoising, peak alignment,
and normalization), differential metabolite selection using univariate statistical analyses, and
the visualization of cellular
subpopulations through multivariate dimensionality reduction methods (such as t-SNE and PCA).
This is followed by metabolite
identification and pathway enrichment analysis to transform high-dimensional, sparse MS data
into robust biological insights.
However, despite advances in sampling, ionization, and data analysis strategies, single-cell
metabolomics inherently captures metabolite
abundance distributions as "snapshots" at specific time points. It remains limited in its
ability to directly resolve the dynamic
directionality and flux changes within metabolic networks. To achieve a deeper understanding of
metabolic function, dynamic tracing of
metabolic fluxes is therefore essential.
Figure 2. A generalized single cell metabolomics data analysis workflow that consists of raw
data pre-processing, univariate analysis and multivariate analysis[1].
Stable Isotope Tracing: Deciphering Metabolic Flux
To elucidate the dynamic activities of metabolic networks — specifically the rates and
pathways of metabolite transformations — researchers
have developed several sophisticated strategies. Among these, SIT has emerged as a core tool
for analyzing the dynamic behavior of
metabolic networks. By tracing the transfer and distribution of labeled atoms across the
network, SIT, when paired with computational
modeling, enables the quantitative inference of reaction rates (i.e., metabolic fluxes).
Principles and Methodologies of SIT
The core principle of SIT involves introducing metabolic precursors (such as glucose or
glutamine) labeled with stable
isotopes (e.g.,
13C,
15N,
2H) into a biological system to trace the transfer and
distribution of these labeled atoms within the metabolic
network over time. By using MS or NMR to detect mass shifts and isotopologue distribution
patterns, and by combining this data with kinetic
modeling, researchers can quantitatively map the activity and flux distribution of metabolic
pathways.
This approach represents a paradigm shift from static "metabolite
concentration" analysis to dynamic "
metabolic flux" resolution.
To use a transportation analogy: just as vehicle density (concentration) does not always
reflect actual traffic flow (flux), a high level
of a metabolite may result from blocked downstream consumption (reduced flux) rather than
active upstream synthesis (increased flux) (Figure 3).
Figure 3. Metabolite levels versus metabolic flux[2].
A typical SIT experimental design involves several key steps. First, the choice of labeled
substrate and labeling position must be guided
by the specific scientific question. For example,
[1,
2-13C] glucose can specifically trace flux through the pentose phosphate pathway,
while
[U-13C] glutamine is widely used to assess
tricarboxylic acid (TCA) cycle activity and reductive carboxylation.
After introducing
the labeled substrates into cell cultures or
in vivo systems, samples are collected at
multiple time points to capture dynamic labeling
kinetics.
Next, isotopologue distributions of extracted metabolites are analyzed via MS to trace the
fate of labeled atoms within metabolic networks.
For instance, the appearance of
13C-labeled carbons from
[U-13C] glucose in TCA cycle intermediates, such as citrate
or succinate, allows for the
quantification of the contribution of glucose oxidation to the TCA cycle. The table below
summarizes recommended tracer selection strategies and
key interpretative metrics for various metabolic pathways-including glycolysis, the pentose
phosphate pathway, fatty acid synthesis, and one-carbon
metabolism-serving as a practical guide for experimental design.
Table 1. Isotopic tracers for measuring pathway activities[2].
| Application |
Tracer |
Metabolite readouts |
Explanation |
| Pentose phosphate pathway
(PPP) |
| PPP overflow |
[1,2-13C] glucose |
Lactate M+1, M+2 |
Flux through the combined oxidative and non-oxidative PPP generates M+1 lactate
from [1,2-13C]glucose, while glycolysis generates only M+2 lactate (Lee et al., 1998). PPP overflow / glycolysis ≈ LacM+1 /
LacM+2. |
| Source of ribose (oxidative vs non-oxidative branch of PPP) |
[1,2-13C] glucose |
Ribose phosphate M+1, M+2 |
The oxPPP make M1 ribose phosphate; the non-oxPPP makes M2. Ratio of M1/M2
depends on the gross flux (net flux + exchange flux) of each branch:
Reversibility of the non-oxPPP can make M2 even if all net ribose production is
by oxPPP. |
| Glycolysis, TCA and
gluconeogenesis |
| Glycolytic rate |
[U-13C] glucose |
Fructose-1,6-bisphosphate (FBP) Dihydroxyacetone phosphate 3-phosphoglycerate
|
Higher flux yields faster labeling. Labeling results should be confirmed by
glucose uptake and lactate excretion measurements. |
| Reversibility of glycolysis |
50% : 50% mix of [U-12C] : [U-13C] glucose |
Glucose-6-phosphate M+3 FBP M+3 |
Feeding a mixture of labeled and unlabeled glucose results in unlabeled and M+3
triose phosphates. Reversibility of aldolase produces M+3 FBP. Fructose
bisphosphatase activity yields M3 glucose-6-phosphate (Park et al., 2016). |
| Gluconeogenesis |
[U-13C] lactate
[U-13C] glutamine |
Glucose-6-phosphate M+2, M+3 3-phosphoglycerate M+2, M+3 |
Lactate and glutamine are major TCA substrates. Flux from TCA to glycolysis
catalyzed by PEPCK results in triose phosphate labeling. Fructose bisphosphatase
activity then makes labeled hexose phosphates. |
| Pyruvate carboxylase contribution to TCA |
[3-13C] glucose
[1-13C] pyruvate |
Aspartate M+3 Malate M+3 |
C1 of pyruvate comes from glucose C3/C4. Pyruvate C1 is lost in making
acetyl-CoA, but can enter TCA via pyruvate carboxylase which makes M1
oxaloacetate and thus M1 aspartate and M1 malate (Sellers et al., 2015). |
| Reductive carboxylation (“backwards” TCA flux) |
[U-13C] glutamine
[1-13C] glutamine |
Citrate M+5, Malate M+3 or Citrate M+1, Malate M+1 |
Reductive carboxylation of α-ketoglutarate (derived from labeled glutamine)
produces M+5 (or M+1) citrate, and subsequent ATP citrate lyase produces M+3 (or
M+1) malate (Yoo et al., 2008). |
| TCA carbon sources |
[U-13C] nutrients |
Succinate Malate Citrate α-ketoglutarate |
Carbon enrichment (number of 13C atoms versus total carbon atoms) reflects
carbon contribution from the nutrient; useful in vivo with correction for
circulating nutrient enrichment (Davidson et al.,
2016; Faubert et al., 2017; Hui et al., 2017). |
In summary, SIT elevates analysis from the static "metabolite concentration" to dynamic "
metabolic flux" resolution. As such, it serves as a
pivotal tool for quantitatively characterizing central carbon metabolism, assessing the
relative contributions of various nutrient substrates,
and deciphering the functional states of specific metabolic pathways.
Tumor Metabolism Study: Case-Based Applications
In tumor biology, metabolic heterogeneity is often tightly coupled with cellular
interactions and spatial organization within the tumor
microenvironment. The following empirical case studies demonstrate how integrating SIT with
single-cell analytical strategies enables the
resolution of these highly dynamic metabolic processes within the spatiotemporal context of
the native tumor microenvironment.
Case Study 1: Synergizing Single-Cell Metabolomics and SIT: Mapping the Dynamic Landscape of
Tumor Heterogeneity
Zhang et al. established a dynamic single-cell metabolomics analysis system to
investigate the metabolic interactions between tumor cells and
other cell types within the tumor microenvironment. SIT and single-cell metabolic analysis
are complementary core technologies for revealing
the metabolic complexity of cancer and its treatment strategies[3].
SIT (using substrates like
13C-glucose and
13C-glutamine) is
the gold standard for the quantitative analysis of metabolic flux. It plays two key
roles in cancer research: first, in systematically elucidating metabolic reprogramming
characteristics, such as precisely tracing aberrant glucose
flows in the Warburg effect or identifying tumor-associated glutamine dependence; and
second, in dynamically revealing treatment response and drug
resistance mechanisms. By comparing metabolic flux changes before and after treatment,
researchers can discover how cancer cells evade therapy
through metabolic adaptations (such as switching energy sources and activating compensatory
pathways), thereby providing direct evidence for the
design of metabolism-targeted combination therapies (e.g., GLS inhibitors combined with
chemotherapy).
Single-cell metabolic analysis (such as MS and single-cell MS) directly addresses the
challenge of tumor metabolic heterogeneity.
Its main functions are: first, to identify key metabolic subpopulations driving tumor
progression at single-cell resolution, such as
pinpointing cancer stem cells or drug-resistant cells with unique metabolic phenotypes; and
second, to analyze intercellular metabolic
interactions within the tumor microenvironment, visually illustrating nutrient competition
and metabolite exchange between tumor cells,
immune cells, and stromal cells. This lays the foundation for developing precision therapies
that target specific malignant subpopulations.
The integration of these two approaches represents a future direction: combining the
high-throughput, quantitative dynamic advantages of SIT
with the resolution of single-cell metabolism holds promise for creating a panoramic,
high-precision tumor metabolic atlas. This will enable a
systemic understanding of metabolic heterogeneity and facilitate the design of personalized
combination therapies that simultaneously address
multiple metabolic vulnerabilities to overcome drug resistance.
Figure 4. Schematic workflow of the dynamic single-cell metabolomics system[3].
Case Study 2: Mapping Tumor Microenvironment Lipogenesis Heterogeneity via Single-Cell
13C-SpaceM
Buglakova
et al. primarily employed
U-13C6-glucose as
a stable isotope tracer. Utilizing the innovative
13C-SpaceM technique
(integrating spatially-resolved MALDI imaging MS with all-ion fragmentation), the
researchers achieved tracing and imaging of
de
novo fatty acid synthesis metabolic activity in tumor tissues at single-cell and
near-single-cell resolution
[4].
Key findings include:
1) In an in vitro liver cancer cell model, the study revealed that the hypoxic
microenvironment significantly reduces the capacity of cells for glucose-dependent fatty
acid synthesis.
2) In an ACLY knockdown model, the team quantified the labeling degree of the cytosolic
lipogenic acetyl-CoA pool at the single-cell level, uncovering substantial spatial
heterogeneity;
3) In an IDH-mutant glioma mouse model, the study achieved near-single-cell resolution
metabolic imaging in tissue for the first time, demonstrating a strong induction of de novo
fatty acid synthesis in tumor regions, increased uptake of monounsaturated and essential
fatty acids, and significant spatial heterogeneity in acetyl-CoA pool labeling.
The 13C-SpaceM method pioneers SIT at near-single-cell spatial resolution, providing a
powerful tool for deciphering metabolic heterogeneity and revealing spatial patterns of
metabolic reprogramming within complex tissue microenvironments. This not only deepens
understanding of tumor metabolic adaptation but also lays the groundwork for discovering new
metabolic dependencies and potential therapeutic targets.
Figure 5. 13C-SpaceM workflow as applied to interrogate de novo fatty acid synthesis[4].
Recommended Stable Isotope-Labeled Compounds Products
Note: MCE can provide products for research use only. We do not sell to patients.
Summary
The integration of single-cell metabolomics with SIT provides a high-resolution research
paradigm for analyzing tumor metabolic
heterogeneity and metabolic reprogramming. By tracing isotope-labeled metabolic fluxes at
single-cell and spatial scales, this strategy
enables the precise characterization of the dynamic metabolic pathway activities. This, in turn,
allows for the identification of key
metabolic features closely associated with tumor initiation, progression, drug resistance, and
therapeutic response.
With continuous advancements in single-cell isolation and detection technologies, as well as
more sensitive isotopic analytical methods,
this integrated approach holds significant potential for exploring tumor metabolic regulatory
mechanisms and advancing precision
metabolic-targeted therapies.
LinkedIn
Researchgate
Facebook
X
Instagram
Mail
Single-Cell Metabolomics: A Comprehensive Overview
Stable Isotope Tracing: Deciphering Metabolic Flux
Tumor Metabolism Study: Case-Based Applications
