Fig. 4

Individualized influence of beauty product application on skin metabolomics profiles over time. a Multivariate statistical analysis (principal coordinate analysis (PCoA)) comparing mass spectrometry data collected over 9 weeks from the skin of 11 individuals, all body parts, combined (first plot from the left) and then displayed separately (arm, armpits, face, feet). Color scale represents volunteer ID. The PCoA was calculated on all samples together, and subsets of the data are shown in this shared space and the other panels. b The molecular profiles collected over 9 weeks from all body parts, combined then separately (arm, armpits, face, feet). c Representative molecular profiles collected over 9 weeks from armpits of 11 individuals (volunteers 1, 2, 3, 4, 5, 6, 7, 9, 10, 11, 12). Color gradient in b and c represents timepoints (time 0 to time 9), ranging from the lightest orange color to the darkest one that represent the earliest (time 0) to the latest (time 9) timepoint, respectively. 0.5 timepoints represent additional timepoints where three selected volunteers were samples (volunteers 4, 9, and 10). PCoA plots were generated using the Bray–Curtis dissimilarity matrix and visualized in Emperor [28]. See also Additional file 1: Figure S3