Fabry disease (FD) is an X-linked lysosomal disease due to a deficiency in the activity of the lysosomal-galactosidase A (GalA), a key enzyme in the glycosphingolipid degradation pathway. FD is a complex disease with a poor genotype–phenotype correlation. In the early stages, FD could involve the peripheral nervous system (acroparesthesias and dysautonomia) and the ski (angiokeratoma), but later kidney, heart or central nervous system impairment may significantly decrease life expectancy. The advent of omics technologies offers the possibility of a global, integrated and systemic approach well-suited for the exploration of this complex disease. In this narrative review, we will focus on the main metabolomic studies, which have underscored the importance of detecting biomarkers for a diagnostic and prognostic purpose in FD. These investigations are potentially useful to explain the wide clinical, biochemical and molecular heterogeneity found in FD patients. Moreover, the quantitative mass spectrometry methods developed to evaluate concentrations of these biomarkers in urine and plasma will be described. Finally, the complex metabolic biomarker profile depicted in FD patients will be reported, which varies according to gender, types of mutations, and therapeutic treatment.