Accumulating evidence demonstrates important roles for metabolism in cell fate determination. However, it is a challenge to assess metabolism at a spatial resolution that acknowledges both heterogeneity and cellular dynamics in its tissue microenvironment. Using a multi-omics platform to study cell-type-specific dynamics in metabolism in complex tissues, we describe the metabolic trajectories during nephrogenesis in the developing human kidney. Exploiting in situ analysis of isotopic labeling, a shift from glycolysis toward fatty acid β-oxidation was observed during the differentiation from the renal vesicle toward the S-shaped body and the proximal tubules. In addition, we show that hiPSC-derived kidney organoids are characterized by a metabolic immature phenotype that fails to use mitochondrial long-chain fatty acids for energy metabolism. Furthermore, supplementation of butyrate enhances tubular epithelial differentiation and maturation in cultured kidney organoids. Our findings highlight the relevance of understanding metabolic trajectories to efficiently guide stem cell differentiation.

A common drawback of metabolic analyses of complex biological samples is the inability to consider cell-to-cell heterogeneity in the context of an organ or tissue. To overcome this limitation, we present an advanced high-spatial-resolution metabolomics approach using matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) combined with isotope tracing. This method allows mapping of cell-type-specific dynamic changes in central carbon metabolism in the context of a complex heterogeneous tissue architecture, such as the kidney. Combined with multiplexed immunofluorescence staining, this method can detect metabolic changes and nutrient partitioning in targeted cell types, as demonstrated in a bilateral renal ischemia–reperfusion injury (bIRI) experimental model. Our approach enables us to identify region-specific metabolic perturbations associated with the lesion and throughout recovery, including unexpected metabolic anomalies in cells with an apparently normal phenotype in the recovery phase. These findings may be relevant to an understanding of the homeostatic capacity of the kidney microenvironment. In sum, this method allows us to achieve resolution at the single-cell level in situ and hence to interpret cell-type-specific metabolic dynamics in the context of structure and metabolism of neighboring cells.