Representative outcomes from the untargeted analysis of VAMS devices, following the procedure described above, are presented here.
The data from the triplicate analysis of a pooled sample prepared from venous blood of 10 individuals, sampled by VAMS, were evaluated, and the characteristics of the resulting metabolic fingerprints are shown. The quality of the data and the captured information from the VAMS devices were assessed using two primary criteria: the total number of detected chemical features and the reproducibility of the features (n = 3), expressed as the percentage relative standard deviation (%RSD). The total number of features was used to assess the breadth of extraction and chemical coverage of the solvent, while the %RSD was monitored to ensure the method's technical precision and robustness; a lower %RSD indicates a highly reproducible extraction process that minimizes analytical noise.
After raw data extraction and processing, a matrix of m/z-retention time (RT) features, along with their intensities, was obtained. Only features with a detectable signal in two out of three replicates were considered, after proper pre-processing and blank-filtering steps. Figure 2A,B displays m/z-over-RT plots for the same sample analyzed according to the described protocols for both LC-MS and GC-MS.
The high number of features eluted early in the RPLC run is attributed to highly hydrophilic metabolites extracted by the applied protocol but poorly retained on the reversed-phase column. Although early-eluting metabolites pose identification challenges, accurate annotations are still possible by using MS/MS fragmentation data alongside comparisons with spectral databases and in-house reference standards. Furthermore, a significant portion of this hydrophilic fraction is expected to be captured via GC-MS analysis. Ultimately, this represents an inherent limitation of the field, as the sheer complexity of the metabolome necessitates the use of multiple analytical platforms.
The multi-platform evaluation of the methanolic extracts clearly demonstrates the highly complementary nature of the LC and GC analytical workflows, with the performance metrics within the expected range for a comprehensive metabolomics profiling. Figure 3A,B depicts the overall feature numbering and the RSD% distribution for each platform. The LC platform provided extensive chemical coverage, resolving 4,169 features, which aligns with its established capacity to capture a vast array of polar, semi-polar, and non-volatile compounds. Concurrently, the GC platform displayed a smaller subset of features, with 187. This was expected, since GC can only detect volatile compounds amenable to derivatization.
Regarding reproducibility, as shown in the violin plot, the overall distribution and spread of RSD% values are remarkably similar between LC-MS and GC-MS. The average %RSD was typical for both platforms, at around 29.7% for LC (2,655 out of 4,169 had an %RSD lower than 30%) and 32.3% for GC (114 out of 187 had an %RSD lower than 30%). While a small number of features did display higher RSD values on both sides, this variation is common and expected in untargeted metabolomics. These higher values typically belong to low-abundance metabolites that are close to the instrument's limit of detection. Overall, the tight distribution profiles confirm satisfactory reproducibility for both analytical setups.
The chemical space of metabolites captured by the protocol can be considered by the metabolite’s annotations. Various compounds were detected by LC-MS and GC-MS, as expected given the inherent complementarity of the two platforms. A complete list of metabolite annotations for VAMS devices is provided in Supplementary Table 1 and Supplementary Table 2. In total, 97 metabolites were detected and annotated (level 1 and 2, as proposed by Sumner et al. (2007)22) by this protocol via the two techniques, in addition to other level 3 putatively characterized compounds. GC annotations were further confirmed by comparing the experimental retention indexes – obtained through comparison with FAMES – with tabulated ones23. Figure 4A,B shows the subclass distribution of the annotated metabolites for both methods, while Figure 4C shows the overlap between the techniques, highlighting their complementarity. Few key intermediary metabolites can be cross-validated by both techniques, while structurally unique or platform-specific features can fill gaps in central carbon and secondary metabolic pathways. Although amino acids and peptides are the major chemical classes in both analyses, only a few metabolites overlap between the techniques. Many metabolite classes are instrument-specific, such as acyl carnitines in LC-MS and hexoses in GC-MS.
To demonstrate the efficiency of methanol as an extraction solvent for both techniques, comparative trials were performed on VAMS samples using various solvent systems, followed by LC-MS and GC-MS analyses. Three distinct solvent systems were compared with methanol in triplicate: pure acetonitrile (ACN), a methanol–water mixture (60:40, v/v), and an acetonitrile–methanol mixture (70:30, v/v). The selection of these solvents and mixtures was informed by preliminary experiments conducted by our group, as well as previously published literature13,14,24,25. The extraction efficiency of each system was evaluated using three key parameters: the total number of detected features, overall feature intensity, and feature reproducibility, expressed as the relative standard deviation (RSD, %).
For the LC-MS workflow, the highest number of features was detected using MeOH-H2O (60:40, v/v), with an average of 4,514 features across triplicates, followed closely by pure MeOH (4,169 features). The ACN–MeOH (70:30, v/v) mixture yielded an average of 3,502 features, while pure ACN yielded the lowest average feature count, 2,836. A graphical representation of these feature counts per replicate is illustrated in Figure 5A.
To examine precision characteristics for each extraction condition, feature intensities and RSDs (%) from triplicate injections were analyzed. A comparison of the total feature abundance and the average %RSD per extract solvent for LC-MS is shown in Figure 5A. The log2 signal distributions were similar for all extracts, with MeOH-H2O 60:40 (v/v), ACN, and MeOH yielding slightly higher abundances than ACN-MeOH 70:30 (v/v). MeOH yielded a lower overall average %RSD, indicating that the extraction was more reproducible with this solvent.
A similar assessment of the efficiency of the extraction solvent system was performed in the GC-MS data. Once more, for the evaluation of the number of features, only those with a signal detected in at least two out of three replicates were considered. The solvent yielding the highest number of features was ACN-MeOH 70:30 (v/v) (216 features) followed closely by MeOH (187 features). A graphical illustration of the number of detected features per extraction solvent and BµS is shown in Figure 5B.
Regarding the overall abundance of the detected features, MeOH appeared to be the most efficient, yielding higher abundance, as shown in Figure 5B. MeOH-H₂O 60:40 (v/v) was the second best, followed by ACN-MeOH 70:30 (v/v). ACN was the least effective solvent. Regarding feature %RSD, MeOH-H₂O 60:40 (v/v) provided the lowest RSDs for VAMS (RSD = 38.57%), closely followed by ACN (42.65%), MeOH (44.14%), and ACN-MeOH 70:30 (v/v) (49.14%). These values are represented in Figure 5B.
The above align well with previously published literature. Couacault et al. (2026)19 reported a comparable number of features in the RP-LC-MS analysis of DBS and VAMS using a similar extraction solvent composition (rehydration with 80 µL of water, followed by the addition of 320 µL of MeOH). Similarly, Cala et al. (2017)11 detected a comparable feature count (153) via GC-MS analysis of VAMS tips. Although most existing studies did not evaluate the precision of these devices' analyses, our overall findings align with expectations for an untargeted analysis. Furthermore, the LC results agree with those reported by Bossi et al. (2025)26 for the HILIC analysis of DBS and qDBS microsamples, in which most detected features had RSDs < 30%.
Overall, these results demonstrate that methanol-based extraction provides a robust balance of feature counts, intensity, and precision (average RSD, %) for both LC-MS and GC-MS analysis of VAMS samples. Therefore, a single methanol extraction protocol can streamline dual sample preparation, facilitate the parallel acquisition of complementary metabolic profiles, and maximize information output.

Figure 1: Protocol schematics showing (A) Sample extraction protocol, (B) LC-MS reconstitution, and (C) GC-MS/MS derivatization. Created in BioRender. Marques, D. (2026) https://BioRender.com/dprdv18, https://biorender.com/xu2qh34, https://biorender.com/7ztz2qv. Please click here to view a larger version of this figure.

Figure 2: m/z and rt comparisons of a representative sample for (A) LC-MS and (B) GC-MS analysis. Each dot represents a single signal before it is grouped into features. The colors represent the intensity of the MS1 signal for a given feature, with intensity increasing from blue to red. Please click here to view a larger version of this figure.

Figure 3: Feature coverage and reproducibility across analytical platforms. (A) Feature Coverage: Total number of unique features detected by LC (teal) and GC (orange), with exact counts annotated above each bar. (B) Feature Reproducibility: Violin plots showing the distribution density of relative standard deviations (RSD%) for all features. Embedded box plots indicate the median and interquartile ranges; individual data points are omitted for visual clarity. Please click here to view a larger version of this figure.

Figure 4: Histograms showing the number of compounds of the different chemical classes detected by (A) LC-MS and (B) GC-MS. (C) Overlap of metabolites annotated by levels 1 and 2 between platforms. Please click here to view a larger version of this figure.

Figure 5: Extraction solvent results from the LC-MS (A) and GC-MS (B) analysis of VAMS devices. The bar plots on the top of the figure represent the number of features yielded from each extraction solvent applied to VAMS devices per replicate. The plots at the bottom of the figure show the sum of intensities and the average RSD (%) for each extraction solvent used with the VAMS device. Circle size represents the total sum intensity of detected features on a log2 scale, reflecting the extraction efficiency. The color gradient indicates analytical precision, expressed as the Relative Standard Deviation (RSD), with yellow representing higher precision and red indicating greater variability. Please click here to view a larger version of this figure.
Table 1: Troubleshooting table. Overview of common sampling problems, reasons, and solutions of sampling using VAMS devices. Please click here to download this file.
Table 2: RP eluents and wash solutions. Overview of the mobile phases, seal wash, purge solvent, and wash solvent compositions utilized for chromatographic analysis. Please click here to download this file.
Supplementary File 1: Additional Experimental Details. LC-MS and GC-MS parameters. Please click here to download this file.
Supplementary Table 1: LC-MS annotations. Compounds identified with the current protocol on the LC-MS platform. Please click here to download this file.
Supplementary Table 2: GC-MS annotations. Compounds identified with the current protocol on the GC-MS platform. Please click here to download this file.