Metabolomics plays a pivotal role in systems biology, and NMR is a central tool with high precision and exceptional resolution of chemical information. Most NMR metabolomic studies are based on 1H 1D spectroscopy, severely limited by peak overlap. 13C NMR benefits from a larger signal dispersion but is barely used in metabolomics due to ca. 6000-fold lower sensitivity. We introduce a new approach, based on hyperpolarized 13C NMR at natural abundance, that circumvents this limitation. A new untargeted NMR-based metabolomic workflow based on dissolution dynamic nuclear polarization (d-DNP) for the first time enabled hyperpolarized natural abundance 13C metabolomics. Statistical analysis of resulting hyperpolarized 13C data distinguishes two groups of plant (tomato) extracts and highlights biomarkers, in full agreement with previous results on the same biological model. We also optimize parameters of the semiautomated d-DNP system suitable for high-throughput studies.