Next Article in Journal
Efficacy and Safety of Insulin Degludec/Insulin Aspart Compared with a Conventional Premixed Insulin or Basal Insulin: A Meta-Analysis
Previous Article in Journal
Perturbation-Based Modeling Unveils the Autophagic Modulation of Chemosensitivity and Immunogenicity in Breast Cancer Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolite Ratios as Quality Indicators for Pre-Analytical Variation in Serum and EDTA Plasma

1
Institute of Clinical Chemistry and Laboratory Diagnostics and Integrated Biobank Jena (IBBJ), University Hospital Jena, Am Klinikum 1, 07747 Jena, Germany
2
Interdisciplinary Bank of Biological Material and Data Würzburg (IBDW), Straubmühlweg 2a, Haus A9, 97078 Würzburg, Germany
3
Institute of Laboratory Medicine, Clinical Chemistry and Molecular Diagnostics, University Hospital Leipzig, 04103 Leipzig, Germany
4
LIFE Leipzig Research Center for Civilization Diseases, University of Leipzig, 04103 Leipzig, Germany
5
Institute of Medical Statistics, Computer and Data Sciences, Jena University Hospital, Bachstrasse 18, 07743 Jena, Germany
*
Authors to whom correspondence should be addressed.
Metabolites 2021, 11(9), 638; https://doi.org/10.3390/metabo11090638
Submission received: 17 August 2021 / Revised: 13 September 2021 / Accepted: 14 September 2021 / Published: 18 September 2021
(This article belongs to the Section Endocrinology and Clinical Metabolic Research)

Abstract

:
In clinical diagnostics and research, blood samples are one of the most frequently used materials. Nevertheless, exploring the chemical composition of human plasma and serum is challenging due to the highly dynamic influence of pre-analytical variation. A prominent example is the variability in pre-centrifugation delay (time-to-centrifugation; TTC). Quality indicators (QI) reflecting sample TTC are of utmost importance in assessing sample history and resulting sample quality, which is essential for accurate diagnostics and conclusive, reproducible research. In the present study, we subjected human blood to varying TTCs at room temperature prior to processing for plasma or serum preparation. Potential sample QIs were identified by Ultra high pressure liquid chromatography tandem mass spectrometry (UHPLC-MS/MS) based metabolite profiling in samples from healthy volunteers (n = 10). Selected QIs were validated by a targeted MS/MS approach in two independent sets of samples from patients (n = 40 and n = 70). In serum, the hypoxanthine/guanosine (HG) and hypoxanthine/inosine (HI) ratios demonstrated high diagnostic performance (Sensitivity/Specificity > 80%) for the discrimination of samples with a TTC > 1 h. We identified several eicosanoids, such as 12-HETE, 15-(S)-HETE, 8-(S)-HETE, 12-oxo-HETE, (±)13-HODE and 12-(S)-HEPE as QIs for a pre-centrifugation delay > 2 h. 12-HETE, 12-oxo-HETE, 8-(S)-HETE, and 12-(S)-HEPE, and the HI- and HG-ratios could be validated in patient samples.

1. Introduction

Access to high-quality, well characterized human biological samples is crucial for accurate diagnostics in the clinical routine laboratory and for reliable, conclusive biomedical research. In particular, omics technologies rely on high-quality samples obtained by standardized pre-analytical workflows [1].
In practice, sampling procedures and pre-analytical conditions vary substantially in different clinical settings. Similarly, relevant information on these pre-analytical conditions is often missing or insufficient for assessing sample quality. Missing information on pre-analytical conditions may be approximated by retrospective analyses of specific biomarkers known to be affected by distinct pre-analytical conditions. Metabolomics appears particularly suited to identifying potential quality indicators (QI), as metabolites are highly prone to pre-analytical variations. Nevertheless, only a few potential QIs have been described for critical steps of sample processing, such as time-to-centrifugation (TTC) [2,3,4,5,6,7], time-to-freeze [8,9,10,11], or temperature variation [3,4,5]. Examples of well-characterized QIs for TTC are (4E,14Z)-sphingadienine-C-18-1-phosphate (S1P-d18:2) [12], glucose, lactate, arginine, ornithine [3,6] and taurine [7,13]. Particularly, S1P-d18:2, tested in more than 1400 randomly selected serum and plasma samples during single- and multicenter trials in 11 biobanks across three countries, has been proven useful as a QI [12]. QIs specifically reflecting prolonged TTF are lysophosphatidylcholines (lysoPCs) [10] and several proteolytic fragments of endogenous peptides such as fibrinogen alpha and beta chains [9,11]. Further studies identified QIs reflecting time delays between blood collection and centrifugation, variations in time to freeze, and ambient temperature. In a previous study, we have shown that a TTC-dependent increase in taurine can be used in both healthy and diseased individuals to distinguish high qualitative biomaterials from samples undergoing a pre-centrifugation delay of more than one hour [13]. TTC related increase of taurin was temperature dependent as also shown for time-dependent changes of other metabolites, for example, 15-hydroxyeicosatetraenoic acid (15-HETE), glutamate, and cysteine [4]. The largest changes in the metabolome were observed before cell separation, depending on duration and temperature [5,14]. Consequently, the TTC should be minimized and elevated ambient temperature avoided.
While QIs usually respond very clearly to a prolonged TTC or TTF, several studies have shown that their interindividual variability is in some cases greater than the changes observed as a result of altered preanalytics [10,15,16]. Hence, general cutoff thresholds for individual metabolites seem to be insufficient in terms of pre-analytical sample quality assessment. Using ratios of metabolites instead of absolute metabolite concentrations might help to increase the robustness of QIs. Metabolite ratios for ornithine/arginine [6], lactate/glucose [5], and total lysophosphatidylcholines/total phosphatidylcholines [10] have thus been suggested for assessing sample quality.
Here, we used a targeted metabolite profiling approach to discover biomarkers suitable as QIs for the pre-analytical phase, in particular the time-to-centrifugation for both human serum and EDTA plasma. In samples from healthy individuals (n = 10), we could identify novel, and confirm already reported, metabolite ratios that might be well suited as pre-analytical QIs. Among the potential candidates are nucleosides (hypoxanthine/inosine, hypoxanthine/guanosine), amino acids (ornithine/arginine), and eicosanoids. The QI candidates were confirmed by targeted MS/MS analysis in two independent reference sample sets from healthy volunteers (validation sample 1—n = 10, validation sample 2—n = 23) and validated in two independent sample sets from patients with rheumatic and cardiovascular diseases (validation sample 3—n = 40, validation sample 4—serum n = 70; EDTA plasma n = 49).

2. Results

2.1. Targeted Metabolite Profiling of Serum and EDTA Plasma Samples from Healthy Individuals for TTC Quality Indicator Discovery

Targeted metabolite profiling was performed in a pilot study with samples from a set of healthy individuals (discovery sample, n = 10) to compare human serum and EDTA plasma processed after 0.5 and 2 h pre-centrifugation delays. Altogether, 752 serum and 714 plasma metabolites were detected. After data pre-processing (see Methods), 674 serum and 632 plasma metabolites were identified that differed between 0.5 h to 2 h TTC. Effects of a prolonged pre-centrifugation delay were more pronounced in serum than in EDTA plasma. In serum and plasma, 10.4% and 3.9% of all metabolites varied between incubation times, respectively.
In serum, 59 metabolites were increased (seven with a fold change >2), and 11 metabolites were decreased (five with a fold change <0.5). The most remarkable changes were observed for several dipeptides, fibrinogen cleavage peptides, nucleosides, and amino acids (Figure 1A, Data S1). The affected metabolites could be assigned to the following metabolic pathways/subpathways: “glutamate metabolism”, “methionine, cysteine, S-adenosylmethionine (SAM) and taurine metabolism”, “purine metabolism”, “urea cycle”, as well as diverse subpathways of the lipid metabolism (Data S1). In EDTA plasma, 21 metabolites were increased (one with a fold change >2), and four metabolites were decreased (Figure 2A). In EDTA plasma, changes were observed for endocannabinoids and isolated metabolites of various subpathways, such as “fatty acid metabolism” or “urea cycle” (Data S2).
Ratios of metabolites that could be affected by enzymatic or non-enzymatic processes involved in the metabolic pathways of respective QIs were calculated. The relationships were visualized by an over-representation analysis of the metabolites significantly affected by TTC, using KEGG (https://www.metaboanalyst.ca, accessed on 30 June 2020) as a pre-defined reference pathway set. Due to the insufficient coverage of the KEGG pathways by the analytes examined, only 42 metabolites in serum and 13 metabolites in EDTA plasma could be included in the analysis. The pathway impact and p-values were calculated and plotted (Figure 1B and Figure 2B). As a result of concentration changes of hub metabolites (such as asparagine, aspartate, arginine, glycine, glutamate, serine, and taurine) in serum, the amino acid metabolism was enriched. Enrichment in the purine and sphingolipid metabolic pathways was observed as well. Sphingosine, sphinganine, and sphingosine-1-phosphate were assigned to the sphingolipid metabolism and showed increased levels after 2 h TTC. Within the purine metabolism, we observed an increase in xanthine, hypoxanthine, and glycine, and a decrease in adenosine, cyclic AMP, guanosine, and inosine (Data S1). Using this information, we calculated metabolite ratios, such as hypoxanthine/inosine (HI-ratio), hypoxanthine/guanosine (HG-ratio,) xanthine/guanosine (XG-ratio), xanthine/inosine (XI-ratio), and ornithine/arginine (OA-ratio) as potential QIs in serum (Figures S1 and S2A). In EDTA plasma, only the “arginine and ornithine metabolism” was enriched (Data S2). The ratio of arginine to ornithine was also selected in EDTA plasma (Figure S2B).
Further potential QIs were sought by ChemRICH [17] analysis to overcome the limitations of pathway enrichment statistics, such as lacking pathway information for metabolites or metabolites common to different pathways. All metabolites with existing PubChem IDs were included in this analysis. Four hundred and nineteen serum metabolites were grouped into 68 metabolite clusters, of which 19 clusters show differences (Figure 1D). Among these metabolic clusters were saturated fatty acids, ethanolamines, linoleic acids, and dipeptides. Additionally, we observed changes in the OH-FA_17_3_1 cluster, comprising eicosanoids, such as 12-HETE and 15-HETE (Supplemental Data S1). Based on these results, metabolites of the eicosanoid pathway were selected for further investigation.
In EDTA plasma, 458 metabolites were grouped in 62 clusters, of which eight clusters showed differences (Figure 2D). The most pronounced changes were observed in saturated lysophosphatidylcholines, ethanolamines, and basic amino acid clusters (Supplemental Data S2).

2.2. Validation of Candidate QIs for TTC in Serum and EDTA Plasma Samples from Healthy Volunteers

Based on the targeted metabolite profiling and pathway enrichment, we first examined eight potential candidate QIs (xanthine, hypoxanthine, inosine, guanosine, dihydroorotate, N-carbamoylaspartate, 7-methylguanine, allantoin) from the purine and pyrimidine pathway in serum and EDTA plasma samples in a second, independent sample set of healthy volunteers (validation sample 1, n = 10). However, we investigated a more extended grid of 0.5, 1, 2 and 4 h TTC (0, only in EDTA-plasma) and used a quantitative MS/MS-based approach. We observed a TTC-dependent increase for xanthine, hypoxanthine, and dihydroorotic acid (Figure 3). Hypoxanthine showed the largest changes after 1, 2 and 4 h compared to 0.5 h pre-centrifugation delay. Of note, dihydroorotic acid showed the strongest increase after 24 h. As expected, the two purine nucleosides guanosine and inosine were sensitive to prolonged pre-processing and decreased after 2 and 4 h relative to the basal concentration at 0.5 h. To achieve robust QIs, we only calculated the ratios of strongly altered metabolites within a specific metabolic pathway. Ratios between hypoxanthine, inosine, xanthine, and guanosine (HI, HG, XI, and XG) were all increased after 2 and 4 h compared to 0.5 h (Figure 3). However, we observed no evidence of alterations in allantoin and 7-methylguanine in serum samples after a 4 h prolonged blood incubation (Figure S3).
In EDTA plasma, only hypoxanthine and dihydroorotic acid increased after 1, 2, and 4 h relative at 0.5 h. Allantoin, 7-methylguanine, inosine, and N-carbamoylaspartate were not altered after 4 h TTC (Figure S4). Guanosine and xanthine concentrations were below the limit of detection in EDTA plasma samples. Unfortunately, the metabolites adenosine and adenosine 3′,5′-cyclic monophosphate [18], which showed a substantial decrease in intensity in the discovery sample, were below the detection limit in the serum and EDTA plasma of the validation samples and could neither be analyzed nor integrated into the analysis of metabolite ratios.
Based on the results of the chemical similarity enrichment analysis, we investigated eicosanoids in the human serum of healthy individuals (validation sample 1). In general, most metabolites were increased after 4 h prolonged blood incubation compared to 0.5 h (Figure 4), except for 12-oxo-ETE, 11,12-DHET, and 14,15-DHET (Figure 4 and Figure S5). The most significant changes were observed for 12-HETE, 15-(S)-HETE, 12-(S)-HEPE, and (±)13-HODE, resulting in fold changes from 1.8 up to 11.2 (Table S1). The highest discrimination between the different time points was observed for 12-HETE. In EDTA plasma, we quantified 16 eicosanoids, of which only (±)9-HODE and (±)13-HODE were elevated (Figure S6A,B). Eicosanoid concentrations were higher in serum compared to EDTA plasma, with elevated concentrations of TxB2 (15.9 fold), 12-(S)-HHT (7.8 fold), 5-HETE (7 fold), and 12-HETE (29.4 fold) up to 29 fold at 0.5 h.
Additionally, we investigated the effect of prolonged blood incubation (0.5, 1, 2, and 4 h) on arginine and ornithine in the third sample set of healthy volunteers (validation sample 2, n = 23). In both serum and EDTA plasma, ornithine concentrations increased and arginine concentrations decreased (Figure 5). While ornithine concentrations were similar, arginine was about ~65% of the basal concentrations in EDTA plasma compared to serum. The OA-ratios were increased in serum and EDTA plasma after 1, 2 and 4 h pre-centrifugation delays.
We calculated candidate-specific cutoff values by analyzing multiple ROCs for all potential QIs (Table 1). Based on these calculations, guanosine (Sen. 83%, Spec. 100%) and hypoxanthine (Sen. 80%, Spec. 86%) were most suitable for discrimination of TTC > 2 h, whereas the HI- (Sen. 86%, Spec. 83%) and HG-ratio (Sen. 80%, Spec. 100%) showed the highest discrimination for TTC > 1 h. Additionally, the eicosanoids ± 13-HODE (Sen. 100%, Spec. 83%), 15(S)-HETE (Sen. 100%, Spec. 100%), 12(S)-HEPE (Sen. 100%, Spec. 100%), 8(S)-HETE (Sen. 100%, Spec. 100%), 12-HETE (Sen. 100%, Spec. 100%) and 12-oxo-ETE (Sen. 100%, Spec. 100%) were shown to predict a prolonged TTC > 2 h in serum with high sensitivity and specificity. Furthermore, when the optimal cutoff for the selected QIs was calculated for EDTA plasma, only the OA-ratio showed good sensitivity and specificity for identification of samples with a prolonged TTC > 1 h (Sen. 73%, Spec. 85%) or TTC > 2 h (Sen. 83%, Spec. 78%).

2.3. Validation of Candidate QIs for TTC in Serum and EDTA Plasma Samples from Patients with Rheumatic and Cardiovascular Diseases

To further validate the accuracy of potential candidate QIs, we applied their cutoffs to samples with known TTCs collected from patients with cardiovascular and rheumatic diseases (validation samples 3 and 4). Best results to predict samples with a prolonged TTC > 1 h in serum were observed using the HI- and HG-ratio. The calculated median optimal cutoffs of 1.41 (HI-ratio) and 2.06 (HG-ratio) achieved a Sen. of 79% and 100% and a Spec. of 85% and 73%, respectively. Interestingly, the combination of both thresholds, where either one or the other must be exceeded, results in a sensitivity of 100% and a specificity of 92%. Furthermore, several eicosanoids such as Tetranor-12(S)-HETE (Sen. 92%, Spec. 100%), 12(S)-HEPE (Sen. 100%, Spec. 86%), 8(S)-HETE (Sen. 100%, Spec. 86%), 12-HETE (Sen. 92%, Spec. 86%) and 12-oxo-ETE (Sen. 100%, Spec. 86%) were shown to be able to predict a prolonged TTC > 2 h in serum. In EDTA plasma samples, only the HI-ratio (Sen. 81%, Spec. 85%) could be identified as a promising QI for a TTC > 1 h. The AO-ratio could not be confirmed as QI in serum or EDTA plasma. A detailed overview of the results of all potential QIs is shown in Table 2. A summary of all potential quality indicators for further investigation is shown in Table 3.

3. Discussion

Delayed sample processing, in particular prolonged time to centrifugation, has a significant impact on the in vivo metabolome in blood derivatives and can substantially affect clinical diagnostics and research [19]. To ensure high quality of human serum and EDTA plasma samples in biobanks, it is essential to minimize TTC and to keep the dwell time of the samples at room temperature as short as possible [14]. However, due to the lack of standardized pre-analytical workflows [1] and quality assurance tools for workflow-validation, sample history is often incomplete or poorly documented and thus, the impact of pre-processing conditions on sample quality is unclear. For this reason, reliable QIs suitable for retrospective assessment of sample history and thus biospecimen quality are needed. While potential QIs reflecting deviations resulting from prolonged TTC have already been identified [5,6,20], their applicability to patients samples has not been sufficiently validated and demonstrated. Consequently, these QIs are not yet routinely used in biobanks workflows. This study describes several potential metabolite ratios or individual metabolites as QIs, identified in healthy cohorts and confirmed in samples from patients with cardiovascular and rheumatic diseases.

3.1. HI- and HG-Ratios as Potential Quality Indicators in Serum

After prolonged TTC at room temperature, we observed either an increase or a decrease in certain metabolites of the purine degradation pathway in serum and EDTA plasma. This variation is consistent with previous studies; for instance, Liu and colleagues [21] described an increase of relative hypoxanthine levels in serum after a delay of 4 h prior to processing.
Moreover, it has been shown that hypoxanthine rapidly increases in whole blood after storage at room temperature before centrifugation [22,23]. Hypoxanthine is endogenously produced in erythrocytes [22], and it is well known that guanosine and inosine are converted by the purine nucleoside phosphorylase (PNP—EC 2.4.2.1), an enzyme that is found in red blood cells (RBCs), to guanine and hypoxanthine [24]. Farthing and colleagues recommended determining inosine and hypoxanthine in plasma samples, as the clotting time of approximately 30 min required for a serum sample would allow significant conversion of inosine to hypoxanthine in the collection tube due to PNP activity from RBCs [25]. RBCs also function as ATP storage vesicles [26], therefore by ATP depletion of RBC [27] and AMP degradation, small polar substances, such as inosine and hypoxanthine, are formed [28] which could passively diffuse into serum or plasma. Interestingly, low PNP activity in plasma [29,30] might explain the stable inosine concentration and the less pronounced increase of hypoxanthine in EDTA plasma compared to serum.
Furthermore, we observed an increase of xanthine in serum samples after prolonged TTC, which might be due to further degradation of hypoxanthine [24] by xanthine oxidase (XO—EC 1.17.3.2). Low activity of XO in plasma might be the reason of inferior xanthine concentrations in EDTA plasma samples [29]. Thus, metabolite ratios of the purine metabolism are particularly suitable to detect TTC-dependent conversion and degradation processes due to the high activity of PNP and XO in serum.
According to the results obtained with samples from healthy individuals and their validation in patients with cardiovascular and rheumatic diseases, it can be concluded that serum samples with pre-centrifugation delays ≤ and >1 h can be reliably discriminated based on the HI- and HG-ratios with high sensitivity and specificity. Additionally, we have observed that the calculation of metabolite ratios may increase the performance of QIs compared to the use of their individual concentrations in cohorts of cardiovascular and rheumatic patients (Table 2). This was also shown by Jain and colleagues who demonstrate that combined QIs increase the prognostic value in samples across multiple clinical centers for a processing delay > 2 h [6].
Even though both HI- and HG-ratios provide promising results as QIs for prolonged TTC, components of the purine metabolism can be affected by multiple individual factors such as exercise [31,32], nutrition, or disease. For example, patients with acute cardiac ischemia [28,33] and multiple sclerosis [34] show altered or elevated levels of hypoxanthine and inosine. Hypoxanthine is also an indicator of hypoxia [35,36] and a marker of energy stress associated with intense physical training [31,32]. Furthermore, patients who received erythrocyte concentrates that are near the end of their shelf life have increased hypoxanthine concentrations of up to 100 µM [37].

3.2. Eicosanoids as Potential Quality Indicators in Serum

Our data indicate that several eicosanoids such as 8-(S)-HETE, 12-HETE, Tetranor-12-(S)-HETE, 15-(S)-HETE, 12-(S)-HEPE, and ±13-HODE increase in a time-dependent manner in human serum after a prolonged TTC. Eicosanoids are biologically active metabolites that are generated by enzymatic oxidation and rearrangement of free polyunsaturated fatty acids (PUFAs) such as arachidonic acid (AA), linolenic acid (LA), or eicosapentaenoic acid (EPA) through concerted action of cyclooxygenases (COX), lipoxygenase (LOX), and cytochrome P-450 (CYP), or via autoxidation [38,39]. These enzymes are strongly activated in platelets during aggregation and coagulation [40,41], induced by clot activators in serum samples immediately after blood collection. The previously observed high concentrations of eicosanoids in serum [42] compared to EDTA plasma results from the increased release and metabolism of PUFAs. For example, Dorow and colleagues found tremendously elevated levels of 11-, 12-, 15-HETE, 12-HHT, and TxB2 of up to 18,803% relative to EDTA plasma [43].
We show that some eicosanoids (8-(S)-HETE, 12-HETE, Tetranor-12-(S)-HETE, 15-(S)-HETE, 12-(S)-HEPE) are particularly suitable to identify samples with a TTC > 2 h with high sensitivity and specificity from both healthy individuals and patients with cardiovascular and rheumatic diseases (Table 1 and Table 2). However, the use of eicosanoids as QIs for delayed sample processing has clear limitations. For example, 12-HEPE and 12-HETE show significant changes (>2.77 × Coefficient of Variation) after two freeze-thaw cycles (FTC). Tetranor-12-HETE was strongly affected by 5 FTCs [43]. Additionally, the storage time at −80 °C can affect eicosanoid concentrations. Dorow and colleagues observed that eicosanoids such as 9-HODE or 13-HODE are unstable during storage for 6 months at −80 °C [43]. For this reason, the use individual eicosanoids as QI for retrospective analysis of already frozen serum samples should be discouraged. Instead, a combination of 8-(S)-HETE, 12-HETE, Tetranor-12-(S)-HETE, 15-(S)-HETE, 12-(S)-HEPE) should be used to estimate the TTC of these samples.

3.3. The Ratio of Ornithine to Arginine

In both serum and EDTA plasma, we observed a TTC-dependent increase in the ornithine/arginine ratio in samples kept at room temperature prior to centrifugation. This might be due to degradation of arginine to ornithine and urea by arginase, an enzyme highly expressed in erythrocytes [44]. Interestingly, the ability of the OA-ratio (Table 1 and Table 2) to discriminate TTC-dependent changes is not as specific and sensitive as predicted by Jain and colleagues [6]. This insufficient discriminatory power might be due to the age- and sex-associated differences between the discovery and the validation samples. While the original analysis was performed only with healthy male volunteers (discovery sample, n = 10), the validation has been done with samples of older female and male volunteers and patients. Saito and colleagues [45] have observed, that ornithine levels are associated with age and sex. In addition, patients with heart failure show increased arginine levels compared to healthy individuals [46].
Our study has several limitations. Dynamic changes in the composition of the serum and plasma metabolome are not exclusively altered by the pre-analytical phase, but are also influenced by several intrinsic and extrinsic factors, such as age and sex [45,47], the physiological rhythm [48], diet [49], exercise [50], fasting [51,52] or smoking status [53].
First, in this proof of concept study, the discovery samples were collected only from male participants with a limited age range compared to the validation cohorts, to focus entirely on the identification of potential quality indicators with as little influence as possible from strong inter-individual differences, for example, hormonal status. However, the successful validation of quality indicators in several heterogeneous validation cohorts points to the robustness of the identified QC indicators. Nevertheless, to avoid age-related bias, cohorts matched for age are recommended and should be used in future independent validation cohorts. Given that, for example, cardiovascular and rheumatic diseases are chronic age-related diseases that predominantly occur in the elderly, this might be challenging [54].
Second, for validation samples 2–4, 12-h fasting before blood collection could not be warranted, since samples were provided from pre-existing clinical biobank collections. Thus, several QC indicators from the group of amino acids, eicosanoids and nucleosides that were identified as quality indicators in the fasting discovery sample, may not have been successfully confirmed in the non-fasting validation samples 2–4.
Third, validation in this study was limited to two clinical entities, for example, patients with cardiovascular and rheumatic diseases. Thus, potential biomarker candidates, identified in this proof of concept study, which may serve as potential Qis, have to be prospectively validated in larger and independent sample cohorts, covering a wide range of clinical entities.

4. Materials and Methods

4.1. Study Design and Pre-Analytical Conditions

Forty-three self-reported healthy volunteers (discovery sample and validation sample 1 and 2), 60 cardiovascular patients and 79 rheumatic patients (validation sample 3 and 4) were included in this study (Figure S7). The demographic and clinical data are shown in Table 4.
Differing quantities (validation sample 1—106 mL, validation sample 2—33.1 mL) of blood were collected from each volunteer by peripheral venipuncture using a 21-gauge Safety-Multifly (Sarstedt AG & Co, KG, Nürmbrecht, Germany, www.sarstedt.com) blood collection system at room temperature. The volunteers were placed in an upright sitting position and the tourniquet released after blood flow started. We used 4.9 mL Serum-gel-S-Monovettes® with a cloth activator (04.1935.001) and 2.7 mL and 4.9 mL K3EDTA-plasma-S-Monovettes® (04.1917.001, 04.1931; Sarstedt AG & Co, KG, Nürmbrecht, Germany, www.sarstedt.com) for blood collection. EDTA plasma tubes were gently mixed by inverting the tubes five times immediately after blood collection. Serum tubes were kept for 30 min in an upright position to ensure a distinct separating layer after centrifugation (i.e., ‘sausage’ pattern). Blood collection tubes were centrifuged at 2500× g for 10 min at 20 °C (discovery sample, validation sample 1, 2, 3 and 4) and at 3309× g for 10 min at 20 °C (cardiovascular samples of validation sample 3 and 4). Blood collection tubes were handled at room temperature according to the conditions in Table 4. Collection time, pre- and post-centrifugation delay were monitored and recorded accordingly. All samples were divided into aliquots and were stored immediately after centrifugation at −80 °C in ultra-low freezers until shipping or further processing. Serum and EDTA plasma samples of all healthy volunteers (discovery sample, validation sample 1 and 2) were collected at the integrated biobank in Jena (IBBJ). Samples from validation sample 3 and 4 were leftover materials from a clinical ward that were divided into aliquots after routine laboratory diagnostics from patients with rheumatic diseases. Samples kindly provided by the Interdisciplinary Bank of Biomaterials and Data in Würzburg (ibdw) were collected from patients suffering from cardiovascular diseases (validation sample 3 and 4). All samples from clinical wards were collected at random collection times.

4.2. Metabolite Profiling in Human Serum and EDTA Plasma

Serum and EDTA plasma samples with a pre-centrifugation delay of 0.5 h and 2 h were profiled by Metabolon Inc. (Durham, NC, USA). Sample handling, sample preparation, and metabolite profiling protocols were as described previously [55,56]. In brief, serum and EDTA plasma samples were prepared using the automated MicroLab STAR® system (Hamilton Germany GmbH, Gräfelfing, Germany, www.hamiltoncompany.com), precipitated with methanol under vigorous shaking for 2 min (GenoGrinder2000, SPEX Sample Prep, Metuchen, NJ, USA, www.spexsampleprep.com) and then centrifuged. The supernatant was divided into five aliquots. These were subjected to four optimized methods for more hydrophilic, hydrophobic, polar, and nonpolar compounds and one aliquot was used as a backup. Chromatographic separation was achieved for all optimized methods using an ACQUITY (Waters Corporation, Milford, MA, USA) ultra-performance liquid chromatography (UPLC) system. The first method was performed using a Waters UPLC BEH C18 column (particle size 1.7 µm, column length 2.1 × 100 mm) and eluted by a gradient using a mixture of water and methanol, containing 0.05% perfluoropentanoic acid (PFPA) and 0.1% formic acid (FA). The second method utilized the same C18 column as the first but differed in the composition of elution solvents that were used for the gradient. Methanol, acetonitrile, and water, containing 0.05% PFPA and 0.1% FA were used at an overall higher organic content. The third method used a separate dedicated C18 column, and the elution was achieved by a gradient using methanol and water, containing 6.5 mM ammonium bicarbonate at pH 8. The fourth method was performed by using a HILIC column (Waters UPLC BEH Amide particle size 1.7 µm, column length 2.1 × 150 mm). A mixture of water and acetonitrile, containing 10 mM ammonium formate at pH 10.8, achieved the elution.
MS detection was performed using a Q-Exactive mass spectrometer (MS) equipped with a heated electrospray ionization source (Thermo Fisher Scientific, Waltham, MA, USA, www.thermofisher.com) operating in positive (first and second method) and negative (third and fourth method) ion mode. The Orbitrap mass analyzer operated at a mass resolution of 35.000 and alternated between MS and data-dependent MSn scans using dynamic exclusion. The scan range covered 70–1000 m/z.
Data extraction and peak identification were performed using Metabolon’s hardware and software [57,58]. Metabolites were identified using a reference library, which is based on authenticated standards that include retention time, accurate m/z ratios (±10 ppm) and associated MS/MS scores. These scores are based on the comparison between reference and experimental spectra of the ions of interest. This process is curated by visual inspection for quality control using software developed at Metabolon [55,56]. Peaks were quantified using the area-under-the-curve.

4.3. Quantitative Analysis of Nucleosides and Related Compounds

Sample analysis was carried out by MS-Omics (MSOmics, Vedbæk, Denmark, www.msomics.com) as follows. Proteins were precipitated using methanol, and the metabolites were extracted using a liquid–liquid extraction using chloroform and water. The aqueous phase was dried under nitrogen flow. Samples were reconstituted in eluent A (ultrapure water with 10 mM ammonium formate + 0.1% formic acid). The analysis was carried out using a UPLC system (UPLC ACQUITY, Waters) coupled with a time of flight mass spectrometer (Xevo G2 ToF, Waters). An electrospray ionization interface was used. The analysis was performed in negative and positive ionization modes. The UPLC was performed using a slightly modified version of the protocol described previously [59]. Data-processing used MZmine 2 [60] followed by curation using a custom made an in-house protocol. Identification of compounds was performed using both peak retention times (compared against authentic standards included in the analytical sequence) and accurate mass (with an acceptable deviation of 0.005 Da).

4.4. Quantitative Analysis of Eicosanoids

Serum and EDTA plasma samples were used to determine polyunsaturated fatty acids and eicosanoids. Sample handling, sample preparation, and metabolite measurements were as described previously [43,61]. In brief, 200 µL of serum or EDTA plasma were mixed with 450 µL precipitation solution containing the deuterium-labeled standard solution (c = 5 ng/mL in methanol-water 50:50 [v/v] and 50 ng/mL for AA-d8) under vigorous shaking for 2 min. Samples were centrifuged at 10,000× g for 5 min; the supernatant was transferred into a new vial and then stored at −80 °C until further analysis.
A Shimadzu UFLC LC-20A Prominence liquid chromatography system (Shimadzu Deutschland GmbH, Duisburg, Germany, www.shimadzu.de) was used. Online solid-phase extraction was achieved by a Strata-X column (particle size 25 µm; column length 20 × 2 mm i.d., Phenomenex, Aschaffenburg, Germany). Chromatographic separation of analytes was performed using a Kinetex C18 column (particle size 2.6 µm, column length 100 × 2.1 mm i.d., Phenomenex).
MS detection was performed using a 5500 QTrap mass spectrometer (AB Sciex LLC, Framingham, MA, USA) equipped with an electrospray ionization source in negative ionization mode. Scheduled multiple reaction monitoring experiments were used for the quantitative analysis of eicosanoids. In serum 14 eicosanoids and in EDTA plasma 16 eicosanoids could be analyzed.

4.5. Quantitative Analysis of Amino Acids

Calibration was performed using a physiological amino acid standard solved in 0.1 N HCl purchased from Laborservice Onken GmbH (product number: 5.403.151, Gründau, Germany). Amino acids had a concentration of 1 mmol/L, except for cystin (500 µmol/L), cystathionine (525 µmol/L), and urea (15 mmol/L). Additionally, alloisoleucin was solved in deionized water to a final concentration of 1 mmol/L. 100 µL of both stock solutions were combined, diluted using 300 µL of loading buffer and frozen at −80 °C until use. ClinCheck® plasma control (level 1) for amino acids was used as internal quality control material (Recipe, Germany).
Next, 400 μL of sample material (serum and EDTA plasma) or 200 µL of controls were mixed with 100 μL/50 µL of a solution containing 10% sulfosalicylic acid to precipitate proteins. Samples were then put in the fridge at 2–8 °C for 0.5 h and then centrifuged at 16,060× g for 15 min. 100 μL of supernatant were diluted with 100 µL lithium loading buffer, at a pH 2.20 (product code: 80-2038-10) and vigorously mixed, and 30 µL were injected and measured.
The amino acids arginine and ornithine were measured using a BioChrom 30plus high-performance liquid chromatography cation-exchange system with ninhydrin detection (Biochrom Ltd., Cambridge, UK, www.biochrom.co.uk). For further details of the widely adopted method, see instructions for Biochrom 30+ Amino Acid Analyzers (2018) (Version 41 56 1783 IVD instruction for use English iss 16, http://www.biochrom.co.uk/user_downloads/?c=17). Amino acid concentrations were calculated offline using EZChrom Elite software (Version 3.3.1, Agilent, Santa Clara, CA, USA).

4.6. Boinformatics and Statistics for QI Discovery by Targeted Metabolite Profiling

Targeted metabolite profiling data were filtered, and missing values imputed using the MetImp 1.2 (https://metabolomics.cc.hawaii.edu/software/MetImp/, accessed on 30 June 2020) web tool. Filtering was applied by using the “modified 80% rule”, meaning that variables can be excluded from the data when the proportion of non-missing elements are accounted for less than 80% among each biological group. Missing values were imputed using the Gibbs sampler left-censored missing value imputation approach [62]. Profiling data were log10 transformed and compared using the paired t-test. To address multiple comparisons, we derived q-values [63] that extended the idea of the false discovery rates proposed by Benjamini and Hochberg [64].
Pathway enrichment analysis was performed using MetaboAnalyst 3.0 [65,66,67]. KEGG was selected as a reference pathway library. Altered metabolites in serum and EDTA plasma were uploaded to MetaboAnalyst 3.0 using their HMDB IDs and were analyzed using over-representation analysis (ORA) [68]. The selected ORA method was Fishers’ exact test. P-Values were adjusted for multiple testing using the Holm–Bonferroni method. We selected the node importance measure: out-degree centrality for topological analysis [69]. Additionally, we performed a chemical similarity enrichment analysis (ChemRICH –chemrich.fiehnlab.ucdavis.edu) and ontology mapping as described previously [17,70]. InChiKeys and SMILES codes were obtained from the PubChem identifier exchange service (https://pubchem.ncbi.nlm.nih.gov/idexchange/idexchange.cgi, accessed on 30 June 2020). Chemical enrichment statistics were calculated using the Kolmogorov–Smirnov test.
Data were also analyzed using Excel (Microsoft, Albuquerque, NM, USA) and OriginPro 2018 (Origin Lab Cooperation, Northampton, MA, USA).

4.7. Statistical Analyses in Validation Samples

In the case of multiple measurements for each healthy volunteer, we applied paired t-tests. Within one panel figure, all reported p-values were Bonferroni corrected for the number of comparisons. We generated receiver operating characteristic curves and areas under the ROC such that the curves contained only one measurement for each healthy volunteer—in the case of multiple measures, we randomly sampled from the measurements without replacement and generated the distribution via random draws. Data were analyzed in R version 3.6.2 [71] using the packages tidyverse [72], caret [73], pROC [74] and ROCR [75].

5. Conclusions

Controlling and assessing the quality and, thus, the fitness-for-purpose of samples is of significant importance for biobanks and research. However, frequently, biobanks are unable to completely control each pre-analytical step outside their direct area of responsibility and might receive samples with inadequately documented history. To overcome these shortcomings, reliable and versatile methods must be in place to ascertain the major pre-analytical factors affecting sample quality in order to achieve biobanking quality objectives. Hence, suitable methods and appropriate biomarkers, serving as QIs, must receive particular attention.
In this study, we have identified potential biomarker candidates, which may serve as QIs for pre-centrifugation delays. In particular, the HI-ratio, HG-ratio and the concentration of several eicosanoids are sensitive to a TTC longer than 1 or 2 h. For some of the selected QIs, we demonstrated high sensitivity and specificity for the discrimination of samples with different TTCs. While the QIs mentioned are promising, further independent prospective studies covering a wider spectrum of diseases are necessary to validate the proposed QIs and establish them for routine quality checks in biobanks.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/metabo11090638/s1, Figure S1: Purine derivatives and their ratios as potential target QIs in human serum, Figure S2: Arginine, ornithine and AO-ratio as potential target QIs in human serum and EDTA plasma, Figure S3: Precentrifugation delay dependent impact on allantoin and 7-methylguanine in serum, Figure S4: Precentrifugation delay dependent impact on selected metabolites in EDTA plasma, Figure S5: Influence of precentrifugation delay on eicosanoids in serum, Figure S6A: Influence of precentrifugation delay on eicosanoids in EDTA plasma, Figure S6B: Influence of precentrifugation delay on eicosanoids in EDTA plasma, Table S1: Fold changes of eicosanoids in serum and EDTA plasma, Data S1 —Serum, Data S2 —EDTA Plasma.

Author Contributions

Conceptualization, S.H., M.K. and A.S.; methodology, S.H., J.W., N.K.; validation, S.H.; formal analysis (especially explorative bioinformatics and pathway analyses), S.H., F.S.; investigation, S.H., M.R., N.K.; resources, M.K., U.C., R.J.; data curation, S.H., J.W.; writing—original draft preparation, S.H.; writing—review and editing, N.K., F.S., J.G., J.W., A.S., M.K.; visualization, S.H.; supervision, M.K.; project administration, S.H.; funding acquisition, M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the German Ministry of Education and Research (BMBF No. 01EK1505B).

Institutional Review Board Statement

The study protocol adheres to the requirements of the Declaration of Helsinki [76] and Taipei [77] and received a positive vote by the ethics committee of the Friedrich Schiller University Jena (registration number 4893-08/16, 3335-01/12 and 3921-11/13). The ibdw kindly provided samples of patients suffering from cardiovascular diseases. The ethics committee of the University of Würzburg (No. 149/11 and 278/17-Z-me) provides positive votes for the usage of the biospecimen and data for this study.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available upon request from the corresponding authors. The data are not publicly available due to the privacy of the patients assisted in the research.

Acknowledgments

We thank S. Bremer-Streck, N. Blaurock, S. Hein and B. Preuschaft for their technical support. We thank the staff of all recruitment sites for their invaluable help with the sample collection and preparation and C. Richert, J. Köhler, K. Stötzer and E. Hanemann for the technical assistance at the integrated biobank in Jena. Furthermore, we thank all participating volunteers and patients for their support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 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] [PubMed] [Green Version]
  2. Yin, P.; Peter, A.; Franken, H.; Zhao, X.; Neukamm, S.S.; Rosenbaum, L.; Lucio, M.; Zell, A.; Haring, H.U.; Xu, G.; et al. Preanalytical aspects and sample quality assessment in metabolomics studies of human blood. Clin. Chem. 2013, 59, 833–845. [Google Scholar] [CrossRef] [PubMed]
  3. Breier, M.; Wahl, S.; Prehn, C.; Fugmann, M.; Ferrari, U.; Weise, M.; Banning, F.; Seissler, J.; Grallert, H.; Adamski, J.; et al. Targeted metabolomics identifies reliable and stable metabolites in human serum and plasma samples. PLoS ONE 2014, 9, e89728. [Google Scholar] [CrossRef] [PubMed]
  4. Kamlage, B.; Maldonado, S.G.; Bethan, B.; Peter, E.; Schmitz, O.; Liebenberg, V.; Schatz, P. Quality markers addressing preanalytical variations of blood and plasma processing identified by broad and targeted metabolite profiling. Clin. Chem. 2014, 60, 399–412. [Google Scholar] [CrossRef] [PubMed]
  5. Jobard, E.; Tredan, O.; Postoly, D.; Andre, F.; Martin, A.L.; Elena-Herrmann, B.; Boyault, S. A systematic evaluation of blood serum and plasma pre-analytics for metabolomics cohort studies. Int. J. Mol. Sci. 2016, 17, 2035. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Jain, M.; Kennedy, A.D.; Elsea, S.H.; Miller, M.J. Analytes related to erythrocyte metabolism are reliable biomarkers for preanalytical error due to delayed plasma processing in metabolomics studies. Clin. Chim. Acta 2017, 466, 105–111. [Google Scholar] [CrossRef] [Green Version]
  7. Kamlage, B.; Neuber, S.; Bethan, B.; Maldonado, S.G.; Wagner-Golbs, A.; Peter, E.; Schmitz, O.; Schatz, P. Impact of prolonged blood incubation and extended serum storage at room temperature on the human serum metabolome. Metabolites 2018, 8, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Gottfried-Blackmore, A.; Rubin, S.J.S.; Bai, L.; Aluko, S.; Yang, Y.; Park, W.; Habtezion, A. Effects of processing conditions on stability of immune analytes in human blood. Sci. Rep. 2020, 10, 17328. [Google Scholar] [CrossRef]
  9. Thumfart, J.O.; Findeisen, P.; Neumaier, M. Systematic identification of endogenous decay-markers for quality assesment of serum specimens. Clin. Chem. Lab. Med. 2011, 49, S221. [Google Scholar]
  10. Anton, G.; Wilson, R.; Yu, Z.H.; Prehn, C.; Zukunft, S.; Adamski, J.; Heier, M.; Meisinger, C.; Romisch-Margl, W.; Wang-Sattler, R.; et al. Pre-analytical sample quality: Metabolite ratios as an intrinsic marker for prolonged room temperature exposure of serum samples. PLoS ONE 2015, 10, e0121495. [Google Scholar] [CrossRef] [PubMed]
  11. Thumfart, J.O.; Abidi, N.; Mindt, S.; Costina, V.; Hofheinz, R.; Klawonn, F.; Neumaier, M.; Findeisen, P. Lc/ms based monitoring of endogenous decay markers for quality assessment of serum specimens. J. Proteomics Bioinform. 2015, 8, 91–97. [Google Scholar]
  12. Liu, X.; Hoene, M.; Yin, P.; Fritsche, L.; Plomgaard, P.; Hansen, J.S.; Nakas, C.T.; Niess, A.M.; Hudemann, J.; Haap, M.; et al. Quality control of serum and plasma by quantification of (4e,14z)-sphingadienine-c18-1-phosphate uncovers common preanalytical errors during handling of whole blood. Clin. Chem. 2018, 64, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Schwarz, N.; Knutti, N.; Rose, M.; Neugebauer, S.; Geiger, J.; Jahns, R.; Klopp, N.; Illig, T.; Mathay, C.; Betsou, F.; et al. Quality assessment of the preanalytical workflow in liquid biobanking: Taurine as a serum-specific quality indicator for preanalytical process variations. Biopreserv. Biobank. 2019, 17, 458–467. [Google Scholar] [CrossRef] [PubMed]
  14. Stevens, V.L.; Hoover, E.; Wang, Y.; Zanetti, K.A. Pre-analytical factors that affect metabolite stability in human urine, plasma, and serum: A review. Metabolites 2019, 9, 156. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Teahan, O.; Gamble, S.; Holmes, E.; Waxman, J.; Nicholson, J.K.; Bevan, C.; Keun, H.C. Impact of analytical bias in metabonomic studies of human blood serum and plasma. Anal. Chem. 2006, 78, 4307–4318. [Google Scholar] [CrossRef] [PubMed]
  16. Hirayama, A.; Sugimoto, M.; Suzuki, A.; Hatakeyama, Y.; Enomoto, A.; Harada, S.; Soga, T.; Tomita, M.; Takebayashi, T. Effects of processing and storage conditions on charged metabolomic profiles in blood. Electrophoresis 2015, 36, 2148–2155. [Google Scholar] [CrossRef] [PubMed]
  17. Barupal, D.K.; Fiehn, O. Chemical similarity enrichment analysis (chemrich) as alternative to biochemical pathway mapping for metabolomic datasets. Sci. Rep. 2017, 7, 1–11. [Google Scholar] [CrossRef]
  18. del Campo, M.L.; Halitschke, R.; Short, S.M.; Lazzaro, B.P.; Kessler, A. Dietary plant phenolic improves survival of bacterial infection in manduca sexta caterpillars. Entomol. Exp. Appl. 2013, 146, 321–331. [Google Scholar] [CrossRef] [Green Version]
  19. Plebani, M.; Sciacovelli, L.; Aita, A.; Padoan, A.; Chiozza, M.L. Quality indicators to detect pre-analytical errors in laboratory testing. Clin. Chim. Acta 2014, 432, 44–48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Trezzi, J.P.; Bulla, A.; Bellora, C.; Rose, M.; Lescuyer, P.; Kiehntopf, M.; Hiller, K.; Betsou, F. Lacascore: A novel plasma sample quality control tool based on ascorbic acid and lactic acid levels. Metabolomics 2016, 12, 96. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Liu, L.S.; Aa, J.Y.; Wang, G.J.; Yan, B.; Zhang, Y.; Wang, X.W.; Zhao, C.Y.; Cao, B.; Shi, J.A.; Li, M.J.; et al. Differences in metabolite profile between blood plasma and serum. Anal. Biochem. 2010, 406, 105–112. [Google Scholar] [CrossRef] [PubMed]
  22. Boulieu, R.; Bory, C.; Baltassat, P.; Gonnet, C. Hypoxanthine and xanthine levels determined by high-performance liquid-chromatography in plasma, erythrocyte, and urine samples from healthy-subjects—The problem of hypoxanthine level evolution as a function of time. Anal. Biochem. 1983, 129, 398–404. [Google Scholar] [CrossRef]
  23. Hiroshige, K.; Takasugi, M.; Yuu, K.; Kuroiwa, A. Studies on the purine-pyrimidine metabolism (2)—Effects of sample preparation and diet on the concentrations of oxypurines in blood and urine. J. UOEH 1994, 16, 167–178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Giuliani, P.; Zuccarini, M.; Buccella, S.; Rossini, M.; D’Alimonte, I.; Ciccarelli, R.; Marzo, M.; Marzo, A.; Di Iorio, P.; Caciagli, F. Development of a new hplc method using fluorescence detection without derivatization for determining purine nucleoside phosphorylase activity in human plasma. J. Chromatogr. B 2016, 1009, 114–121. [Google Scholar] [CrossRef] [PubMed]
  25. Farthing, D.; Sica, D.; Gehr, T.; Wilson, B.; Fakhry, I.; Larus, T.; Farthing, C.; Karnes, H.T. An hplc method for determination of inosine and hypoxanthine in human plasma from healthy volunteers and patients presenting with potential acute cardiac ischemia. J. Chromatogr. B 2007, 854, 158–164. [Google Scholar] [CrossRef] [PubMed]
  26. Xu, Z.S.; Dou, W.K.; Wang, C.; Sun, Y. Stiffness and atp recovery of stored red blood cells in serum. Microsyst. Nanoeng. 2019, 5, 1–9. [Google Scholar] [CrossRef] [PubMed]
  27. Harkness, R.A. Hypoxanthine, xanthine and uridine in body-fluids, indicators of atp depletion. J. Chromatogr.-Biomed. 1988, 429, 255–278. [Google Scholar] [CrossRef]
  28. Farthing, D.E.; Farthing, C.A.; Xi, L. Inosine and hypoxanthine as novel biomarkers for cardiac ischemia: From bench to point-of-care. Exp. Biol. Med. 2015, 240, 821–831. [Google Scholar] [CrossRef] [PubMed]
  29. Yamamoto, T.; Moriwaki, Y.; Takahashi, S.; Tsutsumi, Z.; Yamakita, J.; Nasako, Y.; Hiroishi, K.; Higashino, K. Determination of human plasma xanthine oxidase activity by high-performance liquid chromatography. J. Chromatogr. B 1996, 681, 395–400. [Google Scholar] [CrossRef]
  30. Viegas, T.X.; Omura, G.A.; Stoltz, R.R.; Kisicki, J. Pharmacokinetics and pharmacodynamics of peldesine (bcx-34), a purine nucleoside phosphorylase inhibitor, following single and multiple oral doses in healthy volunteers. J. Clin. Pharmacol. 2000, 40, 410–420, Erratum in 2000, 40, 664–664. [Google Scholar] [CrossRef] [PubMed]
  31. Dudzinska, W.; Lubkowska, A.; Dolegowska, B.; Safranow, K. Blood uridine concentration may be an indicator of the degradation of pyrimidine nucleotides during physical exercise with increasing intensity. J. Physiol. Biochem. 2010, 66, 189–196. [Google Scholar] [CrossRef] [PubMed]
  32. Dudzinska, W.; Suska, M.; Lubkowska, A.; Jakubowska, K.; Olszewska, M.; Safranow, K.; Chlubek, D. Comparison of human erythrocyte purine nucleotide metabolism and blood purine and pyrimidine degradation product concentrations before and after acute exercise in trained and sedentary subjects. J. Physiol. Sci. 2018, 68, 293–305. [Google Scholar] [CrossRef] [PubMed]
  33. Farthing, D.; Xi, L.; Gehr, L.; Sica, D.; Larus, T.; Karnes, H.T. High-performance liquid chromatography (hplc) determination of inosine, a potential biomarker for initial cardiac ischaemia, using isolated mouse hearts. Biomarkers 2006, 11, 449–459. [Google Scholar] [CrossRef] [PubMed]
  34. Tavazzi, B.; Batocchi, A.P.; Amorini, A.M.; Nociti, V.; D’Urso, S.; Longo, S.; Gullotta, S.; Picardi, M.; Lazzarino, G. Serum metabolic profile in multiple sclerosis patients. Mult. Scler. Int. 2011, 2011, 167156. [Google Scholar] [CrossRef] [PubMed]
  35. Sutton, J.; Fox, I.H. Alterations of purine metabolism during strenuous muscular exercise in man. Clin. Res. 1975, 23, A639. [Google Scholar]
  36. Bangsbo, J.; Sjodin, B.; Hellstenwesting, Y. Exchange of hypoxanthine in muscle during intense exercise in man. Acta Physiol. Scand. 1992, 146, 549–550. [Google Scholar] [CrossRef] [PubMed]
  37. Casali, E.; Berni, P.; Spisni, A.; Baricchi, R.; Pertinhez, T.A. Hypoxanthine: A new paradigm to interpret the origin of transfusion toxicity. Blood Transfus. 2016, 14, 555–556. [Google Scholar] [PubMed]
  38. Funk, C.D. Prostaglandins and leukotrienes: Advances in eicosanoid biology. Science 2001, 294, 1871–1875. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Bazan, N.G.; Stark, D.T.; Petasis, N.A. Lipid mediators: Eicosanoids, docosanoids and platelet-activating factor. In Basic Neurochemistry, 8th ed.; Brady, S.T., Siegel, G.J., Albers, R.W., Price, D.L., Eds.; Academic Press: New York, NY, USA, 2012; pp. 643–662. [Google Scholar]
  40. Pasqualini, M.E.; Mohn, C.E.; Petiti, J.P.; Manzo, P.; Eynard, A.R. Cox and lox eicosanoids modulate platelet activation and procoagulation induced by two murine cancer cells. Prostaglandins Leukot. Essent. Fat. Acids (PLEFA) 2000, 63, 377–383. [Google Scholar] [CrossRef] [PubMed]
  41. Dolegowska, B.; Blogowski, W.; Kedzierska, K.; Safranow, K.; Jakubowska, K.; Olszewska, M.; Rac, M.; Chlubek, D.; Ciechanowski, K. Platelets arachidonic acid metabolism in patients with essential hypertension. Platelets 2009, 20, 242–249. [Google Scholar] [CrossRef] [PubMed]
  42. Rund, K.M.; Nolte, F.; Doricic, J.; Greite, R.; Schott, S.; Lichtinghagen, R.; Gueler, F.; Schebb, N.H. Clinical blood sampling for oxylipin analysis—Effect of storage and pneumatic tube transport of blood on free and total oxylipin profile in human plasma and serum. Analyst 2020, 145, 2378–2388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Dorow, J.; Becker, S.; Kortz, L.; Thiery, J.; Hauschildt, S.; Ceglarek, U. Preanalytical investigation of polyunsaturated fatty acids and eicosanoids in human plasma by liquid chromatography-tandem mass spectrometry. Biopreserv. Biobank. 2016, 14, 107–113. [Google Scholar] [CrossRef] [PubMed]
  44. Spector, E.B.; Rice, S.C.H.; Kern, R.M.; Hendrickson, R.; Cederbaum, S.D. Comparison of arginase activity in red-blood-cells of lower mammals, primates, and man—Evolution to high-activity in primates. Am. J. Hum. Genet. 1985, 37, 1138–1145. [Google Scholar] [PubMed]
  45. Saito, K.; Maekawa, K.; Pappan, K.; Urata, M.; Ishikawa, M.; Kumagai, Y.; Saito, Y. The difference in the metabolite profiles between plasma and serum, ages or sexes, and their inter-individual variations in human subjects. Drug Metab. Rev. 2014, 45, 47. [Google Scholar]
  46. Liu, C.; Li, R.H.; Liu, Y.; Li, Z.G.; Sun, Y.J.; Yin, P.Y.; Huang, R.H. Characteristics of blood metabolic profile in coronary heart disease, dilated cardiomyopathy and valvular heart disease induced heart failure. Front. Cardiovasc. Med. 2021, 7, 622236. [Google Scholar] [CrossRef]
  47. Ishikawa, M.; Maekawa, K.; Saito, K.; Senoo, Y.; Urata, M.; Murayama, M.; Tajima, Y.; Kumagai, Y.; Saito, Y. Plasma and serum lipidomics of healthy white adults shows characteristic profiles by subjects’ gender and age. PLoS ONE 2014, 9, e91806. [Google Scholar] [CrossRef] [Green Version]
  48. Minami, Y.; Kasukawa, T.; Kakazu, Y.; Iigo, M.; Sugimoto, M.; Ikeda, S.; Yasui, A.; van der Horst, G.T.J.; Soga, T.; Ueda, H.R. Measurement of internal body time by blood metabolomics. Proc. Natl. Acad. Sci. USA 2009, 106, 9890–9895. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Gibney, M.J.; Walsh, M.; Brennan, L.; Roche, H.M.; German, B.; van Ommen, B. Metabolomics in human nutrition: Opportunities and challenges. Am. J. Clin. Nutr. 2005, 82, 497–503. [Google Scholar] [CrossRef] [PubMed]
  50. Weigert, C.; Lehmann, R.; Hartwig, S.; Lehr, S. The secretome of the working human skeletal muscle-a promising opportunity to combat the metabolic disaster? Proteom. Clin. Appl. 2014, 8, 5–18. [Google Scholar] [CrossRef] [PubMed]
  51. Brauer, R.; Leichtle, A.B.; Fiedler, G.M.; Thiery, J.; Ceglarek, U. Preanalytical standardization of amino acid and acylcarnitine metabolite profiling in human blood using tandem mass spectrometry. Metabolomics 2011, 7, 344–352. [Google Scholar] [CrossRef]
  52. Gillio-Meina, C.; Cepinskas, G.; Cecchini, E.L.; Fraser, D.D. Translational research in pediatrics ii: Blood collection, processing, shipping, and storage. Pediatrics 2013, 131, 754–766. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Beuchel, C.; Becker, S.; Dittrich, J.; Kirsten, H.; Toenjes, A.; Stumvoll, M.; Loeffler, M.; Thiele, H.; Beutner, F.; Thiery, J.; et al. Clinical and lifestyle related factors influencing whole blood metabolite levels—A comparative analysis of three large cohorts. Mol. Metab. 2019, 29, 76–85. [Google Scholar] [CrossRef] [PubMed]
  54. MacNee, W.; Rabinovich, R.A.; Choudhury, G. Ageing and the border between health and disease. Eur. Respir. J. 2014, 44, 1332–1352. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Evans, A.M.; Bridgewater, B.R.; Liu, Q.; Mitchell, M.W.; Robinson, R.J.; Dai, H.; Stewart, S.J.; DeHaven, C.D.; Miller, L.A.D. High resolution mass spectrometry improves data quantity and quality as compared to unit mass resolution mass spectrometry in high-throughput profiling metabolomics. J. Metab. JOM 2014, 4, 1. [Google Scholar]
  56. Evans, A.M.; DeHaven, C.D.; Barrett, T.; Mitchell, M.; Milgram, E. Integrated, nontargeted ultrahigh performance liquid chromatography/electrospray ionization tandem mass spectrometry platform for the identification and relative quantification of the small-molecule complement of biological systems. Anal. Chem. 2009, 81, 6656–6667. [Google Scholar] [CrossRef] [PubMed]
  57. DeHaven, C.; Evans, A.M.; Dai, H.; Lawton, K.A. Software techniques for enabling high-throughput analysis of metabolomic datasets. Metabolomics 2012, 10, 167–192. [Google Scholar] [CrossRef] [Green Version]
  58. Dehaven, C.D.; Evans, A.M.; Dai, H.; Lawton, K.A. Organization of gc/ms and lc/ms metabolomics data into chemical libraries. J. Cheminform. 2010, 2, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Doneanu, C.E.; Chen, W.; Mazzeo, J.R. Uplc/ms monitoring of water-soluble vitamin bs in cell culture media in minutes. Waters Appl. Note 2011, 1–7. Available online: https://www.waters.com/waters/library.htm?locale=en_US&lid=134636355 (accessed on 17 September 2021).
  60. Pluskal, T.; Castillo, S.; Villar-Briones, A.; Oresic, M. Mzmine 2: Modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular profile data. BMC Bioinform. 2010, 11, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Kortz, L.; Dorow, J.; Becker, S.; Thiery, J.; Ceglarek, U. Fast liquid chromatography-quadrupole linear ion trap-mass spectrometry analysis of polyunsaturated fatty acids and eicosanoids in human plasma. J. Chromatogr. B 2013, 927, 209–213. [Google Scholar] [CrossRef] [PubMed]
  62. Wei, R.M.; Wang, J.Y.; Jia, E.; Chen, T.L.; Ni, Y.; Jia, W. Gsimp: A gibbs sampler based left-censored missing value imputation approach for metabolomics studies. PLoS Comput. Biol. 2018, 14, e1005973. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Storey, J.D. The positive false discovery rate: A bayesian interpretation and the q-value. Ann. Stat. 2003, 31, 2013–2035. [Google Scholar] [CrossRef]
  64. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate—A practical and powerful approach to multiple testing. J. R. Stat. Soc. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  65. Chong, J.; Soufan, O.; Li, C.; Caraus, I.; Li, S.Z.; Bourque, G.; Wishart, D.S.; Xia, J.G. Metaboanalyst 4.0: Towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018, 46, W486–W494. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Xia, J.; Wishart, D.S. Using metaboanalyst 3.0 for comprehensive metabolomics data analysis. Curr. Protoc. Bioinform. 2016, 55, 14.10.1–14.10.91. [Google Scholar] [CrossRef] [PubMed]
  67. Xia, J.G.; Sinelnikov, I.V.; Han, B.; Wishart, D.S. Metaboanalyst 3.0-making metabolomics more meaningful. Nucleic Acids Res. 2015, 43, W251–W257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Khatri, P.; Draghici, S. Ontological analysis of gene expression data: Current tools, limitations, and open problems. Bioinformatics 2005, 21, 3587–3595. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Aittokallio, T.; Schwikowski, B. Graph-based methods for analysing networks in cell biology. Brief Bioinform. 2006, 7, 243–255. [Google Scholar] [CrossRef] [PubMed]
  70. Barupal, D.K.; Fan, S.L.; Fiehn, O. Integrating bioinformatics approaches for a comprehensive interpretation of metabolomics datasets. Curr Opin Biotech 2018, 54, 1–9. [Google Scholar] [CrossRef] [PubMed]
  71. R Core Team. R: A language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  72. Wickham, H.; Averick, M.; Bryan, J.; Chang, W.; McGowan, L.D.A.; François, R.; Grolemund, G.; Hayes, A.; Henry, L.; Hester, J.; et al. Welcome to the “tidyverse”. J. Open Source Softw. 2019, 4, 1686. [Google Scholar] [CrossRef]
  73. Kuhn, M. Caret: Classification and Regression Training. 2021. Available online: https://rdrr.io/cran/caret/f/vignettes/caret.Rmd (accessed on 30 June 2020).
  74. Robin, X.; Turck, N.; Hainard, A.; Tiberti, N.; Lisacek, F.; Sanchez, J.C.; Müller, M. Proc: An open-source package for r and s+ to analyze and compare roc curves. BMC Bioinform. 2011, 12, 77. [Google Scholar] [CrossRef] [PubMed]
  75. Sing, T.; Sander, O.; Beerenwinkel, N.; Lengauer, T. Rocr: Visualizing classifier performance in r. Bioinformatics 2005, 21, 7881. [Google Scholar] [CrossRef] [PubMed]
  76. World Medical Association. World medical association declaration of helsinki: Ethical principles for medical research involving human subjects. JAMA 2013, 310, 2191–2194. [Google Scholar]
  77. World Medical Association. WMA Declaration of Taipei on Ethical Considerations Regarding Health Databases and Biobanks. Secondary WMA Declaration of Taipei on Ethical Considerations Regarding Health Databases and Biobanks. 2016. Available online: https://www.wma.net/policies-post/wma-declaration-of-taipei-on-ethical-considerations-regarding-health-databases-and-biobanks/ (accessed on 30 June 2020).
Figure 1. Chemical similarity and pathway enrichment analysis to identify novel QIs for the precentrifugation delay in human serum. (A) Volcano plot of serum metabolites visualizes the effects of prolonged blood incubation after 2 h before centrifugation. Red (increase) and green (decrease) circles indicate altered metabolite levels (discovery sample, n = 10). A paired t-test was performed, and we applied addressed multiple testing by deriving q-values to control the positive false discovery rate. (B) Represents the pathway enrichment analysis using the altered metabolites of (A). Node size represents the impact of the altered metabolites on the pathways and node color visualizes the q-value. (D) ChemRICH set enrichment statistics plot. Each node reflects an altered cluster of metabolites (p-value ≤ 0.05). Node sizes represent the total number of metabolites in each cluster set. Node color scale shows the proportion of increased (red) or decreased (blue) compounds in serum samples with a prolonged blood incubation of 2 h compared to the control at 0.5 h. Purple color nodes include both increased and decreased metabolites. (C) The pathway and (E) cluster occupancy visualize the ratio (in %) between increased and decreased metabolites within each altered pathway/cluster of (A).
Figure 1. Chemical similarity and pathway enrichment analysis to identify novel QIs for the precentrifugation delay in human serum. (A) Volcano plot of serum metabolites visualizes the effects of prolonged blood incubation after 2 h before centrifugation. Red (increase) and green (decrease) circles indicate altered metabolite levels (discovery sample, n = 10). A paired t-test was performed, and we applied addressed multiple testing by deriving q-values to control the positive false discovery rate. (B) Represents the pathway enrichment analysis using the altered metabolites of (A). Node size represents the impact of the altered metabolites on the pathways and node color visualizes the q-value. (D) ChemRICH set enrichment statistics plot. Each node reflects an altered cluster of metabolites (p-value ≤ 0.05). Node sizes represent the total number of metabolites in each cluster set. Node color scale shows the proportion of increased (red) or decreased (blue) compounds in serum samples with a prolonged blood incubation of 2 h compared to the control at 0.5 h. Purple color nodes include both increased and decreased metabolites. (C) The pathway and (E) cluster occupancy visualize the ratio (in %) between increased and decreased metabolites within each altered pathway/cluster of (A).
Metabolites 11 00638 g001
Figure 2. Chemical similarity and pathway enrichment analysis to identify novel QIs for the precentrifugation delay in human EDTA plasma. (A) Volcano plot of serum metabolites visualizes the effects of prolonged blood incubation after 2 h before centrifugation. Red (increase) and green (decrease) circles indicate altered metabolite levels (discovery sample, n = 10). A paired t-test was performed, and we applied addressed multiple testing by deriving q-values to control the positive false discovery rate. (B) Represents the pathway enrichment analysis using the altered metabolites of (A). Node size represents the impact of the altered metabolites on the pathways and node color visualizes the q-value. (D) ChemRICH set enrichment statistics plot. Each node reflects an altered cluster of metabolites (p-value ≤ 0.05). Node sizes represent the total number of metabolites in each cluster set. Node color scale shows the proportion of increased (red) or decreased (blue) compounds in EDTA plasma samples with a prolonged blood incubation of 2 h compared to the control at 0.5 h. Purple color nodes include both increased and decreased metabolites. (C) The pathway and (E) cluster occupancy visualize the ratio (in %) between increased and decreased metabolites within each altered pathway/cluster of (A).
Figure 2. Chemical similarity and pathway enrichment analysis to identify novel QIs for the precentrifugation delay in human EDTA plasma. (A) Volcano plot of serum metabolites visualizes the effects of prolonged blood incubation after 2 h before centrifugation. Red (increase) and green (decrease) circles indicate altered metabolite levels (discovery sample, n = 10). A paired t-test was performed, and we applied addressed multiple testing by deriving q-values to control the positive false discovery rate. (B) Represents the pathway enrichment analysis using the altered metabolites of (A). Node size represents the impact of the altered metabolites on the pathways and node color visualizes the q-value. (D) ChemRICH set enrichment statistics plot. Each node reflects an altered cluster of metabolites (p-value ≤ 0.05). Node sizes represent the total number of metabolites in each cluster set. Node color scale shows the proportion of increased (red) or decreased (blue) compounds in EDTA plasma samples with a prolonged blood incubation of 2 h compared to the control at 0.5 h. Purple color nodes include both increased and decreased metabolites. (C) The pathway and (E) cluster occupancy visualize the ratio (in %) between increased and decreased metabolites within each altered pathway/cluster of (A).
Metabolites 11 00638 g002
Figure 3. Impact of precentrifugation delay on selected purine and pyrimidine derivatives and their corresponding thresholds to verify selected quality indicators in human serum. Shown is the variation of the concentration of selected purine and pyrimidine derivatives in the serum of healthy volunteers after prolonged incubation for 0.5, 1, 2 and 4 h at room temperature (validation sample 1, n = 10). Additionally, diagnostic thresholds were calculated based on the identified pathway connections and the measured concentrations for the selected purine and pyrimidine derivatives. Asterisks indicate differences between the time points as indicated by the respective brackets (* p < 0.05, ** p < 0.01, *** p < 0.001). p-Values were calculated using linear mixed models or paired t-tests and were corrected for multiple testing using Bonferroni post-hoc test.
Figure 3. Impact of precentrifugation delay on selected purine and pyrimidine derivatives and their corresponding thresholds to verify selected quality indicators in human serum. Shown is the variation of the concentration of selected purine and pyrimidine derivatives in the serum of healthy volunteers after prolonged incubation for 0.5, 1, 2 and 4 h at room temperature (validation sample 1, n = 10). Additionally, diagnostic thresholds were calculated based on the identified pathway connections and the measured concentrations for the selected purine and pyrimidine derivatives. Asterisks indicate differences between the time points as indicated by the respective brackets (* p < 0.05, ** p < 0.01, *** p < 0.001). p-Values were calculated using linear mixed models or paired t-tests and were corrected for multiple testing using Bonferroni post-hoc test.
Metabolites 11 00638 g003
Figure 4. Influence of prolonged precentrifugation delay on selected eicosanoids in human serum. Boxplots visualize the variation of the concentration of selected eicosanoids in human serum of healthy volunteers after prolonged incubation for 0.5, 1, 2 and 4 h at room temperature (validation sample 1, n = 10) before centrifugation. Asterisks indicate differences between the time points as indicated by the respective brackets (* p < 0.05, *** p < 0.001). p-Values were calculated using linear mixed models or paired t-tests and were corrected for multiple testing using Bonferroni post-hoc test.
Figure 4. Influence of prolonged precentrifugation delay on selected eicosanoids in human serum. Boxplots visualize the variation of the concentration of selected eicosanoids in human serum of healthy volunteers after prolonged incubation for 0.5, 1, 2 and 4 h at room temperature (validation sample 1, n = 10) before centrifugation. Asterisks indicate differences between the time points as indicated by the respective brackets (* p < 0.05, *** p < 0.001). p-Values were calculated using linear mixed models or paired t-tests and were corrected for multiple testing using Bonferroni post-hoc test.
Metabolites 11 00638 g004
Figure 5. Impact of precentrifugation delay on arginine and ornithine in human serum and EDTA plasma. Boxplots display the variation of the concentration of arginine and ornithine in serum and EDTA plasma (validation sample 2, n = 23) of healthy volunteers after a prolonged incubation for 0, 0.5, 1, 2 and 4 h at room temperature before centrifugation. Additionally, diagnostic thresholds were calculated based on the identified pathway connections and the measured concentrations for both amino acids. Asterisks indicate differences between the time points as indicated by the respective brackets (*** p < 0.001). p-Values were calculated using linear mixed models or paired t-tests and corrected for multiple testing using Bonferroni post-hoc test.
Figure 5. Impact of precentrifugation delay on arginine and ornithine in human serum and EDTA plasma. Boxplots display the variation of the concentration of arginine and ornithine in serum and EDTA plasma (validation sample 2, n = 23) of healthy volunteers after a prolonged incubation for 0, 0.5, 1, 2 and 4 h at room temperature before centrifugation. Additionally, diagnostic thresholds were calculated based on the identified pathway connections and the measured concentrations for both amino acids. Asterisks indicate differences between the time points as indicated by the respective brackets (*** p < 0.001). p-Values were calculated using linear mixed models or paired t-tests and corrected for multiple testing using Bonferroni post-hoc test.
Metabolites 11 00638 g005
Table 1. Accuracy of all potential QIs for discrimination of TTC ≤ 1 h and 2 h in validation samples of healthy volunteers.
Table 1. Accuracy of all potential QIs for discrimination of TTC ≤ 1 h and 2 h in validation samples of healthy volunteers.
Potential QIsMedian Opt. CutoffSpecificity n+/n > 1 hSensitivity 1-(n-/n) < 1 hMedian Opt. CutoffSpecificity n+/n > 2 hSensitivity 1-(n-/n) < 2 h
Purines and Pyrimidines in serum of validation sample 1
Guanosine #0.50 µmol/L100%80%0.28 µmol/L100%83%
Hypoxanthine32 µmol/L75%83%37 µmol/L86%80%
Xanthine1.5 µmol/L75%67%1.6 µmol/L80%67%
Inosine #1.25 µmol/L86%75%1.13 µmol/L100%67%
Dihydroorotic acid32 µmol/L83%75%38 µmol/L94%75%
XG-ratio0.6180%80%0.9389%75%
HI-ratio1.4183%86%1.5778%100%
HG-ratio2.06100%80%2.2588%75%
XI-ratio0.1478%83%0.2473%100%
Eicosanoids in serum of validation sample 1
Tetranor-12(S)-HETE---0.45 ng/mL88%90%
12(S)-HEPE---8.17 ng/mL100%100%
(±)13-HODE---6.83 ng/mL83%100%
(±)9-HODE---5.31 ng/mL62%67%
15(S)-HETE---3.36 ng/mL100%100%
11-HETE---3.84 ng/mL86%75%
8(S)-HETE---2.22 ng/mL100%100%
12-HETE---15.14 ng/mL100%100%
5-HETE---2.63 ng/mL60%60%
12-oxo-ETE---2.62 ng/mL100%100%
Amino acids in serum of validation sample 2
Arginine #86.15 µmol/L64%67%83.4 µmol/L70%65%
Ornithine91.24 µmol/L67%69%106.31 µmol/L74%67%
OA-ratio0.0178%77%0.0575%75%
Purines and Pyrimidines in EDTA plasma of validation sample 1
Hypoxanthine12 µmol/L50%71%11.5 µmol/L80%75%
Inosine #0.5 µmol/L100%33%0.5 µmol/L89%0%
Dihydroorotic acid12.2 µmol/L50%75%12.65 µmol/L67%71%
HI-ratio1.5440%75%1.5167%80%
Eicosanoids in EDTA plasma of validation sample 1
(±)13-HODE---5.45 ng/mL78%67%
(±)9-HODE---2.49 ng/mL88%67%
5-HETE---0.41 ng/mL67%50%
Amino acids in EDTA plasma of validation sample 2
Arginine #51.24 µmol/L73%69%46.62 µmol/L75%65%
Ornithine87.94 µmol/L71%67%110.38 µmol/L88%67%
OA-ratio0.2385%73%0.378%83%
The table displays the calculated median optimal cutoffs for all identified potential QIs in validation sample 1 and 2 in serum and EDTA plasma. # Cutoff values for arginine, inosine and guanosine (concentration decreases over time) are calculated for a TTC ≥ 1 or 2 h, respectively. QIs with the highest sensitivity and specificity are highlighted in bold.
Table 2. Accuracy of all potential QIs for discrimination of TTC ≤ 1 h and 2 h in validation samples of patients with rheumatic and cardiovascular diseases.
Table 2. Accuracy of all potential QIs for discrimination of TTC ≤ 1 h and 2 h in validation samples of patients with rheumatic and cardiovascular diseases.
Potential QIsMedian Opt. CutoffSpecificity n+/n > 1 hSensitivity 1-(n-/n) < 1 hMedian Opt. CutoffSpecificity n+/n > 2 hSensitivity 1-(n-/n) < 2 h
Purines and Pyrimidines in serum of validation sample 3
Guanosine # 0.50 µmol/L43% (6/14)96% (25/26)0.28 µmol/L65% (15/23)71% (12/17)
Hypoxanthine32 µmol/L86% (12/14)65% (17/26)37 µmol/L87% (20/23)47% (8/17)
Xanthine1.5 µmol/L7% (1/14)85% (22/26)1.6 µmol/L22% (5/23)82% (15/26)
Inosine #1.25 µmol/L57% (8/14)88% (23/26)1.13 µmol/L43% (10/23)88% (15/17)
Dihydroorotic acid32 µmol/L79% (11/14)69% (18/26)38 µmol/L91% (21/23)41% (7/17)
XG-ratio0.6136% (5/14)96% (25/26)0.9370% (16/23)71% (12/17)
HI-ratio0.1429% (4/14)92% (24/26)0.2430% (7/23)94% (16/17)
HG-ratio1.4179% (11/14)85% (22/26)1.5770% (16/23)82% (14/17)
XI-ratio2.06100% (14/14)73% (19/26)2.25100% (23/23)35% (6/17)
Eicosanoids in serum of validation sample 3
Tetranor-12(S)-HETE---0.45 ng/mL92% (12/13)100% (7/7)
12(S)-HEPE---8.17 ng/mL100% (13/13)86% (6/7)
(±)13-HODE---6.83 ng/mL62% (8/13)71% (5/7)
(±)9-HODE---5.31 ng/mL46% (6/13)43% (3/7)
15(S)-HETE---3.36 ng/mL85% (11/13)71% (5/7)
11-HETE---3.84 ng/mL77% (10/13)14% (1/7)
8(S)-HETE---2.22 ng/mL100% (13/13)86% (6/7)
12-HETE---15.14 ng/mL92% (12/13)86% (6/7)
5-HETE---2.63 ng/mL69% (9/13)43% (3/7)
12-oxo-ETE---2.62 ng/mL100% (13/13)86% (6/7)
Amino acids in serum of validation sample 4
Arginine #86.15 µmol/L54% (14/26)61% (27/44)83.4 µmol/L51% (24/47)52% (12/23)
Ornithine91.24 µmol/L54% (10/26)61% (27/44)106.31 µmol/L72% (34/47)48% (11/23)
OA-ratio0.0146% (12/26)73% (32/44)0.0551% (24/47)70% (16/23)
Purines and Pyrimidines in EDTA plasma of validation sample 3
Hypoxanthine12 µmol/L62% (8/13)89% (24/27)11.5 µmol/L32% (9/28)83% (10/12)
Inosine #0.5 µmol/L23% (3/13)89% (24/27)0.5 µmol/L11% (3/28)75% (9/12)
Dihydroorotic acid12.2 µmol/L62% (8/13)89% (24/27)12.65 µmol/L36% (10/28)83% (10/12)
HI-ratio1.5485% (11/13)81% (22/27)1.5143% (12/28)75% (9/12)
Eicosanoids in EDTA plasma of validation sample 3
(±)13-HODE---5.45 ng/mL62% (8/13)43% (4/7)
(±)9-HODE---2.49 ng/mL38% (5/13)29% (2/7)
5-HETE---0.41 ng/mL54% (7/13)43% (3/7)
Amino acids in EDTA plasma of validation sample 4
Arginine #51.24 µmol/L60% (3/5)75% (33/44)46.62 µmol/L48% (11/23)81% (21/26)
Ornithine87.94 µmol/L80% (4/5)68% (30/44)110.38 µmol/L70% (16/23)38% (10/26)
OA-ratio0.2360% (3/5)84% (37/44)0.343% (10/23)69% (18/26)
The table displays the calculated median optimal cutoffs for all identified potential QIs in validation sample 1 and 2 in serum and EDTA plasma and applies them in validation sample 3 and 4. # Cutoff values for arginine, inosine and guanosine (concentration decreases over time) are calculated and applied for a TTC ≥ 1 or 2 h, respectively. n+, measurements with a positive test result; n-, measurements with a negative test result; n > 1 h, 2 h measurements with a TTC > 1 or 2 h; n < 1 h, 2 h, measurements with a TTC < 1 or 2 h; QIs with the highest sensitivity and specificity are highlighted in bold.
Table 3. Proposed QIs to be considered for further investigation in serum and EDTA plasma.
Table 3. Proposed QIs to be considered for further investigation in serum and EDTA plasma.
Potential QIs in SerumMedian Opt. CutoffSpecificity
n+/n > 1 h
Sensitivity
1-(n-/n) < 1 h
HG-ratio1.4179% (11/14)85% (22/26)
XI-ratio2.06100% (14/14)73% (19/26)
Specificity
n+/n > 2 h
Sensitivity
1-(n-/n) < 2 h
Tetranor-12(S)-HETE0.45 ng/mL92% (12/13)100% (7/7)
12(S)-HEPE8.17 ng/mL100% (13/13)86% (6/7)
8(S)-HETE2.22 ng/mL100% (13/13)86% (6/7)
12-HETE15.14 ng/mL92% (12/13)86% (6/7)
12-oxo-ETE2.62 ng/mL100% (13/13)86% (6/7)
Potential Qis in EDTA PlasmaMedian Opt. CutoffSpecificity
n+/n > 1 h
Sensitivity
1-(n-/n) < 1 h
HI-ratio1.5485% (11/13)81% (22/27)
The table displays the calculated median optimal cutoffs for all considered QIs in validation sample 1 and 2 in serum and EDTA plasma and applies them in validation sample 3 and 4. n+, measurements with a positive test result; n-, measurements with a negative test result; n > 1 h, 2 h measurements with a TTC > 1 or 2 h; n < 1 h, 2 h, measurements with a TTC < 1 or 2 h.
Table 4. Demographic and clinical data of study cohorts.
Table 4. Demographic and clinical data of study cohorts.
CharacteristicsData
Discovery sample 1
No. of healthy volunteers10
Age (years) median, range33 (28–38)
Sex male/female10/0
Fasting before blood draw (12 h)yes
Non-smokersyes
TimepointsTTC 0.5 and 2 h
Validation sample 1
No. of healthy volunteers10
Age (years) median, range29 (21–54)
Sex male/female7/3
Fasting before blood draw (12 h)yes
Non-smokersyes
TimepointsTTC 0 (only EDTA plasma), 0.5, 1, 2 and 4 h
Validation sample 2
No. of healthy volunteers23
Age (years) median, range43 (25–60)
Sex male/female6/17
Fasting before blood draw (12 h)no
Non-smokersyes
TimepointsTTC 0.5, 1, 2 and 4 h
Validation sample 3
No. of rheumatologic patients20
Age (years) median, range70.5 (21–82)
Sex male/female14/6
Diseaserheumatic
Fasting before blood draw (12 h)no
Validation sample 3
No. of cardilogic patients20
Age (years) median, range64.5 (32–82)
Sex male/female10/10
Diseasecardiovascular
Fasting before blood draw (12 h)no
Validation sample 4
No. of cardiologic patients40
Age (years) median, range73 (48–89)
Sex male/female20/20
Diseasecardiovascular
Fasting before blood draw (12 h)no
Validation sample 4SerumEDTA plasma
No. of rheumatologic patients3049
Age (years) median, range63 (32–88)65 (20–90)
Sex male/female14/1623/26
Diseaserheumaticrheumatic
Fasting before blood draw (12 h)nono
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Heiling, S.; Knutti, N.; Scherr, F.; Geiger, J.; Weikert, J.; Rose, M.; Jahns, R.; Ceglarek, U.; Scherag, A.; Kiehntopf, M. Metabolite Ratios as Quality Indicators for Pre-Analytical Variation in Serum and EDTA Plasma. Metabolites 2021, 11, 638. https://doi.org/10.3390/metabo11090638

AMA Style

Heiling S, Knutti N, Scherr F, Geiger J, Weikert J, Rose M, Jahns R, Ceglarek U, Scherag A, Kiehntopf M. Metabolite Ratios as Quality Indicators for Pre-Analytical Variation in Serum and EDTA Plasma. Metabolites. 2021; 11(9):638. https://doi.org/10.3390/metabo11090638

Chicago/Turabian Style

Heiling, Sven, Nadine Knutti, Franziska Scherr, Jörg Geiger, Juliane Weikert, Michael Rose, Roland Jahns, Uta Ceglarek, André Scherag, and Michael Kiehntopf. 2021. "Metabolite Ratios as Quality Indicators for Pre-Analytical Variation in Serum and EDTA Plasma" Metabolites 11, no. 9: 638. https://doi.org/10.3390/metabo11090638

APA Style

Heiling, S., Knutti, N., Scherr, F., Geiger, J., Weikert, J., Rose, M., Jahns, R., Ceglarek, U., Scherag, A., & Kiehntopf, M. (2021). Metabolite Ratios as Quality Indicators for Pre-Analytical Variation in Serum and EDTA Plasma. Metabolites, 11(9), 638. https://doi.org/10.3390/metabo11090638

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