Metabolic Profiling Methods and Protocols, Second Edition Springer
Metabolic Profiling Methods and Protocols, Second Edition Springer
Metabolic
Profiling
Methods and Protocols
Second Edition
METHODS IN MOLECULAR BIOLOGY
Series Editor
John M. Walker
School of Life and Medical Sciences
University of Hertfordshire
Hatfield, Hertfordshire, UK
Second Edition
Edited by
Ian D. Wilson
Division of Systems Medicine, Imperial College London, London, UK
Editors
Olga Deda Helen G. Gika
Laboratory of Forensic Medicine Laboratory of Forensic Medicine
& Toxicology & Toxicology
Aristotle University of Thessaloniki Aristotle University of Thessaloniki
Thessaloniki, Greece Thessaloniki, Greece
Ian D. Wilson
Division of Systems Medicine
Imperial College London
London, UK
This work is subject to copyright. All rights are solely and exclusively licensed by the Publisher, whether the whole or part
of the material is concerned, specifically the rights of translation, reprinting, reuse of illustrations, recitation,
broadcasting, reproduction on microfilms or in any other physical way, and transmission or information storage and
retrieval, electronic adaptation, computer software, or by similar or dissimilar methodology now known or hereafter
developed.
The use of general descriptive names, registered names, trademarks, service marks, etc. in this publication does not imply,
even in the absence of a specific statement, that such names are exempt from the relevant protective laws and regulations
and therefore free for general use.
The publisher, the authors and the editors are safe to assume that the advice and information in this book are believed to
be true and accurate at the date of publication. Neither the publisher nor the authors or the editors give a warranty,
expressed or implied, with respect to the material contained herein or for any errors or omissions that may have been
made. The publisher remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This Humana imprint is published by the registered company Springer Science+Business Media, LLC, part of Springer
Nature.
The registered company address is: 1 New York Plaza, New York, NY 10004, U.S.A.
This volume is the second edition of Metabolic Profiling: Methods and Protocols, which was
originally published in 2018 as Vol. 1738 in the Methods in Molecular Biology series. The
current edition provides a mixture new and updated methods and protocols from the first
edition. These new contributions are designed to reflect some of the changes in the
development and practice of metabolic phenotyping.
With the chapters in this second edition, we have tried to maintain coverage of impor-
tant topics from the previous edition (e.g., quality control in untargeted metabolic pheno-
typing) while expanding the topics discussed to include other important aspects of this type
of research such as the implementation of quality assurances processes, as discussed by
Kirwan et al (without which any data generated is of doubtful value). Of course, once
good quality data has been acquired, effort has to be put into converting it into useful
information and an important prerequisite for this is careful bio- and chemoinformatic data
analysis as detailed in two contributions by Witting and Rainer and Garcia-Aloy et al.,
respectively.
For those interested in practical aspects of the art of metabolic phenotyping, the volume
contains a number of updated protocols that are focused on untargeted metabolic pheno-
typing (using reversed-phase and ion-pair LC-MS) that can be used for holistic, hypothesis
free, sample analyses. The rapidly maturing field of lipidomics is covered in two contribu-
tions, one on quantitative lipid analysis using supercritical fluid chromatography (SFC-MS)
by Takeda et al and a second, employing UHPLC-TIMS-PASEF-MS, by Denti et al. In
addition there are protocols for metabolite profiling of feces (Deda et al.), the extraction of
plant material for comprehensive metabolome analysis by NMR spectroscopy (Schripsema
and Dagnino), and an LC-ion mobility-MS-based protocol for acquiring metabolome data
on wine and grapes by Sáez et al.
In addition to untargeted methodologies, often used for “biomarker discovery” in the
last few years, we have seen an increasing trend toward the development and application of
targeted metabolite profiling methods. In this volume we have therefore sought to include
protocols that provide access to such methods. The targeted analysis methods described
cover the determination of metanephrines in urine via LC-MS/MS (Thaitumu and Frank), a
“multi-targeted” HILIC- MS/MS method (Virgiliou et al.) and a procedure employing
solid-state high-resolution “Magic Angle Spinning” (HR-MAS) NMR spectroscopy applied
to profiling olive fruit phenolics (Manolopoulou and Spyros).
Metabolic phenotyping (metabonomics/metabolomics) is now well established in all
areas of biological, environmental, and clinical research. In recent years, this has resulted in
an exponential increase in the number of publications, employing both untargeted and
targeted analytical methods, now appearing in journals every year. This enthusiasm for
metabolomics/metabonomics is gratifying to many of those who have watched the field
grow. The rapid expansion in its user base, however, has been much faster than the increase
in the number of trained practitioners of the art who have the necessary expertise in this
metabolomics. Sadly, because of this, many of these newly published articles have suffered
v
vi Preface
from problems such as poor study design, inadequate or unsuitable analytical strategies and
techniques, and subsequent faulty data analysis. Valuable studies are poorly executed,
analyzed, and reported and that is a great shame. The Volume Editors hope that the
contents of the first and second editions of Metabolic Profiling: Methods and Protocols will
provide researchers with pointers toward good practice in metabolic phenotyping.
Preface . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . v
Contributors. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ix
vii
viii Contents
Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 285
Contributors
ix
x Contributors
Abstract
Metabolic profiling continues to develop, and research is now conducted on this topic globally in hundreds
of laboratories, from small groups up to national centers and core facilities. Here we briefly provide a
perspective on the current status and challenges facing metabolic phenotyping (metabonomics/metabo-
lomics) and consider future directions for this important area of biomarker and systems biology research.
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_1,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
1
2 Helen G. Gika et al.
Fig. 1 Publication trends in the areas of genomics/transcriptomics, proteomics, and metabolic phenotyping
from 2003 to 2022 generated from references contained within SCOPUS using the search terms genomics,
proteomics, metabolomics, or metabonomics or metabolic profiling and transcriptomics in the title, abstract,
or keywords. (Search made June 2023)
Metabolic Profiling: A Perspective on the Current Status, Challenges. . . 3
2 Methodology
adjustment). Samples such as cells and tissues are most often ana-
lyzed after extraction into a suitable solvent, drying down, and
redissolution in a 1H NMR compatible solvent. Solid
(or semisolid) samples can be analyzed without extraction using
“magic angle spinning” NMR spectroscopy (although uptake of
this technology in metabolic profiling applications has been lim-
ited) [11]. A noteworthy property of 1H-NMR-spectroscopy-
based methods is that they are generally easily transferable between
laboratories such that, e.g., a one-dimensional 1H NMR spectrum
acquired at 600 MHz in one laboratory should be identical with
one obtained for the same sample, analyzed under the same experi-
mental conditions in another laboratory, irrespective of the
manufacturer.
An illustrative example of the use of 1H NMR spectroscopy to
compare the effects of the liver toxin acetaminophen (APAP) and a
less toxic structural analogue, the meta-isomer AMAP is shown in
Fig. 2. This example shows the 1H NMR spectra obtained for
“aqueous” liver extracts obtained from control (untreated)
AMAP and APAP dosed mice [12]. As these spectra show, it is
not difficult to see the differences in composition of the extracts
resulting from the different dosing regimens, with both drug-
related material and endogenous metabolites clearly present.
MS, like 1H NMR spectroscopy, can also be employed for the
direct analysis of liquid (or gaseous) samples via directly infusing
(DI) suitably prepared samples into the ion source of the mass
spectrometer (DIMS) [13]. DIMS-based analysis, using either
flow injection analysis (FIA), interfaces such as the “nanomate,”
desorption ESI-MS (DESI), or Rapid Evaporative Ionization Mass
Spectrometry (REIMS), can have advantages in terms of simplicity
and speed. However, matrices such as, e.g., tissues, cell extracts,
food, urine, or blood-derived samples can defy simple methods of
characterization because of their sheer complexity and confounding
effects on analyte ionization performance such as ion-suppression/
enhancement. Furthermore, isobaric or isomeric compounds are
difficult to determine confidently if they are present as mixtures
using methods that do not employ a separation step. This has
resulted in the increasing adoption of the use of hyphenated tech-
niques incorporating a chromatographic, electrophoretic, or ion
mobility (IM) separation prior to MS. Thus, LC, GC, IM, and
capillary zone electrophoresis (CZE-MS, or CE-MS) are very
widely used for metabolic phenotyping, with LC-MS (or LC-IM-
MS)-based techniques currently the most popular.
GC-MS-based methods are a priority technology selection for
the analysis of volatiles in samples such as exhaled breath [14] but
can be also used for involatile metabolites following derivatization
[15]. For the GC analysis of involatile analytes, the derivatization
step needs careful optimization as various metabolites can react
with the derivatization reagents at different rates thus resulting in
Metabolic Profiling: A Perspective on the Current Status, Challenges. . . 5
Fig. 2 Representative 1H-NMR spectra of hepatic extract metabolic profiles of the acetaminophen (APAP),
AMAP, and control groups at 1 h. Resonances assigned to drug-related molecules have been colored in
red. Key: APAP/AMAP-G = APAP/AMAP glucuronide; APAP-SG, APAP glutathionyl; APAP-NAC, APAP-N-
acetylcysteinyl; APAP/AMAPNHCOCH3, APAP/AMAP N-acetyl resonance; GSH, reduced glutathione; GSSG,
oxidized glutathione; Phe, phenylalanine; d-3-HB, d-3-hydroxybutyrate; AMP, adenosine monophosphate,
overlapped resonances from glucose/glycogen/maltose labeled. From reference [12] reprinted with
permission
6 Helen G. Gika et al.
Fig. 3 Representation of a 3D mass chromatogram obtained from the reversed-phase analysis of rat urine by
UPLC–TOF-MS. Reprinted from reference [19] with permission
2.2 Data Processing Assuming that the analytical data passes such preliminary scrutiny,
further measures to optimize, e.g., peak alignment, can be per-
formed and the data can be examined for the presence of potential
biomarkers. This part of the process relies heavily on the correct
choice and the correct function of software tools. Some software
may operate as black boxes, not allowing much information on the
theory or the tools used and/or not allowing much freedom in the
selection of parameters. This kind of software also often provides
different levels of multivariate statistics for data analysis, such as
PCA, and offers options for advanced visualization plots. Open-
source software (including web-based data treatment servers) is
10 Helen G. Gika et al.
also available that can be used to perform these tasks, and its use
continues to increase [29]. All such open-source tools necessitate at
least some basic-level knowledge of the use of software language in
R or MATLAB environment, but they offer unparalleled freedom
in optimizing and tailoring the data treatment and data scrutiny
process. They also offer advanced control and functionalities in
developing visualization of the findings and in generating plots,
tables, and illustrations. Such tools can offer impressive outputs;
however, statistical analysis is not yet fool proof, and the famous
Benjamin Disraeli quote concerning statistics is still timely
[30]. Hence, researchers are advised to always verify the validity
of their findings by the use of different statistical analysis tools.
Researchers who are novices to the field are advised to consult
with fellow researchers who are more experienced in the statistical
analysis of metabolomic data. Indeed the concept and the needs of
metabolic phenotyping are different from those applied in statistics
for genomics. Hence, the statistical analysis tools to be selected
and/or their fine-tuning for effective data treatment in metabolo-
mics may vary to a great extent from the treatment of genomic data.
2.3 Metabolite Having completed the data analysis, the structural identification of
Identification peaks that correlate (or appear to be correlated) with whatever
biological problem is being investigated may then be pursued.
However, unequivocal metabolite identification (to level 1 of the
MSI) is not a trivial exercise and is not often achievable [31].
Techniques like 1H NMR spectroscopy or GC-MS provide good
structural information and benefit from large databases of spectral
data that can be search for potential structures. Indeed, in many
cases, positive assignments of identity to the detected metabolites
can be made. However, with LC-MS, it still remains the case that,
while large databases are available that can provide information to
assist in achieving some level of annotation, these do not represent
identifications, and on examination many of the suggestions that
are returned are obviously incorrect having no possible biological
plausibility for the species/sample type being investigated [32]. A
further problem is that even when the putative annotations appear
plausible, it may not be possible to obtain the authentic standards
necessary to unequivocally confirm identity. An increasing con-
cern is that many research groups are outsourcing metabolomic
profiling to central laboratories/commercial providers that do not
disclose information on their proprietary methods and thus provide
only lists of molecules without raw data and methodological pro-
tocols. In such case, the onus is on the recipient of such data to
ensure its validity so that only metabolite identification that is
correct is finally reported [32].
Where a number of candidates for identity are returned from
such searches, the unknowns should be characterized to the extent
possible using the wide range of MS techniques available via
Metabolic Profiling: A Perspective on the Current Status, Challenges. . . 11
3 Conclusions
References
1. Nicholson JK, Lindon JC, Holmes E (1999) 15. Kopka J (2006) Current challenges and devel-
“Metabonomics”: understanding the meta- opments in GC–MS based metabolite profiling
bolic responses of living systems to pathophys- technology. J Biotechnol 124:312–322
iological stimuli via multivariate statistical 16. Moros G, Chatziioannou AC, Gika HG et al
analysis of biological NMR spectroscopic data. (2017) Investigation of the derivatization con-
Xenobiotica 29:1181–1189 ditions for GC–MS metabolomics of biological
2. Fiehn O, Kopka J, Dörmann P et al (2000) samples. Bioanalysis 9:53–65
Metabolite profiling for plant functional geno- 17. Gika HG, Theodoridis GA, Plumb RS et al
mics. Nat Biotechnol 18:1157–1161 (2014) Current practice of liquid
3. Gavaghan CL, Holmes E, Lenz E et al (2000) chromatography–mass spectrometry in meta-
FEBS Lett 484:169–174 bolomics and metabonomics. J Pharm Biomed
4. Dent CE (1952) Lectures on the scientific basis Anal 87:12–25
of medicine, vol 2. Athlone Press, London 18. Rainville PD, Theodoridis G, Plumb RS et al
5. Dalgliesh CE (1956) Two-dimensional paper (2014) Advances in liquid chromatography
chromatography of urinary indoles and related coupled to mass spectrometry for metabolic
substances. Biochem J 64:481–485 phenotyping. TrAC Trends Anal Chem 61:
6. Teranishi R, Mon TR, Robinson AB et al 181–191
(1972) Gas chromatography of volatiles from 19. Theodoridis GA, Gika HG, Want EJ et al
breath and urine. Anal Chem 44:18–20 (2012) Liquid chromatography–mass spec-
7. Pauling L, Robinson AB, Teranishi R et al trometry based global metabolite profiling: a
(1971) Quantitative analysis of urine vapor review. Anal Chim Acta 711:7–16
and breath by gas-liquid partition chromatog- 20. Lioupi A, Marinaki M, Virgiliou C, Begou O,
raphy. Proc Natl Acad Sci USA 68:2374–2376 Gika H, Wilson I, Theodoridis G (2023) Prob-
8. Robinson AB, Pauling L (1974) Techniques of ing the polar metabolome by UHPLC-MS.
orthomolecular diagnosis. Clin Chem 20:961– TrAC Trends Anal Chem. https://doi.org/
965 10.1016/j.trac.2023.117014
9. Scott CD, Chilcote DD, Lee NE (1972) Cou- 21. Plumb RS, Isaac G, Rainville PD, Hill J,
pled anion and cation-exchange chromatogra- Gethings LA, Johnson KA, Lauterbach J, Wil-
phy of complex biochemical mixtures. Anal son ID (2022) High throughput UHPLC-MS-
Chem 44:85–89 based lipidomics using vacuum jacketed col-
umns. J Proteome Res 21:691–701
10. Scott CD, Chilcote DD, Katz S et al (1973)
Advances in the application of high resolution 22. van de Velde B, Guillarme D, Kohler I (2020)
liquid chromatography to the separation of Supercritical fluid chromatography – mass
complex biological mixtures. J Chromatogr spectrometry in metabolomics: past, present,
Sci 11:96–100 and future perspectives. J Chromatogr
B. https://doi.org/10.1016/j.jchromb.2020.
11. Lenz EM, Wilson ID (2007) Analytical strate- 122444
gies in metabonomics. J Proteome Res 6:443–
458 23. Gika HG, Wilson ID, Theodoridis GA (2014)
LC–MS-based holistic metabolic profiling.
12. Kyriakides M, Maitre L, Stamper BD et al Problems, limitations, advantages, and future
(2016) Comparative metabonomic analysis of perspectives. J Chromatogr B 966:1–6
hepatotoxicity induced by acetaminophen and
its less toxic meta-isomer. Arch Toxicol 90: 24. Michopoulos F, Whalley N, Theodoridis G
3073–3085 et al (2014) Targeted profiling of polar intra-
cellular metabolites using ion-pair-high perfor-
13. Han J, Danell RM, Patel JR et al (2008) mance liquid chromatography and-ultra high
Towards high-throughput metabolomics performance liquid chromatography coupled
using ultrahigh-field Fourier transform ion to tandem mass spectrometry: applications to
cyclotron resonance mass spectrometry. Meta- serum, urine and tissue extracts. J Chromatogr
bolomics 4:128–140 A 1349:60–68
14. Theodoridis G, Gika HG, Wilson ID (2011) 25. Mamani-Huanca M, Villaseñor A, Gonzalez-
Mass spectrometry-based holistic analytical Riano C, López-López A, López-
approaches for metabolite profiling in systems Gonzálvez A, Barbas C (2023) Capillary elec-
biology studies. Mass Spectrom Rev 30:884– trophoresis mass spectrometry-based untar-
906 geted metabolomics to approach disease
diagnosis. TrAC Trends Anal Chem:117049
14 Helen G. Gika et al.
26. Zhou Z, Luo M, Chen X et al (2020) Ion 29. Xia J, Wishart DS (2010) MetPA: a web-based
mobility collision cross-section atlas for metabolomics tool for pathway analysis and
known and unknown metabolite annotation visualization. Bioinformatics 26:2342–2344
in untargeted metabolomics. Nat Commun 3 0 . https :// en.w ikiped i a . o r g / w i k i / L i e s , _
11:4334. https://doi.org/10.1038/s41467- damned_lies,_and_statistics
020-18171-8 31. Kodra D, Pousinis P, Vorkas PA et al (2022) Is
27. Gika HG, Theodoridis GA, Wingate JE et al current practice adhering to guidelines pro-
(2007) Within-day reproducibility of an posed for metabolite identification in LC-MS
HPLC- MS-based method for metabonomic untargeted metabolomics? A meta-analysis of
analysis: application to human urine. J Prote- the literature. J Proteome Res 21:590–598
ome Res 6:3291–3303 32. Theodoridis G, Gika H, Raftery D,
28. Gika HG, Zisi C, Theodoridis G, Wilson ID Goodacre R, Plumb RS, Wilson ID (2023)
(2016) Protocol for quality control in meta- Ensuring fact-based metabolite identification
bolic profiling of biological fluids by U in liquid chromatography–mass spectrometry-
(H)PLC-MS. J Chromatogr B 1008:15–25 based metabolomics. Anal Chem 95:3909–
3916
Chapter 2
Abstract
Metabolic profiling (untargeted metabolomics) aims for a global unbiased analysis of metabolites in a cell or
biological system. It remains a highly useful research tool used across various analytical platforms. Incre-
mental improvements across multiple steps in the analytical process may have large consequences for the
end quality of the data. Thus, this chapter concentrates on which aspects of quality assurance can be
implemented by a lab in the (pre-)analytical stages of the analysis to improve the overall end quality of their
data. The scope of this chapter is limited to liquid-chromatography-mass spectrometry (LC-MS)-based
profiling, which is one of the most widely utilized platforms, although the general principles are applicable
to all metabolomics experiments.
Key words Quality assurance, Quality control, Quality management, Validation, Robust method,
Metabolomics, Untargeted, Mass spectrometry, NMR, Data quality
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_2,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
15
16 Jennifer A. Kirwan et al.
2 Background
2.1 Benefits of The major reasons for good quality management in a lab are to
Quality Management improve experimental quality and robustness and to have good
in the Lab traceability of what was done, including any deviations from pro-
tocols or unexpected occurrences. This in turn should improve the
repeatability of experimental work, including in other labs. In
addition, multiple studies across different industries have shown
that implementing a quality standard such as ISO9001 leads to
better financial performance [14]. While the current ISO quality
standards are referred to here for guidance, certification is not
required to implement the quality management guidelines. Good
quality management will be beneficial in any lab, and each lab
should implement practices that best match their size and budget.
In generic terms, the ISO9000 series defines the benefits of a
good quality management system (QMS) [15] as follows:
1. Defining, improving, and controlling processes.
2. Reducing waste.
3. Preventing mistakes.
4. Lowering costs.
5. Facilitating and identifying training opportunities.
6. Engaging staff.
7. Setting organization-wide direction.
Having a defined QMS communicates a readiness to produce
consistent results to both staff and “customers.” The term “cus-
tomer” will be used in the broadest sense throughout this chapter.
Most QMS were developed for industrial objectives where profit
and customer satisfaction are key objectives. In academic, govern-
mental, and nonprofit labs, the customer may be either a collabo-
rator, a program office and their stakeholders, or the individual
themselves requiring high quality data for their research project.
2.2 Defining Needs Before a QMS can be implemented, the requirements of all stake-
and Expectations holders need to be determined. In a metabolomics context, the
requirement of an international company offering pay-for-service
18 Jennifer A. Kirwan et al.
2.3 Risks and A formal QMS communicates a willingness to address both risks
Opportunities and opportunities associated with the organizational or individual
laboratory’s objectives. Risk-based thinking encourages a lab to
reduce risk to the overall results and to seek and implement oppor-
tunities to improve. It therefore involves risk management, i.e., to
think in advance of which steps, in which processes, could go wrong
and to plan in advance to mitigate or avoid these risks and their
effects. In addition, it also uses this assessment process, and the data
monitoring inherent to QA and QC processes, to identify areas
where changes in the process could improve the results in some
way, either through time, cost, or quality improvements. Risk-
based thinking and risk management are central to the aims of
QM and should be a major focus of any metabolomics laboratory.
The two major risks in a metabolomics lab tend to center around
people and instruments. People are our greatest asset, and knowl-
edge loss from the lab is difficult to replace. Instruments are our
most expensive laboratory investments and are central to the qual-
ity of the results, and the risk management of both project time and
quality is ultimately dependent on keeping them operating at their
best. People management starts with good leadership and a clear
vision of how to recruit, retain, and promote the best staff. In
academic environments, the individual lab can do much—even
within an unstable financial environment—to make their key mem-
bers of staff feel valued and appreciated. Readers are referred to the
wide variety of resources on people management and leadership,
details of which are outside the scope of this chapter. Regarding
instruments, thorough research when purchasing new instrumen-
tation, and a structured program of training staff in best practice
maintenance will help to mitigate this risk.
Good risk management includes detailing the critical risks in
each process: running out of the correct consumables, contami-
nated solvents, contamination from plasticizers or other chemicals
from consumables such as Eppendorf tubes, human error (e.g.,
mixing up samples), methods not being robust, etc. Risk manage-
ment strategies to minimize such occurrences include, but are not
limited to, the following:
Quality Assurance in Metabolomics and Metabolic Profiling 19
2.4 Quality Science works on the principle of setting out objectives and finding
Objectives and ways to achieve them. Quality assurance is no different, and a lab
Planning for Success should invest time in defining the quality objectives it is seeking to
meet and how they correspond with the lab’s wider strategic goals,
aims, and budget. Setting quality objectives is similar to setting out
scientific objectives. ISO 9001 particularly emphasizes that quality
objectives should be consistent with the quality policy, measurable,
and relevant to the work and needs of the lab and customers and
consider any applicable requirements.
Implementing a robust QMS requires good leadership, a clear
vision, adequate support and resources to implement the quality
plan, and identifying who the key individuals are for each stage and
how it should be monitored. The timelines and evaluation criteria
also need to be identified. Each individual in the lab should have a
clear understanding of what the lab is trying to achieve and how
their own part plays a role. To this end, formalized protocols should
20 Jennifer A. Kirwan et al.
3 Table of Definitions
3.1 Definitions Definitions are based on ISO/IEC guide 99:2007. Readers are
directed to their websites [18] for a full list. Specific terms listed
here are discussed in ISO 17025 (2017) [3].
Decision rule: Accounts for the measurement uncertainty in a
measurement when trying to show that a test result conforms
with a specified requirement.
Verification: Objective evidence that a given method fulfils speci-
fied requirements. Verification is not synonymous with
validation.
Validation: Where a verification proves that the specified require-
ments are fit for intended use. While both the FDA and EMA
have well recognized guidelines [8, 9] for the validation of
analytical methods for single compounds, validation is more
complex where multiple compounds or untargeted metabolo-
mic methods are concerned. In addition, a method is only
validated for the biological matrix on which the original valida-
tion was conducted, especially for accurate absolute quantifica-
tion [19, 20].
Quality Assurance in Metabolomics and Metabolic Profiling 21
4.1 Organizational The metabolomics analytical pipeline is both long and complex
Context (Fig. 1) and increasingly involves multiple persons, from multiple
scientific disciplines. Technical variation can occur at all stages of
the pipeline, either from biases in sample collection (e.g., smoking,
sex imbalances, different collection clinics, etc.) to
non-standardized protocols or multiple instruments being used or
due to random variation.
While the analytical chemist may be involved in both sample
collection and data processing, their principal role focuses on the
sample handling and chemical analysis which occurs after samples
have been collected, processed, and stored. This chapter will focus
therefore on the laboratory analytical factors which the analyst can
control. While we recognize that multiple instrumentation and
methods are used for metabolomics, we will concentrate here on
mass spectrometry-based analyses. For readers interested in the
sample collection stages, they are referred to the wide variety of
literature on the subject [21–23]. Indeed, sample collection, pro-
cessing, and storage are worth a chapter in their own right, and
these steps need to be done well for subsequent steps to be success-
ful. The bioinformatical and data analysis steps also require strin-
gent quality management. A good understanding of how individual
processes build into a system, and the management of this opera-
tion, contributes to a laboratory’s overall effectiveness and
efficiency.
Fig. 1 The metabolomics analytical pipeline. Steps in blue are those where the analytical chemist tends to be
most involved
22 Jennifer A. Kirwan et al.
Fig. 2 Diagram detailing the sequence of metabolomics QA/QC processes and their interactions. (Adapted
from ISO 9001)
4.3 Monitoring and ISO 9001 and ISO/IEC 17025:2017(E) both specifically mandate
Measuring of that organizations monitor and measure the resources they have or
Resources require to ensure valid results when monitoring. This includes that
they are suitable for purpose, properly maintained to ensure this
suitability, and that they are regularly calibrated. In a metabolomics
lab, this starts with an appropriate laboratory environment, i.e.,
with suitable space, gases, electricity, and other resources to main-
tain a high analytical quality. Critically, maintaining and retaining
good records is seen as an essential part of this process by ISO 9001
and is also an essential part of good scientific practice. To this end,
this chapter provides some essential checklists and templates that
can assist the metabolomics lab in record maintenance.
Quality Assurance in Metabolomics and Metabolic Profiling 23
4.3.2 Maintenance ISO 9001 requires organizations to ensure valid results when mon-
Schedules itoring. In practice, this means regular testing and calibration of all
equipment in the analytical laboratory involved in the monitoring
process. This includes, but is not limited to, pipettes, centrifuges,
balances, centrifugal evaporators, cold rooms, refrigerators, and
freezers, as well as the larger analytical instruments themselves
(chromatography instrumentation, MS, NMR, CE, etc.). Monitor-
ing includes calibration, verification, or both, and this should be
done to international measurement standards whenever possible. If
an international measurement standard is not available, the stan-
dard against which an instrument is being assessed must be formally
documented and retained as part of good quality assurance.
Labs that are already accredited under various national clinical
laboratory accreditation schemes, especially where they conform to
ISO 15189 (2022) (Medical Laboratories–Requirements for qual-
ity and competence), are probably already meeting these require-
ments [27]. However, even for labs that do not seek accreditation,
regular documented maintenance and checking of all systems are a
must for ongoing confidence that results are both precise and
accurate. For more advanced analytical instrumentation, checklists
can serve well (Table 1). At least some of these checks can be carried
out by trained in-house personnel (e.g., calibration accuracy of
pipettes), although we would recommend for safety reasons that
24 Jennifer A. Kirwan et al.
Table 1
A template for a checklist for the regularly scheduled maintenance of a Q-Exactive (QE) mass
spectrometer coupled with an Accela LC system
Monthly maintenance:
Date Initials Date Initials Date Initials Date Initials Date Initials
Date Initials Date Initials Date Initials Date Initials Date Initials
Annual maintenance:
Change fore pump oil (6-12 months) Replace autosampler valve rotor
Check anti-suckback valve on fore pump seals (6-12 months)
Clean fore pump fan guard Replace DLW holding loop
Clean and grease the rolls and rails
of cold stack
(continued)
Quality Assurance in Metabolomics and Metabolic Profiling 25
Table 1
(continued)
As needed:
Date Initials Date Initials Date Initials Date Initials Date Initials
Date Initials Date Initials Date Initials Date Initials Date Initials
Date Initials Date Initials Date Initials Date Initials Date Initials
Date Initials Date Initials Date Initials Date Initials Date Initials
Date Initials Date Initials Date Initials Date Initials Date Initials
This can be adapted to most high-resolution LC-MS instruments with minimal changes
26 Jennifer A. Kirwan et al.
4.3.4 Documentation and There are two reasons to keep good documentation:
Audit Trail
(i) Most importantly, to have good, documented records of every
process procedure, results, and improvements. This provides
an overview of how the lab is performing, both against its own
quality objectives and also in relation to international labs with
similar aims. It is critical for monitoring where processes may
require optimization to maintain or improve quality.
(ii) To comply with good scientific practice and the international
standard laid out in ISO 9001.
All formal documentation created within a lab should have
appropriate identification and description including, but not lim-
ited to, title, date, and author or reference number. In many cases, a
version number may also be appropriate. Control of documenta-
tion, including who has access to it, and any data protection issues
associated with confidential documentation are also important
[2]. Version control of documents is important, and it is strongly
recommended that fully audited software (e.g., GitLab) is used to
maintain the most important documentation such that no changes
can be implemented without them being recorded. Audit trails,
such as those implemented on many vendor proprietary software
programs, enable formal records of who has made which changes at
which timepoints and are easy to implement as a formal record.
Documentation, both external and internal, may be required both
to show that certain conformity standards have been met and for
certification purposes. It is important that such documentation is
stored in a way that does not allow unintentional alterations.
While often considered synonymous, there is a technical differ-
ence between documents and records which determines how they
should be treated. ISO9000:2015(E) [1] defines a document as
both the information and the medium on which it is stored. Exam-
ples of documents include specifications (including the quality
management plan), documented information, and formal records.
It can include written policies, manufacturers’ technical
28 Jennifer A. Kirwan et al.
Table 2
A template for a lab book checklist
(continued)
Quality Assurance in Metabolomics and Metabolic Profiling 31
Table 2
(continued)
□ Were samples dried or otherwise stabilized? What happened to samples after drying (i.e., immediately
analyzed or stored—If so at what temperature)?
Sample analysis
□ Details and results of instrument SST including method and chemical composition of SST. Full details
of how the SST standard was prepared should be included.
□ Exact details of samples used, to include how they have been stored up to now, and for how long, e.g..
the 24 Baker project human serum samples prepared on August 2, 2023, and stored at -80C for 1 week
were taken from freezer for analysis. Full sample details are provided on that date.
□ Which instrument has been used?
□ Which method was used?
□ Which column and batch number?
□ For LC-MS: which buffers were used, how were they made up, and when?
□ For GC-MS: which liner type and gas was used?
□ For other analysis type, please denote specific details.
□ Which chromatography method was used: name, gradient, injection amount, and all details that allow
the experiment to be repeated on another instrument of the same specification?
□ Which MS method was used: all details of method that allow the experiment to be repeated on another
instrument of the same specification?
Data preprocessing
□ Initial check of raw data; were any rejected due to data quality concerns? Which and why?
□ Which software was used (including version number) to carry out which steps in data preprocessing?
Inputs and parameters used should be noted, as is the location of the original data and subsequent data
steps.
□ Are all data preprocessing steps clearly named, and is the order in which they were carried out clear?
Are outputs of each step clearly recorded and easily accessible to authorized individuals?
□ For steps that require a reference dataset, e.g., pqn normalization, is the data that was used as the
reference dataset clear?
□ What filtering steps were carried out on the data; which parameters were used?
□ Which quality control checks have been undertaken on the data? Which samples were excluded as a
result?
Data analysis
□ Is a data analysis plan in place before starting?
(NB: Certain datasets carry a lot of information in what is not there, e.g., there may be biological
information in the missing value patterns seen between genetic knockouts and their wild-type
counterparts; consider this in the data analysis plan)
□ Which dataset is being used for which data analytical step?
□ How have missing values been handled for each data analytical step (e.g., imputed with knn—k = 5 or
statistical analysis set to ignore missing values, etc.)?
(continued)
32 Jennifer A. Kirwan et al.
Table 2
(continued)
Table 3
A template for an MS untargeted measurement suitable for use and adaptation across a range of
mass spectrometers
Sample list
Machine details:
If targeted MSMS is used, include list of targets and collision energies, or link to list of targets.
A text copy of the method should be kept in an appropriate place in the lab folder.
Liquid chromatography
Buffer C N/A
(continued)
34 Jennifer A. Kirwan et al.
Table 3
(continued)
Maintenance:
SST details
Date run
Composition of SST
Method used
Buffer 1
Buffer 2
(SST and calibration results should also be entered on the respective instrument log tables for
ongoing monitoring)
Column: VisionHT C18 Basic; L × I.D. 150 mm × 4.6 mm, 3 μm particle size
(Dr. Maisch) part number 230189
Column: N/A
Chip: N/A
During measurement:
no problems reported:
Injection problems:
Spray problems:
Mass inaccuracy:
Column problems:
Leaks:
Other:
Quality Assurance in Metabolomics and Metabolic Profiling 35
5. Additional Documentation
Additional documentation is any documentation that acts to
improve QA and can include, for example, tracking of projects’
metrics over time to see how interventions improve results,
laboratory management information systems (LMIS) or
freezer records to track where samples are stored and how
they have been processed so far, project management tools to
track project progress, etc. Tracking metrics such as number of
signals, number of missing values, and relative standard devia-
tions of QC samples across multiple projects can be particularly
useful for monitoring how incremental changes over time may
be affecting results. Project management tools have proved
useful in our hands in monitoring progress of projects and
better time management of multiple projects leading to better
efficiency.
4.3.5 Staff Training Training An individual carrying out a process must have the com-
(Including Records) petence to undertake this process. The leader of a lab, and ulti-
mately the organization, needs to monitor this.
has been achieved, and how this will be documented. Who requires
training and how the training lessons will be broken down and
structured are also important. The trainer is ideally both competent
and confident in their own knowledge of the instrument or process
that they are teaching. A checklist is useful, both to ensure that
training is systematic and so that both the trainer and the trainee are
confident that all relevant knowledge has been passed on for that
level (see Table 4). In many cases, it is helpful to have more than one
trainer training an individual since different individuals can stress
different aspects when training. It can also highlight any differences
in practices in the trainers which can then be addressed.
In the area of quality assurance, improving performance in
highly routine (potentially boring) tasks should also not be
neglected, e.g., pipetting, making of standard solutions and buf-
fers, and weighing. The idea of marginal gains is well known in both
sport and surgery and is extremely important in the field of quality
assurance [43–45]. It can often highlight differences in practice
(including health and safety practice) and encourage dialogue.
Indeed, a culture of regularly refreshing these skills for all lab
members in a fun training session encourages both teamwork and
communication. Assessments of such skills should be done in ways
that encourage confidence and for individuals to seek help where
required, e.g., anonymized competitions. Since many labs also have
robots, the robot should also be assessed on the accuracy, precision,
and repeatability of its pipetting, including across different solvent
types.
Table 4
A template for recording training of new lab staff
Table 5
A template for a checklist for an average analytical experiment doing metabolomics
Aims Define wider study aims □ Do all experimentalists understand the wider aims of the study?
□ Have the major objectives and biological questions been defined?
Scope Define scope of □ Are the planned experiment(s) sufficient to answer the biological question(s)?
metabolomics □ Are there sufficient resources available?
□ Are the wider processes well aligned with the metabolomics (e.g., sampling
protocols)?
□ Have appropriate controls been discussed?
□ Are sufficient number of samples planned for the required statistical power?
Collection Sample collection □ Have samples been collected using a metabolomics-specific protocol?
□ Have samples been stored using a metabolomics-specific protocol?
Batch length □ Will the batch length compromise the integrity of the samples?
Pilot dilution series □ Has the sample matrix behavior been tested?
□ Are the samples appropriately diluted for the study?
□ Are different classes tested separately to avoid bias?
Pilot internal standard □ Has the optimal ratio of internal standard to sample been tested?
test □ Are the correct solvents being used to dissolve the internal standards?
□ Are the best concentrations of internal standards being used?
Dilution series of QCs □ Is a dilution series of pooled QCs planned for every analytical batch?
Preanalytical Equipment condition □ Has all equipment been calibrated and properly maintained?
□ Are appropriate vials being used for the experiment (contaminants, size,
batch, etc.)?
□ Are samples kept at constant low temperature (e.g., with ice bath)?
□ What type of pipette is required for accurate pipetting?
□ Which solvents are required and at which grade?
□ Are particular solvent brands preferred?
□ Which sources and grades of water are required?
(continued)
Quality Assurance in Metabolomics and Metabolic Profiling 41
Table 5
(continued)
Analysis Remove samples from - □ Have the samples warmed sufficiently to prevent condensation from
80 °C freezer affecting them?
□ Have the samples sufficiently equilibrated to operating temperature?
Reconstitution □ Have the internal standards been correctly made up and added?
□ Is the reconstituted internal standard sufficiently new?
□ Is an appropriate reconstitution solvent used?
Mobile phases □ Are the mobile phases sufficiently new?
□ Have the mobile phases been appropriately sonicated?
□ Have the mobile phases been appropriately filtered?
□ Is the pH at an appropriate level?
□ Are the solvents sufficiently protected from light?
Laboratory environment □ Is the lab room temperature within an appropriate range?
□ Is the lab relative humidity within an appropriate range?
Autosampler temperature □ Has the autosampler reached the correct set point prior to loading samples?
Chromatographic □ Is a pre-column correctly installed (if being used)?
column □ Has the column(s) been appropriately conditioned?
□ Is the column(s) being used in an appropriate temperature and pH operating
range?
□ Is the column temperature being maintained in an appropriate range?
□ Has the column been monitored for degradation?
Flow rate and gradients □ Have the flow rates and gradients been optimized for reproducibility?
□ Are the correct flow rates and gradients being applied as stipulated in the
method protocol?
Injection volume □ Are injection volumes accurate and reproducible?
□ Are injection volumes within the autosampler’s working specifications?
□ Has the injection volume been optimized for mass loaded on column, sharp
chromatography, etc.?
(continued)
42 Jennifer A. Kirwan et al.
Table 5
(continued)
analytical chemist is most involved (see Fig. 1). Full details of the
steps covered in the checklist, including details of key quality assur-
ance points, the reasoning behind them, and actions to be taken to
improve the steps, are provided below.
Transfer
The aims of quality assurance in transfer are to ensure that the
samples are delivered in a state fit for the intended purpose; that
the number, type, and identity of samples are not in doubt; and that
any concerns about the transfer process are recorded for proper
interpretation of any data collected from these samples. In addition,
there may be a legal necessity when taking possession of physical
samples from other laboratories that tracking and ownership of
these samples is clear. Useful points to specifically check are as
follows:
(i) Number of delivered samples and that they correspond to
what was previously agreed.
(ii) That there are sufficient numbers of samples per statistical
class. Where numbers deviate from the original plan, it is
good to record in writing what the potential implications of
this are for the aims and objectives of the study and to com-
municate this to the important stakeholders before experi-
mental work starts.
(iii) Are all samples useable. Samples sometimes have less volume
than recorded, particularly after long periods of storage where
sublimation can occur. Certain sample types are also suscepti-
ble to other issues, e.g., blood products sometime have large
amounts of fibrous material which can reduce the effective
volume available for use (and may also selectively bind specific
metabolites affecting results).
(iv) Are samples traceable: each sample should ideally come with a
full sample history including temperature, storage, and other
tracking information.
(v) Have samples been entered/updated on the tracking system;
this contributes to good record keeping and ensures the integ-
rity of the tracking system in point (iv).
(vi) Have field blanks been collected; this enables the origin of
contaminants to be properly assessed.
The specific metabolomics aims of the study should have been
defined within the metabolomics scope.
44 Jennifer A. Kirwan et al.
6 Training
7 Preanalytical
8 Instrument Setup
Analysis
Analysis is a multiple step procedure which is highly dependent on
the analytical procedure to be followed. Like sample extraction, all
steps should have been previously validated to ensure that they are
fit for the purpose for the matrix of interest. Table 5 gives some
examples of steps in this procedure including:
(i) Remove samples from -80 freezer or place of storage; allow-
ing to equilibrate to working temperature will reduce the risk
of condensation affecting samples and reduce batch effects.
(ii) Reconstitution: reconstituting in an appropriate solvent mix,
with an appropriate mix of internal standards, at an appropri-
ate dilution is essential to ensuring that metabolites of inter-
est are properly and fully dissolved and can be fully analyzed
as part of the experiment.
(iii) Mobile phases should be made up freshly, filtered, and soni-
cated to reduce air or particulate matter and pH checked.
They should also be protected from light to reduce the risk of
polymerization (a particular issue with acetonitrile). Where
instruments are used for multiple analytical methods, espe-
cially where buffers are used, an extra step to ensure that
crystallization within the instrument lines does not occur if
solvents are changed out can save large amounts of time in
troubleshooting instruments, e.g., the crystallization limit of
ammonium acetate at 95% is 10 mM.
(iv) Lab temperature and humidity monitoring ensures that
instruments are operating at their optimum. Unexpected
changes in these values also enable action to be taken quickly
before too many analyses need to be repeated.
(v) Autosampler temperatures should also be monitored regu-
larly to ensure that the autosampler temperature control is
working and sample integrity can be maintained throughout
analysis. Most modern instruments have an automated mon-
itoring system for this.
(vi) Column integrity: use of a pre-column reduces long-term
costs by increasing the time between full column changes.
Both pre-column and column should be fully conditioned
48 Jennifer A. Kirwan et al.
Disclaimer
The views expressed in this book chapter are those of the authors
and do not necessarily represent the views or policies of the
U.S. Environmental Protection Agency.
Acknowledgments
References
(2020) Dissemination and analysis of the qual- metrology—basic and general concepts and
ity assurance (QA) and quality control associated terms (VIM)
(QC) practices of LC–MS based untargeted 19. Raposo F, Barceló D (2021) Challenges and
metabolomics practitioners. Metabolomics strategies of matrix effects using
16:1–16 chromatography-mass spectrometry: an over-
6. Hibbert DB, Korte E-H, Örnemark U (2021) view from research versus regulatory view-
Metrological and quality concepts in analytical points. TrAC Trends Anal Chem 134:116068
chemistry (IUPAC recommendations 2021). 20. Chamberlain CA, Rubio VY, Garrett TJ (2019)
Pure Appl Chem 93:997–1048 Impact of matrix effects and ionization effi-
7. Broadhurst D, Goodacre R, Reinke SN, ciency in non-quantitative untargeted metabo-
Kuligowski J, Wilson ID, Lewis MR, Dunn lomics. Metabolomics 15:1–9
WB (2018) Guidelines and considerations for 21. Dagher G (2022) Quality matters: interna-
the use of system suitability and quality control tional standards for biobanking. Cell Prolif
samples in mass spectrometry assays applied in 55:e13282
untargeted clinical metabolomic studies. Meta- 22. González-Domı́nguez R, González-Domı́n-
bolomics 14:1–17 guez Á, Sayago A, Fernández-Recamales Á
8. Administration UFaD (2018) Bioanalytical (2020) Recommendations and best practices
method validation guidance for industry. In: for standardizing the pre-analytical processing
Center for Drug Evaluation and Research of blood and urine samples in metabolomics.
CfVM (ed) US Food and Drug Administration Meta 10:229
9. European Medicine Agency (2011). Guideline 23. Biais B, Bernillon S, Deborde C, Cabasson C,
on bioanalytical method validation. In Rolin D, Tadmor Y, Burger J, Schaffer AA,
(CHMP) CfMPfHU (ed) Moing A (2012) Precautions for harvest, sam-
10. Consortium MQAaQC 2023, posting date. pling, storage, and transport of crop plant
[Online] metabolomics samples. In: Plant metabolo-
11. Dudzik D, Barbas-Bernardos C, Garcı́a A, Bar- mics: methods and protocols, pp 51–63
bas C (2018) Quality assurance procedures for 24. Huber L (1996) Quality assurance and instru-
mass spectrometry untargeted metabolomics. mentation. Accred Qual Assur 1:24–34
A review. J Pharm Biomed Anal 147:149–173 25. Lam H (2004) Procurement, qualification, and
12. Dumas-Mallet E, Button KS, Boraud T, calibration of laboratory instruments: an
Gonon F, Munafò MR (2017) Low statistical overview. In: Analytical method validation and
power in biomedical science: a review of three instrument performance verification, pp
human research domains. R Soc Open Sci 4: 139–151
160254 26. Bedson P, Sargent M (1996) The development
13. Barraza F, Arancibia M, Madrid E, Papuzinski and application of guidance on equipment
C (2019) General concepts in biostatistics and qualification of analytical instruments. Accred
clinical epidemiology: random error and sys- Qual Assur 1:265–274
tematic error. Medwave 19 27. Grochau IH, Leal DKB, ten Caten CS (2020)
14. Manders B, de Vries H, Blind K (2013) Acad- European current landscape in laboratory
emy of management proceedings, p 12255 accreditation. Accred Qual Assur 25:303–310
15. Quality AAf 2023, posting date. What is a 28. Hunt CJ (2019) Technical considerations in
Quality Management System (QMS). [Online] the freezing, low-temperature storage and
16. Dunn WB, Broadhurst D, Begley P, Zelena E, thawing of stem cells for cellular therapies.
Francis-McIntyre S, Anderson N, Brown M, Transfus Med Hemother 46:134–150
Knowles JD, Halsall A, Haselden JN (2011) 29. Nussbeck SY, Weil P, Menzel J, Marzec B,
Procedures for large-scale metabolic profiling Lorberg K, Schwappach B (2014) The labora-
of serum and plasma using gas chromatography tory notebook in the 21st century: the elec-
and liquid chromatography coupled to mass tronic laboratory notebook would enhance
spectrometry. Nat Protoc 6:1060–1083 good scientific practice and increase research
17. Anderson HE, Santos IC, Hildenbrand ZL, productivity. EMBO Rep 15:631–634
Schug KA (2019) A review of the analytical 30. Deutsche Forschungsgemeinschaft e.V. (2022)
methods used for beer ingredient and finished 2022, posting date. Guidelines for Safeguard-
product analysis and quality control. Anal ing Good Research Practice (English language
Chim Acta 1085:1–20 translation). Deutsche Forschungsge-
18. Commission IOfSIE (2007) ISO/IEC Guide meinschaft e.V. 1.1. [Online]
99:2007(en) International voacbulary of
Quality Assurance in Metabolomics and Metabolic Profiling 51
31. Wilkinson MD, Dumontier M, Aalbersberg IJ, community questionnaire. Metabolomics 13:
Appleton G, Axton M, Baak A, Blomberg N, 1–6
Boiten J-W, da Silva Santos LB, Bourne PE 39. Charoo NA, Khan MA, Rahman Z (2022)
(2016) The FAIR guiding principles for scien- Data integrity issues in pharmaceutical indus-
tific data management and stewardship. Sci try: common observations, challenges and
Data 3:1–9 mitigations strategies. Int J Pharm:122503
32. Kumar R, Venkatesh K (2022) Proceedings of 40. World Health Organization (2019) Guideline
international conference on deep learning. In on data integrity. WHO Drug Inf 33:773–793
Computing and Intelligence: ICDCI 2021, pp 41. Perkel JM (2019) Ways to avoid a data-storage
687–698 disaster. Nature 568:131–132
33. Hart EM, Barmby P, LeBauer D, 42. Sharma K, Shingatgeri VM, Pal S (2021) Role
Michonneau F, Mount S, Mulrooney P, of data digitization on data integrity. In: Qual-
Poisot T, Woo KH, Zimmerman NB, Hollister ity assurance implementation in research labs,
JW (2016) Ten simple rules for digital data pp 221–245
storage, vol vol. 12. Public Library of Science,
San Francisco, p e1005097 43. Fleming IO, Garratt C, Guha R, Desai J,
Chaubey S, Wang Y, Leonard S, Kunst G
34. Gika HG, Theodoridis GA, Wilson ID (2008) (2016) Aggregation of marginal gains in car-
Liquid chromatography and ultra-performance diac surgery: feasibility of a perioperative care
liquid chromatography–mass spectrometry fin- bundle for enhanced recovery in cardiac surgi-
gerprinting of human urine: sample stability cal patients. J Cardiothorac Vasc Anesth 30:
under different handling and storage condi- 665–670
tions for metabonomics studies. J Chromatogr
A 1189:314–322 44. Tan WS, Tan MY, Lamb BW, Sridhar A,
Mohammed A, Baker H, Nathan S, Briggs T,
35. Gika HG, Theodoridis GA, Wingate JE, Wil- Tan M, Kelly JD (2018) Intracorporeal robot-
son ID (2007) Within-day reproducibility of an assisted radical cystectomy, together with an
HPLC- MS-based method for metabonomic enhanced recovery programme, improves post-
analysis: application to human urine. J Prote- operative outcomes by aggregating marginal
ome Res 6:3291–3303 gains. BJU Int 121:632–639
36. Yin P, Lehmann R, Xu G (2015) Effects of 45. Gilbert D, Henry A, Choudhury A (2022)
pre-analytical processes on blood samples used Marginal gains and clinical trials–improving
in metabolomics studies. Anal Bioanal Chem and influencing practice. Clin Oncol 34:419
407:4879–4892
46. Garcı́a-Rivera MA, Fernández-Ochoa Á,
37. Gil A, Siegel D, Permentier H, Reijngoud DJ, Brüning U, Fritsche-Guenther R, Kirwan JA
Dekker F, Bischoff R (2015) Stability of energy (2022) Identification and validation of small
metabolites—an often overlooked issue in molecule analytes in mouse plasma by liquid
metabolomics studies: a review. Electrophore- chromatography–tandem mass spectrometry:
sis 36:2156–2169 a case study of misidentification of a short-
38. Dunn WB, Broadhurst DI, Edison A, chain fatty acid with a ketone body. Talanta
Guillou C, Viant MR, Bearden DW, Beger 242:123298
RD (2017) Quality assurance and quality con-
trol processes: summary of a metabolomics
Chapter 3
Abstract
Metabolic profiling performed using untargeted metabolomics of different, complex biological samples
aims to apply agnostic/holistic, hypothesis-free, analysis of the small molecules that are present in the
analyzed sample. This approach has been the center of major investments and dedicated efforts from the
research community for many years. However, limitations and challenges remain, particularly with regard to
the validation and the quality of the obtained results. This has led to increasing community engagement,
with the formation of think tanks, the establishment of working groups, and the many seminars on quality
control (QC) in metabolomics. Here we describe a quality control (QC) protocol used to monitor LC-MS-
based metabolomics analysis. A key target is the monitoring of analytical precision. This methodology is
described for the analysis of urine but can be applied to different biological matrices, such as various
biofluids, cell, and tissue extracts.
Key words Quality control, Untargeted metabolomics, Biological samples, Validation, LC-MS
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_3,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
53
54 Olga Begou et al.
2 Materials
2.1 Stock, Working, Stock solutions of all metabolites should be prepared at a concen-
and Calibration tration of 1 mg/mL (or appropriate concentration depending on
Solutions analyte solubility) in methanol or mixtures of methanol with water.
Working standards at concentration of 10 μg/mL are prepared
from the stock solution by appropriate dilution. From the stock
solutions of the standards, prepare a mixture (concentration
1–5 μg/mL) of the compounds of interest by mixing appropriate
volumes.
All solutions should be stored at 20 C.
2.2 Mobile Phase 1. Mobile Phase A: Water +0.1 vol. % formic acid: add 100 μL of
formic acid to 1 L of Millipore water.
2. Mobile Phase B: Acetonitrile +0.1 vol. % formic acid: add
100 μL of formic acid to 1 L of LC-MS grade acetonitrile. In
the case of plasma/serum samples, use methanol in mobile
phase B instead of acetonitrile, or a mixture of methanol with
acetonitrile.
3. Wash solvent: Use a “strong” solvent containing acetonitrile/
water, 80:20 v/v for post injection cleaning cycles and a weak
solvent of water/acetonitrile, 80:20 v/v for pre-injection
washes (see Note 2).
Quality Control and Validation Issues in LC-MS-Based Metabolomics 57
3 Methods
3.1 Analytical 1. In the beginning, analysts must make sure that the chro-
System Preparation matographic system and the mass spectrometer are in suitable
condition for the analysis. Necessary steps: Check the mass
accuracy of the mass spectrometer; if necessary, recalibrate the
mass spectrometer following the appropriate procedure recom-
mended by the manufacturer to achieve maximum mass accu-
racy and resolution. It is advisable that the calibration
procedure should take place every 3–4 months and that the
temperature of the laboratory should not vary significantly.
Detection sensitivity should also be checked comparing analyt-
ical signals with signals obtained from the analysis of a known
sample (system suitability test).
2. Load the required liquid chromatographic method. This would
include adjusting parameters, such as flow rate, column and
autosampler’s temperature, wash cycles, and gradient elution
program.
3. Allow the system to equilibrate and run a no-injection gradient
in order to aid column equilibration.
4. Depending on the matrix of the samples being analyzed, run a
suitable number of pooled QC/bulk matrix samples at the
beginning of the analytical batch to achieve system stability. It
has been observed that for urine analysis 5 QC injections are
needed; for serum or plasma extract this may need to be
increased to 10–20 injections (see Note 3 and 4).
58 Olga Begou et al.
5. The initial 5–10 injections of the pooled QC that are run prior
to the commencement of the analysis of the batch in order to
equilibrate the system also represent a useful “system suitability
test” (SST). The data from these samples should be examined
to confirm that peak shapes are acceptable and that the
expected number of peaks for a sample of that type of matrix
are being seen over the expected retention time range. Evi-
dence of aberrant behavior for the conditioning QCs should
result in the analytical batch not being run until the problem
has been corrected.
3.2 Sample Handling 1. It is advisable to aliquot all test samples on collection into an
and QC Preparation appropriate number of sub-aliquots for later storage and in
order to avoid unnecessary freeze/thaw cycles.
2. Ideally, samples should be stored in freezers as soon as possible
after sampling/collection, and at the lowest available tempera-
ture, at least at 20 C but preferably at 80 C. If possible,
store sub-aliquots in different freezers, in case a freezer mal-
function takes place.
3. Make sure all samples, stock, and working solutions are cor-
rectly and fully labeled. Thaw stock/working solutions and real
samples shortly before use and sample preparation, respectively.
Use appropriate labels that are not affected by freeze-thaw
cycles.
3.3 Sample 1. Vigorously vortex every sample, after thawing, prior to sample
Preparation preparation.
2. Transfer an appropriate volume (e.g., 50–200 μL) of each
sample into a 1.5 mL Eppendorf tube.
3. Depending on the matrix (urine, plasma, serum), a different
sample preparation procedure is recommended (see Note 5 and
6), as explained below in Subheadings 3.3.1, 3.3.2, and 3.3.3.
3.3.1 For Urine Samples 1. For human urine, dilute the sample in ratio 1:3 v/v with
ultrapure water. For animal (e.g., rodent), urine dilution with
methanol rather than water may be required to precipitate
protein.
2. Vortex for 5 min and centrifuge at 12,000 g or higher speeds
for 5–10 min at 4 C.
3. Transfer the clear supernatant into an LC-MS vial (n.b. 96 well
plates or similar can be substituted for LC-MS vials).
4. Place the vial (well plate) into a precooled (4 C) autosampler
for analysis. If analysis is not being undertaken, immediately
store the vial (well plate) in a fridge (0–4 C) if analysis will start
shortly or in a freezer ( 20/80 C) if analysis will be delayed
for longer than a few hours.
Quality Control and Validation Issues in LC-MS-Based Metabolomics 59
3.3.2 For Plasma or 1. Add three times the sample volume of ice-cold methanol or
Serum Samples acetonitrile for protein precipitation. (see Note 7).
2. Vortex for 5 min and centrifuge at 12,000 g or higher speeds
for 5–10 min.
3. Transfer the clear supernatant into Eppendorf tubes.
4. Evaporate to dryness and reconstitute with water/acetonitrile,
95:5 v/v to the initial volume.
5. Repeat step 2.
6. Transfer the clear supernatant into an LC-MS vial (well plate).
7. Place the vial (well plate) into a precooled (4 C) autosampler
for analysis. If analysis is not taking place immediately, follow
the procedure described for urine in Subheading 3.3.1, Step 4.
When taking plasma/serum samples out of the freezer, it is
advisable to centrifuge samples again before analysis.
3.3.3 For Quality Control 1. Create a pool sample by mixing equal volumes of each sample
Samples (QC) analyzed. Make an appropriate number of aliquots (a minimum
of 1 QC for each 10 samples to be analyzed and possibly up to
1 for every 5, see below) and store them at a freezer. If there are
clear groups (e.g., test/patient and control), it may also be
worth preparing “phenotypic” QC s made up of samples
from these groups. Such phenotypic QCs can then be analyzed
alongside the pooled QC to highlight any differences in analyt-
ical properties due to differences in sample composition.
2. Handle each QC as a real sample and follow the sample prepa-
ration procedure developed for the type of matrix to be
analyzed.
3. If the number of the samples being analyzed is large (e.g.,
>1000) and the creation of a pooled sample from all the sam-
ples is impractical, consider making QC samples from only the
samples contained within each analytical run to provide a
within-batch QC. Have a separate “bulk matrix” QC sample
also analyzed within every batch to enable between batch com-
parisons and normalization [16–18].
4. If the acquired samples are only available in small amounts,
such that making a pooled QC of any sort is impractical, use a
bulk sample of the same matrix (ideally from the same species).
3.3.4 For Test Mix 1. Spike a QC sample with an appropriate volume of the metha-
Sample nolic mixture containing appropriate standard metabolites as
mentioned in Subheading 2.1, in order to achieve a final con-
centration of 5–10 μg/mL. It is advisable to perform this
procedure by evaporating the methanolic stock solution and
then reconstituting by dissolution in the QC matrix.
60 Olga Begou et al.
3.4 Sequence 1. Make sure that the order of the samples is randomized, to avoid
Preparation introducing bias due to changes with individual analytical runs
and between run “batch” effects. Randomization of large
sample-sets can be performed using the specific commands in
spreadsheet programs, such as Microsoft excel. QC samples
should not be randomized but inserted regularly in the run
sequence (e.g., one QC sample every 5–10 real samples).
2. The number of QC samples placed in a sequence depends on
the total number of samples analyzed and on the duration of
the analysis. For a small number of samples (<100), one QC
every five sample injections is recommended. For larger batches
of samples (100 to several thousand), the number QC samples
should represent at least 10% of the total real samples analyzed.
3. After column equilibration, in the middle and at the end of
the sequence, insert in the analytical sequence injections of
the spiked test mix solution. Injections of standards solutions
or blanks should be avoided within the batch of test samples,
as the system can be disequilibrated by the injection of
non-matrix solutions (see Notes 3 and 8). Similarly, blank
sample injections should be avoided in the middle of the ana-
lytical batch.
4. A good way to improve the relative quantification aspects in
such analysis is to implement a sequence of injection of serial
dilutions of the QC sample. This enables the detection of
saturating peaks, which would be better detected in dilute
samples (e.g., 1:10 v/v). At the same time, this series of ana-
lyses facilitates the study of the response of non-saturating
peaks in response to dilution. The dilution series can be run
at the end of each batch of samples.
3.5.1 Data Processing From the raw data acquired, observe the peak width, retention
time, peak intensity/signal, mass accuracy, and noise intensity
from indicative chromatographic peaks located at the beginning,
middle, and end of the chromatogram. These can help researchers
to define carefully the optimum parameters for the software
(XCMS, MassLynx etc.) to be used for peak alignment, peak peak-
ing, and integration.
Data from the QC samples are then used to evaluate the quality
of analysis via the following steps: Treat QC samples as a separate
group and process them alone using exactly the same parameters
Quality Control and Validation Issues in LC-MS-Based Metabolomics 61
selected for processing of the whole sample set. Peaks that appear in
less than 70 or 80% of the QC samples should be omitted from the
dataset. Apply Principal Component Analysis (PCA), using
SIMCA-P or other statistical software, to the processed data from
the QC samples in order to observe any trends, such as time-related
drift during the analytical run of these specific samples.
More details on data preprocessing can be found in Chaps. 4
and 5 of this book.
3.5.2 Data Quality For data quality evaluation, there are some steps that should be
Evaluation taken into consideration. It is advisable to subject the data to the
test against preset criteria of quality, as explained below.
Firstly, assess the performance of the analytical run by inspect-
ing an overlay of full scan chromatograms of all the replicates of the
test mix samples, placed at the beginning, middle, and end of the
sequence. Calculate chromatographic data, such as the retention
time (tR), peak area, and height and determine the precision (RSD
%) of each value. The acceptable limits should be less than 2% for
retention time variation and less than 20–30% for peak area varia-
tion. In cases where these conditions are not met, reexamine the
data carefully to find any trend that might indicate system under-
performance or malfunction. MS-derived data should also be
acceptable in terms of mass accuracy.
The next step is to assess the data acquired from QC samples
placed at regular intervals among the samples, in order to
strengthen and corroborate the initial assessment of the perfor-
mance of the analytical run. Evaluate the full scan chromatograms
of all QC samples. It is expected that all QC samples show good
consistency, reproducibility, and small variability with respect to tR,
peak area, and height, with the possible exception of the first
5 to 10 injections (conditioning QCs). An overlay of full scan
chromatograms can provide only an estimation of analytical repeat-
ability and system stability. A thorough check should be done on
extracted ion chromatograms of all peaks along the whole length of
the analysis. Precision of the tR, the mass detected, and the signal
intensity should be unchanged across the chromatogram and
should not show, e.g., time-related effects. These can be assessed
against the preset criteria of, e.g., RSD% should be less than 2% for
tR, less than 25% for the area of abundant peaks, and less than 30%
of low signal peaks), as mentioned above. Peaks that fail such
criteria should be excluded from the dataset (including data from
the test samples). Evidence of poor repeatability for the QC sam-
ples should be taken as a possible reason for not accepting a run. A
recommended scheme of QC in the analytical data acquired by
untargeted profiling is shown in Fig. 1.
If there are any problems concerning the reproducibility of the
data, try to estimate where the inconsistency is appearing. If it is
62 Olga Begou et al.
Fig. 1 Suggested flow chart of QC in untargeted mass spectrometry-based metabolic profiling based on QC
pooled samples analysis [17]
2. Check for any QC outliers or for batch effect due to run order.
3. Examine the QC data for any unexpected trends. If such trends
are found, study additional visualization plots, such as time
series plots, and see for samples that are out of 2 SD or 3 SD
thresholds. Test samples appearing as such “outliers” should
not be discarded outright but rather be thoroughly scrutinized.
QC samples out of 2 SD indicate that the data from the
neighboring test samples should be thoroughly scrutinized
and probably discarded with the samples being reanalyzed.
Thereinafter, in the exported peak table data, look for any
similarities in the pairs (tR and mass) among the peaks reported,
in order to assess the success of peak alignment. Then, using
appropriate statistical software, create graphs with the data
exported from the QC samples to evaluate the analytical perfor-
mance and precision.
Finally, when preparing publications/reports where QCs have
been used, ensure that all relevant information is provided (guide-
lines and proformas for this are given [17]).
Figure 2 presents a schematic of the proposed validation
scheme. It starts with the QC and sequence preparation, by initially
mixing equal volumes of all the samples to be analyzed (formed
form either urine or plasma/serum) in a vial and then injecting the
QC sample at the beginning of the analytical batch (system equili-
bration) and then periodically, e.g., every 5–10 real samples during
analysis. From the resulting dataset, an extracted-ion chromato-
gram (XIC) of all QC samples can be created, via data processing or
data-mining, where data are collected only for specific m/z values
of the compounds of interest. Those results can also be used for the
formation of an analysis chart, showing total sequence order (both
QC and real samples). Data processing is complemented by statis-
tical analysis. Repeatability and run order trends can be assessed by
evaluating %RSD values of QC’s ion intensities. Features presented
in most samples (~70%) with %RSD values <30% represent a good
dataset on which to carry out further work. Finally, PCA score plots
can be created from the examined results, providing information
for the accuracy and precision of the analysis, as well as, for any time
trends according order analysis of both QC and real samples.
4 Notes
Fig. 2 Schematic representation of the acquisition and acceptance procedures followed during and after
analytical batch is analyzed. The process begins with quality control (QC) samples and sequence preparation
(equal volumes of all real samples analyzed are mixed to provide the pooled QC). The QC sample is then
injected 5–10 times at the beginning of the analytical batch, for system equilibration, and typically every 5–10
real samples, throughout the analytical batch(s). From the resultant dataset, an extracted-ion chromatogram
(XIC) can be generated for all QC samples for selected m/z values, where the peak width, retention time, peak
intensity/signal, and mass accuracy can be observed. The final step is the statistical evaluation of the results.
To assess repeatability and trends in the run order, %RSD values of QC’s peak areas are evaluated. A robust
dataset is indicated by features present in most samples (> 70%) with %RSD values lower than 30%.
Additionally, a PCA score plot can be generated, in order to provide insights regarding the robustness of the
analysis. QCs should cluster together
Acknowledgments
References
1. Gika HG, Wilson ID, Theodoridis GA (2019) spectrometry-based metabolic phenotyping
Omics | metabolomics: an analytical (metabolomics/metabonomics) studies in
perspective. In: Encyclopedia of analytical sci- publications. J Chromatogr B 1164:122515
ence. https://doi.org/10.1016/B978-0-12- 7. Gika HG, Theodoridis GA, Earll M et al
409547-2.14003-X (2010) Does the mass spectrometer define the
2. Theodoridis GA, Gika HG, Plumb R, Wilson marker? A comparison of global metabolite
ID (2019) Liquid chromatographic methods profiling data generated simultaneously via
combined with mass spectrometry in UPLC-MS on two different mass spectro-
metabolomics. In: Proteomic and metabolo- meters. Anal Chem 82:8226–8234. https://
mic approaches to biomarker discovery, pp doi.org/10.1021/ac1016612
149–169. https://doi.org/10.1016/B978-0- 8. Rainville PD, Theodoridis G, Plumb RS, Wil-
12-818607-7.00009-8 son ID (2014) Advances in liquid chromatog-
3. Wilson I (2016) Methods and techniques for raphy coupled to mass spectrometry for
metabolic phenotyping. Bioanalysis 9:1–3. metabolic phenotyping. TrAC Trends Anal
https://doi.org/10.4155/bio-2016-4985 Chem 61:181–191. https://doi.org/10.
4. Gika HG, Wilson ID, Theodoridis GA (2014) 1016/j.trac.2014.06.005
LC–MS-based holistic metabolic profiling. 9. Begou O, Gika HG, Wilson ID et al (2017)
Problems, limitations, advantages, and future HyphenatedMS–based targeted approaches in
perspectives. J Chromatogr B 966:1–6. metabolomics. Analyst
https://doi.org/10.1016/j.jchromb.2014. 10. Gika HG, Theodoridis GA, Earll M, Wilson ID
01.054 (2012) A QC approach to the determination of
5. Gika H, Virgiliou C, Theodoridis G, Plumb day-to-day reproducibility and robustness of
RS, Wilson ID (2019) Untargeted LC/MS- LC-MS methods for global metabolite
based metabolic phenotyping metabonomics/ profiling in metabonomics/metabolomics.
metabolomics: the state of the art. J Chroma- Bioanalysis 4:2239–2247. https://doi.org/
togr B 1117:136–147 10.4155/bio.12.212
6. Wilson ID, Theodoridis G, Virgiliou C (2021) 11. Gika HG, Zisi C, Theodoridis G et al (2016)
A perspective on the standards describing mass Protocol for quality control in metabolic
66 Olga Begou et al.
Abstract
Metabolomics data analysis includes, next to the preprocessing, several additional repetitive tasks that can
however be heavily dataset dependent or experiment setup specific due to the vast heterogeneity in
instrumentation, protocols, or also compounds/samples that are being measured. To address this, various
toolboxes and software packages in Python or R have been and are being developed providing researchers
and analysts with bioinformatic/chemoinformatic tools to create their own workflows tailored toward their
specific needs. This chapter presents tools and example workflows for common tasks focusing on the
functionality provided by R packages developed as part of the RforMassSpectrometry initiative. These
tasks include, among others, examples to work with chemical formulae, handle and process mass spectrom-
etry data, or calculate similarities between fragment spectra.
Key words R, RforMassSpectrometry, Formula handling, Mass spectra handling, Spectra similarity
calculation
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_4,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
67
68 Michael Witting and Johannes Rainer
2 Materials
2.1 R Scripts All scripts can be used directly in R or more comfortable with an
IDE-like RStudio, which allows more convenient work without the
need to remember all variables, since they are displayed in a separate
window which also allows manual inspection of their content. The
presented code snippets work with different packages, which are
not part of the base R distribution and have to be therefore installed
manually. A version of R > = 4.3 and Bioconductor > = 3.18 is
required for the R scripts. The required packages, along with their
dependencies, can be installed directly within R issuing the com-
mands below:
install.packages("BiocManager")
library(BiocManager)
BiocManager::install(c("Spectra", "MetaboAnnotation",
Metabolomics Data Analysis 69
"AnnotationHub", "MsBackendMgf",
"msdata", "MetaboCoreUtils",
"pheatmap", "msentropy",
"MobilityTransformR",
"BridgeDbR"))
3 Methods
chem_formula_pasted <-
pasteElements(chem_formula_count)
3.1.2 Calculation of 1. Exact masses can be calculated from chemical formulae without
Exact Masses, Kendrick or with isotopes. Either strings or named lists are accepted as
Masses, or Kendrick Mass input in the function calculateMass().
Defect
calculateMass(chem_formula)
calculateMass(chem_formula_pasted)
calculateMass(iso_formula)
calculateMass(iso_formula_pasted)
3. From the provided formula, the mass and the m/z are calcu-
lated using the function mass2mz() (see below). This returns a
named vector with the m/z values and the input formula as
names.
4. From the m/z value, the Kendrick mass (KM) and Kendrick
mass defect (KMD) are calculated using the functions calcu-
lateKm() and calculateKmd(). By default, the KM and
KMD are calculated for CH2, but also other KMD are possible
(see Note 2).
isRkmd(lipid_rkmd,
rkmdTolerance = 0.02)
3.1.3 Working with 1. From obtained exact masses, the exact m/z can be calculated
Adducts for different adducts. The function mass2mz() performs the
calculation and accepts a vector of exact masses and adduct
definitions as input. Results are returned as a matrix with one
row per input mass and one column per defined adduct.
adductNames()
adductNames(polarity = "positive")
adductNames(polarity = "negative")
mz2mass(c(181.0707, 203.0526,
132.1019, 154.0838),
c("[M+H]+", "[M+Na]+"))
3.1.4 Working with 1. Isotopic pattern can be used for calculation of potential sum
Isotopic Patterns formulae. The package Rdisop allows the decomposition of
masses, m/z, and isotope patterns. Similar to MetaboCoreUtils,
Rdisop also provides a function to calculate masses from chemi-
cal formulae. Additionally, if these formulae represent ions, the
charge can be defined.
getMass(getMolecule("C6H12O6")
getMass(getMolecule("C6H13O6", z = 1)
getMass(getMolecule("C6H12O6")) -
calculateMass("C6H12O6")
getIsotope(getMolecule("C6H13O6", z = 1))
getIsotope(getMolecule("C6H12O6Na", z = 1))
decomposeIsotopes(iso_pattern$masses,
iso_pattern$intensities,
z = 1,
ppm = 5.0)
matchFormula(queries, targets)
matchedData(matchFormula(queries, targets))
3.2 Retention Time Metabolite identification on the highest level requires, according to
Indexing and Mobility the Metabolomics Standard Initiative, information orthogonal to
Transformation MS1, MS2, or MSn. Retention or migration times represent such an
orthogonal information, but separation conditions in LC and CE
are far from being standardized. Even when using nominally the
74 Michael Witting and Johannes Rainer
3.2.1 Retention Time 1. Retention time indexing requires a number of index markers,
Indexing e.g., N-alkyl pyridinium-3-sulfonates (NAPS), as described by
Stoffel et al. In the example below we create a data frame with
retention times and retention indices of these markers. Typi-
cally, the retention index is the number of carbons in the side
chain of NAPS multiplied by 100.
"tauroursodeoxycholic acid"),
rtime = c(0.46, 2.759,
5.222))
3.2.2 Effective Mobility 1. MetaboCoreUtils also offers the possibility to perform effective
Transformation mobility transformation for CE-MS-based metabolomics
experiments. Markers of the electroosmotic flow are required
to perform this conversion.
3.3 Handling and Handling of different kinds of MS spectra, e.g., MS1 isotope pat-
Processing Mass tern or MS2 fragmentation pattern, is an important task for meta-
Spectra Data bolomics data analysis and metabolite identification workflows. The
Spectra package provides a flexible and expandable infrastructure
for loading, handling, and analyzing MS data in R. By default, the
package allows import/export of data from/to files in mzML,
76 Michael Witting and Johannes Rainer
library(Spectra)
library(MsBackendMgf)
map <- spectraVariableMapping(MsBackendMgf())
map <- c(map, c(msLevel = “MSLEVEL”, smiles = “SMILES”))
spectraVariables(gnps)
msLevel(gnps)
gnps$msLevel
5. Peak data, i.e., the m/z and intensity values of the individual
mass peaks per spectrum, can be extracted separately using
either the mz() and intensity() functions or together in a
two-dimensional array using the peaksData() function. For
the former case, a list of numeric values is returned and for the
latter a list of numerical (two-column) matrices.
3.3.2 Processing and 1. Fragment spectra (and also MS1 spectra) can contain a large
Cleaning of Spectra Data amount of low-intensity peaks that might represent in most
cases noise. The filterIntensity() function can be used
to remove such low-signal peaks. In addition to a fixed intensity
cutoff, also a function defining a spectrum-specific cutoff can
be provided. This function is supposed to take intensity values
of one spectrum as an input (parameter x) and to return a
logical vector of length equal to the number of peaks in a
spectrum defining whether a mass peak should be kept or
removed. The function below returns TRUE for all peaks with
an intensity larger than 1% of the maximum intensity of that
spectrum.
median(lengths(gnps))
gnps <- filterIntensity(gnps, intensity = int_cut)
median(lengths(gnps))
78 Michael Witting and Johannes Rainer
3.3.3 Exporting MS Data 1. Data from any Spectra object can be exported to a file in
in MGF Format MGF format using the export() function and defining the
output format by setting backend = MsBackendMgf(). We in
addition pass the mapping between spectra variables and MGF
fields defined above with parameter mapping. Finally, with
optional parameter exportTitle = FALSE, we disable the
default export of a title (name) for each spectrum. The gener-
ated MGF file will then be in GNPS-flavor format. Note that it
would be similarly possible to export data in, e.g., MSP format
by loading the MsBackendMsp package and using back-
end = MsBackendMsp() instead.
3.3.4 Spectra Similarity Calculating similarities between fragment spectra is a core concept
Calculations for most annotation/identification approaches. While a large num-
ber of software and similarity measures exist, we focus here on
functionality provided by the MsCoreUtils, Spectra, and MetaboAn-
notation Bioconductor packages [3]. Spectra similarity calculation
Metabolomics Data Analysis 79
fl <- system.file("TripleTOF-SWATH",
"PestMix1_DDA.mzML",
package = "msdata")
dda <- Spectra(fl, source = MsBackendMzR())
2. We next restrict the data to MS2 spectra and apply the intensity
filter defined in the previous example that removes all peaks
with an intensity below 1% of the maximum peak signal in a
spectrum.
Fig. 1 Heatmap visualizing the pairwise similarity matrix between all fragment spectra of the example dataset.
Spectra are ordered by retention time
library(pheatmap)
pheatmap(sim, cluster_rows = FALSE, cluster_cols = FALSE)
Fig. 2 Mirror plot of two fragment spectra with high similarity. Matching peaks between the two spectra are
indicated with a blue color
10. While a mirror plot is ideal for pairwise comparison, it does not
allow to visualize multiple spectra in a single plot. One possi-
bility for such a visualization would be a simple intensity heat-
map. To create such a plot for the first ten spectra, in our
dataset, we generate next the required two-dimensional repre-
sentation of our data by first binning the spectra along the m/z
dimension with a bin size of 0.1 and then extract and combine
the intensity values for these binned spectra into an intensity
matrix.
Fig. 3 Visualization of multiple spectra. Rows in this plot represent different spectra and columns (binned) of
m/z values with color-coded intensities of fragment peaks in the individual spectra
z <- colSums(imat) == 0
imat <- imat[, !z]
colnames(imat) <- mz(dda_bin)[[1]][!z]
rownames(imat) <- seq_along(dda_bin)
pheatmap(imat, cluster_rows = FALSE, cluster_cols = FALSE)
library(msentropy)
compareSpectra(dda[2], dda[3], MAPFUN = joinPeaksNone,
FUN = msentropy_similarity,
ms2_tolerance_in_ppm = 50,
ms2_tolerance_in_da = 0,
clean_spectra = FALSE)
library(AnnotationHub)
ah <- AnnotationHub()
query(ah, "MassBank")
mb <- ah[["AH107049"]]
Fig. 4 Mirror plot between an experimental fragment spectrum (top) and a reference spectrum from MassBank
(bottom). Matching peaks are indicated by a blue color
Metabolomics Data Analysis 87
3.4 Mapping/ Depending on the further downstream data analysis tool, e.g.,
Conversion of MetaboAnalyst or others, different identifiers are required for path-
Metabolite Identifiers way mapping or similar tasks. Manual search for different identifiers
in database is a tedious task since no harmonized metabolite data-
base exists. Metabolite structures represent their most unique iden-
tifier, but often identifiers from different databases such as KEGG,
HMDB, and ChEBI are preferred. Such identifiers can be con-
verted into each other, for example, using the chemical translation
service. Another possibility that can be easily integrated into R
workflows is BridgeDB via the BridgeDbR package [11, 12].
1. The BridgeDbR package requires the most current bridge
database, which can be found on the project website (see
above). This database is loaded into a variable called mapper
used for the entire conversion process.
getSystemCode("InChIKey")
getSystemCode("ChEBI")
getSystemCode("KEGG Compound")
map(mapper,
"C00002",
source = getSystemCode("KEGG Compound"),
target = getSystemCode("ChEBI"))
if(grepl("^C\crd+$", x)) {
return(getSystemCode("KEGG Compound"))
} else if (grepl("^HMDB\crd+$", x)) {
return(getSystemCode("HMDB"))
} else {
return(NA_character_)
}
4 Notes
References
1. Smith CA, Want EJ, O’Maille G, Abagyan R, 7. Hao J-D, Chen Y-Y, Wang Y-Z, An N, Bai P-R,
Siuzdak G (2006) XCMS: processing mass Zhu Q-F et al (2023) Novel peak shift correc-
spectrometry data for metabolite profiling tion method based on the retention index for
using nonlinear peak alignment, matching, peak alignment in untargeted metabolomics.
and identification. Anal Chem 78(3): Anal Chem 95(35):13330–13337. https://
7 7 9 – 7 8 7 . h t t p s : // d o i . o r g / 1 0 . 1 0 2 1 / doi.org/10.1021/acs.analchem.3c02583
ac051437y 8. Drouin N, Pezzatti J, Gagnebin Y, González-
2. Benton HP, Wong DM, Trauger SA, Siuzdak G Ruiz V, Schappler J, Rudaz S (2018) Effective
(2008) XCMS2: processing tandem mass spec- mobility as a robust criterion for compound
trometry data for metabolite identification and annotation and identification in metabolomics:
structural characterization. Anal Chem 80(16): toward a mobility-based library. Anal Chim
6382–6389 Acta 1032:178–187. https://doi.org/10.
3. Rainer J, Vicini A, Salzer L, Stanstrup J, Badia 1016/j.aca.2018.05.063
JM, Neumann S et al (2022) A modular and 9. González-Ruiz V, Gagnebin Y, Drouin N,
expandable ecosystem for metabolomics data Codesido S, Rudaz S, Schappler J (2018)
annotation in R. Metabolites 12(2):173 ROMANCE: a new software tool to improve
4. Korf A, Vosse C, Schmid R, Helmer PO, data robustness and feature identification in
Jeck V, Hayen H (2018) Three-dimensional CE-MS metabolomics. Electrophoresis
Kendrick mass plots as a tool for graphical 39(9-10):1222–1232. https://doi.org/10.
lipid identification. Rapid Commun Mass 1002/elps.201700427
Spectrom 32(12):981–991. https://doi.org/ 10. Li Y, Kind T, Folz J, Vaniya A, Mehta SS, Fiehn
10.1002/rcm.8117 O (2021) Spectral entropy outperforms
5. Stoffel R, Quilliam MA, Hardt N, Fridstrom A, MS/MS dot product similarity for small-
Witting M (2022) N-Alkylpyridinium sulfo- molecule compound identification. Nat Meth-
nates for retention time indexing in reversed- ods 18(12):1524–1531. https://doi.org/10.
phase-liquid chromatography-mass spectrome- 1038/s41592-021-01331-z
try-based metabolomics. Anal Bioanal Chem 11. van Rijswijk M, Beirnaert C, Caron C,
414(25):7387–7398. https://doi.org/10. Cascante M, Dominguez V, Dunn W et al
1007/s00216-021-03828-0 (2017) The future of metabolomics in ELIXIR
6. Aalizadeh R, Nikolopoulou V, Thomaidis NS [version 2; peer review: 3 approved].
(2022) Development of liquid chro- F1000Research 6:1649. https://doi.org/10.
matographic retention index based on coca- 12688/f1000research.12342.2
mide diethanolamine homologous series (C 12. Batchelor C, Brenninkmeijer CYA,
(n)-DEA). Anal Chem 94(46):15987–15996. Chichester C, Davies M, Digles D, Dunlop I
https://doi.org/10.1021/acs.analchem. et al (2014) Scientific lenses to support multi-
2c02893 ple views over linked chemistry data. Springer,
Cham, pp 98–113
Chapter 5
Abstract
Liquid Chromatography-Mass Spectrometry (LC-MS) untargeted experiments require complex bioinfor-
matic strategies to extract information from the experimental data. Here we discuss the “data preproces-
sing,” the set of procedures performed on the raw data to produce a data matrix which will be the starting
point for the subsequent statistical analysis. Data preprocessing is a crucial step on the path to knowledge
extraction, which should be carefully controlled and optimized in order to maximize the output of any
untargeted metabolomics investigation.
Key words Preprocessing, Peak picking, Retention time correction, Metadata, Quality check, Missing
values
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_5,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
91
92 Mar Garcia-Aloy et al.
Fig. 1 Data analysis workflow of a typical LC-MS untargeted experiment. Only the colored boxes are detailed
in the text
2 Materials
3 Methods
3.1 Metadata The term “experiment metadata” refers to the set of information
Organization describing the samples and the analytical pipeline (see Note 2).
Organizing and storing a detailed record of the sample metadata
(see Note 3) is crucial to perform a reliable interpretation of the
results, evaluating the presence (and the impact) of possible con-
founding factors.
1. Collect and organize (possibly with the help of a spreadsheet)
all the important properties of the samples (e.g., treatment
type, day of collection, origin, etc.). Keep track also of the
sample-specific analytical details such as batch, column, mea-
surement date, and time.
2. It is recommended to define a spreadsheet that contains the file
names of the raw data files along with descriptions of the
94 Mar Garcia-Aloy et al.
3.2 Data Conversion In general, the vendor software of the MS instrument manufacturer
(See Note 4) can be used to perform this task (see Note 5). Alternatively, or if the
vendor software does not support data export in the open mzML
file format, the open source ProteoWizard “msconvert” software
(see Note 6), which is available as a command line tool or as a point-
and-click application for Windows (MSconvertGUI) (see Note 7),
can be used. Also, a msconvert docker image is available that allows
reproducible batch processing and conversion of large sets of files
on, e.g., a linux cluster.
1. Open a terminal or a command window and run “msconvert”
into the raw data directory. As an example, to convert vendor-
specific raw data (here .RAW data files) into mzML files and
store them into the open_data directory, use:
3.3 Preprocessing Each metabolite produces at least one peak in the (m/z, rt) plane.
These are sometimes also referred to as chromatographic peaks to
distinguish them from the mass peaks, which are individual signals
in a single mass spectrum measured by the MS instrument. To
quantify metabolites, a dedicated peak detection algorithm has to
be applied to the raw MS data that is able to distinguish real signal
of an ion from electrical or chemical noise [7, 20]. This process of
chromatographic peak detection (peak picking) has to be per-
formed (separately) on all data files and needs to be customized
depending on the LC-MS setup used for the experiment. Subse-
quently, peaks representing signals from the same type of ions (with
similar m/z and rt) need to be grouped across samples to define the
final list of features (correspondence analysis or grouping). In order
to do that, it is necessary to account for any potential shifts of m/z
and rt between samples (retention time correction) [21]. The
measured m/z values can differ between samples (or even between
scans of the same sample) because even the most advanced mass
spectrometer instrument is not absolutely precise. Rt shifts, instead,
Data Treatment for LC-MS Untargeted Analysis 95
3.4 Preprocessing: First it is necessary to load the MS data, ideally along with the
Data Import corresponding associated metadata. Below we create a data frame
with sample annotation based on the file names of the experiment
and load the MS data into R using the readMSData() function
from the MSnbase R package. With the parameter mode =
"onDisk", only general spectra information is imported into
memory, while m/z and intensity values will be retrieved on-the-
fly from the raw data files allowing also processing of very large MS
experiments on standard computer systems [24].
library(xcms)
cdffiles <- list.files(path = "MTBLS59",
pattern = "CDF",
full.names = TRUE)
myfiles <- grep("neg", cdffiles, value = TRUE)[1:40]
pd <- data.frame(sample_name = basename(myfiles),
do.call(rbind, strsplit(gsub(".CDF", "",
basename(myfiles)), "_")))
colnames(pd) <- c("sample_name", "apple", "group",
"polarity", "sample")
raw_data <- readMSData(
files = myfiles,
pdata = new("NAnnotatedDataFrame", pd),
mode = "onDisk")
3.5 Preprocessing: The main characteristics that distinguish a real peak (representing
Peak Picking signal of an ion) from noise are its intensity, width, and shape.
Different algorithms handle them in different ways [7, 20], but
most algorithms for detection of chromatographic peaks in LC-MS
data assume the peak shape to be nearly gaussian. Most importantly,
not only the algorithm but also the specific parameters of each
algorithm greatly affect the results (see Note 8).
96 Mar Garcia-Aloy et al.
plotChromPeaks(test_peaks)
3.6 Preprocessing: 1. Choose the retention time correction algorithm among the
Retention Time available options (see Note 9). In this example, we use obiwarp
Correction and [26], which is based on Dynamic Time Warping, an algorithm
Grouping (See Note 12) used to find the best stretching of the time dimension of two
time series to make them as similar as possible [27]. LC-MS
data are multidimensional so in obiwarp the similarity is deter-
mined based on the full spectral information.
2. Determine the more relevant parameters of the algorithm. In
obiwarp, the relevant parameters are the choice of the reference
sample for the rt warping (center); the m/z bins used for
retrieving the spectra (profStep); the spectral similarity func-
tion (distFunc); and the penalizations on the warping opti-
mization (gapInit and gapExtend). With the default
parameters (which work reasonably well in most of the cases),
the reference sample is the one containing the largest number
of peaks. The retention time correction can be performed with
the following command:
Fig. 4 Outcome of the retention time correction performed on the MTBLS59 data
as obtained by plotAdjustedRtime (), after minor visual parameters
customization. Each line represents one sample and is colored by sample class.
The retention time correction is different for each sample and it is not linear in
time
plotAdjustedRtime(xdata)
bw(pdp) <- 1
plotChromPeakDensity(myeic, param = pdp)
3.7 Preprocessing: It could happen that peaks belonging to a subset of samples are
Imputation of Missing missing from “their” groups. This is correct if the corresponding
Peaks analyte is not present in a specific sample, but it could also be the
result of a failed identification of the chromatographic peak during
peak detection (e.g., if the signal in that particular sample was too
low or too noisy). These missing peaks result in missing values in
the final feature data matrix, which are usually tricky to handle
during statistical analysis. xcms provides fillChromPeaks(), a
102 Mar Garcia-Aloy et al.
Fig. 6 Peak grouping with reduced bw setting on m/z slice and rt window
containing signals from quercetin-3-glucoside and quercetin-3-galactoside in
negative ionization mode
3.8 Data Matrix Extract the final two-dimensional data matrix with samples on the
Extraction rows and features on the columns, using as intensity the integral of
the area under the peak:
dt <- t(featureValues(xdata,
value = "into",
method = "sum"))
Data Treatment for LC-MS Untargeted Analysis 103
The m/z and rt details about the features can be extracted using
the featureDefinitions() function:
ft <- featureDefinitions(xdata)
3.9 Annotation of At this point, features are only defined by a specific m/z and rt
MS-Based value, which does not provide sufficient information to make sense
Metabolomics Data of the generated data. Only by assigning a name (identity) to the
signals of interest (annotation) enables the biological interpretation
of the obtained results. Several software tools for metabolite anno-
tation and identification are available [28]. Here we use the Meta-
boAnnotation R package [18] because it provides a customizable
infrastructure that can be easily applied to xcms output. MetaboAn-
notation provides a set of functions that allow matching between
query and target entities (such as m/z and/or rt values).
1. Create a reference database containing theoretical values. Such
reference values are usually defined by measuring pure stan-
dards on the same LC-MS setup.
library(MetaboAnnotation)
matched_feat <- matchValues(ft, target_df,
param = MzRtParam(tolerance = 0.005,
ppm = 0, toleranceRt = 5),
rtColname = c("rtmed", "rt"),
mzColname = c("mzmed", "mz"))
matchedData(matched_feat)[whichQuery(matched_feat), ]
104 Mar Garcia-Aloy et al.
4 Notes
Fig. 7 Pick picking outcome on the MTBLS59 dataset: (left) the peak detection is performed using the default
settings for centWave, (right) or the set of parameters defined in Methods Subheading 3.5, step 5
References
1. Ardrey RE (2003) Liquid chromatography–mass 13. Kessner D, Chambers M, Burke R et al (2008)
spectrometry: an introduction. https://online ProteoWizard: open source software for rapid
l i b r a r y. w i l e y. c o m / d o i / b o o k / 1 0 . 1 0 0 2 proteomics tools development. Bioinformatics
/0470867299 24:2534–2536. https://doi.org/10.1093/bio
2. Patti GJ, Yanes O, Siuzdak G (2012) Innova- informatics/btn323
tion: metabolomics: the apogee of the omics 14. Chambers MC, MacLean B, Burke R et al
trilogy. Nat Rev Mol Cell Biol 13:263–269. (2012) A cross-platform toolkit for mass spec-
https://doi.org/10.1038/nrm3314 trometry and proteomics. Nat Biotechnol
3. Alonso A, Marsal S, Julià A (2015) Analytical 3010(30):918–920. https://doi.org/10.
methods in untargeted metabolomics: state of 1038/nbt.2377
the art in 2015. Front Bioeng Biotechnol 3:23. 15. Smith CA, Want EJ, O’Maille G et al
https://doi.org/10.3389/fbioe.2015.00023 (2006) XCMS: processing mass spectrometry
4. Sansone S-A, Rocca-Serra P, Field D et al data for metabolite profiling using nonlinear
(2012) Toward interoperable bioscience data. peak alignment, matching, and identification.
Nat Genet 44:121–126. https://doi.org/10. Anal Chem 78:779–787. https://doi.org/10.
1038/ng.1054 1021/ac051437y
5. R Core Team (2016) R: a language and envi- 16. Tautenhahn R, Bottcher C, Neumann S
ronment for statistical computing. https:// (2008) Highly sensitive feature detection for
www.r-project.org/ high resolution LC/MS. BMC Bioinformatics
6. De Vos RCH, Moco S, Lommen A et al (2007) 9. https://doi.org/10.1186/1471-2105-
Untargeted large-scale plant metabolomics 9-504
using liquid chromatography coupled to mass 17. Benton HP, Want EJ, Ebbels TMD (2010)
spectrometry. Nat Protoc 2. https://doi.org/ Correction of mass calibration gaps in liquid
10.1038/nprot.2007.95 chromatography-mass spectrometry metabolo-
7. Rafiei A, Sleno L (2014) Comparison of peak- mics data. Bioinformatics 26. https://doi.org/
picking workflows for untargeted liquid 10.1093/bioinformatics/btq441
chromatography/high-resolution mass spec- 18. Rainer J, Vicini A, Salzer L et al (2022) A
trometry metabolomics data analysis. Rapid modular and expandable ecosystem for meta-
Commun Mass Spectrom 29. https://doi. bolomics data annotation in R. Metabolites 12:
org/10.1002/rcm.7094 1 7 3 . h t t p s : // d o i . o r g / 1 0 . 3 3 9 0 /
8. Yu T, Park Y, Johnson JM, Jones DP (2009) METABO12020173
apLCMS-adaptive processing of high- 19. Franceschi P, Masuero D, Vrhovsek U et al
resolution LC/MS data. Bioinformatics 25. (2012) A benchmark spike-in data set for bio-
https://doi.org/10.1093/bioinformatics/ marker identification in metabolomics. J Che-
btp291 mom 26. https://doi.org/10.1002/cem.
9. Tauler R, Gorrochategui E, Jaumot J, Tauler R 1420
(2015) A protocol for LC-MS metabolomic 20. Zhang J, Gonzalez E, Hestilow T et al (2009)
data processing using chemometric tools. Pro- Review of peak detection algorithms in liquid-
toc Exch. https://doi.org/10.1038/protex. chromatography-mass spectrometry. Curr
2015.102 Genomics 10. https://doi.org/10.2174/
10. Katajamaa M, Orešič M (2007) Data proces- 138920209789177638
sing for mass spectrometry-based metabolo- 21. Åberg KM, Alm E, Torgrip RJO (2009) The
mics. J Chromatogr A 1158:318–328 correspondence problem for metabonomics
11. Lange E, Tautenhahn R, Neumann S, Gröpl C datasets. Anal Bioanal Chem 394:151–162
(2008) Critical assessment of alignment proce- 22. Smith R, Ventura D, Prince JT (2015) LC-MS
dures for LC-MS proteomics and metabolo- alignment in theory and practice: a comprehen-
mics measurements. BMC Bioinformatics 9: sive algorithmic review. Brief Bioinform 16:
375. https://doi.org/10.1186/1471-2105- 104–117. https://doi.org/10.1093/bib/
9-375 bbt080
12. Rocca-Serra P, Salek RM, Arita M et al (2015) 23. Koch S, Bueschl C, Doppler M et al (2016)
Data standards can boost metabolomics MetMatch: a semi-automated software tool for
research, and if there is a will, there is a way. the comparison and alignment of LC-HRMS
Metabolomics 12:14. https://doi.org/10. data from different metabolomics experiments.
1007/s11306-015-0879-3 Metabolites 6. https://doi.org/10.3390/
metabo6040039
108 Mar Garcia-Aloy et al.
24. Gatto L, Gibb S, Rainer J (2021) MSnbase, Mol Cell Proteomics 14. https://doi.org/10.
efficient and elegant R-based processing and 1074/mcp.O114.039115
visualization of raw mass spectrometry data. J 32. Krzywinski M, Altman N (2014) Points of sig-
Proteome Res 20:1063–1069. https://doi. nificance: designing comparative experiments.
org/10.1021/acs.jproteome.0c00313 Nat Methods 11:597–598
25. Brodsky L, Moussaieff A, Shahaf N et al (2010) 33. Krzywinski M, Altman N (2014) Points of sig-
Evaluation of peak picking quality in LC-MS nificance: analysis of variance and blocking. Nat
metabolomics data. Anal Chem 82:9177– Methods 11:699–700
9187. https://doi.org/10.1021/ac101216e 34. Krzywinski M, Altman N, Blainey P (2014)
26. Forsberg EM, Huan T, Rinehart D et al (2018) Points of significance: nested designs. Nat
Data processing, multi-omic pathway Methods 11. https://doi.org/10.1038/
mapping, and metabolite activity analysis nmeth.3137
using XCMS online. Nat Protoc 13:633–651. 35. Krzywinski M, Altman N (2014) Points of sig-
https://doi.org/10.1038/nprot.2017.151 nificance: two-factor designs. Nat Methods 11:
27. Prince JT, Marcotte EM (2006) Chro- 1187–1188
matographic alignment of ESI-LC-MS proteo- 36. Altman N, Krzywinski M (2014) Points of sig-
mics data sets by ordered bijective interpolated nificance: sources of variation. Nat Methods
warping. Anal Chem 78. https://doi.org/10. 12:5–6
1021/ac0605344
37. Blainey P, Krzywinski M, Altman N (2014)
28. Stanstrup J, Broeckling C, Helmus R et al Points of significance: replication. Nat Meth-
(2019) The metaRbolomics toolbox in Bio- ods 11:879–880
conductor and beyond. Metabolites 9:200.
https://doi.org/10.3390/metabo9100200 38. Altman N, Krzywinski M (2015) Points of sig-
nificance: split plot design. Nat Methods 12:
29. Martens L, Chambers M, Sturm M et al (2011) 165–166
mzML – a community standard for mass spec-
trometry data. Mol Cell Proteomics 10: 39. Haug K, Salek RM, Conesa P et al (2013)
R110.000133 MetaboLights – an open-access general-pur-
pose repository for metabolomics studies and
30. Wilhelm M, Kirchner M, Steen JAJ, Steen H associated meta-data. Nucleic Acids Res 41.
(2012) mz5: space- and time-efficient storage https://doi.org/10.1093/nar/gks1004
of mass spectrometry data sets. Mol Cell Pro-
teomics 11. https://doi.org/10.1074/mcp. 40. Libiseller G, Dvorzak M, Kleb U et al
O111.011379 (2015) IPO: a tool for automated optimization
of XCMS parameters. BMC Bioinformatics 16:
31. Bouyssié D, Dubois M, Nasso S et al 118. https://doi.org/10.1186/s12859-015-
(2015) MzDB: a file format using multiple 0562-8
indexing strategies for the efficient analysis of
large LC-MS/MS and SWATH-MS data sets.
Chapter 6
Abstract
Untargeted analysis by LC-MS is a valuable tool for metabolic profiling (metabonomics/metabolomics),
and applications of this technology have grown rapidly over the past decade. LC-MS offers advantages of
speed, sensitivity, relative ease of sample preparation, and large dynamic range compared to other platforms
in this role. However, like any analytical approach, there are still drawbacks and challenges that have to be
overcome, some of which are being addressed by advances in both column chemistries and instrumentation.
In particular, the combination of LC-MS with ion mobility offers many new possibilities for improved
analyte separation, detection, and structural identification. There are many untargeted LC-MS approaches
which can be applied to metabolic phenotyping, and these usually need to be optimized for the type of
sample, the nature of the study, or the biological question. Some of the main LC-MS approaches for
untargeted metabolic phenotyping are described in detail in the following protocol.
Key words LC-MS, Mass spectrometry, Ion mobility, Liquid chromatography, Untargeted metabolic
profiling
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_6,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
109
110 Ian D. Wilson and Elizabeth Want
Fig. 1 Improved spectral quality due to enabling IMS for two co-eluting lipid species. The elevated collision
energy spectra of these co-eluting lipids (DG [36:5] and SM [32:1]) without the incorporation of IM separation
(a and b) contain the same collection of fragment ions. IM separation prior to CID produces cleaner elevated
collision energy spectra (c and d) improving assignment of fragment ions to the lipid precursor ions increasing
confidence in possible lipid identifications. (From King et al. [32]. With permission)
Fig. 2 LC-MS untargeted workflow showing the main steps of the process. These steps are described in this
protocol
2 Materials
2.6 Data Processing 1. Vendor supplied software or freeware (e.g., XCMS, MZMine
Software 3, MSDial) for pre-processing of raw data. This includes peak
detection, alignment, and normalization (Subheading 3.6).
2. R, Python and associated software packages.
3. Excel or similar.
4. SIMCA, MATLAB, or similar for multivariate analysis.
5. PRISM or JMP or similar for univariate analysis.
114 Ian D. Wilson and Elizabeth Want
3 Methods
Fig. 3 Selected portion of a BPI chromatogram showing serum QC samples overlaid. Inset is an example of
monitoring ion intensity of a specific metabolite over the QC samples
Sample Preparation
1. Urine [24]
(a) RP LC-MS.
– Aliquot an appropriate volume of urine, e.g.,
50–100 μL into a labeled Eppendorf tube.
– Dilute urine sample with appropriate volume of
LC-MS grade water (1:1 for human samples, 1:3 for
rodent and canine samples).
– Centrifuge at 12,000×g for 10 min to remove
precipitate.
– Aliquot sample into clean Eppendorf tube.
– Prepare quality control sample as described in step 6.
– -Aliquot appropriate volume of sample into well plate
or MS vials.
(b) HILIC LC-MS.
– Aliquot appropriate volume of urine, e.g., 50–100 μL
into a labeled Eppendorf tube.
– Dilute urine sample with appropriate volume of aceto-
nitrile (1:1 for human samples, 1:3 for rodent and
canine samples).
– Centrifuge at ~13,000×g for 10 min to remove
precipitate.
– Aliquot sample into clean Eppendorf tube.
– Prepare quality control sample as described in step 6.
– Aliquot appropriate volume of sample into well plate or
MS vials.
2. Serum/Plasma
(a) The key to extracting metabolites from serum/plasma is
the precipitation of proteins. Here, a simple method
involving protein precipitation with methanol is
suggested [25].
Method 1: Preparation for RP LC-MS using methanol
extraction
(b) Aliquot an appropriate volume of serum/plasma into a
labeled Eppendorf tube.
(c) Add cold methanol in a ratio of at least 3:1, e.g., 150 μL
methanol: 50 μL plasma, and vortex for 15 s.
(d) Leave at -20 °C for at least 20 min.
(e) Centrifuge at 14–16,000×g for 10 min.
(f) Remove supernatant to clean labeled Eppendorf tube.
Metabolic Phenotyping by LC-MS 117
3.1 Preparation of Prepare a fresh stock solution of sodium formate; 0.1 mg/mL in
LC-MS Calibration water. Add 1 mL of stock solution to 9 mL of isopropanol (IPA) to
Solution obtain a working solution of 0.01 mg/mL. This solution can be
stored at 4 °C for a number of weeks. Note that an alternative
calibration solution can be used instead, and many instrument
vendors will recommend one.
13. Randomize the study samples, add to the sample analysis list,
and inject a QC sample every 5–10 samples. Aim for at least
10% QC samples per study batch in order to be able to perform
appropriate statistics. End the sample batch with a QC sample,
followed by a blank sample to assess carryover, the QC dilution
series, and a final system suitability/test mix injection to assess
system changes over the run.
14. Injection volume will depend on sample type and instrumenta-
tion, but usually ranges from 2 to 10 μL. Use the QC samples
to assess the appropriate injection volume.
15. When confident that the performance of the system is accept-
able, based on the data from the SST(s) and the conditioning
QC samples, start the analysis of the study samples, including
within-run QC samples.
3.2.1 LC Gradients 1. Example UHPLC gradients for urine, serum, and tissue sam-
ples by both RP and HILIC analysis are shown in Tables 1, 2, 3,
4, 5, 6, and 7. The column dimensions and details are also
provided. Note that these gradients will require optimization
for different column chemistries, lengths, and particle sizes.
3.3 MS Settings Example MS settings, including cone and capillary voltage, are
shown in Table 8 below. Note that these settings are for Waters
Xevo Q-ToF mass spectrometers and will require optimization for
each specific mass spectrometer. However, they can be used as
guidelines, based on an ultraperformance liquid chromatography
(UHPLC) setup with a flow rate of 400–500 μL/min, and a Q-ToF
mass spectrometer.
Table 1
Example gradient for the RP analysis of urine samples
Table 2
Example gradient for the HILIC analysis of urine samples
Table 3
Example gradient for the RP analysis of plasma/serum samples
Table 4
Example gradient for the HILIC analysis of plasma/serum samples and
aqueous tissue extracts
Table 5
Example gradient for the RP analysis of tissue samples—aqueous extract
3.5 Data It is beyond the scope of this protocol to provide a detailed descrip-
Preprocessing tion of data analysis. However, the key data analysis steps pertaining
to LC-MS untargeted analysis for metabolic profiling studies are
described briefly in steps 1–6. There are many different commercial
Metabolic Phenotyping by LC-MS 123
Table 6
Example gradient for the RP analysis of tissue samples—organic extract
Table 7
Example gradient for the HILIC analysis of tissue samples—aqueous
extracts
Table 8
Example MS settings for positive and negative mode ESI analysis using a Q-ToF mass spectrometer
Parameter Setting
Capillary voltage 1–3 kV electrospray (ESI)+, 1–2.5 kV ESI-mode
Cone voltage 30 V
Source temperature e.g. 120 °C
Desolvation temperature e.g. 350 °C
Cone gas flow 25 L/h
Desolvation gas flow 900 L/h
The main parameters are shown. These will need to be optimized for each instrument and ESI mode and are to some
extent dependent on solvent system and mobile phase flow rate
3.6 Data Analysis Once the metabolite table has been produced, this can be exported
into appropriate software, e.g., excel for further analysis, e.g., QC
CV filtering (covered in more detail in Chapter 3).
[28, 29]. However, such databases must be used with care, as many
of the database “hits” that are suggested based on, e.g., molecular
mass alone may, on close inspection, turn out to be biologically
implausible (e.g., pharmaceuticals, pesticides, or unlikely
phytochemicals) [30].
Structural identification is, initially, usually based on MS/MS
experiments that will, hopefully, result in the production of charac-
teristic fragmentation patterns. These fragments will provide high
mass accuracy information in addition to that obtained from the
MS spectrum alone. Such MS/MS spectra can then also be used in
database searches, helping to narrow the list of potential candidates
down from those suggested using MS data alone. In addition,
chromatographic properties and ion mobility-derived collision
cross-section data can be used to support (or disprove) putative
identifications. Untargeted MS/MS data (e.g., DDA experiments)
can be collected on QC conditioning samples or any QC sample
injected throughout the run.
The level of identification required also needs to be considered
against the suggested criteria of, e.g., the Metabolomics Society
Initiative (MSI) which offers a number of levels of increasing cer-
tainty of metabolite identity supported by experimental data [31].
Targeted MS/MS is usually performed on selected samples
containing the metabolite feature(s) of interest in high abundance.
These targeted assays would be performed after univariate or mul-
tivariate analysis has been performed and potential biomarkers
selected. Many of the commonly used databases allow for upload-
ing of MS/MS fragment ion information, e.g., HMDB and
METLIN. HMDB will output predicted MS/MS spectra which
can be visualized against the experimental MS/MS spectra
obtained during a study.
4 Notes
References
1. Holmes E, Loo RL, Stamler J et al (2008) 95(22):8433–8442. https://doi.org/10.
Human metabolic phenotype diversity and its 1021/acs.analchem.2c04343
association with diet and blood pressure. 10. Letertre M, Munjoma MC, Slade SE et al
Nature 453(7193):396–400. https://doi. (2020) Metabolic phenotyping using UPLC–
org/10.1038/nature06882 MS and rapid microbore UPLC–IM–MS:
2. Newgard CB (2017) Metabolomics and meta- determination of the effect of different dietary
bolic diseases: where do we stand? Cell Metab regimes on the urinary metabolome of the rat.
25(1):43–56. https://doi.org/10.1016/j. Chromatographia 83:853–861
cmet.2016.09.018 11. Nassar AF, Wu T, Nassar SF, Wisnewski AV
3. Uthe H, van Dam NM, Hervé MR et al (2021) (2017) UPLC-MS for metabolomics: a giant
Chapter Six – A practical guide to implement- step forward in support of pharmaceutical
ing metabolomics in plant ecology and biodi- research. Drug Discov Today 22(2):463–470.
versity research. In: Pétriacq P, Bouchereau A https://doi.org/10.1016/j.drudis.2016.
(eds) Advances in botanical research, vol 98. 11.020
Elsevier, Amsterdam, pp 163–203 12. Zhao YY, Lin RC (2014) UPLC-MS(E) appli-
4. Saurina J, Sentellas S (2017) Strategies for cation in disease biomarker discovery: the dis-
metabolite profiling based on liquid chroma- coveries in proteomics to metabolomics. Chem
tography. J Chromatogr B 1044–1045:103– Biol Interact 215:7–16. https://doi.org/10.
111. https://doi.org/10.1016/j.jchromb. 1016/j.cbi.2014.02.014
2017.01.011 13. Plumb RS, Gethings LA, Rainville PD et al
5. Haggarty J, Burgess KE (2016) Recent (2023) Advances in high throughput LC/MS
advances in liquid and gas chromatography based metabolomics: a review. TrAC Trends
methodology for extending coverage of the Anal Chem 160:116954
metabolome. Curr Opin Biotechnol 43:77– 14. Siskos AP, Jain P, Römisch-Margl W et al
85. https://doi.org/10.1016/j.copbio.2016. (2016) Interlaboratory reproducibility of a tar-
09.006 geted metabolomics platform for analysis of
6. Markley JL, Brüschweiler R, Edison AS et al human serum and plasma. Anal Chem 89(1):
(2016) The future of NMR-based metabolo- 656–665
mics. Curr Opin Biotechnol 43:34–40. 15. Begou O, Gika HG, Wilson ID et al (2017)
https://doi.org/10.1016/j.copbio.2016. Hyphenated MS-based targeted approaches in
08.001 metabolomics. Analyst 142:3079–3100
7. Chen Y, Xu J, Zhang R, Abliz Z (2016) Meth- 16. Gika H, Virgiliou C, Theodoridis G et al
ods used to increase the comprehensive cover- (2019) Untargeted LC/MS-based metabolic
age of urinary and plasma metabolomes by phenotyping (metabonomics/metabolomics):
MS. Bioanalysis 8(9):981–997. https://doi. the state of the art. J Chromatogr B 1117:
org/10.4155/bio-2015-0010 136–147
8. Wang S, Jiang X, Ding R et al (2022) MS-IDF: 17. Dunn WB, Broadhurst D, Begley P et al (2011)
a software tool for nontargeted identification Human serum metabolome (HUSERMET)
of endogenous metabolites after chemical iso- consortium. Procedures for large-scale meta-
tope labelling based on a narrow mass defect bolic profiling of serum and plasma using gas
filter. Anal Chem 94(7):3194–3202. https:// chromatography and liquid chromatography
doi.org/10.1021/acs.analchem.1c04719 coupled to mass spectrometry. Nat Protocol
9. Liu Y, Zhang Y, Vennekens T et al (2023) 6(7):1060–1083. https://doi.org/10.1038/
MeRgeION: a multifunctional R pipeline for nprot.2011.335
small molecule LC-MS/MS data processing, 18. Gray N, Lewis MR, Plumb RS et al (2015)
searching, and organizing. Anal Chem High-throughput microbore UPLC-MS
Metabolic Phenotyping by LC-MS 129
metabolic phenotyping of urine for large-scale 26. Want EJ, Masson P, Michopoulos F et al
epidemiology studies. J Proteome Res 14(6): (2013) Global metabolic profiling of animal
2714–2721. https://doi.org/10.1021/acs. and human tissues via UPLC-MS. Nat Protoc
jproteome.5b00203 8(1):17–32. https://doi.org/10.1038/nprot.
19. Wilson ID (2015) Metabolic phenotyping by 2012.135
liquid chromatography-mass spectrometry to 27. Veselkov KA, Vingara LK, Masson P et al
study human health and disease. Anal Chem (2011) Optimized preprocessing of ultra-
87(5):2519. https://doi.org/10.1021/acs. performance liquid chromatography/mass
analchem.5b00409. No abstract available spectrometry urinary metabolic profiles for
20. Vorkas PA, Isaac G, Anwar MA et al (2015) improved information recovery. Anal Chem
Untargeted UPLC-MS profiling pipeline to 83(15):5864–5872. https://doi.org/10.
expand tissue metabolome coverage: applica- 1021/ac201065j
tion to cardiovascular disease. Anal Chem 28. Nash WJ, Dunn WB (2019) From mass to
87(8):4184–4193. https://doi.org/10.1021/ metabolite in human untargeted metabolo-
ac503775m mics: recent advances in annotation of metabo-
21. Gika HG, Theodoridis GA, Plumb RS et al lites applying liquid chromatography-mass
(2014) Current practice of liquid spectrometry data. TrAC Trends Anal Chem
chromatography-mass spectrometry in meta- 120:115324. https://doi.org/10.1016/j.
bolomics and metabonomics. J Pharm Biomed trac.2018.11.022
Anal 87:12–25. https://doi.org/10.1016/j. 29. Bocker S (2016) Searching molecular structure
jpba.2013.06.032 databases using tandem MS data: are we
22. Spagou K, Wilson ID, Masson P et al (2011) there yet? Curr Opin Chem Biol 36:1–6.
HILIC-UPLC-MS for exploratory urinary https://doi.org/10.1016/j.cbpa.2016.
metabolic profiling in toxicological studies. 12.010
Anal Chem 83(1):382–390. https://doi.org/ 30. Theodoridis G, Gika H, Raftery D et al (2023)
10.1021/ac102523q Ensuring fact-based metabolite identification
23. Virgiliou C, Sampsonidis I, Gika HG et al in liquid chromatography–mass spectrometry-
(2015) Development and validation of a based metabolomics. Anal Chem 95(8):
HILIC-MS/MS multi-targeted method for 3909–3916
metabolomics applications. Electrophoresis 31. Sumner LW, Amberg A, Barrett D et al (2007)
36:2215–2225 Proposed minimum reporting standards for
24. Want EJ, Wilson ID, Gika H et al (2010) chemical analysis Chemical Analysis Working
Global metabolic profiling procedures for Group (CAWG) Metabolomics Standards Ini-
urine using UPLC-MS. Nat Protocol 5(6): tiative (MSI). Metabolomics 3:211–221.
1005–1018. https://doi.org/10.1038/nprot. https://doi.org/10.1007/s11306-007-
2010.50 0082-2
25. Dunn WB, Broadhurst D, Begley P et al (2011) 32. King AM, Trengove RD, Mullin LG et al
Procedures for large-scale metabolic profiling (2020) Rapid profiling method for the analysis
of serum and plasma using gas chromatography of lipids in human plasma using ion mobility
and liquid chromatography coupled to mass enabled-reversed phase-ultra-high perfor-
spectrometry. Nat Protoc 6(7):1060–1083 mance liquid chromatography/mass spectrom-
etry. J Chromatogr A 1611:460597
Chapter 7
Abstract
Lipidomics has attracted attention in the discovery of unknown biomolecules and for capturing the changes
in metabolism caused by genetic and environmental factors in an unbiased manner. However, obtaining
reliable lipidomics data, including structural diversity and quantification data, is still challenging. Supercrit-
ical fluid chromatography (SFC) is a suitable technique for separating lipid molecules with high throughput
and separation efficiency. Here, we describe a quantitative lipidomics method using SFC coupled with mass
spectrometry. This technique is suitable for characterizing the structural diversity of lipids (e.g., phospho-
lipids, sphingolipids, glycolipids, and glycerolipids) with high quantitative accuracy to understand their
biological functions.
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_7,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
131
132 Hiroaki Takeda et al.
Fig. 1 Chemical diversity of glycerophospholipids and sphingolipids. The abbreviations of lipid classes are
summarized in Table 1
Fig. 2 Characterization of SFC/QqQMS introduced in the chapter. (a) Nomenclature and annotation levels of
each lipid class. (b) Quantification level of each lipid class. (c) Characteristics of the main analytical strategies.
Each performance was determined by three levels. (d) Characteristics of mass spectrometers. The analytical
performance of each instrument was determined by three levels. The characteristic of structural depth is
evaluated when CID is applied
Fig. 3 Features of SFC. (a) Phase diagram. (b) Critical point of each substance. (c) Characteristics of SFC/MS/
MS introduced in the chapter. (d) Schematic diagram of SFC/QqQMS. (e) Lineup of 23 fatty acyl chains
contained in the in-house lipid MRM library [20]
between the polar head groups rather than between the fatty acyl
chains, are widely used. In particular, HILIC has been widely
applied for the separation of lipid classes, including isomers of
glycolipids, such as diastereomers of hexosylceramides (e.g., gluco-
sylceramide (GlcCer) vs. galactosylceramide (GalCer)) [13] and
positional isomers of gangliosides (e.g., GD1a vs. GD1b)
[14]. However, the analytical throughput is low because of the
time-consuming formation of a water-rich layer. Reversed-phase
LC (RPLC) is suitable for untargeted lipidomics because it can
separate lipid molecules based primarily on the differences in the
fatty acyl chains (i.e., the number of carbons and unsaturated
bonds) and slightly polar head groups. However, normalization of
ionization efficiency by ISDs in RPLC/MS is difficult unless lipi-
dome isotope labeling of yeast is used [15], because the degree of
ionization suppression is different for each lipid molecule.
Recently, supercritical fluid chromatography (SFC) has been
used to separate a wide range of lipid molecules [16–18]. In partic-
ular, supercritical carbon dioxide (CO2) is suitable as a mobile
phase for lipid analysis because of its low polarity (similar to that
of n-hexane), low cost, ease of handling, non-toxicity, and
non-flammability (Fig. 3a, b). Because the polarity of the mobile
phase is changed by gradually increasing the composition of the
modifier solvent (e.g., methanol), columns for which the stationary
phase has a high polarity are suitable for SFC. Dissolution solvent
and ionization efficiency result in unique behaviors compared to
that of HILIC because the major mobile phase composition is
supercritical CO2 [16, 19]. Because the weak solvent for HILIC-
based lipid separation has a high content of acetonitrile, it is neces-
sary to dissolve the sample extract in a weak solvent to prevent peak
136 Hiroaki Takeda et al.
Fig. 4 Lipid separation using the DEA column. (a) Separation of lipid classes using the DEA column. The dot
plot represents the gradient condition of the modifier solvent. (b) MS separation of the constituent PC
molecules using the MRM mode. Fragments of fatty acyl chains were selected in the negative-ion mode.
(c) Separation of positional isomers using the DEA column. (The figures are re-used from ref. [20] with minor
modifications)
Fig. 5 Quantification of lipid classes and their constituent molecules. (a) Strategy for lipid quantification. CF
refers to the correction factor determined by calibration curve of charged aerosol detector (CAD). (b)
Separation of lipid molecules by using the DEA column in the rabbit plasma extracts. The number in
parentheses represents the number of lipids detected. (c) Evaluation of quantitative performance by spike-
and-recovery test. (d) Chromatogram of the HepG2 cell extracts by SFC-CAD. Each peak was annotated by
QqQMS. Ceramide phosphoethanolamine (Cer PE) was used as the ISD. (e) Comparison of quantitative levels
between QqQMS and CAD. Lipid extract of HepG2 cells was used for comparison. The figures are re-used from
ref. [20, 21] with minor modifications
2 Materials
2.1 Sample 1. Salt solution (d = 1.006 g/mL): Add sodium chloride (11.4 g)
Preparation and ethylenediaminetetraacetic acid (EDTA) disodium salt
(0.1 g) to a volumetric flask and dissolve them by distilled
2.1.1 Lipoproteins
water (500 mL) and 1 mol/L sodium hydroxide solution
(1 mL). Add distilled water to bring the volume to 1000 mL,
then fine-tune the solution by adding small increments of
distilled water.
2. Salt solution (d = 1.182 g/mL): Add sodium bromide
(24.98 mg) to the salt solution (d = 1.006 g/mL) (100 mL).
3. Polypropylene tube for ultracentrifugation (4.7 mL).
4. Pasteur pipette.
5. Beckman CentriTube Slicer for micro ultracentrifuge tubes.
6. Ultracentrifuge with a fixed-angle rotor.
2.2 Lipid Extraction 1. Methanol for LC/MS: Because methanol is used for lipid
extraction and as a mobile phase for SFC, LC/MS-grade meth-
2.2.1 Liquid-Liquid
anol is recommended to reduce the MS background noise.
Extraction
2. Distilled water for LC/MS: Same precautions as for methanol.
3. Chloroform for high performance LC (HPLC): As chloroform
is converted to toxic phosgene in the presence of light, it
should be stored in the dark.
4. Eppendorf Safe-Lock Tube (2.0 mL): Because Eppendorf
tubes are made of polypropylene, they exhibit chemical resis-
tance and temperature stability. In addition, plasticizers that
increased the MS background were not detected when organic
solvents were used.
5. Micropipettes (20, 200, and 1 mL): Micropipettes with O-ring
materials, which undergo small volumetric expansion when
dispensing organic solvents, and internal springs, which are
resistant to organic solvents, are suitable for accurate dispens-
ing and for avoiding unnecessary contamination.
6. Cooling rack for microtubes: Prior to lipid extraction, the
sample rack should be cooled in a freezer at -30 °C.
140 Hiroaki Takeda et al.
2.3 Lipidome 1. Brown vial for autosampler: Brown vials are useful for protect-
Analysis ing light-sensitive lipids.
2.3.1 Sample 2. Inert glass inserts for sample vials: The silanol groups on the
Preparation glass were converted to hydrophobic groups to prevent the
adsorptive interaction of polar compounds. Because the
amount of affected lipids depends on the contact area with
the glass, the inert glass insert can effectively reduce adsorp-
tion, particularly for the analysis of sample extracts with low
concentrations.
3. Vial cap and slit septum: The slit septum was effective for
inserting the injection needle. When storing the sample vial
after lipid analysis, it is recommended to replace it with a vial
cap without a slit to prevent the volatilization of the sample
elution.
Quantitative Lipidomics Using SFC/MS 141
2.3.2 SFC 1. Nexera UC system (Shimadzu Co., Kyoto, Japan). The back
pressure regulator can be operated without splitting the lines
after passing through the column.
2. CO2 (99.9% grade).
3. Body scale: Because a heavy CO2 bomb is placed on the scale, a
thin scale is recommended.
4. Torus DEA column (3.0 × 100 mm, 1.7 μm, Waters Co.,
Milford, MA): The DEA column can be used for the simulta-
neous analysis of different phospholipids by effectively setting
up scheduled MRM transitions, because it can separate phos-
pholipid classes at the baseline while maintaining good peak
shapes of acidic phospholipids (Fig. 4a) [20].
5. Methanol/water (95/5, v/v) with 0.1% (w/v) ammonium
acetate as a modifier and a makeup pump solvent: LC/MS-
grade methanol and distilled water is recommended to reduce
the MS background noise.
3 Method
3.1 Sample 1. Using a micropipette, transfer plasma or serum (10 μL) into a
Preparation clean tube. Ensure that fibrin is not contaminated (see Notes 1
and 2).
3.1.1 Blood
2. Use the volume to normalize the lipidomics data.
3.1.2 Lipoproteins [23, 1. Add plasma (1.5 mL) to the tubes for ultracentrifugation and
24] spin down at 700 g for 3 min at 4 °C.
2. Add the salt solution (3.2 mL, d = 1.006 g/mL) with a Pasteur
pipette slowly along the tubes and centrifuge the tubes at
375,000 g for 5 h at 4 °C.
3. Collect the upper and lower layers in clean tubes (upper very
low-density lipoprotein (VLDL) fractions, 1.5 mL; lower other
fractions, 3.2 mL). The microslicer is recommended to split the
tubes without contaminating the fractions.
4. If necessary, the fractions are centrifuged again using the salt
solution with the same density (d = 1.006 g/mL) to obtain the
desired fractions.
5. Collect low-density lipoprotein (LDL) and high-density lipo-
protein (HDL) fractions sequentially by repeating the above
steps (VLDL, d < 1.006 g/mL; LDL, 1.019 g/
mL < d < 1.063 g/mL; HDL, d > 1.063 g/mL). Add the
salt solution (d = 1.182 g/mL) until the desired concentration
144 Hiroaki Takeda et al.
3.1.3 Cells 1. Cultivate cells of interest (e.g., cell lines and primary cells) at
37 °C in an incubator. Among dishes, use one dish for cell
counting.
2. Count the number of cells using the cell counter.
3. Wash the cells with PBS.
4. Quench by adding cold methanol and collect the cells in a clean
tube. Care should be taken to not evaporate the methanol by
performing this step promptly.
5. Use the number of cells to normalize the lipidomics data.
3.1.5 Tissues and Organs 1. Dissect the animals as soon as possible and immediately freeze
the collected organs using liquid N2. Be aware that the lipid
metabolism may change depending on the time since death
(see Note 3).
2. If necessary, lyophilize the organs in the dark.
3. Homogenize the dried organs by ball milling at 30 ~ 50 Hz for
2 min after freezing in liquid N2.
4. Collect a sample (5 mg) of the crushed organs into a clean tube.
Attention should be paid to the precision of the electronic
balance.
5. Use dry or wet weight to normalize the lipidomics data.
3.1.6 Tissue Sections 1. Prepare a thin frozen section (approximately 10–20 μm) and
[40] mount it on the glass slide for laser microdissection. Use dry ice
powder to freeze whole brains while preventing sample
cracking.
2. Collect the region of interest (ROI) by laser microdissection on
a cap of either a polymerase chain reaction tube or a 96-well
Quantitative Lipidomics Using SFC/MS 145
3.3 Lipidome 1. Deaerate the mobile phase solvents sufficiently before analysis
Analysis to prevent unexpected pump problems or signal noise caused
by air bubbles. Replace the old mobile phase with a new one
3.3.1 Equilibration of SFC
within 2 weeks to prevent compositional changes caused by
volatilization or microbial contamination.
2. Prepare methanol as a needle wash solvent. Precautions are
same as for modifier and makeup pump solvents (see Note 8).
3. Measure the weight of the CO2 bomb periodically and note
that the remaining amount of CO2 should not be less than
10 kg. An unstable retention time of compounds may be
caused if the correct pressure is not reached.
4. Purge the autosampler and pumps to remove air.
146 Hiroaki Takeda et al.
3.3.2 Equilibration of MS 1. Equilibrate the MS status by using the MS method for lipi-
dome analysis. The QqQMS conditions are as follows: polarity,
both positive- and negative-ion modes; electrospray voltage,
4.0 kV; desolvation line temperature, 250 °C; heat block tem-
perature, 400 °C; nebulizing gas (N2) flow, 3.0 L/min; drying
gas (N2) flow, 10.0 L/min; collision-induced dissociation
(CID) gas (argon) pressure, 0.23 MPa, detector voltage,
2.16 kV, dwell time, 1 ms; pause time, 2 ms; polarity switching
time, 15 ms.
2. Perform autocalibration in both the positive- and negative-ion
modes to check the resolution, mass accuracy, and sensitivity if
problems occur. This will be helpful in troubleshooting if it is
unclear whether the problems originate from the SFC or
the MS.
3.3.3 Sample Analysis 1. Prior to lipid analysis, optimize collision energy for each lipid
class using the internal standards. Flow injection mode is effi-
cient for MRM optimization.
2. Analyze the ISD mixture to check the data quality (e.g., reten-
tion time, peak width, and intensity). Replace the column if the
retention time and peak width are not reproducible. Clean the
Quantitative Lipidomics Using SFC/MS 147
3.4.2 Multi- 1. When using the MS data except for .lcd files, convert them to .
ChromatoAnalysT (for txt files using the MSconvertGUI software.
Targeted Lipidomics) 2. Drag and drop the .txt files to Multi-ChromatoAnalysT and
read the data.
3. Perform peak alignment for every sample and select the start
and end points of the peaks for the calculation of the peak area
(see Notes 15–17).
4. Starting from the second time, insert the retention time infor-
mation determined by the preceding analysis to pick the peaks
automatically. Alternatively, add new data to the previous
data set.
5. Input the limit of quantification (peak intensity) and the con-
centrations of ISDs to quantify lipid molecules automatically.
6. Output the data matrix to an Excel sheet to interpret the data.
148 Hiroaki Takeda et al.
4 Notes
1. Maintain r.t. and time when leaving blood still to prepare the
serum sample.
2. Use the same anticoagulant (e.g., EDTA, citrate, and heparin)
throughout the study. EDTA or heparin is recommended when
analyzing polar metabolites because citrate is a metabolite in
the tricarboxylic acid cycle.
3. Store the dissected biological samples in a freezer at -80 °C
prior to lipid extraction.
4. Maintain the acidic (pH 4) and cold (<4 °C) conditions to
eliminate acyl migration of lysophospholipids during the
extraction process [43]. The extraction of lipids under alkaline
conditions accelerates the acyl migration of lysophospholipids.
5. Non-polar lipid classes such as CEs and TGs should be
extracted using biphasic liquid extraction, because the solubil-
ity of these classes is insufficient when using monophasic liquid
extraction [44].
6. In the analysis of gangliosides, taking into account the dynamic
range of MS, it may be appropriate to use biphasic liquid
extraction, which can further enrich the aqueous layer, rather
than monophasic liquid extraction.
7. Store the dried extract in a freezer at -30 °C until immediately
before lipid analysis.
8. Check carryover levels before lipid analysis. The carryover
levels should be less than 1% to quantify lipids with reliable
accuracy. A mixture of isopropanol, methanol, acetonitrile, and
water (1:1:1:1, v/v/v/v) may be effective for washing the
needle if the carryover levels are serious for some lipids.
Quantitative Lipidomics Using SFC/MS 149
References
1. Han X (2016) Lipidomics for studying metab- spectrometry-based lipidomics. Anal Chim
olism. Nat Rev Endocrinol 12(11):668–679. Acta 1137:156–169. https://doi.org/10.
https://doi.org/10.1038/nrendo.2016.98 1016/j.aca.2020.09.060
2. Xu T, Hu C, Xuan Q, Xu G (2020) Recent 3. Masoodi M, Gastaldelli A, Hyötyl€ainen T,
advances in analytical strategies for mass Arretxe E, Alonso C, Gaggini M, Brosnan J,
150 Hiroaki Takeda et al.
Anstee QM, Millet O, Ortiz P, Mato JM, Yuan M, Zhao XH, Zhou S (2017) Harmoniz-
Dufour JF, Orešič M (2021) Metabolomics ing lipidomics: NIST interlaboratory compari-
and lipidomics in NAFLD: biomarkers and son exercise for lipidomics using SRM 1950-
non-invasive diagnostic tests. Nat Rev Gastro- Metabolites in Frozen Human Plasma. J Lipid
enterol Hepatol 18(12):835–856. https://doi. Res 58(12):2275–2288. https://doi.org/10.
org/10.1038/s41575-021-00502-9 1194/jlr.M079012
4. Harayama T, Riezman H (2018) Understand- 9. Liebisch G, Ahrends R, Arita M, Arita M, Bow-
ing the diversity of membrane lipid composi- den JA, Ejsing CS, Griffiths WJ, Holčapek M,
tion. Nat Rev Mol Cell Biol 19(5):281–296. Köfeler H, Mitchell TW, Wenk MR, Ekroos K
https://doi.org/10.1038/nrm.2017.138 (2019) Lipidomics needs more standardiza-
5. Fahy E, Subramaniam S, Brown HA, Glass CK, tion. Nat Metab 1(8):745–747. https://doi.
Merrill AH Jr, Murphy RC, Raetz CRH, Rus- org/10.1038/s42255-019-0094-z
sell DW, Seyama Y, Shaw W, Shimizu T, 10. Schuhmann K, Moon HK, Thomas H, Acker-
Spener F, van Meer G, VanNieuwenhze MS, man JM, Groessl M, Wagner N, Kellmann M,
White SH, Witztum JL, Dennis EA (2005) A Henry I, Nadler A, Shevchenko A (2019)
comprehensive classification system for lipids. J Quantitative fragmentation model for
Lipid Res 46(5):839–861. https://doi.org/ bottom-up shotgun lipidomics. Anal Chem
10.1194/jlr.E400004-JLR200 91(18):12085–12093. https://doi.org/10.
6. Liebisch G, Fahy E, Aoki J, Dennis EA, 1021/acs.analchem.9b03270
Durand T, Ejsing CS, Fedorova M, 11. O’Donnell VB, Ekroos K, Liebisch G, Wake-
Feussner I, Griffiths WJ, Köfeler H, Merrill lam M (2019) Lipidomics: current state of the
AH Jr, Murphy RC, O’Donnell VB, art in a fast moving field. Wiley Interdiscip Rev
Oskolkova O, Subramaniam S, Wakelam Syst Biol Med 12(1):1466. https://doi.org/
MJO, Spener F (2020) Update on LIPID 10.1002/wsbm.1466
MAPS classification, nomenclature, and short- 12. Han X, Gross RW (2022) The foundations and
hand notation for MS-derived lipid structures. development of lipidomics. J Lipid Res 63(2):
J Lipid Res 61(12):1539–1555. https://doi. 100164. https://doi.org/10.1016/j.jlr.2021.
org/10.1194/jlr.S120001025 100164
7. Cajka T, Fiehn O (2014) Comprehensive anal- 13. von Gerichten J, Schlosser K, Lamprecht D,
ysis of lipids in biological systems by liquid Morace I, Eckhardt M, Wachten D,
chromatography-mass spectrometry. Trends Jennemann R, Gröne HJ, Mack M, Sandhoff
Analyt Chem 61:192–206. https://doi.org/ R (2017) Diastereomer-specific quantification
10.1016/j.trac.2014.04.017 of bioactive hexosylceramides from bacteria
8. Bowden JA, Heckert A, Ulmer CZ, Jones CM, and mammals. J Lipid Res 58(6):1247–1258.
Koelmel JP, Abdullah L, Ahonen L, Alnouti Y, https://doi.org/10.1194/jlr.D076190
Armando AM, Asara JM, Bamba T, Barr JR, 14. Hájek R, Jirásko R, Lı́sa M, Cı́fková E, Holča-
Bergquist J, Borchers CH, Brandsma J, Breit- pek M (2017) Hydrophilic interaction liquid
kopf SB, Cajka T, Cazenave-Gassiot A, chromatography-mass spectrometry character-
Checa A, Cinel MA, Colas RA, Cremers S, ization of gangliosides in biological samples.
Dennis EA, Evans JE, Fauland A, Fiehn O, Anal Chem 89(22):12425–12432. https://
Gardner MS, Garrett TJ, Gotlinger KH, doi.org/10.1021/acs.analchem.7b03523
Han J, Huang Y, Neo AH, Hyötyl€ainen T, 15. Rampler E, Criscuolo A, Zeller M, El Abiead Y,
Izumi Y, Jiang H, Jiang H, Jiang J, Schoeny H, Hermann G, Sokol E, Cook K,
Kachman M, Kiyonami R, Klavins K, Klose C, Peake DA, Delanghe B, Koellensperger G
Köfeler HC, Kolmert J, Koal T, Koster G, (2018) A novel lipidomics workflow for
Kuklenyik Z, Kurland IJ, Leadley M, Lin K, improved human plasma identification and
Maddipati KR, McDougall D, Meikle PJ, Mel- quantification using RPLC-MSn methods and
lett NA, Monnin C, Moseley MA, isotope dilution strategies. Anal Chem 90(11):
Nandakumar R, Oresic M, Patterson RE, 6494–6501. https://doi.org/10.1021/acs.
Peake D, Pierce JS, Post M, Postle AD, analchem.7b05382
Pugh R, Qiu Y, Quehenberger O, Ramrup P,
Rees J, Rembiesa B, Reynaud D, Roth MR, 16. Nováková L, Perrenoud AGG, Francois I,
Sales S, Schuhmann K, Schwartzman ML, Ser- West C, Lesellier E, Guillarme D (2014) Mod-
han CN, Shevchenko A, Somerville SE, ern analytical supercritical fluid chromatogra-
St. John-Williams L, Surma MA, Takeda H, phy using columns packed with sub-2 μm
Thakare R, Thompson JW, Torta F, Triebl A, particles: a tutorial. Anal Chim Acta 824:18–
Trötzmüller M, SJK U, Vuckovic D, Weir JM, 35. https://doi.org/10.1016/j.aca.2014.
Welti R, Wenk MR, Wheelock CE, Yao L, 03.034
Quantitative Lipidomics Using SFC/MS 151
17. Lesellier E, West C (2015) The many faces of Shevchenko A, Torta F, Wakelam MJO, Wheel-
packed column supercritical fluid ock CE, Wenk MR (2018) MS-based lipido-
chromatography—a critical review. J Chroma- mics of human blood plasma: a community-
togr A 1382:2–46. https://doi.org/10.1016/ initiated position paper to develop accepted
j.chroma.2014.12.083 guidelines. J Lipid Res 59(10):2001–2017.
18. Shulaev V, Isaac G (2018) Supercritical fluid https://doi.org/10.1194/jlr.S087163
chromatography coupled to mass 27. Maccarone AT, Duldig J, Mitchell TW,
spectrometry—a metabolomics perspective. J Blanksby SJ, Duchoslav E, Campbell JL
Chromatogr B 1092:499–505. https://doi. (2014) Characterization of acyl chain position
org/10.1016/j.jchromb.2018.06.021 in unsaturated phosphatidylcholines using dif-
19. Fujito Y, Izumi Y, Nakatani K, Takahashi M, ferential mobility-mass spectrometry. J Lipid
Hayakawa Y, Takayama M, Bamba T (2023) Res 55(8):1668–1677. https://doi.org/10.
Understanding the mechanism of CO2-assisted 1194/jlr.M046995
electrospray ionization for parameter optimiza- 28. Šala M, Lı́sa M, Campbell JL, Holčapek M
tion in supercritical fluid chromatography mass (2016) Determination of triacylglycerol regioi-
spectrometry. Anal Chim Acta 1246:340863. somers using differential mobility spectrome-
https://doi.org/10.1016/j.aca.2023.340863 try. Rapid Commun Mass Spectrom 30(2):
20. Takeda H, Izumi Y, Takahashi M, Paxton T, 256–264. https://doi.org/10.1002/rcm.
Tamura S, Koike T, Yu Y, Kato N, Nagase K, 7430
Shiomi M, Bamba T (2018) Widely-targeted 29. Baba T, Campbell JL, Le Blanc JCY, Baker
quantitative lipidomics method by supercritical PRS, Ikeda K (2018) Quantitative structural
fluid chromatography triple quadrupole mass multiclass lipidomics using differential mobil-
spectrometry. J Lipid Res 59(7):1283–1293. ity: electron impact excitation of ions from
https://doi.org/10.1194/jlr.D083014 organics (EIEIO) mass spectrometry. J Lipid
21. Takeda H, Takahashi M, Hara T, Izumi Y, Res 59(5):910–919. https://doi.org/10.
Bamba T (2019) Improved quantitation of 1194/jlr.D083261
lipid classes using supercritical fluid chroma- 30. Zhu Y, Wang W, Yang Z (2020) Combining
tography with a charged aerosol detector. J mass spectrometry with Paternò–Büchi reac-
Lipid Res 60(8):1465–1474. https://doi. tion to determine double-bond positions in
org/10.1194/jlr.D094516 lipids at the single-cell level. Anal Chem
22. Berger TA (2022) The evolution and current 92(16):11380–11387. https://doi.org/10.
state of instrumentation for analytical super- 1021/acs.analchem.0c02245
critical fluid chromatography. J Chromatogr B 31. Takahashi H, Shimabukuro Y, Asakawa D,
1211:123478. https://doi.org/10.1016/j. Yamauchi S, Sekiya S, Iwamoto S, Wada M,
jchromb.2022.123478 Tanaka K (2018) Structural analysis of phos-
23. Takeda H, Koike T, Izumi Y, Yamada T, pholipid using hydrogen abstraction dissocia-
Yoshida M, Shiomi M, Fukusaki E, Bamba T tion and oxygen attachment dissociation in
(2015) Lipidomic analysis of plasma lipopro- tandem mass spectrometry. Anal Chem
tein fractions in myocardial infarction-prone 90(12):7230–7238. https://doi.org/10.
rabbits. J Biosci Bioeng 120(4):476–482. 1021/acs.analchem.8b00322
https://doi.org/10.1016/j.jbiosc.2015. 32. Baba T, Campbell JL, Le Blanc JCY, Baker PRS
02.015 (2017) Distinguishing cis and trans isomers in
24. Takeda H, Izumi Y, Tomita A, Koike T, intact complex lipids using electron impact
Shiomi M, Fukusaki E, Matsuda F, Bamba T excitation of ions from organics mass spec-
(2017) Lipoprotein profiling methodology trometry. Anal Chem 89(14):7307–7315.
based on determination of apolipoprotein con- https://doi.org/10.1021/acs.analchem.
centration. Bioanalysis 9(1):9–19. https://doi. 6b04734
org/10.4155/bio-2016-0234 33. Kanehisa M, Araki M, Goto S, Hattori M,
25. Nishida-Aoki N, Izumi Y, Takeda H, Hirakawa M, Itoh M, Katayama T,
Takahashi M, Ochiya T, Bamba T (2020) Lipi- Kawashima S, Okuda S, Tokimatsu T, Yama-
domic analysis of cells and extracellular vesicles nishi Y (2007) KEGG for linking genomes to
from high- and low-metastatic triple-negative life and the environment. Nucleic Acids Res 36:
breast cancer. Metabolites 10(2):67. https:// 480–484. https://doi.org/10.1093/nar/
doi.org/10.3390/metabo10020067 gkm882
26. Burla B, Arita M, Arita M, Bendt AK, 34. Sud M, Fahy E, Cotter D, Azam K, Vadivelu I,
Cazenave-Gassiot A, Dennis EA, Ekroos K, Burant C, Edison A, Fiehn O, Higashi R, Nair
Han X, Ikeda K, Liebisch G, Lin MK, Loh KS, Sumner S, Subramaniam S (2016) Meta-
TP, Meikle PJ, Orešič M, Quehenberger O, bolomics workbench: an international
152 Hiroaki Takeda et al.
Abstract
The gut’s symbiome, a hidden metabolic organ, has gained scientific interest for its crucial role in human
health. Acting as a biochemical factory, the gut microbiome produces numerous small molecules that
significantly impact host metabolism. Metabolic profiling facilitates the exploration of its influence on
human health and disease through the symbiotic relationship. Fecal metabolomics-based analysis is an
indisputably valuable tool for elucidating the biochemistry of digestion and absorption in the gastrointesti-
nal system, serving as the most suitable specimen to study the symbiotic relationship between the host and
the intestinal microbiota. It is well-established that the balance of the intestinal microbiota changes in
response to various stimuli, both physiological, such as gender, age, diet, and exercise, and pathological,
such as gastrointestinal and hepatic diseases. Fecal samples have been analyzed using widely adopted
analytical techniques, including NMR spectroscopy, GC-MS, and LC-MS/MS. Rat fecal samples are
frequently used and particularly useful substrates for metabolomics-based studies in related fields.
The complexity and diversity of fecal samples necessitate careful and skillful handling to extract metabo-
lites, while avoiding their deterioration, effectively and quantitatively. Several determinative factors, such as
the fecal sample weight to extraction solvent solution volume, the nature and pH value of the extraction
solvent, and the homogenization process, play crucial roles in achieving optimal extraction for obtaining
high-quality metabolic fingerprints, whether for untargeted or targeted metabolomics.
Key words Metabolomics, Metabolic profiling, Metabonomics, Sample preparation, Fecal samples,
Rats, Fecal extract, NMR, GC-MS, LC-MS/MS, Metabolites, Metabolome, Gut microbiome, Gut
microbiota
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_8,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
153
154 Olga Deda et al.
from 2500 years ago that the root of good and bad health lies in the
gut [4].
Gut microbiome is literally a biochemical factory producing
numerous small molecules, derived either from the metabolism of
dietary compounds or from de novo biosynthesis. It contributes to
host metabolism in a great scale since it is estimated that approxi-
mately 50% of fecal and urine metabolites [5, 6] and 10% of the
blood metabolites originated from the gut microbiota [7].
Metabolomics extensively applies to studying how the gut
microbiome affects the host through their symbiotic relationship.
Various analytical methods investigate the metabolic effects result-
ing from gut microbiome disequilibrium (dysbiosis) for biomarker
discovery. Combining metabolomics with high-throughput DNA
and RNA sequencing reveals the link between the metabolome and
gut microbiome in host (human) health and disease [7–21].
While numerous biomarkers associated with the pathophysiol-
ogy of the gastrointestinal tract have been identified through meta-
bolomics analysis [22], there are instances where different
biomarkers were proposed, even when studying similar factors
such as biospecimen, population, or animal models. Several factors
could account for this variation, including experimental design
details, variations in sample collection and preparation, the selec-
tion of analytical parameters, and disparities in data handling and
statistical analysis. Human studies are further complicated by
genetic variations, diverse lifestyles, potential medications, and
other confounding factors, adding to the complexity of disease
biomarker discovery.
lead to reduced levels of short chain fatty acids (SCFAs) and poten-
tially of other volatile compounds [33, 36], thus editing the profile
especially in GC-MS mode.
Among the innumerous metabolite classes, SCFAs are high-
lighted to be of utmost importance and highly correlated to the gut
microbiome symbiosis and its alterations. The determination of
fecal SCFAs by an independent method provides the advantage of
a distinct process, both in sample preparation and analysis, most
suitable for the specific analytes; however, it has the cost of a second
analysis in addition to fecal metabolic profiling.
In addition, more than one extraction cycle could be per-
formed, and combination of extracts could take place [32]. This
practice is not proposed to keep the simplicity of the sample prepa-
ration process.
Finally, depending on the technique by which fecal samples are
analyzed, other conditions and parameters may affect the analysis
performance. For example, in NMR metabolomics-based analysis,
PBS buffer preparation (preparation in D2O, pH value, and salt
addition) appreciably results in collected spectra. The pH value of
PBS noticeably affects the behavior of the chemical shift reference
used. The use of deuterated solvents is not mandatory, but
improves the shimming, mainly in older spectrometers. Further-
more, the oversized signal of a non-deuterated solvent in 1H-NMR
could deteriorate the spectra’s quality. Sample preparation often
includes the addition of sodium azide as preservative and TSP as
chemical shift reference.
Undoubtedly, the derivatization process of the fecal sample is a
determining factor, regarding the GC-MS analysis. Critical factors
such as residual moisture, insufficient derivatization reagents, and
shorter or longer incubation, at lower or higher temperature, could
greatly degrade the obtained chromatograms [37, 38].
Choosing the right extraction solvent, like the mobile phase,
can lead to improved peaks in liquid chromatography. However,
using highly concentrated fecal extracts may cause column blockage
and break down the analytical system.
2 Materials
2.2 Instrumentation 1. 500 MHz (or higher) NMR spectrometer equipped with a
Software 5 mM triple resonance probe at 300 K (or similar). The appro-
priate software to control acquisition of fecal sample (match-
ing, tuning, shimming) to set pulse calibration parameters and
to process.
2. Agilent 7890A GC coupled to a 5975C inert XL EI/CI MSD
with Triple-Axis Detector MS and a CTC-CH 4222 autosam-
pler with an Agilent HP-5 ms (29 m × 250 μm × 0.25 μm) in
split/splitless mode. MSD ChemStation (Agilent Technologies,
California, USA) to acquire and process GC-MS data.
3. ACQUITY UPLC coupled to a Xevo TQD MS System
(Waters, Massachusetts, USA) with an ACQUITY HILIC,
BEH amide column (2.1 × 150 mM, 1.7 μm). Waters Mas-
sLynx® software to collect and process LC-MS/MS data.
3 Methods
3.2.2 Analysis Using GC- 1. 500 mg of smashed stool material are placed in a 2 mL
MS Eppendorf tube.
2. 1 mL of MeOH-CHCl3 1:1 (v/v) is added, in a ratio of 1:
2 fecal sample weight to extraction buffer volume.
3. Vortex-mixing is followed for 2 min.
4. The mixture is homogenized by sonication for 10 min.
5. Organic fecal extract is centrifuged for 20 min at 4 °C
(18,000 × g).
6. 100 μL of the supernatant are placed in a glass GC-MS vial and
evaporated under nitrogen stream.
7. 20 μL methoxyamine hydrochloride in pyridine (40 mg/mL)
and 180 μL MSTFA (1% trimethylchlorosilane [TMCS]) are
added followed by incubation for 90 min at 28 °C and for
30 min at 37 °C, respectively (see Note 8).
8. 5 μL of 1 mM solution of 2-fluorobiphenyl are added, as
internal standard. An alternative option is the addition of
5 μL pentadecan (in pyridine) solution. Before the addition of
internal standard, the vial should be allowed to fall to ambient
temperature (see Note 9).
9. The gas chromatographic analysis (GC-MS) of the fecal
extracts is performed as described in detail, in a previously
published paper from our group [26].
(see Note 10) [26, 34]. Acetonitrile better extracts a great wide
of metabolites, but 1-proponol preferably extracts amino acids
related to bowel disease.
3. Vortex-mixing is followed for 2 min.
4. The mixture is homogenized by sonication for 10 min.
5. Ultracentrifugation for 30 min at 4 °C (20,000 × g) is
followed.
6. Supernatants are filtered through syringe filters PTFE 0.22 μm.
7. The filtrate should be received in an Eppendorf tube and not
immediately in LC-MS vial which has a narrower opening.
8. The required sample volume, depending on the number of
subsequent injections, is added in the LC-MS vial insert prior
to analysis.
9. Analyze by targeted HILIC-MS/MS with a multi-analyte
method (e.g., as the one described in Chap. 10, in the present
book and in a previous published paper [39]) or an untargeted
method.
The applied protocols for analytical techniques used are
presented briefly in Fig. 1.
Fig. 1 Rat fecal sample preparation protocols for metabolomics-based analysis using NMR spectroscopy,
GC-MS, and LC-MS/MS
Rat Fecal Metabolomics-Based Analysis 161
4 Notes
References
1. Guinane CM, Cotter PD (2013) Role of the 5. Zheng X, Xie G, Zhao A et al (2011) The
gut microbiota in health and chronic gastroin- footprints of gut microbial-mammalian co-me-
testinal disease: understanding a hidden meta- tabolism. J Proteome Res 10(12):5512–5522
bolic organ. Ther Adv Gastroenterol 6(4): 6. Wang G, Huang S, Wang Y et al (2019) Bridg-
295–308 ing intestinal immunity and gut microbiota by
2. O’Hara AM, Shanahan F (2006) The gut flora metabolites. Cell Mol Life Sci 76(20):
as a forgotten organ. EMBO Rep 7(7): 3917–3937
688–693 7. Lyu Q, Hsu C-C (2018) Can diet influence our
3. Li X, Liu L, Cao Z et al (2020) Gut microbiota health by altering intestinal microbiota-derived
as an “invisible organ” that modulates the fecal metabolites? mSystems 3(2)
function of drugs. Biomed Pharmacother 8. Nicholson JK, Holmes E, Kinross J et al (2012)
121:109653 Host-gut microbiota metabolic interactions.
4. Sekirov I, Russell SL, Antunes LCM, Finlay BB Science 336(6086):1262–1267
(2010) Gut microbiota in health and disease.
Physiol Rev 90(3):859–904
162 Olga Deda et al.
aging: a longitudinal study in rats. Bioanalysis ethyl chloroformate derivatization. Anal Bio-
9(1):21–36 chem 393(2):163–175
35. Bezabeh T, Somorjai RL, Smith ICP (2009) 38. Gao X, Pujos-Guillot E, Sébédio J-L (2010)
MR metabolomics of fecal extracts: applica- Development of a quantitative metabolomic
tions in the study of bowel diseases. Magn approach to study clinical human fecal water
Reson Chem 47(S1):S54–S61 metabolome based on trimethylsilylation deriv-
36. Monleon D, Garcia-Valles R, Morales JM et al atization and GC/MS analysis. Anal Chem
(2014) Metabolomic analysis of long-term 82(15):6447–6456
spontaneous exercise in mice suggests 39. Virgiliou C, Sampsonidis I, Gika HG,
increased lipolysis and altered glucose metabo- Raikos N, Theodoridis GA (2015) Develop-
lism when animals are at rest. J Appl Physiol ment and validation of a HILIC-MS/MS mul-
1985 117(10):1110–1119 titargeted method for metabolomics
37. Gao X, Pujos-Guillot E, Martin J-F et al (2009) applications. Electrophoresis 36(18):
Metabolite analysis of human fecal water by gas 2215–2225
chromatography/mass spectrometry with
Chapter 9
Abstract
A protocol for the preparation of tissue extracts for the targeted analysis ca. 150 polar metabolites,
including those involved in central carbon metabolism, is described, using a reversed phase ion pair U
(H)PLC-MS method. Data collection enabled in high-resolution mass spectrometry detection provides
highly specific and sensitive acquisition of metabolic intermediates with wide range physicochemical
properties and pathway coverage. Technical aspects are discussed for method transfer along with the
basic principles of sample sequence setup, data analysis, and validation. At last general comments are
given to help the assessment of data quality and system performance.
Key words Polar metabolites, Ion pair chromatography, Mass spectrometry, Metabolite profiling
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_9,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
165
166 Filippos Michopoulos
2 Materials
2.1 Preparation of Water of chromatographic grade purity (18.2 MΩ), organic sol-
Analytical Standards, vents (acetonitrile, methanol, isopropanol), and mobile phase addi-
Test Mixture, and tives (acetic acid, tributylamine) of LC MS grade should be used.
Infusion Solution The same quality criteria must be applied for analytical standards to
obtain the highest available purity.
1. For all analytical standards listed in Table 1, the appropriate
amount is weighed in a 1.5 mL Eppendorf tube to result in a
1 mL 50 mM solution in MeOH/H2O 50/50 v/v (stock
solution). For amino acids and nucleotides, solubility is
improved by the addition of small amounts of HCl, while
other classes of metabolite may require small amounts of
1 M NaOH.
2. Stock solutions must be diluted 1/100 v/v with HPLC grade
water to produce a solution (dilution A) with a nominal con-
centration of 500 μM which will be further diluted 1/50 with
either HPLC water (dilution B) or MeOH/H2O 50/50 v/v
(infusion solution) to result to a final standard concentration of
10 μM.
3. Stock solutions, dilutions A and B should be transferred to
1.8 mL cryovials and stored at -20 °C as reference material
for future use.
4. An aliquot of 30 μL of the appropriate stock solution standards
(see Table 1) can be mixed in a 15 mL Falcon tube and
concentration adjusted to 100 μM with HPLC water before
dividing into 10 μL aliquots (test mixture 1–6) and stored at -
20 °C for in batch validation.
Table 1
List of metabolites measured with the current methodology, mass spectrometer parameters, and text
mixtures composition
(continued)
Ion Pair Chromatography for Endogenous Metabolite LC-MS Analysis in Tissue. . . 169
Table 1
(continued)
(continued)
170 Filippos Michopoulos
Table 1
(continued)
(continued)
Ion Pair Chromatography for Endogenous Metabolite LC-MS Analysis in Tissue. . . 171
Table 1
(continued)
(continued)
172 Filippos Michopoulos
Table 1
(continued)
2.3 Tissue Extraction 1. In a 0.5 L solvent bottle, add 100 mL of HPLC water, 200 mL
Solvent of MeOH, and 200 mL acetonitrile before manual shake to
ensure mixing.
3 Methods
3.1 Tissue Extraction Tissue samples must be processed and extracted from frozen to
and Sample reduce endogenous metabolite degradation. For the present pro-
Preparation tocol, we used a combined extraction homogenization approach
that is performed using a Precellys 24 system with an attached
temperature control unit. For soft tissue samples, we advise the
use of CK14 tubes, while for harder tissue the CK28R format is
more appropriate. Extraction and homogenization is achieved in
the appropriate tube format from a 50 mg of frozen tissue with
1 mL ACN/MeOH/H2O 40/40/20 v/v/v. Tubes must be
shaken at 5000 rpm for 20 sec, and the process must be repeated
three times with an intermittent pause of 30 sec between each
repeat. While the extraction and homogenization is taking place,
the extraction chamber unit must be at the lowest possible temper-
ature and not more than 10 °C. A clear supernatant is obtained after
centrifugation at 10691 g for 5 min at 0 °C, and this is transferred
to cryovial for storage at a minimum -20 °C until analysis. The
extraction-homogenization procedure is repeated with fresh sol-
vent as described above, and the resulted clear supernatant is com-
bined int the same cryovial.
Prior to LC-MS analysis, a minimum of a 1 in 10 dilution with
HPLC water must be performed to avoid excess carryover of highly
abundant metabolites across injections. For different tissues, an
174 Filippos Michopoulos
3.2 U(H)PLC Sample analysis according to the present protocol was performed
Separation on Thermo Ultimate 3000 RS pump combined with an Ultimate
3000 autosampler operating at 4 °C using a Acquity HSS T3 U
(H)PLC column (Waters Corp, 2.1 × 100 mM, 1.8 μm particle
size) with a column temperature maintained at 60 ± 0.5 °C. The
analysis was performed using a 5 μL injection of sample with
gradient elution at flow rate of 400 μL/min with a binary solvent
mixing schedule of 0 min, 0% B 0.5 min, 0% B; 4 min, 5% B; 6 min,
5% B; 6.5 min, 20% B; 8.5 min, 20% B; 14 min, 55% B; 15 min,
100% B; 17 min, 100% B; 18 min, 0% B; and 21 min 0% B.
3.3 Mass Spectrometric data was acquired on ABSCIEX 6600 triple TOF
Spectrometry instrument operating in a negative electrospray ionization acquisi-
tion mode. With high-resolution acquisition, there is no need to
optimize metabolite detection parameters such as declustering
potential, entrance potential, collision energy, and collision exit
potential. The accurate mass information given in Table 1 is calcu-
lated based on the molecular formula for each metabolite. Injection
(5 μL) of an authentic standard at concentration of 5 μM is suffi-
cient to confirm adequate retention on the chromatographic col-
umn and elution time (RT). Interface source parameters were
optimized at the flow rate (0.4 mL/min) of the intended chro-
matographic separation. For this purpose, 10 μL/min of each
metabolite (dilution A) was infused via the syringe pump to a
T-connector where it was combined with 0.39 mL/min 90%
mobile phase A via the UHPLC system delivering to ESI source a
final metabolite concentration of 12.5 μM. The optimal ESI source
parameters were as follows: ion spray voltage -4.5 kV, temperature
500 °C, collision energy -5 V, curtain gas 35, ion source gas(1) 60,
Ion Pair Chromatography for Endogenous Metabolite LC-MS Analysis in Tissue. . . 175
ion source gas(2) 50, and declustering potential -10 V. Gas values
are arbitrary units. For co-eluting and isobaric species where the
TOF MS scan does not provide sufficient resolving capacity, a
multiple reaction monitoring high-resolution scan (MRMHR)
can be used to explore potential differential fragmentation and
resolve the molecules of interest.
3.5 Data Analysis Raw spectrometric data are processed with MuliQuant software to
obtain peak areas for each of the detected metabolites across the
sample set using the TOF MS high-resolution data with 20mDa
extraction window. The first step prior to peak integration is the
visual confirmation of the retention time for each metabolite peak
by comparing the trace obtained for the test mixture, spiked QC,
and the individual samples. Given the inability of obtain biological
matrix free from the endogenous metabolites being determined,
the solvent blank injection data can be used to define the
Ion Pair Chromatography for Endogenous Metabolite LC-MS Analysis in Tissue. . . 177
4 General Remarks
Fig. 1 Typical extracted ion chromatograms of a reference metabolite mixture injection obtained following the
proposed chromatographic separation
References
1. Bales JR, Bell JD, Nicholson JK, Sadler PJ, 10. Lenz EM, Wilson ID (2007) Analytical strate-
Timbrell JA, Hughes RD et al (1988) Meta- gies in Metabonomics. J Proteome Res 6(2):
bolic profiling of body fluids by proton NMR: 443–458
self-poisoning episodes with paracetamol (acet- 11. Lindon JC, Nicholson JK (2008) Analytical
aminophen). Magn Reson Med 6(3):300–306 technologies for metabonomics and metabolo-
2. Nicholson JK, Wilson ID (1989) High resolu- mics, and multi-omic information recovery.
tion proton magnetic resonance spectroscopy TrAC Trends Anal Chem 27(3):194–204
of biological fluids. Prog Nucl Magn Reson 12. Theodoridis G, Gika HG, Wilson ID (2008)
Spectrosc 21(4–5):449–501 LC-MS-based methodology for global metab-
3. Barzilai A, Horowitz A, Geier A, Degani H olite profiling in metabonomics/metabolo-
(1991) Phosphate metabolites and steroid hor- mics. TrAC Trends Anal Chem 27(3):251–260
mone receptors of benign and malignant breast 13. Shulaev V (2006) Metabolomics technology
tumors. A nuclear magnetic resonance study. and bioinformatics. Brief Bioinform 7(2):
Cancer 67(11):2919–2925 128–139
4. Gavaghan CL, Holmes E, Lenz E, Wilson ID, 14. Smith CA, Want EJ, O’Maille G, Abagyan R,
Nicholson JK (2000) An NMR-based metabo- Siuzdak G (2006) XCMS: processing mass
nomic approach to investigate the biochemical spectrometry data for metabolite profiling
consequences of genetic strain differences: using nonlinear peak alignment, matching,
application to the C57BL10J and Alpk: and identification. Anal Chem 78(3):779–787
ApfCD mouse. FEBS Lett 484(3):169–174 15. Xia J, Mandal R, Sinelnikov IV, Broadhurst D,
5. Kurhanewicz J, Vigneron DB, Nelson SJ, Wishart DS (2012) MetaboAnalyst 2.0 – a
Hricak H, MacDonald JM, Konety B et al comprehensive server for metabolomic data
(1995) Citrate as an in vivo marker to discrimi- analysis. Nucleic Acids Res 40(W1):W127–
nate prostate cancer from benign prostatic WW33
hyperplasia and normal prostate 16. Lu W, Bennett BD, Rabinowitz JD (2008)
peripheral zone: detection via localized proton Analytical strategies for LC–MS-based targeted
spectroscopy. Urology 45(3):459–466 metabolomics. J Chromatogr B 871(2):
6. Lynch MJ, Nicholson JK (1997) Proton MRS 236–242
of human prostatic fluid: correlations between 17. Buescher JM, Moco S, Sauer U, Zamboni N
citrate, spermine, and myo-inositol levels and (2010) Ultrahigh performance liquid chroma-
changes with disease. Prostate 30(4):248–255 tography- tandem mass spectrometry method
7. Marx A, de Graaf AA, Wiechert W, Eggeling L, for fast and robust quantification of anionic and
Sahm H (1996) Determination of the fluxes in aromatic metabolites. Anal Chem 82(11):
the central metabolism of Corynebacterium 4403–4412
glutamicum by nuclear magnetic resonance 18. Michopoulos F, Whalley N, Theodoridis G,
spectroscopy combined with metabolite bal- Wilson ID, Dunkley TP, Critchlow SE (2014)
ancing. Biotechnol Bioeng 49(2):111–129 Targeted profiling of polar intracellular meta-
8. Gika HG, Theodoridis GA, Earll M, Snyder bolites using ion-pair-high performance liquid
RW, Sumner SJ, Wilson ID (2010) Does the chromatography and-ultra high performance
mass spectrometer define the marker? A com- liquid chromatography coupled to tandem
parison of global metabolite profiling data gen- mass spectrometry: applications to serum,
erated simultaneously via UPLC-MS on two urine and tissue extracts. J Chromatogr A
different mass spectrometers. Anal Chem 1349:60–68
82(19):8226–8234 19. Gika HG, Theodoridis GA, Vrhovsek U, Mat-
9. Dunn WB, Ellis DI (2005) Metabolomics: cur- tivi F (2012) Quantitative profiling of polar
rent analytical platforms and methodologies. primary metabolites using hydrophilic interac-
TrAC Trends Anal Chem 24(4):285–294 tion ultrahigh performance liquid
180 Filippos Michopoulos
Abstract
Metabolomics aims at identification and quantitation of key end point metabolites, basically polar, in order
to study changes in biochemical activities in response to pathophysiological stimuli or genetic modifica-
tions. Targeted profiling assays enjoying a growing popularity over the last years with LC-MS/MS as a
powerful tool for development of such (semi-)quantitative methods for a large number of metabolites.
Here we describe a method for absolute quantitation of ca. 100 metabolites belonging to key metabolite
classes such as sugars, amino acids, nucleotides, organic acids, and amines with a hydrophilic interaction
liquid chromatography (HILIC) system comprised with ultra (high) performance liquid chromatography
(UHPLC) with detection on a triple quadrupole mass spectrometer operating in both positive and negative
modes.
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_10,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
181
182 Christina Virgiliou et al.
2 Materials
Table 1
Concentration group and concentration (mg/L) of each group in standards 1–9
(continued)
HILIC-MS/MS Multi-targeted Method for Metabolomics Applications 185
Table 1
(continued)
(continued)
186 Christina Virgiliou et al.
Table 1
(continued)
(continued)
HILIC-MS/MS Multi-targeted Method for Metabolomics Applications 187
Table 1
(continued)
Mobile Phase B
Stock Buffer Solution Preparation
In a 50 mL beaker, weight 0.2254 g ammonium formate and add
water (<250 mL) (MilliQ). Sonicate until salt is fully dissolved.
Pour dissolved buffer into a 250 mL volumetric flash and add H2O
MilliQ till the calibration line.
Final ammonium formate buffer concentration: 14.2 mM.
3 Methods
3.1 UPLC Method With regard to UPLC conditions, the parameters for binary solvent
manager are as follows: flow rate is kept constant throughout the
whole analysis 0.5 mL/min and following gradient is programmed:
4 min isocratic step at 100% A, then rising to 40% B linearly over the
next 21 min, and finally reaching 85% B over 5 min.The column
was equilibrated for 10 min in the initial conditions. Regarding the
Sample Manager, flow through needle system is applied in the
present protocol. Injection volume is set at 5 μL and sample tem-
perature at 6 °C. Injection system was subjected to two washing
cycles with a strong solvent and a weak solvent prior to injection
and one cycle of 6 s with strong solvent for post wash. Column
temperature is set to 40 °C.
3.2 MS Method In order to edit the MS method, find the optimum parameters for
MRM transition of each metabolite. For that protocol, manual
optimization was performed in order to find precursor and product
ions and the optimum cone voltage and collision energy. With
regard to capillary voltage, the best possible for most metabolites
detected in positive and negative ionization mode was applied. All
the other parameters were set according to tune page and the linked
calibration file.
Apply multiple reaction monitoring (MRM) mode for the
detection and quantification of all the compounds. Operate elec-
trospray ionization at polarity switching mode. Set capillary voltage
at +3.5 kV or -3.5 kV, block, and desolvation temperatures at 150 °
C and 350 °C, respectively. Set desolvation gas flow rate at 650 L/
h and cone gas at 50 L/h. Optimum cone voltage and collision
energy for each analyte after direct infusion and optimum time
window and dwell times are presented in Table 2.
System start-up and pretests.
Table 2
Analytes and their MRM conditions and retention time
Cone
Monoisotopic Precursor Product Voltage Colission Rt Molecular Dwell
A/A Name Formula Mass Ion Ion (V) Energy (V) Polarity (min) Weight Time
1 2-Hydroxyisobutyric C4H8O3 104.04 103 57 30 10 - 8.0 104.10 0.005
2 2-Hydroxyisovaleric acid C5H10O3 118.06 117 71 30 12 - 6.0 118.13 0.005
3 2-Methylhippuric acid C10H11NO3 193.20 192 148 35 12 - 8.2 193.20 0.005
4 3-Methylhistidine C7H11N3O2 169.09 170 109 30 10 + 19.0 169.18 0.003
5 4-Hydroxyphenyllactate C9H10O4 182.05 181 63 33 12 - 10.0 182.17 0.02
6 a-Ketoglutaric acid C5H6O5 146.01 145 101 20 9 - 16.0 146.11 0.02
7 Acetylcarnitine C9H17NO4 203.12 204 85 30 10 + 14.4 203.23 0.003
8 Adenine C5H5N5 135.05 136 119 40 20 + 3.6 135.13 0.003
9 Adenosine C10H13N5O4 267.10 268 136 20 15 + 4.4 267.24 0.003
10 Adipic Acid C6H10O4 146.06 145 101 25 12 - 16.0 146.14 0.005
12 Alanine C3H7NO2 89.05 90 44 20 10 + 16.0 89.09 0.005
13 Arabitol C5H12O5 152.07 151 89 25 10 - 9.9 152.14 0.02
14 Arginine C6H14N4O2 174.11 175 70 30 19 + 21.9 174.20 0.005
15 Ascorbic acid C6H8O6 176.03 176 70 20 15 + 3.6 176.12 0.02
16 Asparagine C4H8N2O3 132.05 133 74 20 14 + 18.2 132.11 0.02
17 Aspartic acid C4H7NO4 133.04 134 74 18 16 + 21.8 133.11 0.02
18 Benzoic acid C7H6O2 122.04 121 77 25 11 - 1.8 122.12 0.032
HILIC-MS/MS Multi-targeted Method for Metabolomics Applications
(continued)
Table 2
190
(continued)
Cone
Monoisotopic Precursor Product Voltage Colission Rt Molecular Dwell
A/A Name Formula Mass Ion Ion (V) Energy (V) Polarity (min) Weight Time
21 Cadaverine C5H14N2 102.12 103 86 15 8 + 20.4 102.17 0.003
22 Caffeine C8H10N4O2 194.08 195 138 38 18 + 0.9 194.19 0.032
23 Choline C5H14NO 104.11 104 60 40 22 + 7.0 104.17 0.003
Christina Virgiliou et al.
(continued)
Table 2
192
(continued)
Cone
Monoisotopic Precursor Product Voltage Colission Rt Molecular Dwell
A/A Name Formula Mass Ion Ion (V) Energy (V) Polarity (min) Weight Time
67 Monoisoamylamine C5H13N 87.10 88 43 18 11 + 4.6 87.10 0.003
68 N-Acetylaspartate C6H9NO5 175.04 176 134 18 10 + 20.1 175.14 0.005
69 Nicotinamide C6H6N2O 122.05 123 96 40 15 + 1.1 122.12 0.032
Christina Virgiliou et al.
Follow each step here to ensure high quality of your data. For
the system start-up, the system is prepared for analysis and checked
for system suitability.
1. Install mobile phase A, B, wash, purge, and seal wash, and in
the two left channels, install neat acetonitrile and water.
2. Prime eluents for 3 min each.
3. Connect the BEH amide column.
4. Set the column temperature to 40 °C and the temperature of
sample manage to 6.
5. Flush the column with 95% MeCN and 5% water for 15 min at
flow of 0.2 mL/min.
6. Increase flow rate gradually until 0.5 mL/min.
7. Flush column until equilibration (psi delta < 20) (see Note 6).
8. Load UPLC method.
For MS, load the appropriate AcquityDB file with the extension
.ipr and the most recent calibration file.
3.3 Sample Extraction of samples may vary between different matrices. So far,
Preparation the described method or its variants were tested with serum, urine,
amniotic fluid, intra/extra-cellular content, feces, and various types
of animal tissue but also in foods such as honey, muscle tissue, and
flour. The sample preparation procedure for blood serum samples is
presented as follows:
1. Allow samples to thaw in room temperature.
2. Mix 50 μL of sample with 130 μL MeCN, 10 μL H2O, and
10 μL MeOH in an Eppendorf vial of 1 mL by the use of
variable volume pipette 20–200 μL (see Note 7)
3. Vortex for 10 min and centrifuge at 7000rpm for 10 min.
4. Transfer supernatant to an LC-MS with glass insert vial.
5. Immediately transfer to precooled Sample Manager.
For standard addition approach, sample preparation of spiking
samples is as follows (see Note 8):
1. Mix 50 μL of sample with 100 μL of standard calibration
mixture 1 and 50 μL of MeCN in an Eppendorf vial of 1 mL
by the use of variable volume pipette 20–200 μL.
2. Vortex for 10 min and centrifuge at 7000rpm for 10 min.
3. Transfer supernatant to an LC-MS with glass insert vial.
4. Repeat steps 1–3 using increased concentration standard cali-
bration mixtures (at least 3 in total).
5. Immediately transfer to precooled Sample Manager.
HILIC-MS/MS Multi-targeted Method for Metabolomics Applications 195
3.4 UPLC-MS/MS The following steps describe how to setup a sample table for data
Analysis acquisition. Following a sample order has been evaluated to be
optimal for both throughput and quality controls:
1. Run a gradient without injection in order to evaluate column
performance and have a measure of the minimum and maxi-
mum pressure during analysis (see Note 11).
2. Perform six replicate injections of a standard mixture in order
to pretest systems performance (retention time, signal) and to
equilibrate system (see Note 12).
During analysis with external calibration approach
1. Run a before-batch calibration curve.
2. Run injections of a QC sample in order to perform matrix
equilibration of the system. The number of equilibration injec-
tions depends on the analyzed matrix. For cell media and urine,
five injections may be adequate; for blood and tissue samples,
higher numbers may be necessary, depending also on the use of
the column (new columns need more injections to saturate
active sites/equilibrate).
3. Samples are injected in random order in blocks of ten samples.
4. After each block of ten samples, injection of a standard mixture
and a QC sample is performed.
5. Repeat steps 2 and 3 for a maximum of 100 samples (see Note
13).
6. Run an after-batch calibration curve.
During analysis with standard addition approach
1. Run injections of a QC sample in order to perform matrix
equilibration of the system.
2. Samples followed by their paired spiking samples are injected in
order in block of ten samples.
3. After each block of ten samples, injection of a standard mixture
and a QC sample is performed.
4. Repeat steps 2 and 3 for a maximum of 100 samples.
After finishing the sample, batch analysis column and the MS
have to be cleaned for further use.
196 Christina Virgiliou et al.
1. Flush the column with 50% MeCN and 50% H2O for 50 min,
column temperature 50 °C, flow 0.2 mL/min.
2. Flush the column with 50% MeCN and 50% H2O for 30 min,
column temperature 40 °C, flow rate 0.5 mL/min.
3. Clean MS cone as recommended.
3.5 Data Treatment: Quantitation can be performed in both manually and automated
Quantification of ways using vendor software. In the present protocol, TargetLynx
Metabolites (Waters) was used. A Quantify Method must be created before
integration or quantification can be performed. Related software
from other vendors are MutliQuant (Sciex), Xcalibur, and Mas-
sHunter WorkStation - Quantitative Analysis.
TargetLynx data can be saved as .qld files for further manipula-
tion. Complete the summary of the results and .qld files (area,
response, concentration, S/N, SD, measured concentration etc.)
can be exported as .txt files and open with excel Microsoft program
for further treatment.
In case of external calibration curve approach, once you find
the optimum method parameters, calibration of standards and
quantification can be performed directly with only one process
(see Note 14).
When standard addition method is applied, automated quanti-
tation is not available. Integrated results for each sample and
spiked/fortified samples together with corresponding spiking
levels must be imported in a spreadsheet program (excel or similar)
where unknown concentrations of metabolites in samples can be
calculated as intercept/slope.
The described method was adapted and optimized for the absolute
quantitation of 20 AAs (alanine, arginine, asparagine, aspartic acid,
cysteine, glutamic acid, glutamine, glycine, histidine, isoleucine,
leucine, lysine, methionine, phenylalanine, proline, serine, threo-
nine, tryptophan, tyrosine, valine) and 9 AA derivatives (betaine,
choline, hypotaurine, pyroglutamic acid, taurine,
3-methylhistidine, 4-hydroxyproline, sarcosine, creatine).
Preparation of calibration standards includes the division of
AAs into five separate groups according to their concentration in
urine, namely, A, B, C, D and E. Standard working solutions were
freshly prepared before the analysis by dilution with ACN: H2O
(95:5, v/v) at nine concentration levels. Spike the calibrators in
50 μL of synthetic urine (see Note 15) to get a final concentration
range of 0.01–2 μg/mL for group A, 0.02–4 μg/mL for group B,
0.09–18 μg/mL for group C, 0.15–30 μg/mL for group D,
0.425–95 μg/mL for group E, and 1–200 μg/mL for group F
(Table 3). Add, as IS, in all calibrators or QCs an equal volume
HILIC-MS/MS Multi-targeted Method for Metabolomics Applications 197
Table 3
Concentration group and concentration (mg/L) of each group in standards
1–9
A B C D E F
std1 0.02 0.04 0.18 0.3 0.95 2
std2 0.04 0.08 0.36 0.6 1.9 4
std3 0.16 0.32 1.44 2.4 7.6 16
std4 0.4 0.8 3.6 6 19 40
std5 1 2 9 15 47.5 100
std6 2 4 18 30 95 200
std7 2.4 4.8 21.6 36 114 240
std8 3 6 27 45 142.5 300
std9 4 8 36 60 190 400
Metabolites Group concentration
Alanine C
Arginine A
Asparagine D
Aspartate C
Choline B
Creatine E
Cysteine A
Glutamate C
Glutamine F
Glycine F
Histidine E
Hypotaurine A
Isoleucine A
Leucine B
Lysine F
Methionine A
Betaine B
Phenylalanine C
Proline D
Pyroglutamic acid F
Sarcosine A
(continued)
198 Christina Virgiliou et al.
Table 3
(continued)
A B C D E F
Serine D
Taurine F
Threonine E
Tryptophan C
Tyrosine D
Valine C
3-methylhistisine E
4-hydroxyproline A
Table 4
MRM conditions and retention time of amino acids and amino acid derivatives
Precursor
Name ion Product ion Cone voltage (V) Collision energy (V) Rt (min)
Glycine 76 30 35 6 3.8
Glycine 13C 15N 79 32 35 6 3.9
Sarcosine 90 44 20 8 3.1
Alanine 90 44 20 10 3.4
Alanine 13C 15N 94 47 20 10 3.4
Choline 104 59 10 10 1.5
Serine 106 60 20 10 4.3
Hypotaurine 110 92 22 8 3.4
Serine 13C 15N 110 92 20 10 3.3
Proline 116 70 20 25 2.7
Betaine 118 58 10 10 2.2
Valine 118 72 18 9 2.5
Threonine 120 74 20 10 3.6
Proline 13C 15N 122 76 20 25 2.6
Cysteine 122 76 18 12 3.7
Valine 13C 15N 124 77 18 9 2.5
Threonine 13C 15N 125 78 20 10 3.7
Taurine 126 108 25 10 2.7
Pyroglutamic 130 84 30 12 2.3
4OH Proline 132 68 30 25 3.4
Isoleucine 132 69 25 12 2.1
leucine 132 86 20 10 1.9
Creatine 132 90 28 10 3.4
Asparagine 133 74 20 14 4.4
Aspartic acid 134 74 18 16 5.9
Isoleucine 13C 15N 139 92 25 12 2.1
Leucine 13C 15N 139 92 20 10 1.9
Asparagine 13C 15N 139 77 20 14 4.5
Aspartic acid 13C 15N 139 77 18 16 5.9
Glutamine 147 84 20 15 4.2
Lysine 147 84 14 14 6.5
(continued)
200 Christina Virgiliou et al.
Table 4
(continued)
Precursor
Name ion Product ion Cone voltage (V) Collision energy (V) Rt (min)
Glutamic acid 148 130 20 13 4.8
Methionine 150 104 22 9 2.3
Glutamine 13C 15N 154 89 20 15 4.2
Glutamic acid 13C 15N 154 136 20 13 4.8
Lysine 13C 15N 155 90 14 14 6.5
Histidine 156 110 22 13 6.2
Methionine 13C 15N 156 109 22 9 2.4
Histidine 13C 15N 165 118 22 13 6.6
Phenylalanine 166 120 22 12 1.9
3-Methylhistidine 170 109 30 10 6.2
Arginine 175 70 30 19 6.2
Phenylalanine 13C 15N 176 129 22 12 1.9
Tyrosine 182 136 22 13 2.6
Arginine 13C 15N 185 75 30 19 6.1
Tyrosine 13C 15N 192 145 22 13 2.6
Tryptophan 205 146 20 15 1.9
Tryptophan 13C 15N 205 188 20 15 1.8
5 Notes
Fig. 1 Calibration curves comparing parallelism for 3-methylhistidine (a: area on y-axis, b response factor
area 3-methylhistidine/area 13C15N histidine) and pyroglutamic acid (c: area on y-axis, d response factor
area pyroglutamic acid/area 13C15N glutamic acid)
References
1. Nicholson JK, Lindon JC, Holmes E (1999) 6. Cajka T, Fiehn O (2016) Toward merging
“Metabonomics”: understanding the meta- untargeted and targeted methods in mass
bolic responses of living systems to pathophys- spectrometry-based metabolomics and lipido-
iological stimuli via multivariate statistical mics. Anal Chem 88:524–545. https://doi.
analysis of biological NMR spectroscopic data. org/10.1021/acs.analchem.5b04491
Xenobiotica 29:1181–1189. https://doi.org/ 7. Zhou B, Xiao JF, Tuli L, Ressom HW (2012)
10.1080/004982599238047 LC-MS-based metabolomics. Mol BioSyst 8:
2. Nicholson JK, Connelly J, Lindon JC, Holmes 4 7 0 – 4 8 1 . h t t p s : // d o i . o r g / 1 0 . 1 0 3 9 /
E (2002) Metabonomics: a platform for study- c1mb05350g
ing drug toxicity and gene function. Nat Rev 8. Serafimov K, L€ammerhofer M (2022) Meta-
Drug Discov 1:153–161. https://doi.org/10. bolic profiling workflow for cell extracts by
1038/nrd728 targeted hydrophilic interaction liquid
3. Nicholson JK, Lindon JC (2008) Systems biol- chromatography-tandem mass spectrometry. J
ogy: metabonomics. Nature 455:1054–1056. Chromatogr A 1684:463556. https://doi.
https://doi.org/10.1038/4551054a org/10.1016/j.chroma.2022.463556
4. Theodoridis G, Gika HG, Wilson ID (2011) 9. Verpoorte R, Choi YH, Kim HK (2010) Meta-
Mass spectrometry-based holistic analytical bolomics: will it stay? Phytochem Anal 21:2–3.
approaches for metabolite profiling in systems https://doi.org/10.1002/pca.1191
biology studies. Mass Spectrom Rev 30:884– 10. Schrimpe-Rutledge AC, Codreanu SG, Sher-
906. https://doi.org/10.1002/mas.20306 rod SD, McLean JA (2016) Untargeted meta-
5. Gika HG, Theodoridis GA, Plumb RS, Wilson bolomics strategies – challenges and emerging
ID (2014) Current practice of liquid directions. J Am Soc Mass Spectrom 27:1897–
chromatography–mass spectrometry in meta- 1905. https://doi.org/10.1007/s13361-
bolomics and metabonomics. J Pharm Biomed 016-1469-y
Anal 87:12–25. https://doi.org/10.1016/j. 11. Gertsman I, Barshop BA (2018) Promises and
jpba.2013.06.032 pitfalls of untargeted metabolomics. J Inherit
HILIC-MS/MS Multi-targeted Method for Metabolomics Applications 203
Abstract
A simple analytical workflow is described for gas chromatographic-mass spectrometric (GC-MS)-based
chiral profiling of secondary amino acids (AAs) in biological matrices. The sample preparation is carried out
directly in aqueous biological sample extracts and involves in situ heptafluorobutyl chloroformate (HFBCF)
derivatization—liquid-liquid microextraction of nonpolar products into hexane phase followed by
subsequent formation of the corresponding methylamides from the HFB esters by direct treatment with
methylamine reagent solution. The (O, N) HFB-butoxycarbonyl–methylamide AA products (HFBOC-
MA) are separated on a Chirasil-L-Val capillary column and quantitatively measured by GC-MS operated in
selected ion monitoring (SIM) mode. The protocol includes 12 simple pipetting steps and covers the
quantitative analysis of 8 L, D pairs of secondary amino acids, including proline, isomeric 3-, 4-hydro-
xyprolines, pipecolic acid, nipecotic acid, azetidine-2-carboxylic acid, and cis- and trans-5-hydroxy-L-
pipecolic acid using 13C5 –L-proline as an internal standard. The individual analytical steps are commented
on and explained, with emphasis on the chiral GC-MS analysis of secondary amino acids in human urine,
serum, and peptide hydrolysate samples.
Key words Secondary amino acids, Amino acids, GC-MS, Alkyl chloroformate derivatization, Chiral
analysis, Urine, Serum, Peptide hydrolysates, Quantitative analysis
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_11,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
205
206 Stanislav Opekar et al.
Fig. 2 The analytical workflow (from Ref. [1], Springer Nature with permission)
Fig. 3 GC-SIM-MS separation of D and L isomers obtained for the target secondary amino acids; (a) if 4c-Hyp
is not present, (b) if 4 t-Hyp is not present, (c) resolution between 5 t-L-HyPip and 5c-L-HyPip (from Ref. [1],
Springer Nature with permission)
Table 1
The content of secondary amino acids in human serum and urine (n = 5)
2 Materials
2.1 Samples Samples containing no or little protein and cell residues (urine, cell
culture media) or cell and tissue extracts can be analyzed by this
method. Moreover, this method can be applied to precipitated
serum/plasma or peptide hydrolysates.
2.2 Chemicals, 1. The analytical grade chemicals should be used. All solutions,
Solutions, except that containing the HFBCF reagent, should be prepared
Reaction Media in distilled deionized water (DI water, < 1.5 μS.cm1, 25 C)
and stored them at 4 C (unless otherwise indicated).
2. 100 mM NaHCO3 (99,998% purity, Alfa Aesar) solution: Dis-
solve 840 mg in 100 mL of DI water.
3. 1 M NaOH (99.99% purity, Alfa Aesar) solution: Dissolve 4 g
in 100 μL of DI water.
4. Heptafluorobutanol (HFBOH) (99% purity, Chromservis, Pra-
gue, Czech Republic).
5. Heptafluorobutyl chloroformate (HFBCF, 98%, Chromser-
vis), (see Note 2).
6. Solvents: Isooctane (2,2,4-trimethylpentane, 99.5%), pyridine
(p.a., 99.0%), isopropanol (2-propanol, 99.5% purity) (all
Sigma-Aldrich).
7. The organic reaction medium: prepare isooctane, HFBCF, and
HFBOH in a volume ratio 15:4:1 (v/v/v) in a Teflon-capped,
well-tightened 4-mL glass vial. Store in a refrigerator, where
the mixture remains stable for several months.
A Protocol for GC-MS Profiling of Chiral Secondary Amino Acids 211
2.5 Additional 1. The sample preparation glass culture tubes 6 50 mM, mate-
Equipment rial: sodium-potassium silicate.
2. glass (Merci, P/N Z1632000605010), or borosilicate glass
(Kimble-Kontes, P/N 73500-650).
3. Common screw cap Teflon-lined 2 and 4 mL amber vials for
the reagent solutions.
4. An adjustable 50 and 100 μL Transferpettor pipette with a glass
capillary (Brand, P/N 701868 and 701873) for manipulation
with the reagents and their mixtures in isooctane. The pipette
tips with 25 mM capillary (gel-loading type, VWR Int.) for
aspirating the upper organic phase in the 6 50 mM vial.
Table 2
212
HFBOC-MA
HFBOC-HFB derivatives
derivatives (SPS1) (SPS1 + SPS2) Metabolite Database Coding
RT Diagnostic RT Diagnostic
No. Traditional Name IUPAC Name Abbreviation [min] (m/z) [min] (m/z) CAS KEGG HMDB
1 L-Azetidine-2- (2S)-Azetidine-2-carboxylic L-Aze 5.65 282 5.88 282 2133-34-8 C08267 HMDB29615
Stanislav Opekar et al.
3 Methods
3.1 Sampling and 1. A serious attention should be paid to the sample collection,
Storage transport, and storage, because any omission may result in false
results. For urine, collection of the morning second-void sam-
ples is the most common practice and the easiest sampling
method.
2. For urine analysis, store freshly collected samples at 4 C within
2 h. For longer than 48 h storage, freeze the samples and keep
at 20 C (see Note 3).
3.4 GC-MS Analysis 1. The instrument GC–MS conditions are summarized in Table 3.
2. First, analyze the standard mixtures to check the separation
performance, retention times, the analyte peak shape, and
acquisition of the employed fragment ions in the obtained EI
spectra, Table 2.
3. Using a single quadrupole MS analyzer, selected-ion monitor-
ing (SIM) mode is commonly used for quantification. Use the
characteristic m/z ions listed in Table 2.
4. Prepare an appropriate analysis sequence of a sample series
consisting of repeated blank, standard, the QCs (see Subhead-
ing 3.2, step 3), calibration, and real samples. Measure regu-
larly the blanks, QC samples, and standards (at least every
10 sample runs). Analyze samples in a random order to avoid
a systematic error.
5. Change the GC injector liner approximately after 150 samples
depending on the sample matrix. Before use, condition each
new liner by running the following sequence: solvent blank,
standard mixtures, the pooled QC sample extract (twice), and
solvent blank (twice).
216 Stanislav Opekar et al.
Table 3
GC-EI-SIM-MS operating conditions
GC Injector Mode Pulsed splitless (elevated head pressure from 110 to 220 kPa)
Liner Sky® 4 mM I.D. cyclo double taper inlet liner (Restek, P/N 23310)
Temperature 200 C
Injection 0.5 μL
volume
Temperature 200 C
GC Oven_1a Initial 89 C
temperature
Ramp Hold for 1.5 min, 2 C/min to 150 C, hold for 6 min then back
20 C/min to 89 C
Run time 42 min
GC Oven_2 b
Initial 110 C
temperature
Ramp Hold for 1 min, 5 C/min to 180 C, 1 C/min to 200 C hold for
8 min then back 20 C/min to 110 C
Run time 47.5 min
GC column Capillary CP-Chirasil-L-Val (25 m 0.25 mM, 0.12 μm, Agilent
column Technologies, USA, P/N CP7495)
Carrier gas Helium
Flow rate 1.2 mL/min
Mode Constant flow
Outlet pressure Vacuum
GC-MS transfer Temperature 200 C
line
MS quadrupole Ion source 200 C
temperature
Full scan mode m/z 50–900 Da
SIM mode Metabolite SIM m/z ions, see Table 1
Software MSD ChemStation (version E.02, Agilent)
a
HFBOC-HFB analytes after SPS1
b
HFBOC-MA analytes after SPS 1 + 2
3.5 Data Analysis 1. Peak areas for corresponding ions are integrated. Their ratio is
calculated to test for potential interferences.
2. The peaks of particular metabolites are normalized by the peak
area of the internal standard 13C5 proline.
A Protocol for GC-MS Profiling of Chiral Secondary Amino Acids 217
4 Notes
Acknowledgments
References
1. Opekar S, Zahradnickova H, Vodrazka P, chloroformate & methylamine derivatization.
Rimnacova L, Simek P, Moos M (2021) A Amino Acids 53(3):347–358
chiral GC-MS method for analysis of secondary 2. Langrock T, Garcia-Villar N, Hoffmann R
amino acids after heptafluorobutyl (2007) Analysis of hydroxyproline isomers
A Protocol for GC-MS Profiling of Chiral Secondary Amino Acids 219
Abstract
Trapped ion mobility spectrometry (TIMS) using parallel accumulation serial fragmentation (PASEF®) is
an advanced analytical technique that offers several advantages in mass spectrometry (MS)-based lipido-
mics. TIMS provides an additional dimension of separation to mass spectrometry and accurate collision
cross-section (CCS) measurements for ions, aiding in the structural characterization of molecules. This is
especially valuable in lipidomics for identifying and distinguishing isomeric or structurally similar com-
pounds. On the other hand, PASEF technology allows for fast and efficient data acquisition by accumulat-
ing ions in parallel and then serially fragmenting them. This accelerates the analysis process and improves
throughput, making it suitable for high-throughput applications. Moreover, the combination of TIMS and
PASEF reduces co-elution and ion coalescence issues, leading to cleaner and more interpretable mass
spectra. This results in higher data quality and more confident identifications. In this chapter, a data-
dependent TIMS-PASEF® workflow for lipidomics analysis is presented.
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_12,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
221
222 Vanna Denti et al.
1.2 TIMS-PASEF® in IM-MS offers a cutting-edge analytical platform that can separate
Lipidomics ions beyond the mass-to-charge ratio (m/z). This technique per-
mits to resolve isobaric and isomeric compounds using the mobility
of an ion, which is directly related to its size, shape, and charge as
well as to the properties of the buffer gas [7]. In this chapter, we
will focus on trapped ion mobility spectrometry (TIMS) technol-
ogy that was first introduced in 2011 by Fernandez Lima et al.
TIMS reverses the concept of the classical drift tube IMS [8].
In fact, TIMS holds ions in a stationary position against the gas
flow and then releases them depending on their mobility. In partic-
ular, it allows work to be performed in the parallel accumulation
serial fragmentation (PASEF®) operation mode in combination
with conventional LC-MS [9, 10]. In the dual TIMS design, ions
can be accumulated and analyzed simultaneously because the ion
accumulation and ion analysis functions are spatially separated.
The ions are first gathered in the trap; then they are swiftly
delivered to the analyzer for the mobility examination. To prevent
the loss of the ions, the trap is simultaneously loaded with the
following batch of ions.
TIMS-PASEF®, working in data-dependent (DDA) mode, can
extend the dynamic range of mass spectrometry, enabling the
detection of both highly abundant and low-abundance ions in the
same analysis. Therefore, TIMS-DDA-PASEF analysis can be used
to obtain level 2 lipid annotation [11] (see Notes 1 and 2).
In this modality, TIMS and MS/MS precursor selection are
synchronized. This system selects numerous precursors for MS/MS
acquisition in a single TIMS scan by using the temporal separation
of ions from the TIMS instrument.
The precursor ions in TIMS-QTOF equipment pass through
the analytical quadrupole serially rather than in parallel, as they do
in traditional MS/MS. When in PASEF® mode, an ion group is
collected in the upstream trap while another ion group that has
already been accumulated is examined simultaneously in the
UHPLC-TIMS-PASEF®-MS for Lipidomics: From Theory to Practice 223
Fig. 1 Online Parallel Accumulation—Serial Fragmentation (PASEF®) with the timsTOF. (a) Analytes eluting
from the chromatographic column are ionized and enter the mass spectrometer through a glass capillary. (b)
In the dual TIMS analyzer, the first TIMS section traps and stores ion packets, and the second resolves them by
mobility. (c. d) Ion mobility-separated ions are released sequentially from the second TIMS analyzer as a
function of decreasing electrical field strength and yield mobility-resolved mass spectra. (e) In PASEF® MS/MS
scans, the TIMS analyzer and the quadrupole are synchronized, and the quadrupole isolation window switches
within sub-milliseconds between mobility-resolved precursor ions of different m/z. (f) This yields multiple ion
mobility-resolved MS/MS spectra from a single TIMS scan and ensures that multiple trapped precursor ion
species are used for fragmentation. (g) Cleaner MS/MS spectra obtained by the fragmentation of TIMS-
resolved ions
2 Materials
3 Methods
3.1 TIMS-MS All the following steps require the timsControl software.
Settings
1. Turn the instrument into the “Operate state.”
2. Select the appropriate acquisition method.
3. Wait for the instrument set-up (~ 20 min).
In this chapter, the steps for a DDA-PASEF® acquisition are
described. Therefore, the parameters to setup involve both the MS
and IM dimensions. To ensure accurate measurements both in the
MS and the IM dimensions, the instrument must be calibrated
UHPLC-TIMS-PASEF®-MS for Lipidomics: From Theory to Practice 225
3.1.1 Data-Dependent For the purpose of this chapter, a DDA-PASEF® was performed
Acquisition (DDA) Coupled both in positive and negative ion mode. The main parameters of the
with TIMS PASEF® methods applied are reported in Tables 1, 2, 3 and 4. The para-
meters selected for the method must be used for calibration. If the
IM or the m/z range changes after the calibration, the timsControl
software will highlight with a red color the dimension that needs to
be calibrated. For the protocol presented in this chapter, the ana-
lyses were performed with 100–1350 m/z mass range and
0.55–1.90 V*s/cm2 mobility range (see Note 3).
Ensure that the MS method consists of three segments:
1. 0–0.05 min—The instrument operates in full scan and TIMS
on mode. The flow from the UHPLC is infused through a
divert valve to the ESI source.
2. 0.05–0.3 min—The instrument operates in full scan and TIMS
on mode. The UHPLC flow goes to waste and a calibrant
solution (see next paragraph Subheading 4.) is infused through
the divert valve.
3. 0.3–21 min—The instrument operates in PASEF® mode. The
UHPLC flow enters in the ESI source.
Table 1
TIMS setting and source parameters for lipidomics analysis in ESI(+) and (2) for a timsTOF fleX
instrument
TIMS settings
Polarity Mode Ramp time Accumulation time Duty cycle Ramp rate
ESI () Custom 100 ms 100 ms 100% 9.42 HZ
ESI (+) Custom 100 ms 100 ms 100% 9.42 HZ
Principal source parameters
Polarity Transfer
End plate offset Capillary Nebulizer Dry gas Dry temp
ESI () 500 V 3600 V 2.2 Bar 10 l/min 220 C
ESI (+) 500 V 4500 V 2.2 Bar 10 l/min 220 C
226 Vanna Denti et al.
Table 2
Tune parameters for lipidomics analysis in ESI(+) and (2) for a timsTOF fleX instrument
MS tune parameters
Transfer
Polarity Deflection 1 delta Funnel 1RF isCID energy Funnel 2RF Multiple RF
ESI () 80 V 500Vpp 0 eV 250 Vpp 200Vpp
ESI (+) +80 V 500Vpp 0 eV 250 Vpp 200Vpp
Polarity Quadrupole Focus pre TOF Collision cell
Ion energy Transfer time Pre-pulse storage Collision energy Collision RF
ESI () 5 eV 65 μs 5 μs 10 eV 1100 Vpp
ESI (+) 5 eV 65 μs 5 μs 10 eV 1100 Vpp
TIMS tune parameters
Offset
Polarity Deflection transfer Funnel 1In Collision cell in
To capillary exit To deflection discard To deflection transfer
ESI () 20 V 120 V 80 220 V
ESI (+) 20 V 120 V 80 220 V
Polarity Accumulation trap Accumulation exit Ramp start
Funnel 1 in To accumulation transfer To accumulation exit
ESI () 100 V 0 100 V
ESI (+) 100 V 0 100 V
Polarity Ion charge control Advanced parameters
Enable Target intensity Lock accumulation to mobility range
ESI () X 7.50 M X
ESI (+) X 7.50 M X
Make sure to have two Hamilton glass syringes of 500 μL, one
for the first solution and one for the other two.
2. The pre-acquisition calibration is a two-step procedure:
– External m/z calibration is performed by direct infusion
using a divert valve at 180 μL/h of a sodium formate
solution using the high-order polynomial calibration
(HPC) calibration mode.
External mobility calibration is performed by direct
infusion using a divert valve at 180 μL/h of tuning mix
MMI-L using the linear calibration mode.
UHPLC-TIMS-PASEF®-MS for Lipidomics: From Theory to Practice 227
Table 3
PASEF® parameters for lipidomics analysis in ESI(+) and (2) for a timsTOF fleX instrument
PASEF® parameters
Precursor ions
Collision energy settings Active exclusion
Polarity # PASEF ramps Total cycle time Collision energy Release after
ESI () 2 0.32 s 40 eV 0.10 min
ESI (+) 2 0.32 s 40 eV 0.10 min
Polarity Precursor repetition Isolation width
Target intensity Intensity threshold
ESI () 4000 100 2 m/z
ESI (+) 4000 100 2 m/z
Table 4
MetaboScape® import parameters
T-ReX 4D workflow
3.2 UHPLC For the work presented in this chapter, the extracted lipidome was
Preparation separated by LC prior to the TIMS-PASEF® measurement using an
elute UHPLC system. RP chromatographic separation was
achieved using an ACQUITY UPLC CSH C18 Column (Waters).
1. Prepare the following solutions:
– Mobile phase A: 10 mM ammonium acetate:acetonitrile
(40:60 v/v) with 0.1% acetic acid.
– Mobile phase B: Isopropanol:phase A (90:10 v/v).
– Seal wash: Water:isopropanol (80:20 v/v).
2. Mount the column in the column compartment.
3. Install the solutions in the corresponding lines.
4. Insert the wash lines as follows:
(a) Wash 1: in Mobile phase A.
(b) Wash 2: in Mobile phase B.
5. Open the Hystar 6.0 software to control the UHPLC system:
– Prime the seal wash.
– Purge solvent lines A and B.
– Perform the initial wash of the autosampler (Wash
2 + Wash 1).
6. Set the parameters of the chromatography:
– Set the column temperature at 55 C and turn on the
column oven.
– Set the flow rate to 0.250 mL/min.
– Set the injection mode to “μL-pickup.”
– Set the autosampler temperature to 10 C.
– Set the elution gradient timetable:
– 0 min 99% A, 1 min 99% A, 1.10 min 60% A, 5 min 20% A,
11 min 20% A, 12 min 1% A, 18 min 1% A, 18.10 min
60% A, 20 min 99%. Add 1 min of equilibration time after
each run.
7. Turn on the flow and start at least three runs with no injection
to perform column conditioning.
8. Save the method.
The injection volume will be defined when setting the sample
list, as described in the next paragraph.
3.2.1 Sample List Using Hystar 6.0 software, it is possible to set up the sample list
Settings table.
1. Indicate the name and position of the sample.
2. Set the number of injections for each sample.
UHPLC-TIMS-PASEF®-MS for Lipidomics: From Theory to Practice 229
3.2.2 Sample Analysis 1. Perform the lipid extraction. There are numerous techniques
and complete procedures available in the literature (see
Note 5).
2. Prepare a pooled sample using aliquots of the lipid extracts
from each sample. This will act as a quality control (QC)
sample.
3. Your samples should be put into glass vials with screw caps and
PFTE septa. In order to use less sample volume for the analysis,
glass adapters can be used.
4. Place the vials in the autosampler and proceed with the sample
list setting, as described in the previous paragraph.
5. Perform at least four equilibration runs using the QC sample at
the start of each batch of analyses.
6. The order of sample injection should be randomized, and QC
samples should be analyzed during the run. Use these to test
system suitability.
7. Start the data acquisition.
3.2.3 Data Extraction and Once the lipidomic analysis is complete, it is possible to have a quick
Data Processing check on the quality of the data acquisition using DataAnalysis
software (see Note 6).
In Figs. 2 and 3 are illustrated some examples of the output
visible with this software. The base peak chromatogram and the
total ion count of the MS2 sample analyzed in positive and negative
ion modes are shown. Lipidomics analysis obtained with
LC-TIMS-PASEF® mode can be represented as four-dimensional-
data: time, mobility, m/z, and MS/MS.
Checking these results, you can have a qualitative overview of
the analysis. This step is not mandatory but can be a first screening.
To proceed with the lipidomic data analysis, it is necessary to obtain
a list of identified lipids and their intensities in each sample. To do
so, it is possible to use different software: MetaboScape®, (Bruker
Daltonik GmbH, Germany), or Mzmine3 (http://mzmine.github.
io/).
230 Vanna Denti et al.
Fig. 2 The base peak chromatogram of a lipid extract from plasma analyzed as described in the previous
paragraphs, in ESI(+) a and ESI() c. The corresponding total ion chromatogram (TIC) of the MSn ions is
reported in b and d
Data Import
For 4D lipidomic data, we need to define import parameters for all
the data dimensions. In general, it is necessary to take care of the
following steps:
– Set the intensity threshold and the data point expected for
each peak.
– Set RT and m/z ranges. Mobility ranges will be automatically
detected.
– Set the import method for the MS/MS.
– Perform ion deconvolution, setting the type of ion adducts
expected.
– Perform data recalibration considering the calibration segment
at 0–0.3 min.
UHPLC-TIMS-PASEF®-MS for Lipidomics: From Theory to Practice 231
Fig. 3 An example of high quality MS/MS spectra obtained with TIMS-PASEF® analysis in ESI(+) and ESI()
mode (a and b, respectively). The ion ID, m/z of the precursor, RT, and mobility are reported. The MS/MS
fragments can be used to confirm the ion identity of a PC(15:0/18:1)(d7)
232 Vanna Denti et al.
Lipids Annotation
The 4D feature list imported after the previous step can be anno-
tated matching the experimental data with experimental lipid
libraries or with computational algorithms (i.e., Lipid Species
Annotation, LipiDex Library Forge [14]).
MetaboScape® enables the annotation of lipid species consid-
ering accurate m/z ratio, isotopic pattern, and CCS (if available) of
the precursor, as well as MS/MS spectra. An annotation is assigned
if the score for each of these parameters is higher than a defined
threshold. A final total score for the match is assigned for each
identification. If necessary, lipid identification can be curated man-
ually, examining the MS/MS spectra of the feature.
One of the advantages of TIMS-PASEF® analysis is the possi-
bility to obtain clean spectra, as described in the introduction. In
Fig. 3 is possible to observe the MS/MS spectra (A. ESI+, B. ESI)
of the feature annotated as a PC(15:0/18:1)[d7]. The feature has
the same RT and mobility, and the fragments are those specific for
this lipid fragmentation, both in positive and negative mode.
Additionally, an overall view of the lipids identified in an ana-
lytical batch can be offered by plotting a Kendrick Mass-Defect
Plot, and it can be used to further validate the level 2 annotations
assigned by Lipid Species Annotation. In the example in Fig. 4, it is
possible to see the KMD for CH2 defect plotted against the m/z
obtained from the positive and negative (A and B, respectively).
Another advantage of using ion mobility for lipidomics analysis
is that it can potentially discriminate isomeric molecules. In Fig. 5,
it is reported the example of two isomers of a cholesterol ester
(CE) 18:3 putatively annotated with Lipid Species Annotation. In
panel A is possible to see the extracted ion chromatogram (EIC) of
CE(18:3) + NH4+ (664.6031 m/z 0.005 Da) in green. The EICs
of the two isomers (R1 at 1/K0 ¼ 1.394–1.421 and R2
1/K0 ¼ 1.424–1.454. The regions isolated from the heatmap in
panel B allowed to extract the ion mobilograms (EIMs) of R1 and
R2 (panel C), respectively. In panel D, it is possible to observe that
the two species share the same m/z precursor (indicated as a blue
square) and some characteristic peaks (i.e. 369.3516 m/z: -FA 18:3
(+HO) -Cholesterol), despite some differences in the lower mass
range (i.e., 311.2732 m/z and 313.2890 m/z).
Fig. 4 Kendrick Mass-Defect (KMD) Plot showing the level 2 annotated lipids of the two datasets obtained as
described above. The figure above (a) reports the data of the ESI(+) analysis, whereas b reports the data of the
ESI() analysis. As shown in the legend on the right, the colors correspond to the lipid classes identified, while
the radius of the dots is proportional to their CCS. CCS and RT outlier symbols are also reported
234 Vanna Denti et al.
Fig. 5 Example of two isomers of a compound identified as a CE(18:3) by Lipids Species Annotation algorithm.
In panel a, the extracted ion chromatogram (EIC) of the [M + NH4]+ 664.6031 m/z: in green, the overall EIC; in
blue and red, the EICs of the R1 and R2 isomers, respectively. In panel b are indicated the two isomers ions
isolated in the regions of interest R1 and R2. The extracted ions mobilograms (EIMs) of the isolated regions are
reported in panel c. In panel d, the MS/MS spectra of the two isomers of the CE(18:3), in blue the R1 and in red
R2 (see Note 3)
4 Notes
Acknowledgments
References
1. Züllig T, Trötzmüller M, Köfeler HC (2020) 5. Ridgway N, McLeod R (2021) Biochemistry of
Lipidomics from sample preparation to data lipids, lipoproteins and membranes. Elsevier
analysis: a primer. Anal Bioanal Chem 412: 6. Wojcik R, Webb IK, Deng L, Garimella SV,
2191–2209 Prost SA, Ibrahim YM, Baker ES, Smith RD
2. Zhang Z, Singh M, Kindt A, Wegrzyn AB, (2017) Lipid and glycolipid isomer analyses
Pearson MJ, Ali A, Harms AC, Baker P, Han- using ultra-high resolution ion mobility spec-
kemeier T (2023) Development of a targeted trometry separations. Int J Mol Sci
hydrophilic interaction liquid 7. Paglia G, Smith AJ, Astarita G (2022) Ion
chromatography-tandem mass spectrometry mobility mass spectrometry in the omics era:
based lipidomics platform applied to a corona- challenges and opportunities for metabolomics
virus disease severity study. J Chromatogr A and lipidomics. Mass Spectrom Rev 41:722–
1708:464342 765
3. Konjevod M, Sáiz J, Nikolac Perkovic M, 8. Fernandez-Lima F, Kaplan DA, Suetering J
Nedic Erjavec G, Tudor L, Uzun S, et al (2011) Gas-phase separation using a
Kozumplik O, Barbas C, Zarkovic N, Pivac N, trapped ion mobility spectrometer. Int J Ion
Strac DS (2022) Plasma lipidomics in subjects Mobil Spectrom 14
with combat posttraumatic stress disorder. Free 9. Meier F, Park MA, Mann M (2021) Trapped
Radic Biol Med 189:169–177 ion mobility spectrometry and parallel
4. Camunas-Alberca SM, Moran-Garrido M, Sáiz accumulation-serial fragmentation in proteo-
J et al (2023) Integrating the potential of ion mics. Mol Cell Proteomics 20:100138
mobility spectrometry-mass spectrometry in 10. Di Poto C, Tian X, Peng X et al (2021) Meta-
the separation and structural characterisation bolomic profiling of human urine samples
of lipid isomers. Front Mol Biosci 10:1112521
UHPLC-TIMS-PASEF®-MS for Lipidomics: From Theory to Practice 237
Abstract
The final aim of metabolomics is the comprehensive and holistic study of the metabolome in biological
samples. Therefore, the use of instruments that enable the analysis of metabolites belonging to various
chemical classes in a wide range of concentrations is essential, without compromising on robustness,
resolution, sensitivity, specificity, and metabolite annotation. These characteristics are crucial for the analysis
of very complex samples, such as wine, whose metabolome is the result of the sum of metabolites derived
from grapes, yeast(s), bacteria(s), and chemical or physical modification during winemaking. In recent
years, a big advantage, in this direction, was the hardware developments on hyphenated instruments that
enable the integration of liquid chromatography (LC), ion mobility spectrometry (IMS), and mass spec-
trometry (MS). This chapter describes an LC-IMS-MS protocol for the analysis of wine and grape samples
as well as the use of IMS data in metabolite annotation.
Key words Vitis, Metabolomics, Metabolite annotation, Traveling wave ion mobility, CCS, Mass
spectrometry, HDMSE
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_13,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
239
240 Vania Sáez et al.
Fig. 1 Schematic overview of the LC-IMS-MS protocol for wine and grape samples. The UPLC-TWIMS-QToF
protocol has two separations: the first is pre-ionization in the UPLC column, and the second is post-ionization
in the TWIMS device. In the HDMSE DIA method, the trap acts as an ion guide, the TWIMS device achieves the
separation by ion mobility, and the transfer device operates as collision cell. In the HDMSE DIA, the precursor
ion is associated with drift time data, and the fragments are matched with the precursor by retention time
correlation data. This strategy allows to obtain multidimensional data in a single injection
2 Materials
2.2 Devices, Borosilicate glass syringe of 0.5 mL with stainless steel and PTFE
Equipment, and luer lock (see Note 2).
Instruments UPLC Column Acquity TM Premier HSS T3 with VanGuard
Fit, 1.8 μm, 2.1 × 100 mM, Waters™ (SKU:186009471).
Acquity I-class UPLC with a binary solvent manager, sample
manager FTN-1, column manager, and sample organizer (Waters
242 Vania Sáez et al.
3 Methods
3.1.2 CCS Major MIX The CCS calibration mix, referred to as the CCS Major MIX
IMS/ToF Calibration IMS/ToF solution, allows automated or manual mass and CCS
Solution calibration for high-resolution mass spectrometers with ion mobil-
ity. The CCS Major MIX IMS/ToF content (see Note 3) enables
calibration up to 2000 Da in negative and positive acquisition
modes. The solution was prepared using a glass Pasteur pipette
following the detailed procedure described by the manufacturer
[19]. After calibration, the solution should be stored at 4 °C.
3.1.3 Mobile Phase and The chromatographic separation required two mobile phases:
Wash Solution Preparation (A) water with 0.1% formic acid and (B) methanol with 0.1% formic
acid. For mobile phase A in a 2 L volumetric flask, fill half the flask
with LC-MS-grade water, add 2 mL of formic acid, and fill the
volumetric flask with LC-MS-grade water. For mobile phase B in a
2 L volumetric flask, fill 1 L with LC-MS-grade methanol, add
2 mL of formic acid, and fill to its total volume with LC-MS-
grade methanol. Pour the mobile phase into a certified class A
borosilicate bottle.
In addition to mobile phases, this method also requires the
preparation of a syringe wash, seal wash, and purge and wash
solutions. The seal wash solution is composed of water and aceto-
nitrile (90:10). To prepare, mix 450 mL of MS-grade water and
Advanced LC-IMS-MS Protocol for Holistic Metabolite Analysis in Wine. . . 243
Table 1
IMS-MS operating parameters for grape and wine UPLC-TWIMS-QToF protocol
3.2.2 LC Analytical The column manager is set at 40 °C, the injection volume is 2 μL,
Parameters the flow rate is 0.3 mL/min, and the multistep gradient is as
follows: 0.0–1.0 min (0% mobile phase B), 1.0–3.0 min (0–10%
mobile phase B), 3.0–12.0 min (10–40% mobile phase B),
12.0–16.0 min (40–100% mobile phase B), 16.0–18.0 min (100%
mobile phase B), 18.0–18.10 min (100–0% mobile phase B),
18.10–20.0 min (0% mobile phase B). Before inserting the rack
with the sample vials in the autosampler and/or sample organizer,
activate the temperature control and set to 4 °C in the sample
manager of the UPLC system. Place the vials once the temperature
is stable.
3.2.3 TWIMS-QToF The DIA HDMSE method was developed for the resolution acqui-
Method sition mode for both polarities (ESI+ and ESI-), with a normal
dynamic range in continuum acquisition mode ranging from
50 to 2000 Da with a scan time of 0.2 s. The source conditions,
IMS settings, Triwave, and Triwave DC for ESI+ and ESI- para-
meters are described in Table 1 (Subheading 3.2.1).
The collision energies for the functions in the HDMSE method
differ. For ESI+ and ESI-, the low-energy function is disabled (<
5 V). The collision energy was ramped from 20 to 50 V in ESI-
mode acquisition and from 20 to 40 V in ESI+ in the transfer
collision cell for the high energy function. An LM solution (Sub-
heading 3.1.1) was infused throughout the entire data acquisition
process (Table 1). The acquisition interval for LM is 15 s, with a
scan average of 3. The HDMSE method acquires the LM but does
not apply the correction (see Subheading 3.4, Note 5).
3.3 Sample Prior to sample preparation and/or data acquisition using analytical
Preparation and Data methodology, the untargeted metabolomic workflow must include
Acquisition relevant steps, such as experimental design and pre-analytical pro-
cess [20, 21].
The collection and organization of metadata must begin with
experimental design and sampling. The documented metadata
must be clear and accessible to any researcher and must track any
changes or observations made to the samples during the metabo-
lomic workflow (see Note 6). Codify the samples according to a
random sequence for sample preparation, extraction, and analysis
[22]. Label each sample with the relevant information. It is prefer-
able to use printed labels rather than handwritten labels to avoid
unreliable readings and/or ink smears during the metabolomic
workflow.
3.3.3 Build a Sample For data acquisition, injections must follow the previously estab-
Sequence lished random sequence of the sample preparation step. The first
injection of a data batch sample must be a blank, followed by the
injection of 5 pooled QC from the same vial. Begin sample injec-
tions in random order from injection number 6. A pooled QC must
be injected for every 6 sample injections. Add an injection of the
pooled QC at the end of the sample batch.
3.3.4 LC-IMS-MS Start- Prior to initiating data acquisition, clean the source of the IMS-MS
Up and Pretest analyzer (ESI) in accordance with the manufacturer’s instructions.
Load and purge the mobile phases A and B and the wash solutions,
and then prime the mobile phases for 8 min at a ratio of 50% mobile
phase A to 50% mobile phase B. Proceed with the calibration of the
IMS-MS analyzer as described in Subheading 3.2.1. Install the
column, turn on the temperature (40 °C), and gradually bring
the flow rate to 0.3 mL/min. Wait until the pressure difference is
less than 10 psi before starting the data acquisition.
Advanced LC-IMS-MS Protocol for Holistic Metabolite Analysis in Wine. . . 247
3.4 Data Analysis, The LC-IMS-MS data format is supported by MZmine [23],
Processing, and MS-DIAL [24], MetaboScape, Skyline [4, 25], and Progenesis
Visualization QI [8]. This section describes the analysis and processing of data
using Progenesis QI commercial software [26] (see Note 8).
To demonstrate data analysis processing and visualization, a
small training dataset of 18 runs (including 1 blank, 11 samples,
and 7 pooled QC files) of wine samples (4 Cabernet Sauvignon,
4 Barbera, and 3 Timorasso samples) in ESI+ and ESI- was
acquired using the described UPLC-TWIMS-QToF methodology
(see Note 9). The 36 raw data files (ESI+ and ESI-) are accessible
in the MetaboLights (https://www.ebi.ac.uk/metabolights/
MTBLS8329) open repository [27]. Figure 2 shows a multidimen-
sional visualization of the HDMSE data for a wine sample, includ-
ing the chromatograms at low and high energy applied in the
transfer (Fig. 2a, b) and drift time in milliseconds v/s m/z plot
(Fig. 2c).
For the analysis of the training dataset with Progenesis QI,
initially it is necessary to select the data format and ionization
mode in order to import correctly the raw data into the software
(see Note 10) [26]. To start the processing, a pooled QC must be
selected as an alignment reference (see Note 11). If the CCS
calibration is performed (Subheading 3.2.1), the software
Fig. 2 HDMSE DIA chromatogram of Cabernet Sauvignon wine in ESI- (a) low energy, (b) high energy, and (c)
drift time vs. m/z plot for all features with different charges in the range of m/z 50–2000, with a drift time
range of 1–13 milliseconds in ESI- acquisition mode
248 Vania Sáez et al.
Fig. 3 (a) PCA biplot with no filter applied; (b) PCA biplot with a filter p-value <0.01, max abundance >2 for
HDMSE DIA of Cabernet Sauvignon, Timorasso, and Barbera wine samples acquired in ESI+
3.5 Metabolite The strategies for metabolite annotation process of the selected
Annotation features include resources such as in-house libraries (RT, MSE spec-
tra, and TWCCSN2 value) built with analytical standards and man-
aged by UNIFI software (Subheading 3.5.1), spectral public or
published databases/databanks, the Metabolic Profiling CCS
Library, and the curated ion mobility atlas for experimental and
predicted CCS values [10].
3.5.1 In-House Library UNIFI commercial software enables the construction, manage-
ment, and organization of information for the creation of a custo-
mized scientific in-house library for attaining confident annotation
(Rt, m/z, MSE, and CCS value) in real samples.
The first step in creating a customized library is to acquire the
raw data from the reference standard(s). The concentration of
standard compounds can range from 1 ppm to 10 ppm. Import
the HDMSE DIA raw file(s) and create a sample set from the
UNIFI main menu.
The second step is to create the library from the tools tab in the
scientific library section. Each compound in the library is stored as
an “item” and can be created, moved, or copied from the “Manage
Library” section. Each item includes data from ESI+ and ESI-
mode acquisitions. The metabolite structure must be uploaded in
“. mol” format.
The third step is to create an analysis method for the output
HDMSE DIA raw data (select the option to generate a process-only
method for accurate mass screening on IMS data). Through the
menu path Create and Analysis Method, the user can store data in a
library from standard analysis or to analyze samples by comparing a
standard from a library with the selected features. After the creation
of the Analysis Method, a new menu should be displayed including
the tabs: “purpose,” “processing,” “reporting,” and “history.”
Import and filter the ionization technique of the item(s)/com-
pounds from the libraries in the purpose task.
Select the processing parameters; we recommend using a toler-
ance of 0.1 min for the retention time, < 10 ppm for m/z,
and < 2% for the CCS value. The identification of the adducts is
based on the specific setting in the “processing” section, including
the expected adduct and the possible transformations (glucuroni-
dation, glutathione S-conjugation, oxidation, etc.). Configure the
correct LM settings in the “LockMass” tab. Check the calibration
report (Subheading 3.2.1.) for LM and input the reference value
(see Note 5) and the drift time from the calibration report. After
analysis, the settings can be modified for further processing.
The fourth step is to run the analysis method on the sample set
generated in the first step. The processing includes peak and chan-
nel processing, identification, extract ion chromatograms (XIC),
and comparison. The processing step can last from minutes to
250 Vania Sáez et al.
Fig. 4 (a) Component summary of polyphenols identified in a Cabernet Sauvignon wine sample by UNIFI from
raw HDMSE data in ESI-. Four metabolites were annotated at the first level using the in-house library with a Δ
ppm < 5 and Δ CCS% < 1%; (b) HDMSE spectra data for catechin precursor ion; and (c) ion fragments, with
fragment prediction indicated in blue
3.5.2 Reporting an Curated database/libraries, in silico data spectra and CCS predic-
Annotated Feature tion [4, 10, 28], and scoring algorithms are used as resources for
the group of features which are not covered with the in-house
library to support the metabolite annotation.
Advanced LC-IMS-MS Protocol for Holistic Metabolite Analysis in Wine. . . 251
4 Notes
Acknowledgments
References
1. Cumeras R, Figueras E, Davis CE, Baumbach lipid fingerprinting and identification. Anal
JI, Gràcia I (2015) Review on ion mobility Chem 87(2):1137–1144. https://doi.org/10.
spectrometry. Part 2: hyphenated methods 1021/ac503715v
and effects of experimental parameters. Analyst 8. Paglia G, Astarita G (2017) Metabolomics and
140:1391–1410. https://doi.org/10.1039/ lipidomics using traveling-wave ion mobility
c4an01101e mass spectrometry. Nat Protoc 12:797–813.
2. Lapthorn C, Pullen F, Chowdhry BZ (2012) https://doi.org/10.1038/nprot.2017.013
Ion mobility spectrometry-mass spectrometry 9. Righetti L, Dreolin N, Celma A,
(IMS-MS) of small molecules: separating and McCullagh M, Barknowitz G, Sancho JV, Dall
assigning structures to ions. Mass Spectrom ‘Asta C (2020) Travelling wave ion mobility-
Rev 32(1):43–71. https://doi.org/10.1002/ derived collision cross section for mycotoxins:
mas.21349 investigating interlaboratory and interplatform
3. Gabelica V, Shvartsburg AA, Afonso C, reproducibility. J Agric Food Chem 68(39):
Barran P, Benesch JLP, Bleiholder C, Bowers 10937–10943. https://doi.org/10.1021/acs.
MT, Bilbao A, Bush MF, Campbell JL, Cam- jafc.0c04498
puzano IDG, Causon T, Clowers BH, Creaser 10. Zhou Z, Luo M, Chen X, Yin Y, Xiong X,
CS, De Pauw E, Far J, Fernandez-Lima F, Wang R, Zhu ZJ (2020) Ion mobility collision
Fjeldsted JC, Giles K, Groessl M, Hogan CJ, cross-section atlas for known and unknown
Hann S, Kim HI, Kurulugama RT, May JC, metabolite annotation in untargeted metabo-
McLean JA, Pagel K, Richardson K, Ridgeway lomics. Nat Commun 11. https://doi.org/10.
ME, Rosu F, Sobott F, Thalassinos K, Valen- 1038/s41467-020-18171-8
tine SJ, Wyttenbach T (2019) Recommenda- 11. Blaženović I, Kind T, Ji J, Fiehn O (2018)
tions for reporting ion mobility mass Software tools and approaches for compound
spectrometry measurements. Mass Spectrom identification of LC-MS/MS data in metabo-
Rev 38(3):291–320. https://doi.org/10. lomics. Meta 8:31. https://doi.org/10.3390/
1002/mas.21585 metabo8020031
4. Paglia G, Smith AJ, Astarita G (2022) Ion 12. Diamantidou D, Zotou A, Theodoridis G
mobility mass spectrometry in the omics era: (2018) Wine and grape marc spirits metabolo-
challenges and opportunities for metabolomics mics. Metabolomics 14. https://doi.org/10.
and lipidomics. Mass Spectrom Rev 41:722– 1007/S11306-018-1458-1
765. https://doi.org/10.1002/mas.21686
13. Arapitsas P, Ugliano M, Marangon M,
5. May JC, McLean JA (2015) Ion mobility-mass Piombino P, Rolle L, Gerbi V, Versari A, Mat-
spectrometry: time-dispersive instrumentation. tivi F (2020) Use of untargeted liquid
Anal Chem 87(3):1422–1436. https://doi. chromatography-mass spectrometry metabo-
org/10.1021/ac504720m lome to discriminate italian monovarietal red
6. Richardson K, Giles K (2021) Chapter 4. wines, produced in their different terroirs. J
Travelling wave ion mobility separation: basics Agric Food Chem 68(47):13353–13366.
and calibration. In: Ion mobility-mass spec- https://doi.org/10.1021/acs.jafc.0c00879
trometry, fundamentals and applications, pp 14. Sáez V, Schober D, González Á, Arapitsas P
8 3 – 1 0 4 . h t t p s : // d o i . o r g / 1 0 . 1 0 3 9 / (2021) LC–MS-based metabolomics discrimi-
9781839162886-00083 nates premium from standard chilean
7. Paglia G, Angel P, Williams JP, Richardson K, cv. Cabernet sauvignon wines from different
Olivos HJ, Thompson JW, Menikarachchi L, valleys. Meta:11. https://doi.org/10.3390/
Lai S, Walsh C, Moseley A, Plumb RS, Grant metabo11120829
DF, Palsson BO, Langridge J, Geromanos S, 15. Arapitsas P, Della CA, Gika H, Narduzzi L,
Astarita G (2015) Ion mobility-derived colli- Mattivi F, Theodoridis G (2016) Studying the
sion cross section as an additional measure for effect of storage conditions on the metabolite
Advanced LC-IMS-MS Protocol for Holistic Metabolite Analysis in Wine. . . 255
guidelines for fair data and metadata manage- analysis: chemical analysis working group
ment. Meta 11. https://doi.org/10.3390/ (CAWG) metabolomics standards initiative
metabo11110757 (MSI). Metabolomics 3:211–221. https://
31. Sumner LW, Amberg A, Barrett D, Beale MH, doi.org/10.1007/s11306-007-0082-2
Beger R, Daykin CA, Fan TWM, Fiehn O, 32. Schymanski EL, Jeon J, Gulde R, Fenner K,
Goodacre R, Griffin JL, Hankemeier T, Ruff M, Singer HP, Hollender J (2014) Iden-
Hardy N, Harnly J, Higashi R, Kopka J, Lane tifying small molecules via high resolution mass
AN, Lindon JC, Marriott P, Nicholls AW, Reily spectrometry: communicating confidence.
MD, Thaden JJ, Viant MR (2007) Proposed Environ Sci Technol 48(4):2097–2098.
minimum reporting standards for chemical https://doi.org/10.1021/es5002105
Chapter 14
Abstract
Metanephrines (metanephrine [MN] and normetanephrine [NMN]) are O-methylated metabolites derived
from the catecholamines, epinephrine, and norepinephrine, respectively. High concentrations of metane-
phrines have been observed in individuals with pheochromocytoma, a neuroendocrine tumor. Measure-
ment of metanephrines in urine is used to screen for the tumor. Analysis using liquid chromatography-
tandem mass spectrometry (LC-MS/MS) is recommended due to the high sensitivity, specificity, and
throughput of the technique. Herein, we describe an optimized LC-MS/MS assay for the analysis of
urinary metanephrines.
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_14,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
257
258 Marlene N. Thaitumu and Elizabeth L. Frank
2 Materials
2.1 Samples Timed (24-h) or random urine. Timed specimens should be refri-
gerated during collection. Specimens are stable for 2 weeks refri-
gerated or 1 month frozen at or below 20 C (see Note 2).
2.3.1 Standards: 0.5 mM 1. Prepare 0.1 L of 0.5 mM combined stock solution of MN (()
MN and NMN Combined metanephrine hydrochloride—Toronto Research Chemicals,
Stock Solution North York, Ontario, Canada) and NMN (DL-normetanephr-
ine hydrochloride—Sigma-Aldrich, St. Louis, Missouri, USA).
The metanephrine standards are in the form of hydrochloride
salts. Prepare stock solutions using information from Table 1.
2. Determine the metanephrines weights for a 0.5 mmoL solu-
tion using the following formulas:
(a) Weight ¼ Concentration ðmMÞ Volume ðLÞ
mg
Molecularweight mmoL =Purity mg
233:69
(b) MN : Weight ¼ 0:500 mM 0:100 L mmoL ¼ 11:92 mg
0:98
mg
219:67
(c) NMN : Weight ¼ 0:500 mM 0:100 L mmoL ¼ 11:21 mg
0:98
3. Using an analytical balance accurate to 5 decimal places, accu-
rately weigh MN and NMN in separate 10 mL glass screw cap
tubes.
4. Transfer both standards to the same 100 mL volumetric flask
by dissolving and washing out the glass tubes with the equiva-
lent of three tube volumes using deionized water to ensure
complete transfer.
5. Using a volumetric pipette, add 10 mL of 0.05 mol/L hydro-
chloric acid (HCL) prepared in house to the flask. Final HCL
concentration in the stock is 5 mM.
Table 1
Values for calculation of weights for a combined MN and NMN 0.5 mM stock standard
Table 2
Volume of 0.5 mM stock standard (MN and NMN) solution required for
250 mL calibrators preparation
2.3.2 Standards: 1. Dilute the stock to five combined MN and NMN calibration
Combined MN and NMN standards: 25, 100, 500, 2000, and 7000 nmol/L using the
Calibrator Standards (25, information in Table 2.
100, 500, 2000, and 2. Using a graduated cylinder, add 200 mL of deionized water to
7000 nmol/L) five 250 mL volumetric flasks correctly labeled with the calibra-
tor concentration.
2. Using a volumetric pipette, add 25 mL of 0.05 mol/L HCL
(prepared in house) to the flasks. Final concentration of HCL
in each flask is 5 mM.
3. Using a pipette, add the appropriate volume of the combined
stock solution for each calibrator per information in Table 2.
Bring the calibrators to 250 mL volume with deionized water.
4. Mix the calibrator solutions thoroughly by inverting the volu-
metric flask (capped with a stopper) three times (see Note 4).
Table 3
Information for preparation of a combined MN and NMN 0.5 mM ISTD stock solution
Compound Conc (mM) Vol (L) ISTD HCL salt MW (g/moL) Purity
rac-Metanephrine-d3·HCl (α1, β2) 0.0125 0.5 236.71 0.98
rac-Normetanephrine-d3·HCl (α1, β2) 0.0125 0.5 222.68 0.98
mg
236:71
(b) MN d3 weight ¼ 0:125 mM 0:5 L mmoL ¼ 1:51 mg
0:98
mg
222:68
(c) NMN d3 weight ¼ 0:0125 mM 0:5 L mmoL ¼ 1:42 mg
0:98
2. Using an analytical balance accurate to 5 decimal places, accu-
rately weigh MN and NMN ISTDs in a separate 10 mL glass
screw cap tubes.
3. Transfer both MN and NMN ISTDs to the same 0.5 L volu-
metric flask by dissolving and washing out the glass tubes with
the equivalent of three tube volumes using deionized water to
ensure complete transfer.
4. Using a graduated cylinder, add 50 mL of 0.05 mol/L HCL
solution (prepared in house) to the volumetric flask (final
concentration of HCL in solution is 5 mM).
5. Bring the combined ISTD working solution to 0.5 L volume
with deionized water. Mix the ISTD working solutions thor-
oughly by inverting the volumetric flask (capped with a stop-
per) three times (see Note 5).
3 Methods
Table 4
LC gradient for separation of urinary metanephrines
3.4 Data Analysis 1. Data is analyzed using appropriate software, e.g., SCIEX Ana-
lyst software.
2. Calibration for each analyte with 1/x weighting should pro-
duce a correlation coefficient > 0.99.
3. Evaluate chromatography of calibrators and controls for reten-
tion time, peak shape, IS area counts, and ion ratios (see Fig. 2,
from Clarke and Frank [19]). Minimum IS peak area counts
must be established during method validation.
4. QC concentrations must fall within 2 standard deviations of
the set mean (see Note 7).
5. Data can be analyzed using an alternative set of transitions for
MN, if necessary due to interference (see Note 8).
Analysis of Urinary Metanephrines Using Liquid Chromatography Tandem Mass. . . 265
Table 5
Metanephrines MRMs and DPs
4 Notes
References
1. Joyce CM, Melvin A, O’Shea PM, Costelloe recommendations for clinical practice from
SJ, O’Halloran DJ (2020) Case report of a the first international symposium. Nat Clin
phantom pheochromocytoma. Biochem Med Pract Endocrinol Metab 3(2):92–102
30(2):325–330 9. Pamporaki C, Därr R, Bursztyn M, Glöckner S,
2. Zhang W, Li X, Li W, Zhang Y, Cai J, Feng S, Bornstein SR, Lenders JWM, Pacak K,
Sun Z (2024) Clinical diagnosis of pheochro- Krinner A, Eisenhofer G (2013) Plasma-free
mocytoma and paraganglioma-induced sec- vs deconjugated metanephrines for diagnosis
ondary hypertension through UPLC-MS/MS of phaeochromocytoma. Clin Endocrinol
analysis of plasma catecholamines and their 79(4):476–483
metabolites. J Clin Hypertens 26(4):416–424 10. Eisenhofer G, Prejbisz A, Peitzsch M,
3. Jain A, Baracco R, Kapur G (2019) Pheochro- Pamporaki C, Masjkur J, Rogowski-Lehmann-
mocytoma and paraganglioma—an update on N, Langton K, Tsourdi E, Pe˛czkowska M,
diagnosis, evaluation, and management. Fliedner S et al (2018) Biochemical diagnosis
Pediatr Nephrol 35(4):581–594 of chromaffin cell tumors in patients at high
4. Ceccato F, Mantero F (2019) Monogenic and low risk of disease: plasma versus urinary
forms of hypertension. Endocrinol Metab free or deconjugated O-methylated catechol-
Clin N Am 48(4):795–810 amine metabolites. Clin Chem 64(11):
5. Lenders JWM, Duh Q-Y, Eisenhofer G, 1646–1656
Gimenez-Roqueplo A-P, Grebe SKG, Murad 11. Eisenhofer G, Åneman A, Hooper D,
MH, Naruse M, Pacak K, Young WF (2014) Rundqvist B, Friberg P (2002) Mesenteric
Pheochromocytoma and paraganglioma: an organ production, hepatic metabolism, and
endocrine society clinical practice guideline. J renal elimination of norepinephrine and its
Clin Endocrinol Metabol 99(6):1915–1942 metabolites in humans. J Neurochem 66(4):
6. Eisenhofer G, Peitzsch M, McWhinney BC 1565–1573
(2016) Impact of LC-MS/MS on the labora- 12. Eisenhofer G, Coughteree MW, Goldstein DS
tory diagnosis of catecholamine-producing (1999) Dopamine sulphate: an enigma
tumors. TrAC Trends Anal Chem 84 resolved. Clin Exp Pharmacol Physiol Suppl
(PB):106–116 26:S41–S53
7. Kim HJ, Lee JI, Cho YY, Lee SY, Kim JH, Jung 13. Eisenhofer G, Lenders J (2012) Rapid circula-
BC, Kim SW, Chung JH, Min Y-K, Lee M-S tory clearances and half-lives of plasma free
et al (2015) Diagnostic accuracy of plasma free metanephrines. Clin Endocrinol 77(3):
metanephrines in a seated position compared 484–485
with 24-hour urinary metanephrines in the 14. Grouzmann E, Drouard-Troalen L, Baudin E,
investigation of pheochromocytoma. Endocr J Pierre-François P, Müller B, Grand D, Buclin T
62(3):243–250 (2010) Diagnostic accuracy of free and total
8. Pacak K, Eisenhofer G, Ahlman H, Bornstein metanephrines in plasma and fractionated
SR, Gimenez-Roqueplo A-P, Grossman AB, metanephrines in urine of patients with pheo-
Kimura N, Mannelli M, McNicol AM, Tischler chromocytoma. Eur J Endocrinol 162(5):
AS (2007) Pheochromocytoma: 951–960
268 Marlene N. Thaitumu and Elizabeth L. Frank
Abstract
Metabolomics is the area of research, which strives to obtain complete metabolic fingerprints, to detect
differences between them and to provide hypothesis to explain those differences (Schripsema J, Dagnino D,
Handbook of chemical and biological plant analytical methods. Wiley, New York, 2015). However,
obtaining complete metabolic fingerprints is not an easy task. Metabolite extraction is a key step during
this process, and much research has been devoted to finding the best solvent mixture to extract as much
metabolites as possible.
Here a procedure is described for analysis of both polar and apolar metabolites using a two-phase
extraction system. D2O and CDCl3 are the solvents of choice, and their major advantage is that, for the
identification of the compounds, standard databases can be used because D2O and CDCl3 are the solvents
most commonly used for pure compound NMR spectra. The procedure enables the absolute quantification
of components due to the addition of suitable internal standards. The extracts are also suitable for further
analysis with other systems like LC-MS or GC-MS.
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_15,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
269
270 Jan Schripsema and Denise Dagnino
Fig. 1 1H NMR spectra from the (a) D2O and (b) CDCl3 extracts of coffee powder obtained with the two-phase
extraction. In the CDCl3 phase, signals from caffeine (marked with C) and pyridine (marked with P) are
indicated. In the D2O phase. Signals from TMSP and trigonelline (marked with T) are indicated
2 Materials
Table 1
Calculations with the NMR data shown in Fig. 1. For each spectrum, the
calculation of the absolute quantity of a specific compound is illustrated
2.2 Preparation of Deuterated solvents should be at least 99.8% deuterated. The inter-
the Internal Standard nal standards should be of analytical grade or higher. Weighing of
Solution (ISS) standards and solvents should be carried out with maximum
precision.
Internal standard Solution 1 (ISS1): Weigh 150 mg of pyridine
(MW: 79.10 g/moL) in an empty bottle and add 100 g
CDCl3.
Internal standard Solution 2 (ISS2): Weigh 50 mg of TMSP (MW:
172.27 g/moL) in an empty bottle and add 100 g D2O.
Internal standard solutions can be kept for up to 1 year at room
temperature, if properly handled. Flasks should be kept tightly
closed when not in use. In case of handling a large number of
samples, distribute the prepared ISS in smaller bottles so as to
minimize the time each flask is open.
3 Methods
3.1 Extraction and The extraction described here has been used for dried plant material
Analysis (freeze dried) and also for the analysis of dried bacteria and food
stuff such as teas, coffee, butter, and cheese. The procedure is
274 Jan Schripsema and Denise Dagnino
Fig. 2 Schematic representation of the two-phase extraction for comprehensive analysis of the plant
metabolome by NMR
4 Notes
References
1. Schripsema J, Dagnino D (2015) 7. Beltran A, Suarez M, Rodrı́guez MA et al
Metabolomics. In: Hostettmann K, (2012) Assessment of compatibility between
Stuppner H, Marston A, Chen S (eds) Hand- extraction methods for NMR and LC/MS-
book of chemical and biological plant analytical based metabolomics. Anal Chem 84:5838–
methods, 1st edn. Wiley, New York 5844. https://doi.org/10.1021/ac3005567
2. Schripsema J (2010) Application of NMR in 8. Kim HK, Choi YH, Verpoorte R (2010)
plant metabolomics: techniques, problems and NMR-based metabolomic analysis of plants.
prospects. Phytochem Anal 21:14–21. https:// Nat Protoc 5:536–549. https://doi.org/10.
doi.org/10.1002/pca.1185 1038/nprot.2009.237
3. Deda O, Gika HG, Wilson IA, Theodoridis GA 9. Schripsema J (2019) Similarity and differential
(2015) An overview of fecal preparation for NMR spectroscopy in metabolomics: applica-
global metabolic profiling. J Pharm Biomed tion to the analysis of vegetable oils with H-1
Anal 113:137–150. https://doi.org/10. and C-13 NMR. Metabolomics 15:39.
1016/j.jpba.2015.02.006 https://doi.org/10.1007/s11306-019-
4. Kim H-S, Park SJ, Hyun S-H et al (2011) 1502-9
Biochemical monitoring of black raspberry 10. Schripsema J, Merlim RD, Parvan LG et al
(Rubus coreanus Miquel) fruits according to (2022) Towards a holistic view of tablet quality,
maturation stage by 1H NMR using multiple an extensive study on paracetamol tablets with
solvent systems. Food Res Int 44:1977–1987. nuclear magnetic resonance using similarity
https://doi.org/10.1016/j.foodres.2011. calculations, differential NMR and hierarchical
01.023 cluster analysis. J Pharm Biomed Anal 215:
5. Heyman HM, Meyer JJM (2012) NMR-based 114773. https://doi.org/10.1016/j.jpba.
metabolomics as a quality control tool for 2022.114773
herbal products. S Afr J Bot 82:21–32. 11. Schripsema J (2008) Comprehensive analysis
https://doi.org/10.1016/j.sajb.2012.04.001 of polar and apolar constituents of butter and
6. Kim HK, Verpoorte R (2010) Sample prepara- margarine by nuclear magnetic resonance,
tion for plant metabolomics. Phytochem Anal reflecting quality and production processes. J
21:4–13. https://doi.org/10.1002/pca.1188 Agric Food Chem 56:2547–2552. https://doi.
org/10.1021/jf072274c
Chapter 16
Abstract
High Resolution-Magic Angle Spinning (HR-MAS) solid-state NMR spectroscopy is finding increasing
application in the analysis of solid foods, bypassing the need for complicated solvent extraction procedures.
In the present protocol, we report a simple analytical approach based on HR-MAS NMR spectroscopy for
the phenolic profiling of olive fruits, flesh, or skin. This approach allows the facile characterization of
phenolic compounds in olive fruits cultivated for extra-virgin olive oil production as a function of matura-
tion and variety, in addition to processing technology for table olives.
Key words NMR spectroscopy, Solid state, Phenolic profiling, Food analysis, Authentication
1 Introduction
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1_16,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2025
277
278 Efstathia Manolopoulou and Apostolos Spyros
2 Materials
3 Methods
3.1 Sample 1. Store 1–2 g of olive fruit, flesh, or skin at -18 °C for 24 h (see
Preparation Note 1).
2. Cut frozen sample to small pieces and weigh.
3. Place sample in 50 mL glass flask, connect to freeze-drier, and
freeze-dry for at least 16 h.
4. Weigh sample again to calculate the amount of moisture.
5. Grind sample with a pestle and mortar under liquid nitrogen
(see Note 2).
6. Pre-weigh the empty rotor and introduce 5–10 mg of frozen
olive fruit/flesh/skin powder using the rotor loading funnel
accordingly (see Note 3).
7. Add carefully 35–40 μL of deuterated methanol into the rotor/
insert via pipette (see Note 4).
8. Completely seal the rotor by placing the insert screw and the
rotor cap using the provided rotor tools.
9. Insert the rotor in the probe either directly or using an auto-
sampler, if available.
10. Run NMR experiment protocol.
3.2 1H NMR 1. Set the probe temperature to 298 K, spin the sample at the
Spectroscopy magic angle at 4 kHz, and wait (3–5 min) until the sample
Experimental Protocol temperature is equilibrated.
with Water 2. Lock, tune, and shim the sample according to standard NMR
Presaturation spectrometer procedures.
3. Load a standard solvent presaturation (Bruker: zgpr or zgcppr)
pulse program with default spectrometer parameters, and set
the frequency of the residual water signal exactly on resonance.
4. Record a water-suppressed 1H NMR spectrum with para-
meters: SW = 20 ppm, NS = 160 scans, DS = 16 dummy
scans, AQ = 4.0 s, D1 = 1 s.
5. Perform Fourier transformation, phase correction, and base-
line correction according to standard NMR spectrometer
(or processing software) procedures.
280 Efstathia Manolopoulou and Apostolos Spyros
3.3 NMR Data Initial phenolic profiling may be achieved by comparing the 1H
Analysis NMR spectra of olives with those of standard olive phenols
obtained in the same solvent, and it should be noted that liquid-
state NMR spectra of phenol standards may also be used for identi-
fication, since chemical shift differences between liquid-state and
HR-MAS NMR spectra of olives, if present, are usually minimal.
Furthermore, the standard arsenal of homonuclear and heteronuc-
lear 2D NMR experiments for structure elucidation (gCOSY,
gHSQC, gHMBC, etc.) is available also for the HR-MAS NMR
spectroscopy platform and may be used for compound identifica-
tion and verification. Other means of phenol identification include
the use of publicly available NMR spectral databases, such as
FoodDB [26] and the Biological Magnetic Resonance Bank [27].
Figure 1 depicts the aromatic region of the HR-MAS solid-
state 1H NMR spectra of olive flesh obtained from four different
olive varieties, along with the NMR spectra of four standard phe-
nolic compounds (verbascoside 1, oleocanthal 2, oleacin 3, oleur-
opein 4) that have been reported as constituents of olive oil and
olives [18, 28–30]. Koroneiki, Chondrolia, and Lianolia samples
Fig. 1 High Resolution-Magic Angle Spinning (HR-MAS) solid-state 1H NMR spectra (aromatic region) of model
compounds (verbascoside 1, oleocanthal 2, oleacin 3, oleuropein 4) and of olive fruit (flesh) of different
varieties (Koroneiki, Chondrolia, Kavala, Lianolia) in MeOD-d4, obtained by applying the method developed in
this study at a proton frequency of 400.2 MHz. Numbers 5 and 6 correspond to the compounds ligstroside and
cornoside, respectively
Olive Fruit Phenolic Profiling Using High Resolution-Magic Angle. . . 281
were obtained from Crete, Greece, while the Kavala sample was a
local olive variety obtained from Kavala, Greece. The phenolic
composition of the olive flesh samples of Fig. 1 differs significantly,
as evident by comparison with the NMR spectra of the standard
compounds. Oleacin 3 is an important and well-known polyphenol
also present in high-quality extra-virgin olive oil and is found to be
present in Chondrolia, Kavala, and Lianolia olive fruit. (see Note 5)
Koroneiki variety on the other hand is found to contain mainly
oleuropein 4, ligstroside 5, and verbascoside 1. Chondrolia and
Kavala olives also contain another interesting phenolic compound,
cornoside 6, [31] and smaller amounts of hydroxytyrosol and
tyrosol, while the doublet at δ 5.95 in the NMR spectrum of
Chondrolia is assigned to halleridone, which is a well-known deg-
radation product of cornoside [29, 31].
As a second example, Fig. 2 depicts the HR-MAS solid-state
1 1
H- H gCOSY 2D NMR spectrum (aromatic region) of olive skin
obtained from processed table olives of Chalkidiki variety. Corno-
side and halleridone are identified as components of the olive skin
4 Notes
1. If only the flesh or skin are the focus of analysis, they must be
separated carefully before freezing.
2. Grinding was not used for olive skin samples, instead these
were just cut in very small pieces after freeze drying.
3. It is important to complete the rotor loading procedure as fast
as possible, since partial thawing of the olive sample will make
the use of the funnel ineffective. In this respect, the loading
funnel may also be cooled right before use.
4. If the sample quantity is less than 5 mg, the solvent volume may
be slightly increased to 45 μL, in order to make sure that the
active volume of the rotor after sealing does not contain any air
pockets that will make shimming inefficient.
5. Oleacin and oleocanthal are present in their hemiacetal form in
CD3OD solutions, as already reported [32].
Acknowledgments
References
1. Valentini M, Ritota M, Cafiero C, Cozzolino S, 2. Spyros A, Dais P (2012) NMR spectroscopy in
Leita L, Sequi P (2011) The HRMAS-NMR food analysis. RSC food analysis monographs
tool in foodstuff characterisation. Magn Reson 3. Santos ADC, Fonseca FA, Lião LM, Alcantara
Chem 49(Suppl 1):S121–S1S5. https://doi. GB, Barison A (2015) High-resolution magic
org/10.1002/mrc.2826 angle spinning nuclear magnetic resonance in
Olive Fruit Phenolic Profiling Using High Resolution-Magic Angle. . . 283
foodstuff analysis. TrAC Trends Anal Chem 13. Vermathen M, Marzorati M, Diserens G,
73:10–18. https://doi.org/10.1016/j.trac. Baumgartner D, Good C, Gasser F et al
2015.05.003 (2017) Metabolic profiling of apples from dif-
4. Corsaro C, Cicero N, Mallamace D, Vasi S, ferent production systems before and after con-
Naccari C, Salvo A et al (2016) HR-MAS and trolled atmosphere (CA) storage studied by 1H
NMR towards foodomics. Food Res Int 89: high resolution-magic angle spinning
1085–1094. https://doi.org/10.1016/j. (HR-MAS) NMR. Food Chem 233:391–400.
foodres.2016.09.033 https://doi.org/10.1016/j.foodchem.2017.
5. Sundekilde UK, Eggers N, Bertram HC 04.089
(2019) NMR-based metabolomics of 14. Song EH, Kim HJ, Jeong J, Chung HJ, Kim
food. In: Methods in molecular biology. HY, Bang E et al (2016) A 1H HR-MAS
Humana Press, pp 335–344 NMR-based metabolomic study for metabolic
6. Tang F, Vasas M, Hatzakis E, Spyros A (2019) characterization of rice grain from various
Magnetic resonance applications in food Oryza sativa L. Cultivars. J Agric Food Chem
analysis. In: Webb GA (ed) Annual reports on 64(15):3009–3016. https://doi.org/10.
NMR spectroscopy. Academic Press, pp 1021/acs.jafc.5b05667
239–306. https://doi.org/10.1016/bs. 15. Augustijn D, de Groot HJM, Alia A (2021)
arnmr.2019.04.005 HR-MAS NMR applications in plant metabo-
7. Brescia MA, Jambrenghi AC, Di Martino V, lomics. Molecules 26(4). https://doi.org/10.
Sacco D, Giannico F, Vonghia G et al (2002) 3390/molecules26040931
High resolution nuclear magnetic resonance 16. Larsen FH, Blennow A, Engelsen SB (2008)
spectroscopy (NMR) studies on meat compo- Starch granule hydration-a MAS NMR investi-
nents: potentialities and prospects. Ital J Anim gation. Food Biophys 3(1):25–32. https://doi.
Sci 1(2):151–158. https://doi.org/10.4081/ org/10.1007/s11483-007-9045-4
ijas.2002.151 17. Manolopoulou E, Spinella A, Caponetti E,
8. Shintu L, Caldarelli S, Franke BM (2007) Dais P, Spyros A (2013) High resolution-
Pre-selection of potential molecular markers magic angle spinning NMR study of olive
for the geographic origin of dried beef by leaves. In: van Duynhoven J, Belton PS, Webb
HR-MAS NMR spectroscopy. Meat Sci 76(4): GA, van As H (eds) Magnetic resonance in
700–707. https://doi.org/10.1016/j.mea food science: food for thought. The Royal
tsci.2007.02.010 Society of Chemistry, pp 164–170
9. Panthi RR, Sundekilde UK, Kelly AL, 18. Beteinakis S, Papachristodoulou A, Gogou G,
Hennessy D, Kilcawley KN, Mannion DT Katsikis S, Mikros E, Halabalaki M (2020)
et al (2019) Influence of herd diet on the NMR-based metabolic profiling of edible
metabolome of Maasdam cheeses. Food Res olives – determination of quality parameters.
Int 123:722–731. https://doi.org/10.1016/ Molecules 25(15):3339. https://doi.org/10.
j.foodres.2019.05.026 3390/molecules25153339
10. Vermathen M, Marzorati M, Vermathen P 19. Irmak Ş, Dıraman H (2022) Chemometric
(2012) Exploring high-resolution magic angle characterization of raw olives from important
spinning (HR-MAS) NMR spectroscopy for Turkish table olive cultivars Cvs. Using HPLC–
metabonomic analysis of apples. Chimia DAD method based on their biophenolic pro-
66(10):747–751. https://doi.org/10.2533/ files. J Food Sci Technol 59(5):1792–1800.
chimia.2012.747 https://doi.org/10.1007/s13197-021-
11. Corsaro C, Mallamace D, Dugo G (2014) A 05190-6
quantitative 1H HR-MAS study of the meta- 20. Rocha J, Borges N, Pinho O (2020)
bolic profile of the PGI cherry tomato of Table olives and health: a review. J Nutr Sci 9:
pachino. In: Nuclear magnetic resonance e57. https://doi.org/10.1017/jns.2020.50
(NMR): theory, applications and technology. 21. Infante R, Infante M, Pastore D, Pacifici F,
Nova Science, pp 163–178 Chiereghin F, Malatesta G et al (2023) An
12. Cicero N, Corsaro C, Salvo A, Vasi S, Giofré appraisal of the oleocanthal-rich extra virgin
SV, Ferrantelli V et al (2015) The metabolic olive oil (EVOO) and its potential anticancer
profile of lemon juice by proton and neuroprotective properties. Int J Mol Sci
HR-MAS NMR: the case of the PGI Interdo- 2 4 ( 2 4 ) . h t t p s : // d o i . o r g / 1 0 . 3 3 9 0 /
nato lemon of Messina. Nat Prod Res 29(20): ijms242417323
1894–1902. https://doi.org/10.1080/ 22. El Qarnifa S, El Antari A, Hafidi A (2019)
14786419.2015.1012166 Effect of maturity and environmental
284 Efstathia Manolopoulou and Apostolos Spyros
A H
Alkyl chloroformate derivatization............................... 206 HDMSE ...............................................240, 245, 247, 252
Amino acids (AAs) ................................... 6, 11, 160, 167, HILIC-MS/MS ...........................................160, 181–202
182, 183, 198, 199, 205–218
Authentication............................................................... 277 I
Identifications.......................................3, 6, 8, 10, 11, 20,
B
27, 32, 46, 56, 60, 68, 70, 73, 75, 78, 84, 99, 103,
Biological samples .............................................6, 56, 109, 110, 111, 114, 121, 123, 125, 166, 174, 175,
112, 114, 118, 124, 131–149 181, 182, 205, 232, 249, 253, 270, 277, 280
Biomarker discovery...................................................... 154 Ion mobility (IM) ........................................ 4, 8, 26, 110,
113, 114, 119, 222, 224, 225, 227, 232, 235,
C 239–242, 249, 253
Cells .............................................1, 3, 4, 11, 54, 55, 109, Ion pair chromatography..................................... 166, 182
112, 114, 131, 138, 142, 144, 181, 195, 198,
L
210, 226, 235, 241, 244, 245, 257
Chiral analysis ....................................................... 205–207 Lipidomics .......................................................6, 112, 127,
Collision cross section (CCS)............................. 8, 11, 26, 131–149, 221–236
119, 232, 233, 235, 240, 242–244, 247–253 Lipoproteins ................................................ 138, 143, 144
Liquid chromatography (LC).........................1, 110, 134,
D 156, 167–173, 183–188
Liquid chromatography-mass spectrometry
Data qualities....................................31, 61–63, 111, 115,
146, 175, 178, 179, 234 (LC-MS) ....................................... 3, 4, 10, 19, 25,
47, 53–65, 91–106, 109–128, 157, 160,
E 165–179, 194, 195, 198, 240–242, 270, 271
Liquid chromatography-tandem mass spectrometry
Extracellular vesicles............................................. 139, 144 (LC-MS/MS) .......................................... 110, 155,
157, 158, 160, 161, 187, 258, 259, 266
F
Liquid-liquid extraction.............................. 139, 140, 145
Fecal extracts ..............................155, 156, 158, 159, 161
Fecal samples ..................................................39, 154–161 M
Food analysis ................................................................. 277 Mass spectra handling ........................................ 68, 75–77
Formula handling............................................................ 69 Mass spectrometry (MS)..................................... 2, 26, 29,
39, 53, 55, 75–77, 79, 83, 85, 87, 109, 121,
G 131–149, 165, 174, 181, 239
Gas chromatography-mass spectrometry Metabolic fingerprints................................................... 270
(GC-MS)................................ 6, 10, 31, 109, 156, Metabolic phenotyping....................................... 1–4, 6, 8,
157, 159–161, 205–218, 258, 270 11, 109–128
Glycolipids ..................................................................... 135 Metabolic profiling..................................... 2, 4, 8, 15–49,
Gut microbiome................................................... 154, 156 54, 56, 62, 109–112, 119, 122–124, 155, 156,
Gut microbiota.............................................................. 154 166, 249
Olga Deda et al. (eds.), Metabolic Profiling: Methods and Protocols, Methods in Molecular Biology, vol. 2891,
https://doi.org/10.1007/978-1-0716-4334-1,
© The Editor(s) (if applicable) and The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer
Nature 2025
285
METABOLIC PROFILING: METHODS AND PROTOCOLS
286 Index
Metabolite annotation ....................................67, 68, 103, Quality management (QM).................. 15–18, 21, 27, 38
240, 249–251, 253 Quantification .............................................. 3, 12, 16, 20,
Metabolite profiling .................................. 1–11, 181, 218 38, 54, 60, 110, 114, 132, 134, 137, 138, 141,
Metabolites ............................................. 1, 19, 53, 67–87, 143, 147, 166, 181, 182, 187, 196, 202, 208,
94, 109, 142, 154, 165, 181, 214, 240, 257 214, 215, 217, 264, 266, 270–272, 277
Metabolome ......................................................3, 6, 8, 11, Quantitative analysis ................................... 141, 196, 235
109–111, 114, 142, 154, 165, 240, 269–276
Metabolomics ................................................1, 2, 6, 9–12, R
15–49, 53–65, 67–87, 103, 104, 123, 125, 142,
Rats ...................................................................7, 160, 161
154–156, 181–202, 218, 240, 245, 246, 252, Retention time correction .............................93–101, 106
253, 269–272 RforMassSpectrometry ................................................... 68
Metabonomics..........................1–3, 9, 11, 165, 166, 181
Robust method ........................................... 15, 16, 18, 45
Metadata ..................................................... 32, 76, 93–95,
104, 105, 245, 252 S
Metanephrines (MNs) ........................ 257–261, 263–266
Missing values............................................ 31, 35, 99, 252 Sample preparation ...........................................3, 6, 7, 19,
24, 30, 44, 58–60, 64, 111, 114, 116, 117, 126,
N 138–140, 143–145, 154–156, 158–160,
172–174, 181, 182, 194, 195, 198, 207, 211,
NMR spectroscopy................................... 3, 4, 10, 11, 53, 213, 214, 218, 221, 241, 245–246, 252, 258, 279
109, 155–158, 160, 165, 166, 277–281 Secondary AAs ..................................................... 205–218
Nuclear magnetic resonance (NMR) ................. 2–4, 155,
Serum............................................. 1, 6, 7, 31, 55–59, 63,
181, 269–276 109, 112, 114–116, 119, 121, 122, 143, 148,
166, 182, 194, 200, 202, 206–210, 214, 217
O
Solid-state ............................................................. 277–282
Organs ................................... 1, 131, 139, 142, 144, 153 Spectra similarity calculations ...................................78–83
Sphingolipids ........................................................ 132, 145
P Supercritical fluid chromatography (SFC)............. 6, 8, 9,
Parallel accumulation serial fragmentation 131–149
(PASEF) ..........................222, 223, 225, 227, 235
T
Peak picking ......................................... 57, 62, 67, 93–97,
105, 106, 123, 143, 147, 235, 248, 252 Targeted metabolomics .................................54, 111, 182
Peptide hydrolysates...................207, 208, 210, 214, 218 Tissues.....................................................3, 4, 6, 7, 46, 54,
Phenolic profiling................................................. 277–282 55, 109–112, 114, 117–119, 122, 123, 126, 139,
Pheochromocytoma (PHEO) ............................. 257, 258 142, 144, 155, 165–179, 181, 194, 195, 210
Phospholipids .............................................. 141, 145, 149 Trapped ion mobility spectrometry (TIMS) ..... 113, 222,
Plants ......................................6, 109, 269–276, 278, 282 223, 225, 226, 235, 240
Plasma ...............................................1, 3, 6, 7, 30, 55–59, Traveling wave ion mobility ......................................... 240
63, 112, 114, 116, 119, 121, 122, 127, 138, 140, Two-phase extraction.................................. 270, 271, 274
142, 143, 166, 182, 210, 214, 217, 230, 235,
257, 258, 265 U
Polar analytes........................................................ 111, 182 Untargeted metabolic profiling........................... 110, 111
Polar metabolites.............................................6, 111, 114,
Untargeted metabolomics ........................... 9, 16, 20, 44,
148, 166, 182, 205 53–56, 240, 245
Preprocessing........................................31, 61, 65, 67, 92, Urine......................................1, 3, 4, 6, 7, 54, 55, 57–59,
94–103, 105, 106, 122, 123, 252, 270
63, 109, 112, 114, 116, 118–121, 124, 154, 166,
182, 194–200, 202, 207–210, 214, 217, 218,
Q
258, 259, 262, 263, 265, 266
Quality assurance (QA).............................. 15, 16, 18, 19,
21–43, 45, 49, 54, 55 V
Quality checks ............................93, 95, 96, 99, 101, 234 Validation.............................................. 20, 37, 44, 53–65,
Quality control (QC)......................................8, 9, 15, 18,
167, 174, 175, 177, 252, 264
22, 29–31, 35, 36, 44, 46, 48, 53–65, 96, 105, Vitis ................................................................................ 239
111, 114–120, 123–125, 127, 128, 147, 149,
174–176, 178, 179, 195, 214, 215, 217, 218,
229, 245–248, 251–253, 262