Background Primary hyperparathyroidism (PHPT) may be related to a single gland disease or multiglandular disease, which requires specific treatments. At present, an operation is the only curative treatment for PHPT. Currently, there are no biomarkers available to identify these 2 entities (single vs. multiple gland disease). The aims of the present study were to compare (1) the tissue metabolomics profiles between PHPT and renal hyperparathyroidism (secondary and tertiary) and (2) single gland disease with multiglandular disease in PHPT using metabolomics analysis. Methods The method used was 1H high-resolution magic angle spinning nuclear magnetic resonance spectroscopy. Forty-three samples from 32 patients suffering from hyperparathyroidism were included in this study. Results Significant differences in the metabolomics profile were assessed according to PHPT and renal hyperparathyroidism. A bicomponent orthogonal partial least square-discriminant analysis showed a clear distinction between PHPT and renal hyperparathyroidism (R2Y = 0.85, Q2 = 0.63). Interestingly, the model also distinguished single gland disease from multiglandular disease (R2Y = 0.96, Q2 = 0.55). A network analysis was also performed using the Algorithm to Determine Expected Metabolite Level Alterations Using Mutual Information (ADEMA). Single gland disease was accurately predicted by ADEMA and was associated with higher levels of phosphorylcholine, choline, glycerophosphocholine, fumarate, succinate, lactate, glucose, glutamine, and ascorbate compared with multiglandular disease. Conclusion This study shows for the first time that 1H high-resolution magic angle spinning nuclear magnetic resonance spectroscopy is a reliable and fast technique to distinguish single gland disease from multiglandular disease in patients with PHPT. The potential use of this method as an intraoperative tool requires specific further studies.