Next Article in Journal
Chemical Fingerprints of Emotional Body Odor
Previous Article in Journal
Integrative Metabolomics for Assessing the Effect of Insect (Hermetia illucens) Protein Extract on Rainbow Trout Metabolism
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Dried Blood Spot-Based Metabolomics: A Targeted, Broad-Spectrum, Single-Injection Method

1
The Mitochondrial and Metabolic Disease Center, School of Medicine, University of California, San Diego, CA 92103, USA
2
Department of Medicine, School of Medicine, University of California, San Diego, CA 92103, USA
3
Department of Neuroscience, School of Medicine, University of California, San Diego, CA 92103, USA
4
Department of Pediatrics, School of Medicine, University of California, San Diego, CA 92103, USA
*
Authors to whom correspondence should be addressed.
Metabolites 2020, 10(3), 82; https://doi.org/10.3390/metabo10030082
Submission received: 2 January 2020 / Revised: 19 February 2020 / Accepted: 24 February 2020 / Published: 27 February 2020

Abstract

:
Dried blood spots (DBS) have proven to be a powerful sampling and storage method for newborn screening and many other applications. However, DBS methods have not yet been optimized for broad-spectrum targeted metabolomic analysis. In this study, we developed a robust, DBS-based, broad-spectrum, targeted metabolomic method that was able to measure over 400 metabolites from a 6.3 mm punch from standard Whatman 903TM filter paper cards. The effects of blood spot volumes, hematocrit, vacutainer chemistry, extraction methods, carryover, and comparability with plasma and fingerstick capillary blood samples were analyzed. The stability of over 400 metabolites stored under varying conditions over one year was also tested. No significant impacts of blood volume and hematocrit variations were observed when the spotted blood volume was over 60 µL and the hematocrit was between 31% and 50%. The median area under the curve (AUC) of metabolites in the DBS metabolome declined by 40% in the first 3 months and then did not decline further for at least 1 year. All originally detectable metabolites remained within detectable limits. The optimal storage conditions for metabolomic analysis were −80 °C with desiccants and without an O2 scavenger. The method was clinically validated for its potential utility in the diagnosis of the mitochondrial disease mitochondrial encephalomyopathy, lactic acidosis, and stroke-like episodes (MELAS). Our method provides a convenient alternative to freezing, storing, and shipping liquid blood samples for comparative metabolomic studies.

1. Introduction

Metabolomics is emerging as an important tool in clinical chemistry to identify biomarkers, find new roles of metabolism in disease, and discover novel metabolic drug targets [1,2,3]. Blood-based samples, including serum and plasma, are probably the best studied biofluids for metabolomics [4]. However, serum and plasma preparation requires laboratory healthcare professionals and stringent storage conditions to avoid metabolite degradation. Dried blood spots have been used for newborn screening in millions of children since the 1960s [5]. In this paper, we explore the utility of dried blood spot (DBS) analysis in metabolomic studies.
Samples prepared and stored as dried blood spots on standardized filter paper cards may be an attractive alternative to storage, shipment, and analysis of conventional liquid venous blood samples for many studies. One advantage of dried blood spots is that quantitative results for many targeted molecules can still be obtained within relatively large margins of variation in the original amount of blood and hematocrit spotted. This occurs because of reproducible spreading of the sample by the fibers of standardized filter paper. For example, a 1/8th inch (3.175 mm) diameter punch from a standard dried blood spot on Whatman 903TM paper contains the residue from 3.1 µL of blood [6], without the need for great volumetric precision in the amount of blood applied as long as the amount is sufficient (50–100 µL; 1–2 drops) to fill at least 90% of the 12-mm dotted circles inscribed on the cards [7].
Fingerstick lancet and pipetting from vacutainer tubes of collected venous blood are two common methods for DBS preparation. With regard to DBS prepared by fingerstick, it is minimally invasive, requires a small blood volume, and the samples can be collected by the patients themselves and sent to the assigned labs by regular mail with low cost [8]. With regard to DBS prepared by health professionals by pipetting from blood collected in vacutainer tubes, samples stored as DBS may be useful for many applications. The benefits of dried blood spots are of particular interest for comparative metabolomics studies conducted in remote locations, and those having large cohorts from multiple collaborating centers [9]. In recent years, due to the above benefits, DBS have also been increasingly used in various other fields, such as therapeutic drug and abused substances monitoring [8,10], pharmacokinetics [11], proteomics [12], lipidomics [13,14], and metabolomics [15,16,17,18].
Although metabolomic analysis of DBS samples has been reported in a few studies, there is still a lack of systematic optimization of the specific protocols for DBS preparation and storage for metabolomics. In addition, previous DBS-based metabolomics methods often require multiple injections and platforms to increase the metabolome coverage. These methods have only measured a small number of metabolites [13,16]. A single liquid chromatography-tandem mass spectrometry (LC-MS/MS) platform targeting a large number of metabolites with a broad range of chemical classes is highly desirable to fulfill different application requirements, reduce the cost, and increase the productivity of DBS-based metabolomics. Moreover, the stability of the whole blood metabolome in DBS under different storage conditions remains largely unknown except for some routinely analyzed metabolites in newborn screening, such as acylcarnitines and amino acids [4].
MELAS (mitochondrial encephalomyopathy, lactic acidosis and stroke-like episodes) is a mitochondrial disorder with symptoms that can also include headache, fatigue, diabetes, and hearing loss [19]. Nearly 80% of MELAS patients have an adenine-to-guanine transition mutation at position 3243 of their mitochondrial DNA (A3243G) located in the mitochondrial tRNA-leucine (UUR) gene [20]. The metabolomic analysis of disease-specific metabolic fingerprints in the DBS of MELAS patients might provide useful information for the diagnosis and potential treatment targets of MELAS.
In this study, we developed a DBS-based broad-spectrum, targeted metabolomic method and performed the systematic optimization of critical steps in the workflow. The effect of variations in blood spot volumes, extraction methods, carryover, the comparability with other types of blood samples, and the potential influence of hematocrit variations were analyzed. Additionally, the stability of over 400 metabolites stored as DBS samples under three different storage conditions was analyzed periodically for up to one year. The analytical reproducibility and potential clinical utility were demonstrated with DBS samples prepared from 12 pediatric patients with MELAS and 12 age- and sex-matched controls.

2. Materials and Methods

2.1. Informed Consent

The study was conducted with signed informed consent under the approved (2/24/2014) University of California, San Diego (UCSD) institutional review board (IRB) protocol #140072. The elements and design of this study are illustrated in Figure 1A–D.

2.2. DBS and DPS (Dried Plasma Spot) Preparation

Whole venous blood was collected using Li-heparin vacutainer tubes (Becton-Dickinson #367884, Franklin Lakes, NJ, USA) or ethylenediaminetetraacetic acid (EDTA) vacutainer tubes (Beckton-Dickinson #367861) mixed by inversion 10 times. Plasma was obtained by centrifugation of the whole blood at 900× g for 10 min at room temperature (RT). The DBS samples were prepared by pipetting 70 µL of Li-heparin whole blood (Li-hep DBS) or EDTA whole blood (EDTA DBS) on to Whatman 903TM paper cards (GE Healthcare, #10534612, Chicago, IL, USA). Dried plasma spot (DPS) samples were prepared by spotting 40 µL of Li-hep plasma onto Whatman 903TM paper cards. Additionally, blood via fingerstick was taken using a sterile lancet and directly spotted onto the circles of a paper card (Fingerstick DBS). One to two drops was sufficient to fill the 12-mm inner diameter of the dotted circle when applied in the center. The DBS cards were air-dried at room temperature (RT) for 24 h. The disks were punched out from the DBS using a ¼ inch diameter (6.3 mm) one-hole puncher (ACCO, Catalog: 74005). Punching from an empty area of the sampling card (Blank disks) was performed in between every sample.

2.3. Method Optimization

Blood spotting volumes on Whatman 903TM paper cards were optimized. Five blood spotting volumes (40, 50, 60, 70, and 80 µL) were tested. The volumes of extraction buffer (100, 200, 300, 400, 500 µL) and the extraction methods (20 min on ice, overnight at RT, and sonication for 5 min on ice) were also optimized. The potential cross-contamination was assessed by comparing the level of a metabolite in a DBS sample with the levels in the following 2 blank paper punches. Each condition was measured in triplicate.

2.4. Sample Type Comparisons

Samples were collected from 3 subjects with duplicates for each condition (Figure 1B). Sample comparisons included Li-hep DBS versus liquid whole blood from lithium-heparin venous blood, DBS versus DPS, DBS obtained from Li-heparin versus EDTA vacutainer tubes, and fingerstick DBS versus Li-heparin DBS.

2.5. Stability Analysis

DBS samples from 3 individuals were collected at the baseline (after drying for 24 h), 1, 2, 3, 7, 9, and 12 months (Figure 1C). The following storage conditions were tested: (1) Room temperature (RT) with desiccant and an oxygen scavenger (RT + desic + O2 scavenger), (2) −80 °C with desiccant (−80 °C + desic); and (3) −80 °C with desiccant and an oxygen scavenger (−80 °C + desic + O2 scavenger). In each case, the DBS was placed in a light-tight, metalized Mylar bag (PackFreshUSA, Cat. #MGB50, Ontario, CA, USA) along with two desiccant packs (1.5 g/pack), with or without one bag of oxygen scavenger (300 cc, PackFreshUSA, Cat. #F3C1C, ferrous iron) and was then heat sealed.

2.6. Hematocrit Effects

Plasma and whole blood samples from the same individual were mixed to obtain hematocrit levels ranging from 0% to 50%. A total of 70 µL of the mixture was spotted on Whatman 903TM paper cards and dried for 48 h before analysis. Three punches per hematocrit level were collected for further analysis.

2.7. Metabolite Extraction

The 6.3 mm disk punch was transferred to a 2 mL tube and 200 µL of extraction buffer consisting of prechilled MEOH-ACN-H2O (40:40:20, v/v/v), and custom-synthesized stable isotope internal standards were added. The list of internal standards has been previously described [21]. After incubation on ice for 20 min, it was centrifuged at 16,000 g for 10 min at 4 °C. The supernatant was transferred to a new tube and stored at −80 °C until analysis (Figure 1E).

2.8. Metabolomic Analysis

Samples were analyzed on a Shimadzu LC-20AD high-performance liquid chromatography (HPLC) system coupled with a SCIEX Qtrap 5500 mass spectrometer (LC-ESI-MS/MS) using a broad-spectrum targeted metabolomic method as previously described [21]. Briefly, metabolites were separated on a polymer-based NH2 column (250 × 2 mm, 4 µm, Asahipak NH2 P-40 2E, Showa Denko America, Inc., New York, NY, USA) in hydrophilic interaction liquid chromatography (HILIC) mode. A total of 622 endogenous metabolites and 78 stable isotope internal standards were targeted in a single injection [21]. The peaks were manually checked in Multiquant 3.0 (Sciex) and the peak area was exported for further analysis. The spiked stable isotope internal standards were used to monitor the metabolite extraction and MS detection.

2.9. Clinical Application in MELAS

Whole venous blood was collected from 12 MELAS patients, ranging in age from 0.3 to 11.6 years, and 12 age and sex-matched controls using Li-Heparin vacutainer tubes. The DBS samples were prepared by pipetting 70 µL of whole blood on to Whatman 903TM paper cards (Figure 1D). The DBS samples were stored at RT with desiccant and were analyzed within a week after DBS preparation using the metabolomic method described above. The patients’ characteristics are listed in Table S1. The reproducibility of the assay was evaluated using 6 replicates of the pooled MELAS and control DBS extract analyzed on each of the 5 days [22]. The mean intra-assay, inter-assay, and total relative standard deviations (RSDs) for each pool were calculated [22].

2.10. Bioinformatic and Statistical Analysis

The data were expressed as mean ± SD unless otherwise indicated in the figure legends. A Student’s t-test or one-way ANOVA followed by Tukey’s test was used to compare the group differences. p < 0.05 was considered as statistically significant. The correlations were analyzed by Spearman’s rank correlation. Frequency distribution analysis was applied to show the range of RSDs for all detected metabolites.
Metabolomic analysis was performed in MetaboAnalyst 4.0 (https://www.metaboanalyst.ca/). Briefly, the mass chromatographic areas under the curve (AUCs) were log2-transformed, scaled by the SD of each metabolite, and analyzed by partial least squares discriminant analysis (PLS-DA). Metabolites with the highest separation power in PLS-DA were ranked by variable importance in projection (VIP). VIP scores > 1.5 were considered significant. Sets of 3–6 metabolites were selected manually from top 50 significant metabolites as candidate diagnostic classifiers. The diagnostic performance of the selected classifiers was then visualized and quantified by the receiver operator characteristic (ROC) curve analysis with a random forest (RF) algorithm. Classifier robustness was estimated by permutation test (permutation times = 1000). Sensitivity, specificity, and accuracy were estimated by 2 × 2 contingency table analysis in Prism 8.0 (GraphPad Software, La Jolla, CA, USA, www.graphpad.com).

3. Results

3.1. Method Optimization

We first optimized the blood spot volumes on GE Whatman 903TM protein saver cards. A subset of metabolites (47 compounds), selected to cover the entire LC gradient and all the chemical classes we targeted, was analyzed in more detail. Figure S1 shows the correlation between the spot volumes and the diameters of DBSs. The peak area for 47 representative metabolites increased from 40 to 60 µL of spotted volume. (Figure 2A and Table S2). We found that blood volumes of 60–80 µL spotted on Whatman 903TM cards resulted in an even spreading of the sample and equivalent yields in the punch. For convenience, we used 70 µL for all subsequent experiments. This produced a blood spot diameter of 13.5 ± 0.55 mm. It was calculated that a 6.3 mm diameter hole punch of the dried blood spot contained the dried residue from 13.4 µL of blood (Figure S1B).
We used methanol: acetonitrile: water (MEOH-ACN-H2O; 40:40:20) as the buffer for metabolite extraction [21] and optimized the extraction volume for the DBS discs. The yield remained similar with ≥ 200 µL of extraction buffer (Figure 2B). Three different extraction protocols were tested including incubation on ice for 20 min, incubation at 4 °C overnight, and sonication in an ice bath for 5 min. No significant differences between the three extraction methods were identified (Figure 2C).
Cross-contamination during the processing of DBS for metabolomic analysis was minimal. The average peak intensity of the first blank disc following DBS sample was 217-fold lower than DBS sample, which was similar to the second punch of the blank disc and the solvent blank (Figure 2D).

3.2. Sample Type Comparisons

We next analyzed the similarities and differences between different sample types. A comparison of Li-heparin DBS to Li-heparin liquid whole blood revealed no significant difference in the total number of metabolites detected (S/N > 3) (431 ± 4 vs 433 ± 2, p > 0.05). Four hundred and twenty-four (424) detected metabolites overlapped, nine metabolites were only detectable in DBS, and 11 metabolites were unique in the liquid whole blood (Table S3). The AUCs correlated well between DBS and liquid whole blood samples (Spearman r = 0.874, p < 0.01, Figure 3A). The similarities and differences between DBS and DPS are described in the supplemental results (Figure S2 and Table S4).
The anticoagulant used in the blood collection tube affected the resulting metabolome of DBS. Overall, the metabolites from Li-heparin DBS were highly correlated with those from EDTA DBS (Spearman r = 0.981, p < 0.01, Figure 3B). However, there were also some notable differences: EDTA DBS yielded higher AUCs for nucleotides and organic acids, while, adenosine, glutamine, L-lysine, and p-hydroxyphenylacetic acid were dramatically lower (p < 0.05) (Figure 3C and Table S5).
Li-heparin DBS also correlated well with fingerstick DBS (Spearman r = 0.985, p < 0.01, Figure 3D). Samples collected as lithium heparin venous blood and fingerstick capillary blood samples collected by lancet and stored as DBS were indistinguishable by principal component analysis (PCA). Volcano plot analysis showed that out of the 424 detected metabolites in our study, only 10 metabolites were significantly different between Li-heparin DBS and fingerstick DBS (Table S6).

3.3. One-Year Stability Analysis

We first examined the overall changes in the geometric means of the AUCs for all detectable metabolites over a 1-year period. Based on the analytical variances of our method for DBS (Figure S3), metabolites were considered as stable when the changes were ≤ 15%. Regardless of storage condition, the geometric mean was stable for the first month, rapidly declined to 58%–60% of the baseline by 3 months, and then remained unchanged for the subsequent 9 months (Figure 4A–C). A more detailed analysis, however, revealed differences in individual metabolite stability between the 3 storage conditions (Figure 4D). The best performing condition overall was −80 °C with desiccant and without an O2 scavenger (Figure 4D–G).
Surprisingly, the presence of an O2 scavenger showed a negative impact on the stability of metabolites stored at −80 °C. The inclusion of an O2 scavenger resulted in fewer unchanged metabolites (50 vs 95) (Table S7) and more metabolites with a > 50% decrease (136 vs 103) (Figure 4D and Table S8). The presence of an O2 scavenger also caused a subset of metabolites, primarily hydroxyeicosatetraenoic acids (HETEs) and eicosanoids, to have increased AUCs (Figure S4A,B and Table S9). This increase was noted at 1 month and persisted throughout the year of analysis. This increase was most notable at room temperature, but also occurred in the samples stored at −80 °C with an O2 scavenger (66 and 26 metabolites, respectively).
Phospholipids, sphingolipids, acylcarnitines, amino acids, and steroids were the main chemical classes with ≥ 30% reduction (Figure 4E). A few examples from each chemical class demonstrating the changes in metabolite abundance within a year at RT with desiccant and an O2 scavenger are shown in Figure S4C–L. Interestingly, we noticed that both the number of decreased phospholipids and sphingolipids and the percentage of decline from the baseline were similar in all three storage conditions (Figure 4E,F). In terms of phospholipid subclasses, phosphatidylserine (PS), phosphatidylinositol (PI), phosphatidylethanolamine (PE), Bis(monoacylglycero)phosphate (BMP) and phosphatidylglycerol (PG) had a more significant reduction than cardiolipins (CL), phosphatidic acid (PA), and phosphatidylcholine (PC) (Figure S5A). The saturation of the phospholipid acyl chain had no effect on the stability of phospholipids on DBS cards, and the levels of reduction for sphingolipids were not correlated with acyl chain length (Figure S5B,C).
In contrast, the decline of acylcarnitine and amino acid AUCs from the baseline was decreased at −80 °C compared to room temperature storage (Figure 4E). Out of the 31 acylcarnitines measured, 26 were reduced by 54.5% ± 16.2% at room temperature, but only three acylcarnitines were reduced at −80 °C with desiccant, with an average reduction of 36.9% ± 6.8% (Figure 4F). Similarly, the number of amino acids with ≥ 30% reduction dropped from 24 (out of 55) when the cards were stored at room temperature to 10 for the DBS samples stored at −80 °C with desiccant.
Out of 404 detected metabolites, the number of metabolites that was not significantly altered after one year of storage ranged from 29 metabolites at RT with desiccant and an O2 scavenger to 50 at −80 °C with desiccant and an O2 scavenger, and 95 at −80 °C with desiccant only (Figure 4D and Table S7). A Venn diagram of the overlapping metabolites between the 3 conditions is shown in Figure 4G. Interestingly, 10 metabolites remained stable at all three conditions. These included creatine, hydroxyproline, leucine, acetoacetic acid, and O-acetylserine (Figure 4G and Table S7).

3.4. Clinical Application in MELAS

The average intra-day and inter-day RSDs were 9.7% and 10.9% for controls, and 9.1% and 13.9% for the pooled MELAS samples (Figure S3A–D), well within the limit of 20% for RSDs listed in the Food and Drug Administration (FDA)’s guideline for bio-analytical method validation using LC-MS [23].
The raw DBS metabolomic data for MELAS samples and controls are in Table S10. PLS-DA revealed a clear separation of the blood metabolome of MELAS patients from the controls (Figure 5A). VIP analysis listed the top 20 metabolites with the highest separation power in PLS-DA (Figure 5B). Fatty acid oxidation and synthesis, plasmalogens, phospholipids, and sphingolipids were the major metabolic pathways disturbed in MELAS patients (Figure 5C). Figure 5D shows the metabolites that were significantly correlated with heteroplasmy.
The highest area under the ROC curve was 0.91 (95% CI: 0.72–1) using three metabolites including 3-hydroxydodecanoylcarnitine, plasmalogen (p18:0/22:6), and glutathione (Figure 5E). The sensitivity was 84.6% (95% CI: 54.6–98.1%) and the specificity was 90.9% (95% CI: 58.7–99.8%) (Table S11). Similar diagnostic accuracy, sensitivity, and specificity were obtained in the liquid whole blood from lithium-heparin samples using the same set of metabolites (Figure 5F and Table S11). Univariate statistical analysis showed that 3-hydroxydodecanoyl carnitine and plasmalogen (p18:0/22:6) were significantly higher in MELAS samples than the controls (p = 0.001 for both metabolites) (Figure 5G–I). In contrast, reduced glutathione (GSH) was dramatically lower in MELAS compared to the controls (p = 0.008).

4. Discussion

The objective of this study was the development and optimization of a dried blood spot-based, broad-spectrum targeted metabolomic method. We performed the first systematic optimization of DBS preparation, metabolite extraction and detection, and one-year storage stability in different conditions for over 400 metabolites in DBS. In addition, we examined the similarities and differences of the metabolome from DBS as compared to other types of blood samples.
In DBS preparation, blood volume and hematocrit variations are common challenges. Several sophisticated approaches have been proposed for the estimation and correction of blood volume variations [24,25]. Here, we found that the differences in blood spot volume did not affect our metabolomic analysis if the spot volume was over 60 µL on FDA-approved Whatman 903TM protein saver cards. The mass of metabolites in the 6.3 mm punch remained constant as long as at least 90% of the 12 mm dotted circle on the card was filled. This was true if 60 µL, 70 µL, or 80 µL of blood was applied. The hematocrit effect is another potential obstacle to DBS application [26]. Our results showed that hematocrits between 31% and 50% did not play a significant role as a source of variation in broad-spectrum metabolomics results (Figure S6 and Table S12). Previous studies have similarly shown no apparent impact of hematocrit on the quantification of 25-hydroxyvitamin D3 and cyclosporine A with hematocrits ranging from 20% to 60% [27,28].
Drolet et al. reported the detection of 300 polar metabolites from DBS using the combination of HILIC-MS/MS and GC-MS [13]. Our method routinely detected about 430 metabolites from one DBS sample using a single injection. Moreover, 11 chemical classes were covered in our assay, such as lipids and lipid-like molecules, amino acids and derivatives, organic acids and derivatives, nucleotides and derivatives, and acylcarnitines (Figure S7). To our knowledge, this is the first method that measures such a large number of metabolites with a broad range of chemical classes from DBS using a targeted approach in a single injection.
Our results highlight the importance of using a consistent sample type and collection method for metabolomic studies. DBS samples prepared by Li-heparin or EDTA collection tubes have similar levels of precision and reproducibility compared to liquid whole blood, but show clear differences related to the specific class of metabolites detected. Our lab has previously reported similar metabolomic differences between Li-heparin plasma and EDTA plasma samples [21]. The most striking differences were the significant decrease of phosphorylated nucleotides and an increase in nucleosides such as adenosine in Li-heparin samples compared with EDTA samples.
Regarding differences between capillary and venous blood, little has been described beyond some basic biochemical parameters, such as complete blood counts and hemoglobin [29]. Our analysis showed that samples collected by fingerstick lancet or as Li-heparin venous blood and stored as dried blood spots were superimposable by PCA (Figure 3E). This was true despite the comparison between a precisely measured 70 µL of heparinized whole blood spotted with a pipetman, and fingerstick lancet samples that contained variable volumes of blood. A small number (10 of 424) of the metabolites were found to be different in samples collected by fingerstick and heparinized venous blood (Table S6), but 97.6% (414 of 424) of the measured metabolites were indistinguishable by these two sample collection methods. The choice of any specific sampling method for dried blood spot preparation, e.g., fingerstick lancet, EDTA-whole blood, or Li-heparin whole blood, should be standardized at the beginning and not mixed within any given metabolomic study.
It is critical to evaluate the stability of metabolites in DBS to enable the optimization of processes and the proper interpretation of data. Our results demonstrate that most metabolites were stable for the first month of storage, regardless of condition. Similar results have been reported for the stability of amino acids, acylcarnitines, and vitamin D in DBS stored for one month at room temperature with desiccants [30,31]. Besides these routinely analyzed metabolites, in our study, we expanded the stability analysis to other key metabolites in DBS, and for up to one year of storage. We showed that the decline primarily occurred between 1 and 3 months of storage, then did not decrease further for at least 1 year. Additionally, metabolites in DBS were more stable at −80 °C than those stored at room temperature even though the decline cannot be completely prevented. Ten metabolites were found to be unchanged for the entire year, regardless of the storage condition. Further studies are needed to determine if any of these analytes can be used to normalize for storage time. Our results highlight the need for case and control samples used in comparative metabolomic studies to be matched for sample collection method, storage time, and storage conditions. This study focused on the relative changes in the AUCs of the targeted metabolites compared to the baseline sample and did not attempt to determine the absolute concentration of any single metabolite.
In addition to the decrease in the stability of many metabolites in the presence of an O2 scavenger, the finding of an apparent increase in abundance of some metabolites was unexpected. This increase occurred rapidly in the first month, slowly approached a plateau after 3 months of storage, then stabilized (Figure S3A,B). The affected metabolites included arachidonic acid, and several of its metabolites, such as 5-HETE and 5, 6-DHET (Table S9). Since the number of metabolites that increased in the presence of more hypoxic environment was less in the cold (−80 °C) than at room temperature, these results may reflect the residual activity of the enzymes needed to synthesize these molecules. For example, lower oxygen tensions may have resulted in improved preservation of several blood proteins, including enzymes like phospholipase A2 (PLA2) needed to release arachidonate from phospholipids, arachidonate 5-lipoxygenase (ALOX), and other lipoxygenases that convert arachidonate to 5-HETE and related oxylipins. After 3 months of storage, no further changes were noted (Figure S3).
In this study, we used our dried blood spot method to examine the metabolome in 12 MELAS patients and 12 matched controls. We identified 3-hydroxydodecanoylcarnitine (C12-OH), plasmalogen (p18:0/22:6), and glutathione (GSH) in DBS metabolome as the potential biomarkers for MELAS. Acylcarnitines are involved in the mitochondrial metabolism of fatty acids. Defects of fatty acid oxidation and the carnitine shuttle system had been observed in inborn errors of metabolism with myopathies including MELAS [32]. Similar to our findings, the circulating hydroxylated long-chain acylcarnitines were shown to be significantly elevated in the plasma of patients with mitochondrial myopathy [33]. Plasmalogens are a subtype of phospholipids characterized by the presence of a vinyl–ether bond at the sn-1 position of the phosphoglycerol moiety. The initial steps of plasmalogen biosynthesis are confined to peroxisomes. Due to their unique structure and composition, plasmalogens are enriched in highly aerobic tissues dependent on mitochondria for energy like brain, heart, muscle, and nerve. Plasmalogens play many crucial roles in cellular functions, including reservoirs for second messengers, working as endogenous antioxidants, as well as the structural components of cellular membranes [34,35]. The aberrant metabolism of plasmalogens had been found in Alzheimer’s disease [36], Parkinson’s disease [37], Niemann-Pick disease [38], and cancers [39]. The role of plasmalogens in mitochondrial disorders remains largely unknown. Mitochondrial disorders are associated with redox imbalance, and GSH is the low molecular weight thiol involved in mitochondrial redox signaling. Decreased levels of peripheral whole blood GSH had been previously reported in mitochondrial disease patients, including Leigh syndrome, electron transport chain abnormalities, MELAS, and mtDNA depletion syndrome [40].
Our study had some limitations. First, we only performed the stability analysis for up to one year under three different storage conditions. The stability beyond one year of storage remains unknown. Second, this study did not include a room temperature condition without desiccant or an oxygen scavenger.

5. Conclusions

We describe a reproducible dried blood spot-based, targeted metabolomic workflow that is able to reliably detect about 430 metabolites from a diverse array of chemical classes in a single injection and LC-MS/MS analysis. The DBS metabolome had excellent correlations with other blood sample types. The blood volume and hematocrit variations had no significant impact on the assay for a spot volume > 60 µL and hematocrit > 30%. We found that DBS samples were most stable when stored in −80 °C with desiccants and without an O2 scavenger. The median AUC for the extracted metabolites stabilized at about 60% of the baseline value after 3 months and did not decline further for at least 1 year of storage. Additionally, all metabolites measured at the baseline remained quantifiable for at least 1 year. The proof-of-concept study with DBS samples from MELAS patients and contemporaneously collected disease controls demonstrates the practical clinical value of the described methods. The case and control samples used in comparative metabolomic studies need to be matched for the sample collection method, storage time, and storage conditions.

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-1989/10/3/82/s1, Figure S1: The correlation between the whole blood spot volumes and the diameters of the dried blood spots (DBSs). Figure S2: The similarities and differences between DBS and DPS. Figure S3: Frequency distribution analysis of intra-day and inter-day reproducibility for all detected metabolites in DBS of control and MELAS patients. Figure S4: The changes of the representative metabolites on DBS stored at room temperature with an O2 scavenger and desiccant for a year. Figure S5: The effect of the head group, acyl chain length, and saturation on the stability of phospholipids and sphingolipids on DBS. Figure S6: The impact of the varying hematocrit values on the metabolomic profiling of DBS. Figure S7: Chemical classes interrogated by DBS Metabolomics. Table S1: The subjects’ characteristics for MELAS patients and controls. Table S2: The effect of blood spot volumes on the AUCs of 47 representative metabolites. Table S3: Metabolites only detected in either Li-heparin DBS or Li-heparin liquid whole blood. Table S4: The full list of differential metabolites between DBS and DPS identified by volcano plot analysis. Table S5: The full list of differential metabolites between Li-hep DBS and EDTA DBS identified by volcano plot analysis. Table S6: The full list of differential metabolites between Li-heparin DBS and fingerstick DBS identified by volcano plot analysis. Table S7: The list of stable metabolites on DBS after one year of storage at three different conditions. Table S8: The list of metabolites with decreased AUCs after one-year of storage. Table S9: The list of metabolites with increased AUCs after one-year of storage. Table S10: The raw DBS metabolomic data for MELAS study. Table S11: Broad-spectrum targeted metabolomic profiling of DBS and liquid whole blood yields a similar diagnostic accuracy for MELAS. Table S12: The list of metabolites that were not significantly associated with hematocrit changes.

Author Contributions

Conceptualization, K.L. and R.K.N.; methodology, K.L., J.C.N., R.K.N., J.M.M. and L.W.; software, J.M.M.; validation, K.L., J.C.N. and R.K.N.; formal analysis, K.L., J.C.N. and R.K.N.; investigation, K.L. and J.C.N.; resources, J.M.M.; data curation, K.L. and J.C.N.; writing—original draft preparation, K.L.; J.C.N.; J.M.M.; L.W. and R.K.N.; writing—review and editing, K.L., J.C.N., J.M.M., L.W. and R.K.N.; visualization, K.L., and J.M.M.; supervision, R.K.N.; project administration, R.K.N.; funding acquisition, R.K.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the UCSD Christini Fund, The Lennox Foundation, and the UCSD Mitochondrial Disease Research Fund. Funding for the mass spectrometers used in this study was provided by a gift from the Jane Botsford-Johnson Foundation.

Acknowledgments

The authors thank the patients and families who donated their time and effort in helping to make this study possible.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Naviaux, R.K.; Naviaux, J.C.; Li, K.; Bright, A.T.; Alaynick, W.A.; Wang, L.; Baxter, A.; Nathan, N.; Anderson, W.; Gordon, E. Metabolic features of chronic fatigue syndrome. Proc. Natl. Acad. Sci. USA 2016, 113, E5472–E5480. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Naviaux, R.K.; Naviaux, J.C.; Li, K.; Wang, L.; Monk, J.M.; Bright, A.T.; Koslik, H.J.; Ritchie, J.B.; Golomb, B.A. Metabolic features of Gulf War illness. PLoS ONE 2019, 14, e0219531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Johnson, C.H.; Ivanisevic, J.; Siuzdak, G. Metabolomics: Beyond biomarkers and towards mechanisms. Nat. Rev. Mol. Cell Biol. 2016, 17, 451–459. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Kirwan, J.A.; Brennan, L.; Broadhurst, D.; Fiehn, O.; Cascante, M.; Dunn, W.B.; Schmidt, M.A.; Velagapudi, V. Preanalytical Processing and Biobanking Procedures of Biological Samples for Metabolomics Research: A White Paper, Community Perspective (for “Precision Medicine and Pharmacometabolomics Task Group”—The Metabolomics Society Initiative). Clin. Chem. 2018, 64, 1158–1182. [Google Scholar] [CrossRef] [Green Version]
  5. Guthrie, R.; Susi, A. A Simple Phenylalanine Method for Detecting Phenylketonuria in Large Populations of Newborn Infants. Pediatrics 1963, 32, 338–343. [Google Scholar]
  6. Zytkovicz, T.H.; Fitzgerald, E.F.; Marsden, D.; Larson, C.A.; Shih, V.E.; Johnson, D.M.; Strauss, A.W.; Comeau, A.M.; Eaton, R.B.; Grady, G.F. Tandem mass spectrometric analysis for amino, organic, and fatty acid disorders in newborn dried blood spots: A two-year summary from the New England Newborn Screening Program. Clin. Chem. 2001, 47, 1945–1955. [Google Scholar] [CrossRef]
  7. Kvaskoff, D.; Ko, P.; Simila, H.A.; Eyles, D.W. Distribution of 25-hydroxyvitamin D3 in dried blood spots and implications for its quantitation by tandem mass spectrometry. J. Chromatogr. B 2012, 901, 47–52. [Google Scholar] [CrossRef]
  8. Martial, L.C.; Aarnoutse, R.E.; Schreuder, M.F.; Henriet, S.S.; Bruggemann, R.J.; Joore, M.A. Cost Evaluation of Dried Blood Spot Home Sampling as Compared to Conventional Sampling for Therapeutic Drug Monitoring in Children. PLoS ONE 2016, 11, e0167433. [Google Scholar] [CrossRef]
  9. Albani, V.; Celis-Morales, C.; Marsaux, C.F.; Forster, H.; O’Donovan, C.B.; Woolhead, C.; Macready, A.L.; Fallaize, R.; Navas-Carretero, S.; San-Cristobal, R.; et al. Exploring the association of dairy product intake with the fatty acids C15:0 and C17:0 measured from dried blood spots in a multipopulation cohort: Findings from the Food4Me study. Mol. Nutr. Food Res. 2016, 60, 834–845. [Google Scholar] [CrossRef] [Green Version]
  10. Mercolini, L.; Mandrioli, R.; Sorella, V.; Somaini, L.; Giocondi, D.; Serpelloni, G.; Raggi, M.A. Dried blood spots: Liquid chromatography-mass spectrometry analysis of Delta9-tetrahydrocannabinol and its main metabolites. J. Chromatogr. A 2013, 1271, 33–40. [Google Scholar] [CrossRef]
  11. Martial, L.C.; Kerkhoff, J.; Martinez, N.; Rodriguez, M.; Coronel, R.; Molinas, G.; Roman, M.; Gomez, R.; Aguirre, S.; Jongedijk, E.; et al. Evaluation of dried blood spot sampling for pharmacokinetic research and therapeutic drug monitoring of anti-tuberculosis drugs in children. Int. J. Antimicrob. Agents 2018, 52, 109–113. [Google Scholar] [CrossRef] [PubMed]
  12. Bjorkesten, J.; Enroth, S.; Shen, Q.; Wik, L.; Hougaard, D.M.; Cohen, A.S.; Sorensen, L.; Giedraitis, V.; Ingelsson, M.; Larsson, A.; et al. Stability of Proteins in Dried Blood Spot Biobanks. Mol. Cell. Proteom. 2017, 16, 1286–1296. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Gao, F.; McDaniel, J.; Chen, E.Y.; Rockwell, H.E.; Drolet, J.; Vishnudas, V.K.; Tolstikov, V.; Sarangarajan, R.; Narain, N.R.; Kiebish, M.A. Dynamic and temporal assessment of human dried blood spot MS/MS(ALL) shotgun lipidomics analysis. Nutr. Metab. Lond 2017, 14, 28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Koulman, A.; Prentice, P.; Wong, M.C.Y.; Matthews, L.; Bond, N.J.; Eiden, M.; Griffin, J.L.; Dunger, D.B. The development and validation of a fast and robust dried blood spot based lipid profiling method to study infant metabolism. Metabolomics 2014, 10, 1018–1025. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Haijes, H.A.; Willemsen, M.; Van der Ham, M.; Gerrits, J.; Pras-Raves, M.L.; Prinsen, H.; Van Hasselt, P.M.; De Sain-van der Velden, M.G.M.; Verhoeven-Duif, N.M.; Jans, J.J.M. Direct Infusion Based Metabolomics Identifies Metabolic Disease in Patients’ Dried Blood Spots and Plasma. Metabolites 2019, 9, 12. [Google Scholar] [CrossRef] [Green Version]
  16. Petrick, L.; Edmands, W.; Schiffman, C.; Grigoryan, H.; Perttula, K.; Yano, Y.; Dudoit, S.; Whitehead, T.; Metayer, C.; Rappaport, S. An untargeted metabolomics method for archived newborn dried blood spots in epidemiologic studies. Metabolomics 2017, 13, 27. [Google Scholar] [CrossRef] [Green Version]
  17. Wang, Q.; Sun, T.; Cao, Y.; Gao, P.; Dong, J.; Fang, Y.; Fang, Z.; Sun, X.; Zhu, Z. A dried blood spot mass spectrometry metabolomic approach for rapid breast cancer detection. Onco Targets Ther. 2016, 9, 1389–1398. [Google Scholar] [CrossRef] [Green Version]
  18. Jing, Y.; Wu, X.; Gao, P.; Fang, Z.Z.; Wu, J.J.; Wang, Q.J.; Li, C.; Zhu, Z.T.; Cao, Y.F. Rapid Differentiating Colorectal Cancer and Colorectal Polyp Using Dried Blood Spot Mass Spectrometry Metabolomic Approach. IUBMB Life 2017, 69, 347–354. [Google Scholar] [CrossRef] [Green Version]
  19. Gorman, G.S.; Chinnery, P.F.; DiMauro, S.; Hirano, M.; Koga, Y.; McFarland, R.; Suomalainen, A.; Thorburn, D.R.; Zeviani, M.; Turnbull, D.M. Mitochondrial diseases. Nat. Rev. Dis. Primers 2016, 2, 16080. [Google Scholar] [CrossRef]
  20. Koga, Y.; Povalko, N.; Nishioka, J.; Katayama, K.; Yatsuga, S.; Matsuishi, T. Molecular pathology of MELAS and L-arginine effects. Biochim. Biophys. Acta 2012, 1820, 608–614. [Google Scholar] [CrossRef]
  21. Li, K.; Naviaux, J.C.; Bright, A.T.; Wang, L.; Naviaux, R.K. A robust, single-injection method for targeted, broad-spectrum plasma metabolomics. Metabolomics 2017, 13, 122. [Google Scholar] [CrossRef] [PubMed]
  22. Grant, R.P.; Hoofnagle, A.N. From lost in translation to paradise found: Enabling protein biomarker method transfer by mass spectrometry. Clin. Chem. 2014, 60, 941–944. [Google Scholar] [CrossRef] [PubMed]
  23. Kaza, M.; Karazniewicz-Lada, M.; Kosicka, K.; Siemiatkowska, A.; Rudzki, P.J. Bioanalytical method validation: New FDA guidance vs. EMA guideline. Better or worse? J. Pharm. Biomed. Anal. 2019, 165, 381–385. [Google Scholar] [CrossRef] [PubMed]
  24. Liao, H.W.; Lin, S.W.; Chen, G.Y.; Kuo, C.H. Estimation and Correction of the Blood Volume Variations of Dried Blood Spots Using a Postcolumn Infused-Internal Standard Strategy with LC-Electrospray Ionization-MS. Anal. Chem. 2016, 88, 6457–6464. [Google Scholar] [CrossRef]
  25. Zheng, N.; Yuan, L.; Ji, Q.C.; Mangus, H.; Song, Y.; Frost, C.; Zeng, J.; Aubry, A.F.; Arnold, M.E. “Center punch” and “whole spot” bioanalysis of apixaban in human dried blood spot samples by UHPLC-MS/MS. J. Chromatogr. B 2015, 988, 66–74. [Google Scholar] [CrossRef] [PubMed]
  26. Capiau, S.; Wilk, L.S.; Aalders, M.C.; Stove, C.P. A Novel, Nondestructive, Dried Blood Spot-Based Hematocrit Prediction Method Using Noncontact Diffuse Reflectance Spectroscopy. Anal. Chem. 2016, 88, 6538–6546. [Google Scholar] [CrossRef] [Green Version]
  27. Newman, M.S.; Brandon, T.R.; Groves, M.N.; Gregory, W.L.; Kapur, S.; Zava, D.T. A liquid chromatography/tandem mass spectrometry method for determination of 25-hydroxy vitamin D2 and 25-hydroxy vitamin D3 in dried blood spots: A potential adjunct to diabetes and cardiometabolic risk screening. J. Diabetes Sci. Technol. 2009, 3, 156–162. [Google Scholar] [CrossRef] [Green Version]
  28. Wilhelm, A.J.; den Burger, J.C.; Vos, R.M.; Chahbouni, A.; Sinjewel, A. Analysis of cyclosporin A in dried blood spots using liquid chromatography tandem mass spectrometry. J. Chromatogr. B 2009, 877, 1595–1598. [Google Scholar] [CrossRef]
  29. Kayiran, S.M.; Ozbek, N.; Turan, M.; Gurakan, B. Significant differences between capillary and venous complete blood counts in the neonatal period. Clin. Lab. Haematol. 2003, 25, 9–16. [Google Scholar] [CrossRef]
  30. Han, J.; Higgins, R.; Lim, M.D.; Lin, K.; Yang, J.; Borchers, C.H. Short-Term Stabilities of 21 Amino Acids in Dried Blood Spots. Clin. Chem. 2018, 64, 400–402. [Google Scholar] [CrossRef] [Green Version]
  31. Rola, R.; Kowalski, K.; Bienkowski, T.; Kolodynska-Goworek, A.; Studzinska, S. Development of a method for multiple vitamin D metabolite measurements by liquid chromatography coupled with tandem mass spectrometry in dried blood spots. Analyst 2018, 144, 299–309. [Google Scholar] [CrossRef] [PubMed]
  32. El-Gharbawy, A.; Vockley, J. Inborn Errors of Metabolism with Myopathy: Defects of Fatty Acid Oxidation and the Carnitine Shuttle System. Pediatr. Clin. N. Am. 2018, 65, 317–335. [Google Scholar] [CrossRef] [PubMed]
  33. Vissing, C.R.; Duno, M.; Wibrand, F.; Christensen, M.; Vissing, J. Hydroxylated Long-Chain Acylcarnitines are Biomarkers of Mitochondrial Myopathy. J. Clin. Endocrinol. Metab. 2019, 104, 5968–5976. [Google Scholar] [CrossRef] [PubMed]
  34. Braverman, N.E.; Moser, A.B. Functions of plasmalogen lipids in health and disease. Biochim. Biophys. Acta 2012, 1822, 1442–1452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Nagan, N.; Zoeller, R.A. Plasmalogens: Biosynthesis and functions. Prog. Lipid Res. 2001, 40, 199–229. [Google Scholar] [CrossRef]
  36. Su, X.Q.; Wang, J.; Sinclair, A.J. Plasmalogens and Alzheimer’s disease: A review. Lipids Health Dis. 2019, 18, 100. [Google Scholar] [CrossRef] [Green Version]
  37. Calvano, C.D.; Ventura, G.; Sardanelli, A.M.M.; Savino, L.; Losito, I.; Michele, G.; Palmisano, F.; Cataldi, T.R.I. Searching for Potential Lipid Biomarkers of Parkinson’s Disease in Parkin-Mutant Human Skin Fibroblasts by HILIC-ESI-MS/MS: Preliminary Findings. Int. J. Mol. Sci. 2019, 20, 3341. [Google Scholar] [CrossRef] [Green Version]
  38. Schedin, S.; Sindelar, P.J.; Pentchev, P.; Brunk, U.; Dallner, G. Peroxisomal impairment in Niemann-Pick type C disease. J. Biol. Chem. 1997, 272, 6245–6251. [Google Scholar] [CrossRef] [Green Version]
  39. Messias, M.C.F.; Mecatti, G.C.; Priolli, D.G.; de Oliveira Carvalho, P. Plasmalogen lipids: Functional mechanism and their involvement in gastrointestinal cancer. Lipids Health Dis. 2018, 17, 41. [Google Scholar] [CrossRef] [Green Version]
  40. Enns, G.M.; Moore, T.; Le, A.; Atkuri, K.; Shah, M.K.; Cusmano-Ozog, K.; Niemi, A.K.; Cowan, T.M. Degree of glutathione deficiency and redox imbalance depend on subtype of mitochondrial disease and clinical status. PLoS ONE 2014, 9, e100001. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Overview of the study design and clinical cohorts. (A) Extraction and carryover. (B) Sample types. (C) Storage stability. (D) MELAS study. (E) DBS extraction and analysis.
Figure 1. Overview of the study design and clinical cohorts. (A) Extraction and carryover. (B) Sample types. (C) Storage stability. (D) MELAS study. (E) DBS extraction and analysis.
Metabolites 10 00082 g001
Figure 2. Dried blood spot (DBS) workflow optimization. (A) The effect of blood spot volumes. (B) The effect of the volumes of the extraction buffer. (C) The different extraction methods. (D) The effect of cross-contamination between punches. A total of 47 representative metabolites covering the whole gradient were selected for the analysis. Data were mean ± SD. One-way ANOVA followed by Tukey’s test was used. **p < 0.01 and ns indicated non-significant.
Figure 2. Dried blood spot (DBS) workflow optimization. (A) The effect of blood spot volumes. (B) The effect of the volumes of the extraction buffer. (C) The different extraction methods. (D) The effect of cross-contamination between punches. A total of 47 representative metabolites covering the whole gradient were selected for the analysis. Data were mean ± SD. One-way ANOVA followed by Tukey’s test was used. **p < 0.01 and ns indicated non-significant.
Metabolites 10 00082 g002
Figure 3. Sample type comparisons. (A) The correlations of areas under the curve (AUCs) of metabolites between wet blood and DBS. All 430 detected metabolites were used and the AUCs were log 2 transformed. The Spearman rank correlation was performed. The other correlation analyses in this figure were similarily performed. (B) The correlations of metabolite AUCs between Li-hep DBS and EDTA DBS. (C) The differentially abundant metabolites between Li-hep DBS and EDTA DBS visualized on the heatmap. The top 25 significant metabolites ranked by the student’s t-test were listed. p < 0.01. (D) The correlations of metabolite AUCs between Li-hep DBS and fingerstick DBS. (E) A Principal component analysis (PCA) of the five different sample collection methods showed that samples collected by fingerstick lancet and heparinized venous blood were superimposable.
Figure 3. Sample type comparisons. (A) The correlations of areas under the curve (AUCs) of metabolites between wet blood and DBS. All 430 detected metabolites were used and the AUCs were log 2 transformed. The Spearman rank correlation was performed. The other correlation analyses in this figure were similarily performed. (B) The correlations of metabolite AUCs between Li-hep DBS and EDTA DBS. (C) The differentially abundant metabolites between Li-hep DBS and EDTA DBS visualized on the heatmap. The top 25 significant metabolites ranked by the student’s t-test were listed. p < 0.01. (D) The correlations of metabolite AUCs between Li-hep DBS and fingerstick DBS. (E) A Principal component analysis (PCA) of the five different sample collection methods showed that samples collected by fingerstick lancet and heparinized venous blood were superimposable.
Metabolites 10 00082 g003
Figure 4. Stability analysis. (AC) The fractional change over a year from the baseline for the geomean of all detected metabolites. (A) DBS stored at room temperature (RT) with desiccant and an O2 scavenger (RT with desiccant and an O2 scavenger). (B) DBS stored in −80 °C with desiccant and an O2 scavenger. (C) DBS stored in −80 °C with desiccant. The shaded area represents the standard deviations (SD). (D) The number of metabolites altered in a year of storage in different conditions. (E) The chemical classes for the metabolites with ≥ 30% decline in one year. The number of metabolites changed was listed in the graph and the total number of targeted metabolites in each chemical class was shown in the legend. (F) The percentage of reduction (%) from the baseline for the metabolites in the major chemical classes. The value in the bar graph was the total number of metabolites in each chemical class. A one-way ANOVA followed by Tukey’s test was conducted. Data indexed by different letters (“a” and “b”) were significantly different from each other (p < 0.05) in each storage condition. (G) The number of stable metabolites over one year of storage in each storage condition. The names of the stable metabolites are listed in Table S7.
Figure 4. Stability analysis. (AC) The fractional change over a year from the baseline for the geomean of all detected metabolites. (A) DBS stored at room temperature (RT) with desiccant and an O2 scavenger (RT with desiccant and an O2 scavenger). (B) DBS stored in −80 °C with desiccant and an O2 scavenger. (C) DBS stored in −80 °C with desiccant. The shaded area represents the standard deviations (SD). (D) The number of metabolites altered in a year of storage in different conditions. (E) The chemical classes for the metabolites with ≥ 30% decline in one year. The number of metabolites changed was listed in the graph and the total number of targeted metabolites in each chemical class was shown in the legend. (F) The percentage of reduction (%) from the baseline for the metabolites in the major chemical classes. The value in the bar graph was the total number of metabolites in each chemical class. A one-way ANOVA followed by Tukey’s test was conducted. Data indexed by different letters (“a” and “b”) were significantly different from each other (p < 0.05) in each storage condition. (G) The number of stable metabolites over one year of storage in each storage condition. The names of the stable metabolites are listed in Table S7.
Metabolites 10 00082 g004
Figure 5. Metabolomic fingerprints and biomarkers in DBS of MELAS patients. (A) Partial least squares discriminant analysis (PLS-DA) revealed the complete separation of MELAS patients from the controls (N = 12/group). (B) variable importance in projection (VIP) plot of the top 20 metabolites with the highest separation power between the controls and MELAS patients in PLS-DA. A VIP score > 1.5 was considered as statistically significant. (C) A pathway analysis of the significantly altered metabolites in the DBS of MELAS patients compared to the controls. The fractional impact of a pathway is defined as the sum of the absolute values of z-scores for the significant metabolites (PLS-DA VIP score > 1.5) divided by the sum of the absolute values of z-scores for all the detected metabolites in that pathway. (D) The top 20 metabolites correlated with heteroplasmy. A Spearman correlation analysis was performed. (E) A receiver operator characteristic (ROC) analysis to evaluate the diagnostic performance of three selected metabolites in DBS for MELAS. (F) The diagnostic performance of the same three metabolites in Li-heparin liquid whole blood for MELAS. Data were log 2 transformed before analysis. (GI) Univariate statistical analysis of three selected biomarkers in the DBS of MELAS and controls. A Student’s t-test was performed after the log 2 transformation of the data.
Figure 5. Metabolomic fingerprints and biomarkers in DBS of MELAS patients. (A) Partial least squares discriminant analysis (PLS-DA) revealed the complete separation of MELAS patients from the controls (N = 12/group). (B) variable importance in projection (VIP) plot of the top 20 metabolites with the highest separation power between the controls and MELAS patients in PLS-DA. A VIP score > 1.5 was considered as statistically significant. (C) A pathway analysis of the significantly altered metabolites in the DBS of MELAS patients compared to the controls. The fractional impact of a pathway is defined as the sum of the absolute values of z-scores for the significant metabolites (PLS-DA VIP score > 1.5) divided by the sum of the absolute values of z-scores for all the detected metabolites in that pathway. (D) The top 20 metabolites correlated with heteroplasmy. A Spearman correlation analysis was performed. (E) A receiver operator characteristic (ROC) analysis to evaluate the diagnostic performance of three selected metabolites in DBS for MELAS. (F) The diagnostic performance of the same three metabolites in Li-heparin liquid whole blood for MELAS. Data were log 2 transformed before analysis. (GI) Univariate statistical analysis of three selected biomarkers in the DBS of MELAS and controls. A Student’s t-test was performed after the log 2 transformation of the data.
Metabolites 10 00082 g005

Share and Cite

MDPI and ACS Style

Li, K.; Naviaux, J.C.; Monk, J.M.; Wang, L.; Naviaux, R.K. Improved Dried Blood Spot-Based Metabolomics: A Targeted, Broad-Spectrum, Single-Injection Method. Metabolites 2020, 10, 82. https://doi.org/10.3390/metabo10030082

AMA Style

Li K, Naviaux JC, Monk JM, Wang L, Naviaux RK. Improved Dried Blood Spot-Based Metabolomics: A Targeted, Broad-Spectrum, Single-Injection Method. Metabolites. 2020; 10(3):82. https://doi.org/10.3390/metabo10030082

Chicago/Turabian Style

Li, Kefeng, Jane C. Naviaux, Jonathan M. Monk, Lin Wang, and Robert K. Naviaux. 2020. "Improved Dried Blood Spot-Based Metabolomics: A Targeted, Broad-Spectrum, Single-Injection Method" Metabolites 10, no. 3: 82. https://doi.org/10.3390/metabo10030082

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop