Next Article in Journal
Night Shift Work Affects Urine Metabolite Profiles of Nurses with Early Chronotype
Previous Article in Journal
A Multi-Methodological Protocol to Characterize PDO Olive Oils
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Validation and Automation of a High-Throughput Multitargeted Method for Semiquantification of Endogenous Metabolites from Different Biological Matrices Using Tandem Mass Spectrometry

1
Metabolomics Unit, Institute for Molecular Medicine Finland (FIMM), University of Helsinki, HiLIFE, Tukholmankatu 8, Biomedicum 2U, 00290 Helsinki, Finland
2
Roche Diagnostics GmbH, Nonnenwald 2, 82377 Penzberg, Germany
3
Computational Systems Medicine group, University of Helsinki, 00290 Helsinki, Finland
4
Technical Research Center of Finland, P.O. Box 1000, 02044 Espoo, Finland
5
Network Pharmacology for Precision Medicine Group, University of Helsinki, 00290 Helsinki, Finland
6
Finnish Customs Laboratory, Tekniikantie 13, 02150 Espoo, Finland
*
Author to whom correspondence should be addressed.
Metabolites 2018, 8(3), 44; https://doi.org/10.3390/metabo8030044
Submission received: 20 June 2018 / Revised: 27 July 2018 / Accepted: 3 August 2018 / Published: 5 August 2018

Abstract

:
The use of metabolomics profiling to understand the metabolism under different physiological states has increased in recent years, which created the need for robust analytical platforms. Here, we present a validated method for targeted and semiquantitative analysis of 102 polar metabolites that cover major metabolic pathways from 24 classes in a single 17.5-min assay. The method has been optimized for a wide range of biological matrices from various organisms, and involves automated sample preparation and data processing using an inhouse developed R-package. To ensure reliability, the method was validated for accuracy, precision, selectivity, specificity, linearity, recovery, and stability according to European Medicines Agency guidelines. We demonstrated an excellent repeatability of retention times (CV < 4%), calibration curves (R2 ≥ 0.980) in their respective wide dynamic concentration ranges (CV < 3%), and concentrations (CV < 25%) of quality control samples interspersed within 25 batches analyzed over a period of one year. The robustness was demonstrated through a high correlation between metabolite concentrations measured using our method and the NIST reference values (R2 = 0.967), including cross-platform comparability against the BIOCRATES AbsoluteIDQp180 kit (R2 = 0.975) and NMR analyses (R2 = 0.884). We have shown that our method can be successfully applied in many biomedical research fields and clinical trials, including epidemiological studies for biomarker discovery. In summary, a thorough validation demonstrated that our method is reproducible, robust, reliable, and suitable for metabolomics studies.

1. Introduction

Metabolomics has a great influence on many disciplines, as metabolites are intermediates or end products of cellular functions. Hence, metabolomics can be used as a powerful tool in understanding, diagnosing, and managing different pathophysiological conditions. It is therefore essential to be able to identify and measure metabolites from different biological matrices [1]. Although global metabolomics has been widely used in discovery studies for understanding cellular responses to normal and abnormal biological conditions, targeted metabolomics has more advantages for addressing biological questions in a more hypothesis-driven manner than global untargeted metabolomics [2]. Furthermore, targeted metabolomics can quantify metabolites that are low in abundance, which are difficult to assess using an untargeted approach.
An appropriate sample pretreatment is required to obtain reproducible and high-quality quantitative data in targeted metabolomics. However, metabolites are present in a wide dynamic range with great diversity in physicochemical properties in the biological matrices [3]. Recent advancements in extraction techniques and automated approaches for sample preparation have partially satisfied the demands of targeted metabolomics. However, there are still many outstanding challenges, such as the matrix effect and laboratory-to-laboratory variations associated with sample preparation. Hence, standardization of sample preparation is a fundamental requirement in metabolomics studies [4].
Secondly, a robust analytical methodology is required for an accurate quantification of metabolites with good reproducibility over an extended period of time [5]. Molecular diversity is a major problem that hinders the separation of all preselected metabolites in a single chromatographic run and the detection of all separated metabolites with a minimum technical variation [6]. Tandem mass spectrometry (MS) is a technique used predominantly due to its high sensitivity and high throughput for the detection of metabolites. A combination of MS and separation techniques is used to increase the sensitivity and reliability of analytical methods for the analysis of metabolites from complex biological matrices [7]. In addition, the latest developments in triple quadrupole instrumentation strengthened the possibilities to develop multianalyte methods in a single injection that yield reliable and quantitative data [8]. Several analytical methods have also been developed for semiquantitative measurement of large number of metabolites in a single run [9,10,11,12,13,14,15].
Even though liquid chromatography-MS (LC-MS) is a method of choice in targeted metabolomics, obtaining an accurate quantification and long-term data reproducibility remains an analytical challenge. This is due to limitations such as matrix effect, MS performance drift, and LC column contamination and aging [16]. Finally, although instrument vendor software for data processing provides some crucial functions, such as peak integration (rendering the data from high-throughput metabolomics experiments into numerical values that represent metabolite concentrations), the lack of automation in downstream data processing and quality control remains a major bottleneck in high-throughput analyses. Thus, an efficient pipeline is necessary to enable rapid, accurate, and standardized processing of these data. Such a pipeline should facilitate automated analyses, while at the same time allowing the user to fine-tune the parameters for accurate data processing. Taken together, there is a need for development and validation of standardized, robust, and quantitative methods for large-scale targeted metabolomics studies in a high-throughput manner to minimize the bias associated with sample preparation and the analytical technique used.
We have previously developed a robust, reproducible, and high-throughput targeted method for measuring 102 polar metabolites from various biological classes semiquantitatively in a single injection [17,18]. Furthermore, the method has been used in many biomedical and clinical studies for biomarker discovery [17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39]. The metabolites were selected according to the following criteria: have important roles in many biological processes, are known biomarkers in several diseases and technical feasibility for developing an analytical method that covers all the metabolites in a single assay. The selected metabolites come from 24 different classes (Table 1), covering a wide range of metabolic pathways. We created an inhouse metabolite database by manually curating all the available information (i.e., names, HMDB, PUBCHEM, KEGG Ids, chemical properties, reported normal and abnormal concentration ranges, links to their structures) from the Human Metabolome Database (HMDB). The selected metabolites were separated by Hydrophilic Interaction Liquid Chromatography (HILIC) and measured using triple quadrupole tandem MS. A detailed method description is provided in the supplementary text.
The primary objective of this work is to show the robustness of our previously developed analytical method through a thorough validation according to European Medicines Agency (EMA) guidelines. We also demonstrate the automation of tedious and manual data-processing tasks in high-throughput metabolomics analyses using an inhouse developed R-package. The R-package automates various corrections and normalization steps to convert the raw peak area data to molecular concentrations for each compound in each sample and also provides quality evaluation of the data and reduces the manual workload significantly.

2. Materials and Methods

2.1. Chemicals and Reagents

All metabolite standards were purchased from Sigma-Aldrich (St. Louis, MO, USA). Internal standards were ordered from Cambridge Isotope Laboratory. Inc. (Tewksbury, MA, USA). LC-MS-grade solvents, 2-proponol, acetonitrile, and methanol (HiPerSolv) were obtained from VWR International (Helsinki, Finland). Analytical-grade chemicals (formic acid, ammonium formate, and ammonium hydroxide) were obtained from Sigma-Aldrich. Deionized water (18 MΩ·cm at 25 °C) used for solution preparation was made using a Milli-Q water purification system (Bamstead EASYpure RoDi ultrapure water purification system, Thermo scientific, Waltham, OH, USA). Mouse tissues, including heart, liver, brain, spleen, and muscles were obtained from Innovative Research Laboratory (Novi, MI, USA). The whole blood, from which serum was prepared during method optimization and validation, was obtained from the Finnish Red Cross blood service (Helsinki, Finland). Cell samples were provided by our research collaborators. NIST Standard reference material (SRM) 1950 plasma was purchased from Sigma-Aldrich (Gillingham, UK).

2.2. Metabolite Extraction Protocol and Instrumentation

All metabolites were extracted, separated with HILIC (Acquity BEH amide, 2.1 × 100 mm, 1.7 µ), and analyzed with a Waters Xevo TQ-S triple quadrupole mass spectrometer using our previously published protocol [18]. The protocol for tissues and adherent cells was optimized for better recovery and chromatography and to cover a wide range of tissue and cell types with a single protocol. For tissue sample extractions, 90/10% ACN/H2O + 1% formic acid was used instead of 80/20% ACN/H2O + 1% formic acid during the second step of extraction. Additionally, during cell pellet sample extraction, 80/20% ACN/H2O + 1% formic acid was replaced with 90/10% ACN/H2O + 1% formic acid. After optimization, we used the tissue protocol for analysis of various biological matrices, such as heart, liver, placenta, brain, muscles, spleen tissues, dental carries, and Drosophila larvae (weight 10–20 ± 5 mg), C. elegans (2000–4000 worms/sample), dried blood spots (12–15 punched spots), and fecal samples (weight 20–50 ± 5 mg). The cell pellet protocol was used for all types of adherent cells (around 1 million cells/sample) and E. coli (3-day old bacterial lawns) and S. cerevisiae (OD600nm ~ 2.2) samples. The biofluid protocol was used for all types of biofluids, such as blood, plasma, serum, cell culture supernatant, CSF, and urine (sample volume 50–100 µL).

2.3. Method Validation

Validation of the method was performed to verify various parameters and the reliability of the developed method for the analysis of a large number of samples. The method was validated according to EMA guidelines for bioanalytical method validation in terms of selectivity, specificity, linearity, accuracy, precision, extraction recovery, matrix effect, and stability [40]. In addition, we used pooled healthy human serum samples as internal quality control (QC) samples in all studies to correctly signal drift during sample runs and to improve confidence in the statistical data. QC samples at high, medium, and low (for serum) or high and low concentration levels (for tissues) were prepared by spiking a mixed standard solution in their respective homogenized biological matrices to perform all the method validation experiments. We performed the validation for commonly used biological samples in metabolomics analyses, such as biofluid (serum), tissue (liver, brain and spleen), and cell samples. An aqueous calibration curve was used to calculate the concentration values during the method validation. The instrument performance for response reproducibility and sensitivity was always verified by 6 consecutive injections of a medium concentration solution at the start of any experiment.

2.3.1. Selectivity and Specificity

The selectivity and specificity for each metabolite were investigated using serum-spiked samples (N = 6) with known amount of standard. Chromatographic interferences from other endogenous compounds of the biological matrix at the retention time of the target analyte for a particular metabolite were verified. The chromatographic peaks from spiked samples were compared with the standards by the retention times and if required from their respective MRM spectra.

2.3.2. Linearity, Accuracy, and Precision

To assess the linearity, accuracy, and precision, six replicates of spiked QC samples at high, medium, and low concentrations along with calibration curve standards were injected on three separate days. Calibration curve standards of 11 points were prepared via serial dilution and each calibration curve has over 1000-fold dynamic concentration range. The curve was plotted by using the peak area response ratios (standard/labeled standard) versus the concentrations of the individual metabolites. Each calibration curve was statistically evaluated and constructed using appropriate regression models, weighing factors and transformations. The accuracy was calculated as the measured value divided by the nominal value at each concentration level of the calibration curve standards in all three batches. Inter- and intrabatch variability was calculated by measuring coefficient of variation (%CV) at each QC concentration level.

2.3.3. Recovery and Matrix Effect

The recovery efficiencies for each analyzed metabolite were determined by comparing analytical results from QC samples spiked with a standard mixture before and after extraction using different concentrations. The spike concentrations covered the calibration range. The matrix effect (percentage of ion suppression or enhancement of the MS signal) was determined by comparing the analytical response of the QC samples that were spiked after extraction with the analytical response of aqueous spiked samples (diluent spiked with respective concentrations of QCs). Since there were endogenous metabolites, we subtracted the endogenous concentrations from the samples that were spiked. This experiment was performed using 6 QC replicates.

2.3.4. Stability of the Metabolites

Wet extract, freeze-thaw, and stock solution stability for all metabolites were determined to check the integrity of the analytes in solvents and in QC samples at different conditions. To determine the wet extract stability, six replicates of extracted QC samples were kept in the autosampler at 5 °C. The same samples in the same sequence were reinjected with freshly extracted QC samples and the results were compared.
Freeze-thaw stability was evaluated up to three cycles by freezing and thawing the spiked QC samples stored at –80 °C (for 12–16 h) and comparing the concentrations against the freshly thawed and spiked QC samples.
Long-term stock solution stability for metabolite stock solutions and intermediate solutions were checked by comparing the mean peak area of freshly prepared solutions with stored solutions at 4 °C. All stability experiments were performed with 6 replicates of QCs.

2.3.5. Carryover

The carryover was evaluated by injecting the highest standard concentration (ULOQ) of the metabolites in the calibration curve followed by a series of blank injections and lowest standard concentration (LLOQ). The blank samples were evaluated for any signal at the retention time of particular metabolites and signal intensities of the blank samples were compared with the LLOQ samples. The acceptance criteria for carryover was set at 20% of the peak area corresponding to the LLOQ level as per the EMA guidelines for bioanalysis.

2.3.6. QC Samples

Internal QC samples were prepared after separating serum from pooled healthy human blood samples. A volume of 350 µL of serum was aliquoted and stored at –80 °C after providing a lot number and QC number. The concentration of QC samples that were incorporated in batches during the metabolomics studies was calculated for all the metabolites along with the experimental samples. Average concentrations (μmol/L) and %CV of the QC samples were calculated for each metabolite. The data were saved along with QC lot numbers, batch name, and run date. The QC data were collected from six different lots for a period of 5.5 years (N = 539 replicates). An internal QC database has been maintained and used for quality checks.

2.3.7. Comparison with Reference Material

To evaluate the performance of our semiquantitative method, commercially available standard-reference plasma (NIST SRM 1950) [41] was analyzed using our method (N = 8 replicates). The concentration values from the matched 17 metabolites were compared with the given standard reference values.

2.3.8. Cross-Platform Comparison

To further evaluate the robustness and performance of our method, we performed a cross-platform comparison using two completely different analytical platforms: (1) the commercially available AbsoluteIDQ p180 targeted metabolomics assay kit using LC-MS/MS and (2) a nuclear magnetic resonance (NMR) platform. We sent our internal QC samples to the BIOCRATES Life Sciences AG (Innsbruck, Austria) (N = 3 replicates) and to the NMR Metabolomics Laboratory, School of Pharmacy, University of Eastern Finland (Kuopio, Finland) (N = 3 replicates). Our QC samples were extracted and analyzed as described previously for the AbsoluteIDQ p180 kit [42] and for the NMR analysis of small molecules [43]. We compared these results with the results obtained from our method (N = 4–5 replicates).

2.4. Statistical Analyses

To estimate the median concentration values (μmol/L) of the metabolites from the QC samples (N = 539), we fitted a linear mixed model with the MCMCglmm R-package [44] using an expanded parameter formulation and default settings. Observed data were assumed to be log-normally distributed and corrected for the six different QC lots. Credibility intervals (95%) of the median concentration values were computed from 20,000 samples of the posterior distribution. Error bars shown in the scatter plots are 95% confidence intervals, while the coefficient of determination R2 was obtained from the linear regression between the variables. Coefficient of variation (CV) percentages was calculated as a measure of variability. To automate the downstream processing of data produced by the instrument vendor software (TargetLynx, v4.1), we built a data-processing package called “Unlynx” in R statistical programming language. The R-package is available upon request.

2.5. Automated Data Processing

The “Unlynx” package parses the output of TargetLynx software (i.e., raw data containing the concentration values in PPB units) and produces a processed dataset in an Excel spreadsheet after performing a series of preprocessing operations.
The preprocessing steps included the following:
(i)
Molecular weight normalization, in which the ppb values are normalized by the molecular weight of each compound, thereby converting the data from ppb units to µmoles.
(ii)
Process efficiency correction for the semiquantification of metabolites without internal standards.
(iii)
Normalization using dilution factor for specific sample type if dilution was needed.
(iv)
Cell number normalization (for cell samples) to convert the concentration values per million cells.
(v)
Calculation of mean, standard deviation, and relative standard deviation (RSD) of molecular concentrations (resulting from the previous steps) for each phenotypic group.
(vi)
Outlier detection in each phenotypic group; if the concentration value of a compound in a sample is more than one or two standard deviations (SD) away from the mean of the phenotypic group, then it is marked as an outlier in the Excel data set in two different colors.
(vii)
QC check by comparing the RSD of QC samples in the current dataset against the internal database of QC sample RSDs (based on interday RSDs recorded over one year).

3. Results and Discussions

3.1. Extraction Method Optimization

The primary objective of this work was to optimize and validate our previously published protocol for different types of biological matrices. For tissue samples (placenta, liver, heart, brain, spleen, and muscles), the sample volumes of the tissues and extraction solvent volumes were optimized to fit the concentrations of most of the metabolites within the linearity of calibration curve for reliable results. We observed that most of the metabolites could be semiquantified within the calibration curve range with 20 ± 5 mg of sample weight.
Furthermore, we optimized the protocol with extraction solvent for tissues and adherent cells. Some of the metabolites (in particular inositol, GABA, asymmetric dimethylarginine, symmetric dimethylarginine, spermidine, ribose-5-phosphate, and orotic acid) had poor separation and irreproducible chromatography. Interference of isobaric compounds with other metabolites was also observed due to poor separation. Thus, different compositions of the extraction solvent were assessed to achieve the acceptable chromatography. We observed that modification of acetonitrile content from 80% to 90% and applying longer equilibration time for the HILIC column yielded acceptable chromatography and also good separation for most of the metabolites.
We also optimized the extraction protocol for different sample types from various organisms (human, mouse, rat, dog, drosophila, nematode, yeast, bacteria), such as tissue types (spleen, pancreas, muscles, adipose, endometrium, testicles, lung), biofluids (plasma, blood, urine, cyst fluid, bile, CSF, saliva), cell types (adherent cells, cell suspension, bone marrow cells, extracellular vesicles, mitochondrial isolates, E. coli, and S. cerevesiae) and other sample types such as fecal samples, dried blood samples, dental carries, biofilm, drosophila (larvae and whole flies), and C. elegans.
A flowchart describing various stages of our workflow (preanalytical stage, sample analysis, analytical stage, data processing, and automation with quality checks) is shown in Figure 1.

3.2. Method Validation

3.2.1. Selectivity and Specificity

There were no significant interference peaks from the matrix components in their respective retention-time windows, indicating the selectivity of the metabolites in our method. We repeated the injections for five times from all different serum samples and confirmed every peak to have eluted only from the target analyte, indicating that they are specific to their corresponding MRM transitions. The chromatograms for all the 102 metabolites of QC sample were given in Supplementary Figure S1.

3.2.2. Linearity, Accuracy, and Precision

To cover a broad concentration range, when fitting the heteroscedastic calibration data, linear or quadratic regression models with appropriate weighing factors (1/x) and transformations (log-log) were used [45]. The coefficient of determination (R2) value for each metabolite was greater than 0.980 at their respective concentration range, except for some metabolites such as aspartate, uracil, 2-deoxyuridine sucrose, and chenodeoxycholic acid (likely due to their broad peak shapes and poor recovery at lower concentration, Table S1).
The concentration precision for QC samples was calculated by measuring %CV at high, medium, and low concentration level of QCs (N = 6 replicates). In general, intra- and interday precision (CV) values were within 15% for 94 metabolites except acetoacetic acid, folic acid, sucrose, homoserine, 2-deoxyuridine, and cholic acid at high concentration (Figure 2). However, at low concentrations, more than 20% CV was observed for NAD and myo-inositol. This might be due to low recoveries of these compounds. Hence, these compounds cannot be measured reliably with this method even semiquantitatively.

3.2.3. Recovery and Matrix Effect

For 85–90 metabolites, recoveries were found to be between 50 to 120% with good repeatability at all three concentrations levels (low, medium, and high) in both biofluid (serum) and tissues (brain, liver, and spleen). However, compounds such as UDP-glucose, IMP, cGMP, d-ribose 5-phosphate, NAD, AMP, homocysteine, carnosine, and glutathione had recoveries less than 30% in serum. However, CV of recoveries at low, medium, and high concentration levels was within 25% except for cGMP in serum. Metabolites such as 1-methylhistamine, aspartate, glutamine, adenosine, and glutathione had recoveries over 120%. Histidine, ornithine, cystathionine, 3-OH-dl-kynurenine, carnosine, AMP, NAD, cGMP, IMP, and UDP-glucose had less than 30% recovery in some tissues types, indicating matrix effect or degradation (Figure 3). However, CV of repeatability at every concentration level was within 15% except for 1-methylhistamine, aspartate, glutamine, glutathione, NAD, and UDP-glucose in some tissues.
The matrix effect values were observed to be within the range of 0.6 to 1.8 (below 1 indicates ion suppression and above 1 indicates ion enhancement) for the metabolites in serum and tissues. The challenge of the matrix effect can be overcome by having individual isotope-labeled internal standards for each individual compound for true quantification. However, this is not practically possible for high-throughput metabolomics analyses. This is due to high costs and also because not all internal standards are commercially available. In our method, we selected 12 labeled internal standards (Table S1), which represent chemically similar classes for optimal correction. This is because the matrix effect was expected to be the same for an analyte and its labeled isotope analogue. The process efficiency percentages were calculated for the metabolites without internal standards. The analyte concentrations determined through the external calibration were divided with the total process efficiency values to correct the concentration values of the analyte in the given biological sample. The repeatability of the matrix effect in terms of CV was also less than 25% for most of the compounds. Reliable measurements are accordingly possible.

3.2.4. Stability

Some of the endogenous metabolites are not stable due to degradation or conversion reactions. Hence, the stabilities of all metabolites were assessed under different conditions. For wet extract stability, approximately 90% of the metabolites were stable (stabilities range between 85 and 115%) for 35 h at 5 °C in the autosampler (Figure 4A).
For freeze and thaw cycle stability, most of the metabolites were stable even after three freeze and thaw cycles, with the exception of cGMP, succinate, glutathione, and homocysteine (stability below 30%, Figure 4B). This information is particularly important for clinical studies, where samples are often thawed once or twice.
To determine the stability of working solutions, we started evaluating the stability from intermediate solutions for all the metabolites. Most of the metabolites were stable for 245 days at intermediate concentration when stored at 4 °C. However, 16 metabolites had low stability at intermediate concentration; stabilities were thus determined at stock-level concentration for these compounds. We observed that the stock solutions were stable for 56 days except taurocholic acid, sucrose, UDP-glucose, and glutamine (Figure 5). Hence, these stock solutions were freshly prepared during the analysis. The stability of internal standards solutions was also assessed and they were stable for one year.

3.2.5. Carryover

In general, for the majority of the analyte MRM channels neither a peak nor any interference in the blank samples was detected after injection of the metabolite standard with high concentration. For compounds such as spermidine, succinate, AMP, and IMP, carryover eluted constantly even after washing but was not significantly high. Other than these compounds, we can conclude that the column, needle, syringe, and seal washes were sufficient to avoid any intersample carryover.

3.2.6. Reproducibility

To ensure good quality of the data, internal QC samples were incorporated to a batch of samples and run after every tenth experimental sample. QC data were collected from 25 different batches that were performed during various metabolomics studies over a period of 1 year. Mean concentrations and %CV values of QC replicates within each batch were calculated for all the 25 batches. We observed that approximately 80 to 85% of the metabolites were always present within 25% of CV values (Figure 6). The higher %CV values for the remaining metabolites could be partially explained by low abundance in human serum, low recovery, or poor chromatography; these were consistently found to be below LLOQ within the 25 batches.
In addition, %CV values for retention times and R2 values of calibration curves for each metabolite in all the 25 batches were calculated to verify the reproducibility. Based on these results, the repeatability was excellent except for a few compounds over a period of 1 year. No drifting effect for the retention times (CV < 4%) was observed, and excellent reproducibility was observed for R2 values of calibration curves (CV < 3%) (Figure 7). On the basis of these results, our method can be considered accurate, reliable, and reproducible.

3.2.7. Quality Management

To obtain reproducible and accurate data, we set up a strict quality management and electronic lab notebook system. To reduce the bias from sample analysis, we always double-randomized the samples (i.e., one before the sample extraction step and one before injecting into the LCMS system across different phenotypes of the samples). For stabilization of response and retention time, we always verified a few runs of highest calibration level 11 before injecting the experimental samples. During the stabilization process, we also verified the chromatography, including peak shape, retention time, and response of all the metabolites. Any significant changes in the intensity, peak shape, retention time, and system pressure were thoroughly investigated and corrected by resolving the problems before injection of experimental samples.
To ensure the integrity of LCMS runs, QC samples were run at every tenth experimental sample and a blank sample at every fifth run during all the metabolomics studies within a batch. Furthermore, chromatography and response of QC samples (including chromatography of some metabolites and IS response variation) and blank runs were always verified after completion of the runs and before starting data processing. In case of any abnormality observed for particular samples, those samples were reinjected or reanalyzed. Only after passing these quality checks we proceeded further to process the data. This included verifying the accuracy of calibration curve standards, chromatography peak integrations, IS response variation, and verifying LLOQ and ULOQ for each sample for all metabolites within a batch. The high-throughput targeted metabolomics workflow is shown in Figure 8.
We collected concentration values (μmol/L) for our QC samples within metabolomics studies conducted over a period of 5.5 years from six different lots (N = 539 replicates). The median values of each metabolite together with a 95% credibility interval are presented in Table 1. These represent a reference level for a population of healthy adult individuals.

3.2.8. Robustness and Cross-Platform Comparison

To verify the performance of our method, we analyzed the NIST standard reference material SRM 1950 plasma. The correlation coefficient for 17 matched metabolites between the given reference values and from our semiquantitative method was 0.967, indicating the high performance of our method (Figure 9A).
Furthermore, we verified cross-platform comparability. This was achieved by comparing metabolite concentrations analyzed using our method against two completely different analytical platforms (BIOCRATES AbsoluteIDQ p180 kit and NMR) in our QC samples. We obtained a high correlation coefficient for matched 38 metabolites measured using BIOCRATES AbsoluteIDQ p180 kit and our method (R2 = 0.975) (Figure 9B) and for matched 22 metabolites measured using NMR and our method (R2 = 0.884) (Figure 9C). These results demonstrate the robustness of our method.

3.3. Automated Data Processing

After the raw data processing using the instrument-coupled software (TargetLynx), an unstructured flat text file containing information such as sample ids, file name, and concentrations (PPB) for all metabolites is generated in a complex format. For example, in a single-batch run of 85 samples, after the data processing, the concentration values for all samples are separately obtained under each metabolite. This means that, if a data matrix of 85 samples × 100 metabolites is desired, each individual concentration value must be copied and pasted in another sheet in a tabular format. Then, the PPB values are to be converted into µm and corrected for process efficiencies and dilution factors (if any) and normalized with tissue weight or cell number depending on the sample type. Apart from these steps, as the samples are analyzed in a randomized manner, the next steps would include rearranging the experimental samples according to phenotypic group and separating the QC samples from the study samples. This manual data processing creates a ready-to-use data matrix for visualization and for downstream statistical analyses, but is tedious, time-consuming, and, more importantly, prone to errors.
To automate the manual data processing, we have implemented a software package “Unlynx” in R statistical language. This package takes the raw data produced by TargetLynx software as input and produces processed data into ready-to-use spreadsheets. We use Unlynx to covert the PPB values to µmol/L, µmol/g, or µmol/million cells by dividing with molecular weight of the respective metabolites and other correction factors (e.g., weight of tissues, number of cells, and dilution factor) appropriate for the sample type. In addition, mean concentrations and %CV of QC samples for all the metabolites are calculated and used to evaluate quality by comparing with the inhouse QC database constructed based on inter-day %CV. Furthermore, Unlynx retrieves LLOQ, ULOQ, and outlier values in each phenotypic group according to one and two SDs, retention time values, and R2 values for each metabolite. With the automated data processing, we reduced a two-day manual workload to a few minutes.
Thus, after processing the data using TargetLynx, we routinely exported all results for the automated data processing using Unlynx. Typically, the data resulting from such automated data processing undergoes more specialized data analyses (such as statistical hypotheses testing, classification, regression, and clustering) aimed at answering specific scientific questions related to the study design.

3.4. Applicability of the Method

We have applied our fully validated analytical methodology in various international and national biomedical research projects, epidemiological studies, clinical studies including dietary interventions, and clinical trials. We have successfully implemented our technology in the following research fields, including but not limited to: mitochondrial metabolism/disorders [17,19,20,21,22,23,24,25], cancer [26,27], bone metabolism [18], endocrinology [28,29], psychiatric disorders [30], inflammatory bowel disease [31], viral infections [32,33,34,35], allergies [36,37], circadian rhythms [38], and pain research [39].

4. Conclusions

We validated our high-throughput targeted and semiquantitative analytical method according to the EMA guidelines for bioanalytical methods. The key features of our method include,
(i)
Optimization: well-characterized protocols for various biological matrices from different organisms enabled to study wide variety of research projects.
(ii)
Accuracy/Precision: the targeted and semiquantitative analysis using 102 external 11-point calibration curves, including 12 labeled internal standards in every analysis, made it possible to compare the data within and between the studies.
(iii)
Quality management: standard operating protocols, good laboratory practices, strict quality-management system, and proper documentation using electronic laboratory notebook enabled to check/retrieve very old data.
(iv)
High-throughput: automated sample preparation and short analysis time (17.5 min), enabled high-throughput capabilities, which is the most desired feature for large-scale analyses.
(v)
Stability: long-term stability studies in stock and intermediate solutions, wet extract, and freeze-thaw stability studies, critical for projects based on clinical/biobank samples.
(vi)
Automation: downstream data processing steps in an automated manner reduces the interbatch variation and human errors, valuable for analyzing population cohorts and in epidemiological studies.
(vii)
Reproducibility: low %CV of concentrations, retention times, correlation coefficient of calibration curves for an extended period of time, very important in metabolomics studies when comparing the data produced at different points of time.
(viii)
Reliability: as shown by the excellent correlation between metabolite concentrations measured using our method and the NIST SRM plasma reference values, our method serves a standardised and reliable platform for metabolomics studies.
(ix)
Robustness: Our results demonstrate an excellent cross-platform comparability with two completely different analytical platforms, a highly desirable criterion in multicenter studies when comparing the data across different laboratories using different instrumentation, protocols and analytical platforms.
(x)
Data sharing: the huge QC sample database of healthy adults (N = 539) collected for six years and shared with the scientific community, provide normal reference values such as those provided in the HMDB database.
The above-mentioned aspects are critical in the metabolomics analysis, as the International Metabolomics Community has been putting a lot of efforts in data standardization, QC, reproducibility, robustness, and data sharing with the goal of moving towards applicability and integration of metabolomics data in precision/personalised medicine. In this context, our method provides a very timely and important contribution to the field. Furthermore, we have widely applied this method in many biomedical research projects, clinical trials, epidemiological studies, and in biomarker discovery.

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-1989/8/3/44/s1, Table S1: List of metabolites showing molecular weight, retention time, linearity of calibration, and compound dependent MS parameters and catalogue numbers. Supplementary text: A detailed method description. Figure S1: Chromatograms for all the 102 metabolites of QC sample.

Author Contributions

J.N. performed all the validation experiments, processed the data, and wrote the manuscript. M.K. and J.N. were involved in extraction protocol optimization. G.P. developed the R-package (Unlynx) for automated data processing. A.P. performed the QC data analysis, prepared the related table, and wrote the corresponding section. V.V. conceived, designed, and supervised the study, and wrote and edited the manuscript. All authors edited the manuscript.

Funding

This research was funded by Biocenter Finland, HiLIFE, University of Helsinki, Finland.

Acknowledgments

The authors would like to thank all the collaborators, especially Anu Wartiovaara for providing the biological samples; Jean-Christophe Yorke, Vasudev Kantae, Kenneth Nazir, and Bharat Gajera for technical assistance; and Janne Wallenius for reviewing the manuscript.

Conflicts of Interest

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

References

  1. Gowda, G.A.N.; Djukovic, D. Overview of mass spectrometry-based metabolomics: Opportunities and challenges. Methods Mol. Biol. 2014, 1198, 3–12. [Google Scholar] [CrossRef] [PubMed]
  2. Dudley, E.; Yousef, M.; Wang, Y.; Griffiths, W.J. Targeted metabolomics and mass spectrometry. Adv. Protein Chem. Struct. Biol. 2010, 80, 45–83. [Google Scholar] [CrossRef] [PubMed]
  3. Raterink, R.-J.; Lindenburg, P.W.; Vreeken, R.J.; Ramautar, R.; Hankemeier, T. Recent developments in sample-pretreatment techniques for mass spectrometry-based metabolomics. TrAC Trends Anal. Chem. 2014, 61, 157–167. [Google Scholar] [CrossRef]
  4. Gong, Z.-G.; Hu, J.; Wu, X.; Xu, Y.-J. The recent developments in sample preparation for mass spectrometry-based metabolomics. Crit. Rev. Anal. Chem. 2017, 47, 325–331. [Google Scholar] [CrossRef] [PubMed]
  5. Dettmer, K.; Aronov, P.A.; Hammock, B.D. Mass spectrometry-based metabolomics. Mass Spectrom. Rev. 2007, 26, 51–78. [Google Scholar] [CrossRef] [PubMed]
  6. Abdel Rahman, A.M.; Pawling, J.; Ryczko, M.; Caudy, A.A.; Dennis, J.W. Targeted metabolomics in cultured cells and tissues by mass spectrometry: Method development and validation. Anal. Chim. Acta 2014, 845, 53–61. [Google Scholar] [CrossRef] [PubMed]
  7. Villas-Bôas, S.G.; Mas, S.; Åkesson, M.; Smedsgaard, J.; Nielsen, J. Mass spectrometry in metabolome analysis. Mass Spectrom. Rev. 2005, 24, 613–646. [Google Scholar] [CrossRef] [PubMed]
  8. Virgiliou, C.; Sampsonidis, I.; Gika, H.G.; Raikos, N.; Theodoridis, G.A. Development and validation of a hilic-ms/ms multitargeted method for metabolomics applications. Electrophoresis 2015, 36, 2215–2225. [Google Scholar] [CrossRef] [PubMed]
  9. 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]
  10. Liu, P.; Qi, C.-B.; Zhu, Q.-F.; Yuan, B.-F.; Feng, Y.-Q. Determination of thiol metabolites in human urine by stable isotope labeling in combination with pseudo-targeted mass spectrometry analysis. Sci. Rep. 2016, 6, 21433. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Rosebrock, A.P. Targeted full-scan lc-ms metabolomics: Simultaneous quantitation of knowns and feature discovery provide the best of both worlds. Bioanalysis 2017, 9, 5–8. [Google Scholar] [CrossRef] [PubMed]
  12. Wei, R.; Li, G.; Seymour, A.B. High-throughput and multiplexed lc/ms/mrm method for targeted metabolomics. Anal. Chem. 2010, 82, 5527–5533. [Google Scholar] [CrossRef] [PubMed]
  13. Yuan, M.; Breitkopf, S.B.; Yang, X.; Asara, J.M. A positive/negative ion–switching, targeted mass spectrometry–based metabolomics platform for bodily fluids, cells, and fresh and fixed tissue. Nat. Protoc. 2012, 7, 872–881. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Gu, H.; Zhang, P.; Zhu, J.; Raftery, D. Globally Optimized Targeted Mass Spectrometry: Reliable Metabolomics Analysis with Broad Coverage. Anal. Chem. 2015, 87, 12355–12362. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Roberts, L.D.; Souza, A.L.; Gerszten, R.E.; Clish, C.B. Targeted metabolomics. Curr. Protoc. Mol. Biol. 2012, 1–24. [Google Scholar] [CrossRef] [PubMed]
  16. Peng, J.; Chen, Y.-T.; Chen, C.-L.; Li, L. Development of a universal metabolome-standard method for long-term lc-ms metabolome profiling and its application for bladder cancer urine-metabolite-biomarker discovery. Anal. Chem. 2014, 86, 6540–6547. [Google Scholar] [CrossRef] [PubMed]
  17. Khan, N.A.; Auranen, M.; Paetau, I.; Pirinen, E.; Euro, L.; Forsström, S.; Pasila, L.; Velagapudi, V.; Carroll, C.J.; Auwerx, J.; et al. Effective treatment of mitochondrial myopathy by nicotinamide riboside, a vitamin B3. EMBO Mol. Med. 2014, 6, 721–731. [Google Scholar] [CrossRef] [PubMed]
  18. Roman-Garcia, P.; Quiros-Gonzalez, I.; Mottram, L.; Lieben, L.; Sharan, K.; Wangwiwatsin, A.; Tubio, J.; Lewis, K.; Wilkinson, D.; Santhanam, B.; et al. Vitamin b12–dependent taurine synthesis regulates growth and bone mass. J. Clin. Investig. 2014, 124, 2988–3002. [Google Scholar] [CrossRef] [PubMed]
  19. Nikkanen, J.; Forsström, S.; Euro, L.; Paetau, I.; Kohnz, R.A.; Wang, L.; Chilov, D.; Viinamäki, J.; Roivainen, A.; Marjamäki, P.; et al. Mitochondrial DNA Replication Defects Disturb Cellular dNTP Pools And Remodel One-Carbon Metabolism. Cell Metab. 2016, 23, 635–648. [Google Scholar] [CrossRef] [PubMed]
  20. Ahola, S.; Auranen, M.; Isohanni, P.; Niemisalo, S.; Urho, N.; Buzkova, J.; Velagapudi, V.; Lundbom, N.; Hakkarainen, A.; Muurinen, T.; et al. Modified Atkins diet induces subacute selective ragged red fiber lysis in mitochondrial myopathy patients. EMBO Mol. Med. 2016, 11, 1234–1247. [Google Scholar] [CrossRef] [PubMed]
  21. Rajendran, J.; Tomašić, N.; Kotarsky, H.; Hansson, E.; Velagapudi, V.; Kallijärvi, J.; Fellman, V. Effect of high-carbohydrate diet on plasma metabolomics in mice with mitochondrial respiratory chain complex III deficiency due to a Bcs1l mutation. Int. J. Mol. Sci. 2016, 17, 1824. [Google Scholar] [CrossRef] [PubMed]
  22. Szibor, M.; Dhandapani, P.K.; Dufour, E.; Holmström, K.M.; Zhuang, Y.; Salwig, I.; Wittig, I.; Heidler, J.; Gizatullina, Z.; Gainutdinov, T.; et al. Broad AOX expression in a genetically tractable mouse model does not disturb normal physiology. Dis. Model. Mech. 2017, 2, 163–171. [Google Scholar] [CrossRef] [PubMed]
  23. Schatton, D.; Pla-Martin, D.; Marx, M.C.; Hansen, H.; Mourier, A.; Nemazanyy, I.; Pessia, A.; Zentis, P.; Corona, T.; Kondylis, V.; et al. CLUH regulates mitochondrial metabolism by controlling translation and decay of target mRNAs. J. Cell Biol. 2017, 216, 675–693. [Google Scholar] [CrossRef] [PubMed]
  24. Purhonen, J.; Rajendran, J.; Mörgelin, M.; Uusi-Rauva, K.; Katayama, S.; Krjutskov, K.; Einarsdottir, E.; Velagapudi, V.; Kere, J.; Jauhiainen, M.; et al. Ketogenic diet attenuates hepatopathy in mouse model of respiratory chain complex III deficiency caused by a Bcs1l mutation. Sci. Rep. 2017, 7, 957. [Google Scholar] [CrossRef] [PubMed]
  25. Khan, N.A.; Nikkanen, J.; Yatsuga, S.; Wang, L.; Jackson, C.; Pessia, A.; Riikka, K.; Velagapudi, V.; Suomalainen, A. mTORC1 Regulates Mitochondrial Folate Cycle and the Integrated Stress Response in Mitochondrial Disease. Cell Metab. 2017, 26, 419–428. [Google Scholar] [CrossRef] [PubMed]
  26. Scott, T.A.; Quintaneiro, L.M.; Norvaisas, P.; Lui, P.P.; Wilson, M.P.; Leung, K.Y.; Herrera-Dominguez, L.; Sudiwala, S.; Pessia, A.; Clayton, P.T.; et al. Host-Microbe Co-metabolism Dictates Cancer Drug Efficacy in C. elegans. Cell 2017, 169, 442–456. [Google Scholar] [CrossRef] [PubMed]
  27. Puhka, M.; Takatalo, M.; Nordberg, M.E.; Valkonen, S.; Nandania, J.; Aatonen, M.; Yliperttula, M.; Laitinen, S.; Velgapudi, V.; Mirtti, T.; et al. Metabolomic profiling of extracellular vesicles and alternative normalisation methodsreveal enriched metabolites and strategies to study prostate cancer related changes. Theranostics 2017, 7, 3824–3841. [Google Scholar] [CrossRef] [PubMed]
  28. Schrade, A.; Kyrönlahti, A.; Akinrinade, O.; Pihlajoki, M.; Häkkinen, M.; Fischer, S.; Alastalo, T.-P.; Velagapudi, V.; Toppari, J.; Wilson, D.B.; et al. GATA4 is a key regulator of steroidogenesis and glycolysis in mouse Leydig cells. Endocrinology 2015, 156, 1860–1872. [Google Scholar] [CrossRef] [PubMed]
  29. Schrade, A.; Kyrönlahti, A.; Akinrinade, O.; Pihlajoki, M.; Fischer, S.; Rodriguez, V.M.; Velagapudi, V.; Toppari, J.; Wilson, D.B.; Heikinheimo, M. GATA4 regulates blood-testis barrier function and lactate metabolism in mouse Sertoli cells. Endocrinology 2016, 157, 2416–2431. [Google Scholar] [CrossRef] [PubMed]
  30. Ali-Sisto, T.; Tolmunen, T.; Toffol, E.; Viinamäki, H.; Mäntyselkä, P.; Valkonen-Korhonen, M.; Honkalampi, K.; Ruusunen, A.; Velagapudi, V.; Lehto, S.M. Purine metabolism is dysregulated in patients with major depressive disorder. Psychoneuroendocrinology 2016, 70, 25–32. [Google Scholar] [CrossRef] [PubMed]
  31. Kolho, K.-L.; Pessia, A.; Jaakkola, T.; de Vos, W.M.; Velagapudi, V. Fecal and serum metabolomics in pediatric inflammatory bowel disease. J. Crohn’s Colitis 2017, 11, 321–334. [Google Scholar] [CrossRef]
  32. Fu, Y.; Gaelings, L.; Söderholm, S.; Belanov, S.; Nandania, J.; Nyman, T.A.; Matikainen, S.; Anders, S.; Velagapudi, V.; Kainov, D.E. NJ872 inhibits influenza A virus replication without altering cellular antiviral responses. Antiviral Res. 2016, 133, 23–31. [Google Scholar] [CrossRef] [PubMed]
  33. Gaelings, L.; Söderholm, S.; Bugai, A.; Fu, Y.; Nandania, J.; Schepens, B.; Lorey, M.B.; Tynell, J.; Ginste, L.V.; Goffic, R.L.; et al. Regulation of Kynurenine Biosynthesis during Influenza Virus Infection. FEBS J. 2017, 284, 222–236. [Google Scholar] [CrossRef] [PubMed]
  34. Kuivanen, S.; Bespalov, M.M.; Nandania, J.; Ianevski, A.; Velagapudi, V.; De Brabander, J.K.; Kainov, D.E.; Vapalahti, O. Obatoclax, saliphenylhalamide and gemcitabine inhibit Zika virus infection in vitro and differentially affect cellular signaling, transcription and metabolism. Antiviral Res. 2017, 139, 117–128. [Google Scholar] [CrossRef] [PubMed]
  35. Bulanova, D.; Ianevski, A.; Bugai, A.; Akimov, Y.; Kuivanen, S.; Paavilainen, H.; Kakkola, L.; Nandania, J.; Turunen, L.; Ohman, T.; et al. Antiviral Properties of Chemical Inhibitors of Cellular Anti-Apoptotic Bcl-2 Proteins. Viruses 2017, 9, 271. [Google Scholar] [CrossRef] [PubMed]
  36. Zinkevičienė, A.; Kainov, D.; Girkontaitė, I.; Lastauskienė, E.; Kvedarienė, V.; Fu, Y.; Anders, S.; Velagapudi, V. Activation of Tryptophan and Phenylalanine Catabolism in the Remission Phase of Allergic Contact Dermatitis: A Pilot Study. Int. Arch. Allergy Immunol. 2016, 170, 262–268. [Google Scholar] [CrossRef] [PubMed]
  37. Haapaniemi, E.M.; Fogarty, C.L.; Katayama, S.; Vihinen, H.; Keskitalo, S.; Ilander, M.; Krjutškov, K.; Mustjoki, S.; Lehto, M.; Hautala, T.; et al. Combined immunodeficiency with hypoglycemia caused by mutations in hypoxia up-regulated 1. J. Allergy Clin. Immunol. 2017, 139, 1391–1393. [Google Scholar] [CrossRef] [PubMed]
  38. Rey, G.; Valekunja, U.K.; Feeney, K.A.; Wulund, L.; Milev, N.B.; Stangherlin, A.; Bollepalli, L.; Velagapudi, V.; O’Neill, J.S.; Reddy, A.B. The Pentose Phosphate Pathway Regulates the Circadian Clock. Cell Metab. 2016, 24, 462–473. [Google Scholar] [CrossRef] [PubMed]
  39. Mäntyselkä, P.; Ali-Sisto, T.; Kautiainen, H.; Savolainen, J.; Niskanen, L.; Viinamäki, H.; Velagapudi, V.; Lehto, S.M. The association between musculoskeletal pain and circulating ornithine: A population-based study. Pain Med. 2017, 18, 1145–1151. [Google Scholar] [CrossRef] [PubMed]
  40. European Medicines Agency. Guideline on Bioanalytical Method Validation; European Medicines Agency: London, UK, 2011. [Google Scholar]
  41. Yamil, S.-M.; Mark, S.L.; Lisa, E.K.; Maureen, L.S.; Kelly, H.T.; Paul, A.R.; William, G.M.; Daniel, W.B.; Tracey, B.S.; Dmitrii, V.T.; et al. Metabolite Profiling of a NIST Standard Reference Material for Human Plasma (SRM 1950): GC/MS, LC/MS, NMR and Clinical Laboratory Analyses, Libraries and Web-based resources. Anal. Chem. 2013, 24, 11725–11731. [Google Scholar] [CrossRef]
  42. Siskos, A.P.; Jain, P.; Römisch-Margl, W.; Bennett, M.; Achaintre, D.; Asad, Y.; Marney, L.; Richardson, L.; Koulman, A.; Griffin, J.L.; et al. Interlaboratory Reproducibility of a Targeted Metabolomics Platform for Analysis of Human Serum and Plasma. Anal. Chem. 2017, 89, 656–665. [Google Scholar] [CrossRef] [PubMed]
  43. Mihaleva, V.V.; Korhonen, S.P.; van Duynhoven, J.; Niemitz, M.; Vervoort, J.; Jacobs, D.M. Automated quantum mechanical total line shape fitting model for quantitative NMR-based profiling of human serum metabolites. Anal. Bioanal. Chem. 2014, 406, 3091–3102. [Google Scholar] [CrossRef] [PubMed]
  44. Hadfield, J.D. Mcmc methods for multi-response generalized linear mixed models: The mcmcglmm R package. J. Stat. Softw. 2010, 33, 22. [Google Scholar] [CrossRef]
  45. Singtoroj, T.; Tarning, J.; Annerberg, A.; Ashton, M.; Bergqvist, Y.; White, N.J.; Lindegardh, N.; Day, N.P. A new approach to evaluate regression models during validation of bioanalytical assays. J. Pharm. Biomed. Anal. 2006, 41, 219–227. [Google Scholar] [CrossRef] [PubMed]
Figure 1. A flowchart describing preanalytical stage, sample treatment and analysis, analytical stage, automation, and quality checks.
Figure 1. A flowchart describing preanalytical stage, sample treatment and analysis, analytical stage, automation, and quality checks.
Metabolites 08 00044 g001
Figure 2. Interbatch variation of high concentrations of spiked quality control (QC) serum samples (N = 6) analyzed on 3 different days. Data not shown for g-glutamylcysteine and orotic acid (coefficient of variation (CV) > 30%).
Figure 2. Interbatch variation of high concentrations of spiked quality control (QC) serum samples (N = 6) analyzed on 3 different days. Data not shown for g-glutamylcysteine and orotic acid (coefficient of variation (CV) > 30%).
Metabolites 08 00044 g002
Figure 3. Percentage mean recoveries for all metabolites at low, medium, and high concentration levels of QCs (N = 6 at each level) spiked in serum, brain, spleen, liver, and adherent cells. Data not shown for metabolites with poor chromatography and irreproducible results at different concentration levels.
Figure 3. Percentage mean recoveries for all metabolites at low, medium, and high concentration levels of QCs (N = 6 at each level) spiked in serum, brain, spleen, liver, and adherent cells. Data not shown for metabolites with poor chromatography and irreproducible results at different concentration levels.
Metabolites 08 00044 g003
Figure 4. Percentage mean stability was calculated for three different concentrations (low, medium, and high) spiked in QC serum (N = 6 at each level) during (A) wet extract and (B) freeze-thaw stability for all the metabolites.
Figure 4. Percentage mean stability was calculated for three different concentrations (low, medium, and high) spiked in QC serum (N = 6 at each level) during (A) wet extract and (B) freeze-thaw stability for all the metabolites.
Metabolites 08 00044 g004
Figure 5. Bar graph representing the stability of stock solutions and intermediate solutions of metabolites that were used to prepare calibration curve standards during the analysis. First, the stability was determined at the intermediate concentration levels (green bar) and then at the stock level (blue bar) only for those metabolites, which were not stable at the intermediate concentration level.
Figure 5. Bar graph representing the stability of stock solutions and intermediate solutions of metabolites that were used to prepare calibration curve standards during the analysis. First, the stability was determined at the intermediate concentration levels (green bar) and then at the stock level (blue bar) only for those metabolites, which were not stable at the intermediate concentration level.
Metabolites 08 00044 g005
Figure 6. Percentage of metabolites with less than 25% CV values of QC concentrations in 25 different batches analyzed over a period of 1 year.
Figure 6. Percentage of metabolites with less than 25% CV values of QC concentrations in 25 different batches analyzed over a period of 1 year.
Metabolites 08 00044 g006
Figure 7. Reproducibility of (A) retention time of respective chromatographic peaks and (B) regression coefficient values (R2) from external calibration curve standards for each metabolite analyzed in 25 different batches over a period of 1 year.
Figure 7. Reproducibility of (A) retention time of respective chromatographic peaks and (B) regression coefficient values (R2) from external calibration curve standards for each metabolite analyzed in 25 different batches over a period of 1 year.
Metabolites 08 00044 g007
Figure 8. A typical workflow for high-throughput targeted metabolomics analysis.
Figure 8. A typical workflow for high-throughput targeted metabolomics analysis.
Metabolites 08 00044 g008
Figure 9. (A) Comparison of metabolite concentration reference values given in the NIST SRM 1950 plasma against our method. Comparison of metabolite concentrations in our QC samples measured between our method and (B) BIOCRATES AbsoluteIDQ p180 kit, and (C) nuclear magnetic resonance (NMR) analyses. Error bars represent 95% confidence intervals.
Figure 9. (A) Comparison of metabolite concentration reference values given in the NIST SRM 1950 plasma against our method. Comparison of metabolite concentrations in our QC samples measured between our method and (B) BIOCRATES AbsoluteIDQ p180 kit, and (C) nuclear magnetic resonance (NMR) analyses. Error bars represent 95% confidence intervals.
Metabolites 08 00044 g009
Table 1. Median concentration levels (μmol/L) measured in pooled healthy adult serum samples.
Table 1. Median concentration levels (μmol/L) measured in pooled healthy adult serum samples.
Population Median (µmol/L)
95% Credibility Interval
Class and Metabolite NameHMDB IdEstimateLowerUpper
1. Alpha Amino Acids and Derivatives
2-AminoisobutyrateHMDB00019061.1280.8091.509
4-l-HydroxyprolineHMDB000072515.89511.30021.184
5-HydroxytryptophanHMDB00004720.0430.0310.058
ADMAHMDB00015390.9630.2482.230
AlanineHMDB0000161477.946339.667635.383
AminoadipateHMDB00005102.1211.5022.818
ArginineHMDB000051784.90260.795113.172
AsparagineHMDB000016847.20433.63562.727
AspartateHMDB000019126.93218.78635.424
BetaineHMDB0000043100.04572.251134.758
CitrullineHMDB000090427.81519.81236.972
CreatineHMDB000006461.18044.48182.161
CreatinineHMDB000056266.08547.09588.066
CystathionineHMDB00000990.1310.0930.174
DimethylglycineHMDB00000923.5412.5214.733
GABAHMDB00001120.1950.1380.259
G-GlutamylcysteineHMDB00010492.9662.1163.946
GlutamateHMDB000014853.44637.69170.613
GlutamineHMDB0000641791.142560.9991050.605
GlutathioneHMDB00001250.0210.0150.028
GlycineHMDB0000123243.717172.929325.417
GuanidoacetateHMDB00001282.6351.8803.498
HistidineHMDB000017788.61662.718117.848
HomocysteineHMDB00007420.4800.1241.153
HomoserineHMDB00007190.3370.2400.448
HydroxykynurenineHMDB00007320.0960.0680.128
IsoleucineHMDB000017283.31658.352110.232
KynurenineHMDB00006841.1460.8131.520
LeucineHMDB0000687126.07290.139167.353
LysineHMDB0000182176.519127.474236.393
MethionineHMDB000069629.20121.08139.249
OrnithineHMDB000021490.17765.043120.796
PhenylalanineHMDB000015988.32663.916117.977
ProlineHMDB0000162251.689180.529335.079
SDMAHMDB00033342.8622.0603.834
SerineHMDB0000187148.482107.012198.507
ThreonineHMDB0000167152.013108.810203.651
TryptophanHMDB000092933.96124.10445.261
TyrosineHMDB000015865.55845.77886.262
ValineHMDB0000883394.932282.460526.776
2. Benzoic Acids and Derivatives
3-HydroxanthranilateHMDB00014760.1880.1340.251
HippurateHMDB00007146.1014.3628.124
3. Beta Amino Acids and Derivatives
CarnosineHMDB00000330.0140.0100.019
PantothenateHMDB00002100.3100.2200.411
4. Bile Acids, Alcohols and Derivatives
ChenodeoxycholateHMDB000051853.66138.58871.987
CholateHMDB00006190.6760.4830.900
GlycocholateHMDB00001380.3730.2670.498
TaurochenodesoxycholateHMDB00009510.5070.3590.673
5. Carbohydrates and Carbohydrate Conjugates
d-Ribose-5-PHMDB00015481.2730.9191.710
GlyceraldehydeHMDB0001051239.946172.630322.782
SucroseHMDB00002581.4171.0151.886
6. Dialkylamines
SpermidineHMDB000125733.60123.95844.759
7. Dicarboxylic Acids and Derivatives
SuccinateHMDB00002547.9125.59710.491
8. Fatty Acyls
AcetylcarnitineHMDB00002019.7092.32322.264
DecanoylcarnitineHMDB00006510.3050.0870.720
HexanoylcarnitineHMDB00007050.0550.0150.129
IsobutyrylcarnitineHMDB00007360.2480.0610.578
IsovalerylcarnitineHMDB00006880.1020.0270.240
OctanoylcarnitineHMDB00007910.2970.0760.691
PropionylcarnitineHMDB00008240.4230.1190.998
9. Folates
FolateHMDB00001210.0110.0030.027
10. Glucuronic Acid and Derivatives
GlucuronateHMDB00001271.9601.4042.627
11. Imidazoles
1-MethylhistamineHMDB00008980.0060.0040.008
AllantoinHMDB00004622.4471.7413.254
12. Indoles and Derivatives
5-HydroxyindoleacetateHMDB00007630.0740.0530.099
13. Keto Acids and Derivatives
AcetoacetateHMDB00000606.7134.7988.988
14. Organic Phosphoric Acids and Derivatives
PhosphoethanolamineHMDB00002243.3162.3564.420
15. Organosulfonic Acids
TaurineHMDB0000251221.359156.971294.936
TaurocholateHMDB00000360.0830.0600.112
16. Oxides
Trimethylamine N-oxideHMDB00009251.4771.0551.973
17. Phenols
HomogentisateHMDB00001300.1150.0820.153
NormetanephrineHMDB00008190.00100.00070.0014
18. Pteridines and Derivatives
NeopterinHMDB00008450.0050.0040.007
19. Purines and Derivatives
AdenineHMDB00000340.0070.0050.009
AdenosineHMDB00000500.0080.0060.011
AMPHMDB00000450.1060.0750.140
cAMPHMDB00000580.0050.0030.006
cGMPHMDB00013140.0090.0020.026
GuanosineHMDB00001330.4450.3150.593
HypoxanthineHMDB000015758.83841.90878.087
IMPHMDB00001750.2120.1510.282
InosineHMDB000019536.06125.79848.008
XanthineHMDB00002923.9262.7905.211
XanthosineHMDB00002990.3520.2490.466
20. Pyridines and Derivatives
4-PyridoxateHMDB00000170.0570.0420.077
CotinineHMDB00010460.5310.3820.710
NADHMDB00009020.0150.0110.020
NiacinamideHMDB00014060.3950.2790.524
NicotinateHMDB00014880.0120.0090.017
PyridoxineHMDB00002390.00070.00050.0009
21. Pyrimidines and Derivatives
CytidineHMDB00000890.0030.0020.004
CytosineHMDB00006300.0800.0560.106
DeoxycytidineHMDB00000140.8710.4401.461
DeoxyuridineHMDB00000120.5430.3910.726
OrotateHMDB00002260.0360.0060.096
UDP GlucoseHMDB00002860.2320.1650.308
UracilHMDB00003000.0580.0420.078
22. Quaternary Ammonium Salts
CarnitineHMDB000006285.66961.210114.060
CholineHMDB000009795.17067.565126.279
23. Quinolines and Derivatives
KynurenateHMDB00007150.0440.0310.058
24. Sugar Alcohols
Myo-inositolHMDB000021116.98912.07522.606
SorbitolHMDB00002473.6922.6174.904

Share and Cite

MDPI and ACS Style

Nandania, J.; Peddinti, G.; Pessia, A.; Kokkonen, M.; Velagapudi, V. Validation and Automation of a High-Throughput Multitargeted Method for Semiquantification of Endogenous Metabolites from Different Biological Matrices Using Tandem Mass Spectrometry. Metabolites 2018, 8, 44. https://doi.org/10.3390/metabo8030044

AMA Style

Nandania J, Peddinti G, Pessia A, Kokkonen M, Velagapudi V. Validation and Automation of a High-Throughput Multitargeted Method for Semiquantification of Endogenous Metabolites from Different Biological Matrices Using Tandem Mass Spectrometry. Metabolites. 2018; 8(3):44. https://doi.org/10.3390/metabo8030044

Chicago/Turabian Style

Nandania, Jatin, Gopal Peddinti, Alberto Pessia, Meri Kokkonen, and Vidya Velagapudi. 2018. "Validation and Automation of a High-Throughput Multitargeted Method for Semiquantification of Endogenous Metabolites from Different Biological Matrices Using Tandem Mass Spectrometry" Metabolites 8, no. 3: 44. https://doi.org/10.3390/metabo8030044

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