Next Article in Journal
Metabolomic Analysis of Response to Nitrogen-Limiting Conditions in Yarrowia spp.
Previous Article in Journal
SeMPI 2.0—A Web Server for PKS and NRPS Predictions Combined with Metabolite Screening in Natural Product Databases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Holistic Metabolomic Laboratory-Developed Test (LDT): Development and Use for the Diagnosis of Early-Stage Parkinson’s Disease

by
Petr G. Lokhov
1,*,
Dmitry L. Maslov
1,
Steven Lichtenberg
2,
Oxana P. Trifonova
1 and
Elena E. Balashova
1
1
Institute of Biomedical Chemistry, 10 Building 8, Pogodinskaya Street, 119121 Moscow, Russia
2
Metabometrics, Inc, 651 N Broad St, Suite 205 #1370, Middletown, DE 19709, USA
*
Author to whom correspondence should be addressed.
Metabolites 2021, 11(1), 14; https://doi.org/10.3390/metabo11010014
Submission received: 9 October 2020 / Revised: 24 December 2020 / Accepted: 24 December 2020 / Published: 29 December 2020
(This article belongs to the Section Advances in Metabolomics)

Abstract

:
A laboratory-developed test (LDT) is a type of in vitro diagnostic test that is developed and used within a single laboratory. The holistic metabolomic LDT integrating the currently available data on human metabolic pathways, changes in the concentrations of low-molecular-weight compounds in the human blood during diseases and other conditions, and their prevalent location in the body was developed. That is, the LDT uses all of the accumulated metabolic data relevant for disease diagnosis and high-resolution mass spectrometry with data processing by in-house software. In this study, the LDT was applied to diagnose early-stage Parkinson’s disease (PD), which currently lacks available laboratory tests. The use of the LDT for blood plasma samples confirmed its ability for such diagnostics with 73% accuracy. The diagnosis was based on relevant data, such as the detection of overrepresented metabolite sets associated with PD and other neurodegenerative diseases. Additionally, the ability of the LDT to detect normal composition of low-molecular-weight compounds in blood was demonstrated, thus providing a definition of healthy at the molecular level. This LDT approach as a screening tool can be used for the further widespread testing for other diseases, since ‘omics’ tests, to which the metabolomic LDT belongs, cover a variety of them.

Graphical Abstract

1. Introduction

Metabolomics studies have demonstrated the possibility of using the identification of metabolites for the successful diagnosis of many diseases [1]. These data give hope for the successful application of metabolomics methods in medicine. Therefore, researchers are trying to create such omics tests for disease diagnosis, risk assessment of their development, and determination of the patient’s response to treatment [2]. However, omics tests use in clinical practice is very challenging due to the complexity of most omics technologies, thus making their standardization for acceptance in clinics extremely difficult [3]. Consequently, considering omics tests as in-house laboratory-developed tests (LDTs) is the most evident solution to this problem. An LDT is a type of in vitro diagnostic test that is developed and used within the same laboratory [4] and is used to measure a wide range of substances, including nucleic acids, proteins, and low-molecular-weight compounds in different biological samples. To date, numerous LDTs have been developed for the diagnosis of various diseases, including cancers, infections, genetic disorders, and other pathologies [5,6,7,8,9,10,11]. In this work, an LDT comprising the latest advancements in metabolomics was developed.
Recently, we established an LDT that can reveal overrepresentation of pathways, providing a basis for diagnosis [12]. The LDT workflow included blood sample preparation, high-resolution mass spectrometry, pre-processing of mass spectra, a compound annotation algorithm, and a statistical model testing the retrieved metabolic data against human pathways. The diagnostic performance of this LDT was demonstrated for the diagnosis of Parkinson’s disease (PD). Herein, the LDT was further developed by introducing information about metabolites, such as their known connection with diseases, conditions of the body, and location in the body. This improvement completes the development of the LDT, making it a holistic metabolomic test whereby the currently available data in the human metabolome database (HMDB [13]) are applied through panoramic measurement of metabolites in the biomaterial. This updated LDT uses disease-associated metabolite sets (631 diseases), pathway-associated metabolite sets (808 human metabolic pathways), metabolite sets associated with abnormal concentrations of metabolites (352 conditions), and location-based metabolite sets (110 sets based on their location in organs, tissues, and subcellular localization). Updated LDT was also tested as a diagnostic method for early-stage PD.
PD is a widespread neurodegenerative disease. Its incidence has increased dramatically in recent years due to the aging population. Due to the etiologic and pathogenic heterogeneity of PD [14,15], the discovery of biomarkers for diagnosis of PD is complicated and has not yet been successful [16]. Metabolomics technologies that make it possible to measure the entire collection of low-molecular-weight compounds of a sample may help in this situation. Thus, it is reasonable to test the updated LDT on early-stage PD. In addition, the success of such diagnostics can more clearly confirm the validity of the metabolomic LDT.

2. Results

2.1. Mass Spectrometric Analysis of Compounds in Blood

Mass spectrometric analysis, as the first analytical block of the LDT (Figure 1), generated high-resolution spectra of approximately 10,000 mass peaks of low-molecular-weight compounds in the blood plasma samples. The measured masses were submitted to the bioinformatic treatment block of the LDT that resulted in the annotation of 709 compounds (Table 1).

2.2. Metabolite Set Overrepresentation Patterns in the LDT Output

A case-control comparison revealed patterns of PD in the LDT output. The LDT output was generated as a metabolite set names cloud for the controls and cases. Among the top overrepresented disease-associated metabolite sets, the diseases semantically similar to PD were presented (Alzheimer’s disease, Lewy body disease, and frontotemporal dementia; Figure 1a). The PD-relevant patterns were also found in the pathway-associated metabolite sets (Figure 2b) and the metabolite sets associated with abnormal concentrations of metabolites (Figure 2c). With the use of the location-based metabolite sets, excluding a neuron-associated set, it was difficult to associate the top overrepresented sets with PD directly (Figure 2d).
The metabolite set representation scores for the control samples, samples from PD patients, and the top 20 overrepresented metabolite sets for the PD patients are listed in Table 2.

2.3. Diagnostic Performance of the LDT

Using the LDT, the diagnosis of PD based on diagnostic score reached an accuracy of 73% (Table 2). It is noteworthy that the high diagnostic performance practically did not depend on the type of metabolite sets used, that is, whether the disease-associated sets, pathway-associated sets, abnormal concentration sets, or localization-based sets were used did not matter. This is explained by the fact that the same metabolites, indicating PD, were applied to these groups.
Upon analysis of the diagnostic performance based on overrepresented metabolite sets, it was noted that a neurodegenerative disease of the central nervous system could be suspected in patients. Thus, the metabolite sets of Alzheimer’s disease, Lewy body disease, and frontotemporal dementia were overrepresented. Almost all of the overrepresented pathways (Table 2) were also relevant to PD. Thus, dysregulation of transcription and translation is described in PD [18,19]. Moreover, it is a known fact that dopamine level is connected with PD [20]. The role of lipid metabolism and mitochondria is also described for PD [21]. Furthermore, overrepresented pterin synthesis is directly connected with neurotransmitters; Segawa syndrome together with guanosine triphosphate cyclohydrolase deficiency is related to Dopa-responsive dystonia; and 6-pyruvoyltetrahydropterin synthase deficiency is a neurodegenerative disease that, similar to dihydropteridine reductase (DHPR) deficiency, is treated by levodopa. Overexpression was also observed in sets of metabolites associated with internal organs, for example, high overexpression in a set of metabolites associated with colorectal cancer. Synucleinopathy can explain such overexpression. Synucleinopathy develops in various parts of the nervous system at PD and leads to denervation of the heart, disturbances in the large intestine, esophagus, kidneys, etc. A metabolite set associated directly with PD was not presented in the top list of overrepresented disease-related sets. Perhaps this is due to the lack of such a set that is applicable to the LDT.
The overrepresented metabolite sets associated with abnormal concentrations also contributed to the diagnostic performance of the LDT. Regarding the top overexpressed sets (schizophrenia, alcohol intoxication and drunk driver, pellagra, etc.; see Table 2), it is known that dopamine and dopaminergic neurons play an important role in schizophrenia, as well as in PD [22,23]. Moreover, it has been argued that a functional excess of dopamine or oversensitivity of certain dopamine receptors is one of the causal factors in schizophrenia. In schizophrenia, the antipsychotic effects of traditional ‘neuroleptic’ drugs, such as chlorpromazine, are highly correlated with their ability to block dopamine receptors and reduce the effects of dopamine. The overrepresentation of alcohol intoxication and drunk driver sets may be explained by the fact that alcohol has a powerful effect on dopamine activity in the brain, which has been revealed in animals [24] and human studies [25]. Pellagra is due to a diet that does not contain enough niacin and tryptophan, which, in turn, can be converted into serotonin and are altered in PD [26].

2.4. Diagnosis of PD by the LDT

Figure 3 shows the metabolite set representation scores for each participant in this study. This figure shows that diagnostics based on the metabolite set representation score is possible. Unlike analysis of the case-control sets, which reveals the common group patterns, personal data analysis is more complicated. Aged patients have a whole range of diseases leading to the overrepresentation of different metabolite sets, thus making diagnosis very difficult. However, for some individuals, it is quite possible. Figure 4 presents the LDT output as a metabolite set names cloud for one person.
The overrepresentation of metabolite sets associated with neurodegenerative diseases (Figure 4a), together with the overrepresentation of the Dopa-responsive dystonia pathway (Figure 4b), which is a movement disorder characterized by muscle tone and Parkinsonian features, the change in the concentration of substances matching alcohol intoxication (Figure 4c), and priority localization of such metabolites (Figure 4c) with an abnormal concentration in neurons, would allow health professionals to suspect neurodegenerative disease and motivate them to go deeper and to obtain additional support for a PD diagnosis. An overrepresentation of the metabolic sets associated with Fabry disease [27], pellagra, which, as mentioned above, has a connection with PD, as well as many Dopa-related pathways (Segawa syndrome, pterine biosynthesis, sepiapterin reductase deficiency, guanosine triphosphate cyclohydrolase deficiency, 6-pyruvoyltetrahydropterin synthase deficiency (PTPS), and DHPR-deficiency) was observed. In other words, the LDT indicates which direction to search, thus helping the clinician to select of confirmatory, targeted tests. For this patient, an appointment for a single-photon emission computerized tomography scan is fully justified, and PD will be diagnosed.

2.5. LDT Output for a ‘Healthy’ Individual

In the control group, there were individuals for whom the LDT did not reveal any overrepresentation in the metabolite sets. Thus, from the point of view of metabolic processes, these people can be considered ‘healthy,’ since the composition of low-molecular-weight compounds in their blood corresponds to the age norm. An example of the LDT results for such a person is presented in Figure 5.

3. Discussion

The Food and Drug Administration (USA) considers LDTs as tests that are designed, manufactured, and used inside the same laboratory [4]. It simplifies the implementation of metabolomics-based tests, bringing protocols and standardization activities to single laboratory routines. The developed LDT is based on direct mass spectrometry of blood plasma, which has been widely used in metabolomics and, in particular, in the laboratory where LDTs were developed for the study of cancers [28,29,30,31,32], diabetes [33], and PD [34]. This type of mass spectrometry is characterized by a high processing speed and a relatively high reproducibility [35,36,37], which are important for the use of mass spectrometry in the clinic. The mass spectrometry data processing, like peak alignment and data standardization, was specially developed for high-resolution mass spectra and successfully used for many years in studies of blood plasma [29,38] and now is implemented in the LDT.
Generally, mass spectrometry allows the detection of hundreds of compounds in metabolomics studies, which is crucial for obtaining biochemical information [39]. Unfortunately, the vast majority of compounds in the sample remain unknown [40]; current annotation methods require a clear mass spectrometric picture of compounds or its fragments, which can be obtained only for well-separated and abundant metabolites. In the LDT described in this study, a recently developed biochemical context-driven annotation is realized for annotation of compounds, which uses the knowledge of their biotransformation in metabolic pathways. This approach was introduced by Rogers and coworkers [41] and further updated by Silva and coworkers [42]. Later, the suitability of this approach for blood plasma samples was also demonstrated [17]. This updated algorithm was implemented in the LDT, which allowed the annotation of more than 700 metabolites per sample. The obtained metabolite annotations were classified as putatively annotated compounds (level 2 of metabolite identification), according to the Metabolomics Standards Initiative standard [43], because two independent orthogonal features of each metabolite were used for annotation (accurate mass tag and biochemical context). Thus, annotation results do not include the most robust identifications at level 1, which is acceptable for medical purposes and requires a chemical standard for identification. Obviously, for big data, to which the metabolomics data relates, a level 1 often is impossible, thus making it reasonable to use the described approach as a screening technology, which helps the clinician to optimize the selection of confirmatory, secondary tests [44].
The selection of PD to test this LDT was not an accident. Previously, it has been shown that an LDT can reveal the pathway overrepresentation efficient for the diagnosis of PD [12], thus making it reasonable to complete LDT development by introducing metabolic data about diseases, different organism conditions, and metabolite location. Figure 3 shows the representation scores for each participant in this study and each metabolite set. This figure confirms that each metabolite set contributes to diagnostics, and the accuracy of such kind diagnostics is 73% (Table 3). It should be noted, the early laboratory diagnosis of PD is currently unavailable and urgently needed for effective therapy [45,46,47]. However, the multifactorial nature of PD complicates the development of conventional biomarker-based tests. The clinical application of ‘panoramic’ methods, to which metabolomic methods are related, have the complex workflow that makes their standardization and following registration illusive. The usage of the LDT in such a situation overcomes this obstacle because all of the LDT-related routines are located in a single laboratory.
The developed LDT uses currently available data on the concentrations of metabolites in humans and uses them to analyze panoramically measured blood composition data. Along with the use of modern data processing algorithms, it can be argued that the LDT is an omics test that demonstrates the current diagnostic capabilities of metabolomics, the most obvious of which are as follows:
  • Confirmation of a person’s healthy state. This option of the LDT is the most obvious; the output of the LDT in this case is self-explanatory and comprehensively confirms human health at the molecular level. The LDT shows that the detected deviations in the blood composition do not form any patterns specific to a disease or pathology. So, the LDT is ready for use to determine wellness and longevity. It is expected that the healthy state can be confirmed by the LDT and that any abnormalities that will appear at the molecular level can be detected in a timely manner, which lays the foundation for a long and quality life.
  • Score-based diagnostics. Score-based diagnostics requires control samples and samples from a cohort of patients with disease. The advantage of such diagnostics is the absence of human error in diagnosis and possible full automation.
  • Disease diagnosis based on metabolite set overrepresentation (i.e., without diagnostic scoring). This option of the LDT is ready to use (i.e., cohorts are not required) for the diagnosis of a wide diversity of diseases. The metabolite set names cloud allows visualization of the LDT output data that a physician can interpret. An example of this is demonstrated in this paper for the diagnosis of PD, although, among the LDT outputs, there were also results that were difficult to interpret. It is possible that the effectiveness of the LDT output interpretation will increase as further LDT output data are accumulated. Most importantly, the LDT is panoramic in terms of measuring substances and untargeted in terms of diagnosing diseases, which in the end makes it especially valuable.
The main disadvantage of the LDT is its complexity, which leads to the lack of strictly standardized protocols and direct dependence on the equipment used and the experience of the staff working on it. In general, therefore, the LDT is implemented in a single laboratory and, as a rule, is not translated to other laboratories. But this problem can be substantially ignored by the compatibility of the mass spectrometric measurement in the LDT with a dried blood spot [48] because a dried blood spot can be obtained without assistance at home and transported on a blood sample card at room temperature by mail, thus making the LDT available almost everywhere. It seems that such a type of laboratory diagnostic method that does not require direct contact with people during blood sampling will be relevant in light of the pandemic that has swept the world.

4. Materials and Methods

4.1. Mass Spectra of Blood Plasma

Samples of blood plasma used in this study were taken from the previously published study [34]. Table 3 presents the clinical characteristics of the cohort.
Table 3. Study cohort characteristics.
Table 3. Study cohort characteristics.
CharacteristicsValues
Subjects with PDControl Subjects
Number2828
Gender (male/female)14/1414/14
Age (years; mean ± s.d. (range))62.6 ± 8.6 (37–77)62.8 ± 8.7 (45–77)
PD stages (1/1.5/2/2.5) 16/6/12/4-
1 PD stages are according to Hoehn and Yahr scale [49].
Samples were analyzed with a maXis hybrid quadrupole time-of-flight mass spectrometer with an electrospray ionization source as described in the previously published study [12]. Normalization of mass peak intensities was performed as described previously [38]. The alignment of the m/z values of the mass peaks to the different mass spectra was performed as described previously [29]. The resulting m/z values with a nonzero mass peak intensity for more than nine samples (removes noise and suspect data) were submitted to the metabolite search block of the in-house software.

4.2. Compound Annotation

The search for correspondence of each mass peak to metabolite identifiers was done by the metabolite search block of the in-house software as described in the previous publication [12]. The HMDB (www.hmdb.ca) was used as the source of the m/z values and identifiers from the Kyoto Encyclopedia of Genes and Genomes database (KEGG IDs). A compound annotation algorithm was recently developed and described in detail [42]. This algorithm uses metabolic pathway data and allows for the effective annotation of low-molecular-weight blood components (metabolome) with relatively high speed. In the list of compound names, many candidates, on average, were associated with one mass. The main task of the algorithm is to compare the obtained experimental data, i.e., mass spectra, with the available information on biochemical pathways and to decline all false candidates. It is known that the concentrations of compounds involved in the same pathways correlate [50]. Thus, if mass spectrometry data for a set of samples are available, the correlation between the mass of interest and other mass peaks can be found. The masses of these correlating peaks can also be associated with a set of compounds in which their locations in the metabolic pathway must be bunched around the compound with the true annotation. The details of the application of this algorithm for blood plasma samples have been described previously [17]. The next release of this algorithm with updated source code to make it more suitable for blood samples was recently used in the LDT [12] and applied in this work.

4.3. Metabolite Set Overrepresentation Analysis

Analyzed HMDB data were used to compile metabolite sets: disease-associated metabolite sets (631 diseases), pathway-associated metabolite sets (808 human metabolic pathways), metabolite sets associated with the abnormal concentrations of metabolites (352 conditions), and location-based metabolite sets (110 sets based on location in organs, tissues, and subcellular localization). Selected metabolites in the case-control or individual sample studies were projected on these metabolite sets, and the obtained results were compared with the projections, which were performed 30,000 times with the same number of randomly selected metabolites. The obtained results were normalized to produce a metabolite set representation score for each set as described previously [12]. This score also was used to visualize the results as a metabolite set names cloud in which the font size is related. The scores for the top 20 overrepresented metabolite sets were summarized to produce the final diagnostic score for each person who participated in the study.
To reveal the PD-associated pattern in metabolite sets, the mean value for representation scores for cases were compared with those of the controls.
To estimate diagnostic performance of LDT, as well as to estimate the metabolic data projection on the separate metabolite sets, the perfcurve function of the MATLAB program was used. This function presents the accuracy, sensitivity, and specificity for each point of receiver operating characteristic (ROC) curve and selects the optimal values.

4.4. Analysis of Individual Samples by the LDT

To reveal mass peaks belonged to metabolites with abnormal concentration, the Z-score calculation and the leave-one-out approach were applied for mass peak intensities as described previously [51]. This method involves the one-by-one removal of each data point (sample) from the dataset and recalculation of the model parameters based on the remaining data. The model is then tested by the excluded sample.
The LDT workflow used to analyze blood plasma samples is presented in Figure 1. The in-house software, as a part of the LDT, was implemented in MATLAB and was used for data pre-processing, database searching, and overrepresentation analysis. To perform all calculations, a Lenovo (Intel® Xeon® E-2176M CPU 2.70 GHz, Windows 10 Pro) computer was used.

5. Conclusions

LDTs or, more familiarly, ‘home brew’ tests, have been around for decades. Traditionally, their scope is small for low-risk diagnostic applications. Today, more complex LDTs can be used. Being a ‘direct-to-customer’ test, LDTs provide clinical results to a wide range of customers: physicians and their patients, researchers, citizen scientists, and simply educated people. While the current work demonstrated a metabolomic LDT for supporting the diagnosis of PD, the omics nature of the LDT suggests that it can be used for a variety of diseases. Information, in the name cloud form, in the LDT output about the state of the organism can be easy and quick to read by a wide range of customers although such an output is formed by metabolic big data. Diagnostics through the use of the metabolomic LDT for a wide variety of diseases on a single dried blood spot obtained at home lays the groundwork for improvements in terms of accessibility, price, and versatility of laboratory diagnostics, which can lead to an improvement in people’s living standards.

Author Contributions

Project administration, conceptualization, funding acquisition, P.G.L.; methodology, O.P.T. and D.L.M.; resources and samples curation, O.P.T. and D.L.M.; investigation, validation, writing—original draft preparation, review and editing, E.E.B. and P.G.L.; visualization, S.L. and E.E.B.; formal analysis, O.P.T. and P.G.L.; data curation, software and programming, S.L. and P.G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Russian Foundation for Basic Research (grant № 19-29-01125 “A metabolomic analysis of biochemical disorders in the prodromal stage of Parkinson’s disease for a systematic assessment of the risk of disease and personalization of pharmacotherapy”). In-house software development was supported by Metabometrics, Inc.

Institutional Review Board Statement

Ethical review and approval were waived for this study, due to samples of blood plasma were taken from the previously published study. Institutional Review Board Statement with the approval number were provided in this original study [34].

Informed Consent Statement

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

Data Available Statement

The data presented in this study are openly available in FigShare at doi 10.6084/m9.figshare.13621277.

Conflicts of Interest

P.G.L. is a shareholder of Metabometrics, Inc. The other authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Trifonova, O.; Lokhov, P.; Archakov, A. Postgenomics diagnostics: Metabolomics approaches to human blood profiling. OMICS 2013, 17, 550–559. [Google Scholar] [CrossRef]
  2. Omenn, G.S.; DeAngelis, C.D.; DeMets, D.L.; Fleming, T.R.; Geller, G.; Gray, J.; Hayes, D.F.; Henderson, G.I.; Kessler, L.; Lapidus, S.; et al. Evolution of Translational Omics: Lessons Learned and the Path Forward—Institute of Medicine; National Academies Press: Washington, DC, USA, 2012; ISBN 0309224185. [Google Scholar]
  3. Nass, S.J.; Moses, H.L. Cancer Biomarkers: The Promises and Challenges of Improving Detection and Treatment; National Academies Press: Washington, DC, USA, 2007; ISBN 0309667119. [Google Scholar]
  4. FDA. Laboratory Developed Tests. Available online: https://www.fda.gov/medical-devices/vitro-diagnostics/laboratory-developed-tests (accessed on 5 October 2018).
  5. Burchard, P.R.; Abou, T.A.N.; Lefferts, J.A.; Lewis, L.D.; Tsongalis, G.J.; Cervinski, M.A. Development of a rapid clinical TPMT genotyping assay. Clin. Biochem. 2014, 47, 126–129. [Google Scholar] [CrossRef] [PubMed]
  6. Munari, E.; Zamboni, G.; Lunardi, G.; Marconi, M.; Brunelli, M.; Martignoni, G.; Netto, G.J.; Quatrini, L.; Vacca, P.; Moretta, L.; et al. PD-L1 expression in non–small cell lung cancer: Evaluation of the diagnostic accuracy of a laboratory-developed test using clone E1L3N in comparison with 22C3 and SP263 assays. Hum. Pathol. 2019, 90, 54–59. [Google Scholar] [CrossRef] [PubMed]
  7. Fiset, P.O.; Labbé, C.; Young, K.; Craddock, K.J.; Smith, A.C.; Tanguay, J.; Pintilie, M.; Wang, R.; Torlakovic, E.; Cheung, C.; et al. Anaplastic lymphoma kinase 5A4 immunohistochemistry as a diagnostic assay in lung cancer: A Canadian reference testing center’s results in population-based reflex testing. Cancer 2019, 125, 4043–4051. [Google Scholar] [CrossRef] [PubMed]
  8. Tinawi-Aljundi, R.; King, L.; Knuth, S.T.; Gildea, M.; Ng, C.; Kahl, J.; Dion, J.; Young, C.; Schervish, E.W.; Frontera, J.R.; et al. One-year monitoring of an oligonucleotide fluorescence in situ hybridization probe panel laboratory-developed test for bladder cancer detection. Res. Rep. Urol. 2015, 7, 49–55. [Google Scholar]
  9. Brukner, I.; Eintracht, S.; Forgetta, V.; Papadakis, A.I.; Spatz, A.; Oughton, M. Laboratory-developed test for detection of acute Clostridium difficile infections with the capacity for quantitative sample normalization. Diagn. Microbiol. Infect. Dis. 2019, 95, 113–118. [Google Scholar] [CrossRef]
  10. Kulis-Horn, R.K.; Tiemann, C. Evaluation of a laboratory-developed test for simultaneous detection of norovirus and rotavirus by real-time RT-PCR on the Panther Fusion® system. Eur. J. Clin. Microbiol. Infect. Dis. 2020, 39, 103–112. [Google Scholar] [CrossRef] [Green Version]
  11. Calvert, J.; Saber, N.; Hoffman, J.; Das, R. Machine-learning-based laboratory developed test for the diagnosis of sepsis in high-risk patients. Diagnostics 2019, 9, 20. [Google Scholar] [CrossRef] [Green Version]
  12. Lokhov, P.G.; Trifonova, O.P.; Maslov, D.L.; Lichtenberg, S.; Balashova, E.E. Diagnosis of Parkinson’s disease by a metabolomics-based laboratory-developed test (LDT). Diagnostics 2020, 10, 332. [Google Scholar] [CrossRef]
  13. Wishart, D.S.; Feunang, Y.D.; Marcu, A.; Guo, A.C.; Liang, K.; Vázquez-Fresno, R.; Sajed, T.; Johnson, D.; Li, C.; Karu, N.; et al. HMDB 4.0: The human metabolome database for 2018. Nucleic Acids Res. 2018, 46, D608–D617. [Google Scholar] [CrossRef]
  14. Sommer, A.; Winner, B.; Prots, I. The Trojan horse—Neuroinflammatory impact of T cells in neurodegenerative diseases. Mol. Neurodegener. 2017, 12, 78. [Google Scholar] [CrossRef] [Green Version]
  15. Lei, S.; Powers, R. NMR Metabolomics Analysis of Parkinson’s Disease. Curr. Metab. 2013, 1, 191–209. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Chen-Plotkin, A.S.; Albin, R.; Alcalay, R.; Babcock, D.; Bajaj, V.; Bowman, D.; Buko, A.; Cedarbaum, J.; Chelsky, D.; Cookson, M.R.; et al. Finding useful biomarkers for Parkinson s disease. Sci. Transl. Med. 2018, 10, 454. [Google Scholar] [CrossRef] [PubMed]
  17. Lokhov, P.G.; Balashova, E.E.; Trifonova, O.P.; Maslov, D.L.; Ponomarenko, E.A.; Archakov, A.I. Mass Spectrometry-Based Metabolomics Analysis of Obese Patients’ Blood Plasma. Int. J. Mol. Sci. 2020, 21, 568. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Dorval, V.; Hébert, S.S. LRRK2 in transcription and translation regulation: Relevance for Parkinson’s disease. Front. Neurol. 2012, 3, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Parlato, R.; Liss, B. How Parkinson’s disease meets nucleolar stress. Biochim. Biophys. Acta Mol. Basis Dis. 2014, 1842, 791–797. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Kalia, L.V.; Lang, A.E. Parkinson’s disease. Lancet 2015, 386, 896–912. [Google Scholar] [CrossRef]
  21. Eckert, G.P.; Lipka, U.; Muller, W.E. Omega-3 fatty acids in neurodegenerative diseases: Focus on mitochondria. Prostaglandins Leukot. Essent. Fat. Acids 2013, 88, 105–114. [Google Scholar] [CrossRef]
  22. Birtwistle, J.; Baldwin, D. Role of dopamine in schizophrenia and Parkinson’s disease. Br. J. Nurs. 1998, 7, 832–834. [Google Scholar] [CrossRef]
  23. Brisch, R.; Saniotis, A.; Wolf, R.; Bielau, H.; Bernstein, H.G.; Steiner, J.; Bogerts, B.; Braun, K.; Kumaratilake, J.; Henneberg, M.; et al. The role of dopamine in schizophrenia from a neurobiological and evolutionary perspective: Old fashioned, but still in vogue. Front. Psychiatry 2014, 19, 47. [Google Scholar]
  24. Engleman, E.A.; Keen, E.J.; Tilford, S.S.; Thielen, R.J.; Morzorati, S.L. Ethanol drinking reduces extracellular dopamine levels in the posterior ventral tegmental area of nondependent alcohol-preferring rats. Alcohol 2011, 45, 549–557. [Google Scholar] [CrossRef] [Green Version]
  25. Urban, N.B.L.; Kegeles, L.S.; Slifstein, M.; Xu, X.; Martinez, D.; Sakr, E.; Castillo, F.; Moadel, T.; O’Malley, S.S.; Krystal, J.H.; et al. Sex differences in striatal dopamine release in young adults after oral alcohol challenge: A positron emission tomography imaging study with [11C]raclopride. Biol. Psychiatry 2010, 68, 689–696. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Szabó, N.; Kincses, Z.T.; Toldi, J.; Vécsei, L. Altered tryptophan metabolism in Parkinson’s disease: A possible novel therapeutic approach. J. Neurol. Sci. 2011, 310, 256–260. [Google Scholar] [CrossRef] [PubMed]
  27. Wise, A.H.; Yang, A.; Naik, H.; Stauffer, C.; Zeid, N.; Liong, C.; Balwani, M.; Desnick, R.J.; Alcalay, R.N. Parkinson’s disease prevalence in Fabry disease: A survey study. Mol. Genet. Metab. Rep. 2018, 14, 27–30. [Google Scholar] [CrossRef] [PubMed]
  28. Lokhov, P.G.; Dashtiev, M.I.; Moshkovskii, S.A.; Archakov, A.I. Metabolite profiling of blood plasma of patients with prostate cancer. Metabolomics 2010, 6, 156–163. [Google Scholar] [CrossRef]
  29. Lokhov, P.G.; Kharybin, O.N.; Archakov, A.I. Diagnosis of lung cancer based on direct-infusion electrospray mass spectrometry of blood plasma metabolites. Int. J. Mass Spectrom. 2011, 309, 200–205. [Google Scholar] [CrossRef]
  30. Lokhov, P.G.; Trifonova, O.P.; Maslov, D.L.; Archakov, A.I. Blood plasma metabolites and the risk of developing lung cancer in Russia. Eur. J. Cancer Prev. 2013, 22, 335–341. [Google Scholar] [CrossRef]
  31. Lokhov, P.G.; Balashova, E.E.; Voskresenskaya, A.A.; Trifonova, O.P.; Maslov, D.L.; Archakov, A.I. Mass spectrometric signatures of blood plasma metabolome for disease diagnostics. Biomed. Rep. 2016, 4, 122–126. [Google Scholar] [CrossRef]
  32. Lokhov, P.G.; Dashtiev, M.I.; Bondartsov, L.V.; Lisitsa, A.V.; Moshkovskiĭ, S.A.; Archakov, A.I. Metabolic fingerprinting of blood plasma for patients with prostate cancer. Biomed. Khim 2009, 55, 247–254. [Google Scholar] [CrossRef]
  33. Lokhov, P.G.; Trifonova, O.P.; Maslov, D.L.; Balashova, E.E.; Archakov, A.I.; Shestakova, E.A.; Shestakova, M.V.; Dedov, I.I. Diagnosing impaired glucose tolerance using direct infusion mass spectrometry of blood plasma. PLoS ONE 2014, 9, e105343. [Google Scholar] [CrossRef] [Green Version]
  34. Balashova, E.E.; Lokhov, P.G.; Maslov, D.L.; Trifonova, O.P.; Khasanova, D.M.; Zalyalova, Z.A.; Nigmatullina, R.R.; Archakov, A.I.; Ugrumov, M.V. Plasma Metabolome Signature in Patients with Early-stage Parkinson Disease. Curr. Metab. 2018, 6, 75–82. [Google Scholar] [CrossRef]
  35. Trifonova, O.P.; Lokhov, P.G.; Archakov, A.I. Metabolic profiling of human blood. Biochem. Suppl. Ser. B Biomed. Chem. 2013, 7, 179–186. [Google Scholar] [CrossRef]
  36. Lokhov, P.G.; Archakov, A.I. Mass spectrometry methods in metabolomics. Biochem. Suppl. Ser. B Biomed. Chem. 2009, 3, 1–9. [Google Scholar] [CrossRef]
  37. Lokhov, P.G.; Arckakov, A.I. Mass spectrometry methods in metabolomics. Biomeditsinskaya Khimiya 2008, 54, 497–511. [Google Scholar] [CrossRef]
  38. Lokhov, P.G.; Maslov, D.L.; Kharibin, O.N.; Balashova, E.E.; Archakov, A.I. Label-free data standardization for clinical metabolomics. BioData Min. 2017, 10, 10. [Google Scholar] [CrossRef] [Green Version]
  39. Viant, M.R.; Kurland, I.J.; Jones, M.R.; Dunn, W.B. How close are we to complete annotation of metabolomes? Curr. Opin. Chem. Biol. 2017, 36, 64–69. [Google Scholar] [CrossRef]
  40. de Jong, F.; Beecher, C.; Raskind, A.; Chamberlain, C.; Guingab, J.; Garrett, T. MetaboNews. Available online: http://www.metabonews.ca/Aug2017/MetaboNews_Aug2017.htm (accessed on 10 August 2017).
  41. Rogers, S.; Scheltema, R.A.; Girolami, M.; Breitling, R. Probabilistic assignment of formulas to mass peaks in metabolomics experiments. Bioinformatics 2009, 25, 512–518. [Google Scholar] [CrossRef] [Green Version]
  42. Silva, R.R.; Jourdan, F.; Salvanha, D.M.; Letisse, F.; Jamin, E.L.; Guidetti-Gonzalez, S.; Labate, C.A.; Vêncio, R.Z.N. ProbMetab: An R package for Bayesian probabilistic annotation of LC-MS-based metabolomics. Bioinformatics 2014, 30, 1336–1337. [Google Scholar] [CrossRef] [Green Version]
  43. Sumner, L.W.; Amberg, A.; Barrett, D.; Beale, M.H.; Beger, R.; Daykin, C.A.; Fan, T.W.-M.; Fiehn, O.; Goodacre, R.; Griffin, J.L.; et al. Proposed minimum reporting standards for chemical analysis Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). Metabolomics 2007, 3, 211–221. [Google Scholar] [CrossRef] [Green Version]
  44. Kennedy, A.D.; Wittmann, B.M.; Evans, A.M.; Miller, L.A.D.; Toal, D.R.; Lonergan, S.; Elsea, S.H.; Pappan, K.L. Metabolomics in the clinic: A review of the shared and unique features of untargeted metabolomics for clinical research and clinical testing. J. Mass Spectrom. 2018, 53, 1143–1154. [Google Scholar] [CrossRef]
  45. Blochberger, A.; Jones, S. Parkinson’s disease clinical features and diagnosis. Clin. Pharm. 2011, 3, 361–366. [Google Scholar]
  46. Gelb, D.J.; Oliver, E.; Gilman, S. Diagnostic criteria for Parkinson disease. Arch. Neurol. 1999, 56, 33–39. [Google Scholar] [CrossRef] [PubMed]
  47. Trifonova, O.P.; Maslov, D.L.; Balashova, E.E.; Urazgildeeva, G.R.; Abaimov, D.A.; Fedotova, E.Y.; Poleschuk, V.V.; Llarioshkin, S.N.; Lokhov, P.G. Parkinson’s Disease: Available Clinical and Promising Omics Tests for Diagnostics, Disease Risk Assessment, and Pharmacotherapy Personalization. Diagnostics 2020, 10, 339. [Google Scholar] [CrossRef] [PubMed]
  48. Trifonova, O.P.; Maslov, D.L.; Balashova, E.E.; Lokhov, P.G. Evaluation of dried blood spot sampling for clinical metabolomics: Effects of different papers and sample storage stability. Metabolites 2019, 9, 277. [Google Scholar] [CrossRef] [Green Version]
  49. Goetz, C.G.; Poewe, W.; Rascol, O.; Sampaio, C.; Stebbins, G.T.; Counsell, C.; Giladi, N.; Holloway, R.G.; Moore, C.G.; Wenning, G.K.; et al. Movement Disorder Society Task Force report on the Hoehn and Yahr staging scale: Status and recommendations. Mov. Disord. 2004, 19, 1020–1028. [Google Scholar] [CrossRef]
  50. Steuer, R.; Kurths, J.; Fiehn, O.; Weckwerth, W. Observing and interpreting correlations in metabolomic networks. Bioinformatics 2003, 19, 1019–1026. [Google Scholar] [CrossRef] [Green Version]
  51. Martens, H.A.; Dardenne, P. Validation and verification of regression in small data sets. Chemom. Intell. Lab. Syst. 1998, 44, 99–121. [Google Scholar] [CrossRef]
Figure 1. The laboratory-developed test (LDT) workflow. Blood plasma samples (or dried blood spots) are collected (A) and transported to the laboratory. In the laboratory, after sample preparation (B) and high-resolution direct mass spectrometry (C), the mass spectra of the blood plasma samples are obtained. The obtained masses of compounds after preprocessing (D) are submitted to the metabolite search block (E) to find metabolite identifiers from Kyoto Encyclopedia of Genes and Genomes database (KEGG) database matching the m/z values. Matched KEGG IDs is submitted to a compound annotation algorithm (F) [17], and the retrieved results are used for the overrepresented metabolite set analysis (G). Finally, overrepresented metabolite sets from an individual are visualized as a metabolite set names cloud, where the font size corresponds to the representation value (score) (H).
Figure 1. The laboratory-developed test (LDT) workflow. Blood plasma samples (or dried blood spots) are collected (A) and transported to the laboratory. In the laboratory, after sample preparation (B) and high-resolution direct mass spectrometry (C), the mass spectra of the blood plasma samples are obtained. The obtained masses of compounds after preprocessing (D) are submitted to the metabolite search block (E) to find metabolite identifiers from Kyoto Encyclopedia of Genes and Genomes database (KEGG) database matching the m/z values. Matched KEGG IDs is submitted to a compound annotation algorithm (F) [17], and the retrieved results are used for the overrepresented metabolite set analysis (G). Finally, overrepresented metabolite sets from an individual are visualized as a metabolite set names cloud, where the font size corresponds to the representation value (score) (H).
Metabolites 11 00014 g001
Figure 2. LDT output for the ‘case-control’ study of Parkinson’s disease (PD). The output is presented for disease-associated metabolite sets (a), pathway-associated metabolite sets (b), metabolite sets associated with abnormal concentrations of metabolites (c), and location-based metabolite sets (d). The top five overrepresented metabolite sets are colored in red.
Figure 2. LDT output for the ‘case-control’ study of Parkinson’s disease (PD). The output is presented for disease-associated metabolite sets (a), pathway-associated metabolite sets (b), metabolite sets associated with abnormal concentrations of metabolites (c), and location-based metabolite sets (d). The top five overrepresented metabolite sets are colored in red.
Metabolites 11 00014 g002aMetabolites 11 00014 g002b
Figure 3. Metabolite set representation scores produced by the laboratory-developed test (LDT) for control subjects and patients with Parkinson’s disease (PD). The rows correspond to the metabolite sets presented in Table 2. The widest bands correspond to a metabolite set representation score of 75.
Figure 3. Metabolite set representation scores produced by the laboratory-developed test (LDT) for control subjects and patients with Parkinson’s disease (PD). The rows correspond to the metabolite sets presented in Table 2. The widest bands correspond to a metabolite set representation score of 75.
Metabolites 11 00014 g003
Figure 4. The laboratory-developed test (LDT) outputs for a person with Parkinson’s disease (PD). The LDT outputs for the overrepresentation pattern in disease-associated metabolite sets (a), pathway-associated metabolite sets (b), metabolite sets associated with abnormal concentrations (c), and location-based metabolite sets (d). [[X]] is a marker corresponding to a representation score of 50. The overrepresented metabolite sets with the top five scores are colored in red.
Figure 4. The laboratory-developed test (LDT) outputs for a person with Parkinson’s disease (PD). The LDT outputs for the overrepresentation pattern in disease-associated metabolite sets (a), pathway-associated metabolite sets (b), metabolite sets associated with abnormal concentrations (c), and location-based metabolite sets (d). [[X]] is a marker corresponding to a representation score of 50. The overrepresented metabolite sets with the top five scores are colored in red.
Metabolites 11 00014 g004aMetabolites 11 00014 g004b
Figure 5. The laboratory-developed test (LDT) outputs for a ‘healthy’ individual. The LDT outputs for the overrepresentation pattern in disease-associated metabolite sets (a), pathway-associated metabolite sets (b), metabolite sets associated with abnormal concentrations (c), and location-based metabolite sets (d). No overrepresented metabolite sets were detected, thus confirming that the individual does not have any deviation in the metabolite composition of their blood. [[X]] is a marker corresponding to a representation score of 50.
Figure 5. The laboratory-developed test (LDT) outputs for a ‘healthy’ individual. The LDT outputs for the overrepresentation pattern in disease-associated metabolite sets (a), pathway-associated metabolite sets (b), metabolite sets associated with abnormal concentrations (c), and location-based metabolite sets (d). No overrepresented metabolite sets were detected, thus confirming that the individual does not have any deviation in the metabolite composition of their blood. [[X]] is a marker corresponding to a representation score of 50.
Metabolites 11 00014 g005
Table 1. Variables associated with this study.
Table 1. Variables associated with this study.
ParameterValue
Detection mass range (m/z)45–900
Detected compound mass peaks (mean ± s.d.)9664 ± 620 1
Masses submitted to metabolite search block14,857
‘Mass peak/metabolite name’ pairs submitted to the annotation algorithm31,724
Mass peaks with annotated compound(s)2741
Unique compound annotations709
1 average ± standard deviation.
Table 2. Criteria for laboratory-developed test diagnostics for Parkinson’s disease.
Table 2. Criteria for laboratory-developed test diagnostics for Parkinson’s disease.
Metabolite SetRepresentation Score Over-RepresentationDiagnostic Performance
ControlsCasesSensitivitySpecificityAccuracy
Disease-associated metabolite sets
Alzheimer’s disease21.953.831.975.071.473.2
Frontotemporal dementia24.453.028.675.067.971.4
Lewy body disease24.453.028.675.067.971.4
Early preeclampsia15.042.527.567.971.469.6
Autosomal dominant polycystic kidney disease12.640.127.460.778.669.6
Pregnancy16.843.626.875.067.971.4
Ulcerative colitis27.553.025.578.660.769.6
Colorectal cancer26.852.025.189.357.173.2
Periodontal disease20.845.825.060.778.669.6
Pancreatic cancer23.247.724.564.371.467.9
Late-onset preeclampsia12.136.724.560.775.067.9
Crohn’s disease25.950.024.271.464.367.9
Schizophrenia25.048.323.264.371.467.9
Eosinophilic esophagitis28.551.122.664.371.467.9
Lipoyltransferase 1 deficiency10.232.722.571.467.969.6
Leukemia9.031.322.375.075.075.0
Maple syrup urine disease16.337.521.167.967.967.9
Perillyl alcohol administration for cancer treatment12.933.921.167.971.469.6
Heart failure3.324.020.753.685.769.6
Rheumatoid arthritis1.620.819.246.489.367.9
Parameters for the whole group of metabolite sets: 78.660.769.6
Pathway-associated metabolite sets
Transcription/translation17.443.225.771.471.471.4
Dopa-responsive dystonia13.434.821.446.482.164.3
Fatty acid elongation in mitochondria13.434.821.446.482.164.3
Hyperphenylalaniemia due to guanosine triphosphate cyclohydrolase deficiency13.434.821.446.482.164.3
Hyperphenylalaninemia due to 6-pyruvoyltetrahydropterin synthase deficiency (PTPS)13.434.821.446.482.164.3
Hyperphenylalaninemia due to DHPR-deficiency13.434.821.446.482.164.3
Long-chain-3-hydroxyacyl-coa dehydrogenase deficiency (LCHAD)13.434.821.446.482.164.3
Pterine biosynthesis13.434.821.446.482.164.3
Segawa syndrome13.434.821.446.482.164.3
Sepiapterin reductase deficiency13.434.821.446.482.164.3
Glutaminolysis and cancer4.122.017.857.175.066.1
Ubiquinone biosynthesis0.116.116.028.692.960.7
Aspartate metabolism4.219.415.257.178.667.9
Canavan disease4.219.415.257.178.667.9
Hypoacetylaspartia4.219.415.257.178.667.9
2-Hydroxyglutric aciduria (D and L Form)1.016.115.125.0100.062.5
4-Hydroxybutyric Aciduria/succinic semialdehyde Dehydrogenase deficiency1.016.115.125.0100.062.5
Glutamate metabolism1.016.115.125.0100.062.5
Homocarnosinosis1.016.115.125.0100.062.5
Hyperinsulinism-hyperammonemia syndrome1.016.115.125.0100.062.5
Parameters for the whole group of metabolite sets: 82.164.373.2
Abnormal concentration-associated metabolite sets
Schizophrenia10.933.622.764.375.069.6
Alcohol intoxication2.724.121.432.196.464.3
Drunk driver2.724.121.432.196.464.3
Pellagra13.434.821.446.482.164.3
Heart failure3.324.020.653.689.371.4
Fabry disease8.018.810.725.089.357.1
Epilepsy0.56.86.321.492.957.1
Heart transplant5.011.06.050.067.958.9
Lesch-Nyhan syndrome0.75.85.175.057.166.1
Dimethylglycinuria4.79.85.132.185.758.9
Menstrual cycle (follicular phase)0.02.72.742.982.162.5
Menstrual cycle (luteal phase)0.02.72.742.982.162.5
Menstrual cycle (midcycle)0.02.72.742.982.162.5
ACTH deficiency, isolated0.02.72.739.385.762.5
Small intestinal bacterial overgrowth0.02.72.728.689.358.9
Crohn’s disease0.02.72.728.689.358.9
HIV and diarrhea0.02.72.728.689.358.9
Glucocorticoid resistance0.02.72.735.789.362.5
Tic disorder0.02.72.735.789.362.5
Thymidine phosphorylase deficiency0.02.72.739.382.160.7
Parameters for the whole group of metabolite sets: 82.164.373.2
Location-based metabolite sets
Testes11.039.428.457.185.771.4
Prostate22.851.228.467.975.071.4
Kidney15.844.128.367.975.071.4
Fibroblasts19.547.828.278.667.973.2
Placenta14.642.127.564.378.671.4
Spleen12.840.127.360.782.171.4
Intestine19.746.426.767.971.469.6
Bladder16.442.726.471.475.073.2
Neuron11.136.825.678.671.475.0
Pancreas15.540.625.171.471.471.4
Gut8.231.623.460.778.669.6
Platelet7.027.620.560.775.067.9
Liver15.935.719.864.371.467.9
Muscle13.432.819.471.478.675.0
Skeletal muscle10.828.017.260.775.067.9
All tissues13.630.416.875.060.767.9
Skin3.520.116.664.364.364.3
Myelin5.318.413.160.785.773.2
Adipose tissue5.917.711.864.371.467.9
Stratum corneum1.112.511.450.089.369.6
Parameters for the whole group of metabolite sets: 71.471.471.4
Parameters for the whole groups of metabolite sets: 89.357.173.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lokhov, P.G.; Maslov, D.L.; Lichtenberg, S.; Trifonova, O.P.; Balashova, E.E. Holistic Metabolomic Laboratory-Developed Test (LDT): Development and Use for the Diagnosis of Early-Stage Parkinson’s Disease. Metabolites 2021, 11, 14. https://doi.org/10.3390/metabo11010014

AMA Style

Lokhov PG, Maslov DL, Lichtenberg S, Trifonova OP, Balashova EE. Holistic Metabolomic Laboratory-Developed Test (LDT): Development and Use for the Diagnosis of Early-Stage Parkinson’s Disease. Metabolites. 2021; 11(1):14. https://doi.org/10.3390/metabo11010014

Chicago/Turabian Style

Lokhov, Petr G., Dmitry L. Maslov, Steven Lichtenberg, Oxana P. Trifonova, and Elena E. Balashova. 2021. "Holistic Metabolomic Laboratory-Developed Test (LDT): Development and Use for the Diagnosis of Early-Stage Parkinson’s Disease" Metabolites 11, no. 1: 14. https://doi.org/10.3390/metabo11010014

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