Next Article in Journal
A Case of Suspected Hyperphenylalaninemia at Newborn Screening by Tandem Mass Spectrometry during Total Parenteral Nutrition
Next Article in Special Issue
General Guidelines for Sample Preparation Strategies in HR-µMAS NMR-based Metabolomics of Microscopic Specimens
Previous Article in Journal
Dexamethasone-Induced Perturbations in Tissue Metabolomics Revealed by Chemical Isotope Labeling LC-MS Analysis
Previous Article in Special Issue
Metabolomics in the Context of Plant Natural Products Research: From Sample Preparation to Metabolite Analysis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 10
Issue 2
10.3390/metabo10020043
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Current Practice in Untargeted Human Milk Metabolomics

by
Isabel Ten-Doménech
1,†,
Victoria Ramos-Garcia
1,†,
José David Piñeiro-Ramos
1,
María Gormaz
1,2,
Anna Parra-Llorca
1,
Máximo Vento
1,2,
Julia Kuligowski
1,* and
Guillermo Quintás
3,4
1
Neonatal Research Unit, Health Research Institute La Fe, Avenida Fernando Abril Martorell 106, 46026 Valencia, Spain
2
Division of Neonatology, University & Polytechnic Hospital La Fe, Avenida Fernando Abril Martorell 106, 46026 Valencia, Spain
3
Health and Biomedicine, Leitat Technological Center, Carrer de la Innovació, 2, 08225 Terrassa, Spain
4
Unidad Analítica, Health Research Institute La Fe, Avenida Fernando Abril Martorell 106, 46026 Valencia, Spain
*
Author to whom correspondence should be addressed.
†
Both authors contributed equally.
Metabolites 2020, 10(2), 43; https://doi.org/10.3390/metabo10020043
Submission received: 29 November 2019 / Revised: 10 January 2020 / Accepted: 19 January 2020 / Published: 22 January 2020
(This article belongs to the Special Issue Sample Preparation in Metabolomics)
Browse Figures
Versions Notes

Abstract

:
Human milk (HM) is considered the gold standard for infant nutrition. HM contains macro- and micronutrients, as well as a range of bioactive compounds (hormones, growth factors, cell debris, etc.). The analysis of the complex and dynamic composition of HM has been a permanent challenge for researchers. The use of novel, cutting-edge techniques involving different metabolomics platforms has permitted to expand knowledge on the variable composition of HM. This review aims to present the state-of-the-art in untargeted metabolomic studies of HM, with emphasis on sampling, extraction and analysis steps. Workflows available from the literature have been critically revised and compared, including a comprehensive assessment of the achievable metabolome coverage. Based on the scientific evidence available, recommendations for future untargeted HM metabolomics studies are included.

1. Introduction

Human milk (HM) has been markedly established as the optimal way of providing infants with the necessary nutrients and bioactive factors for their early development. Many health associations and organisms, including World Health Organization, recommend exclusive breastfeeding for the first six months of life [1]. Health benefits of HM for infants include reduced mortality and morbidity, including sepsis, respiratory diseases, otitis media, gastroenteritis, and urinary tract infections, among others [2]. In addition, studies reporting on long-term benefits of HM consumption such as lower risk of suffering from type 1 diabetes and inflammatory bowel disease or overweight in adulthood emerged [3]. HM may also be associated with a slightly improved neurological outcome as cohort studies report [4], especially in preterm infants [5], although potential confounders must be accounted for [6].
HM composition is dynamic and influenced by several factors including genetics, gestational and infant’s age, circadian rhythm, maternal nutrition, or ethnicity. It provides a series of nutrients such as lipids, proteins, carbohydrates, and vitamins, jointly with a number of bioactive factors that contribute to several physiological activities in the newborn infant as well as to short- and long-term outcomes [7,8]. Living cells including stem cells, hormones, growth factors, enzymes, microbiota, and even genetic material are part of this vast array of HM components with impact in early development, particularly the immune system [9]. In addition, HM appears to be one of the richest sources of microRNAs [10]. On the other hand, because of the maternal environmental exposure and lifestyle, the presence of some contaminants such as persistent organic pollutants or pharmacologically active substances in HM has been described [11,12].
Due to its complex composition, the analysis of HM is not straightforward. While the advent of “omics” approaches has offered valuable insights into the composition of this unique biofluid, untargeted metabolomic and lipidomic studies have only recently been applied to HM [13]. The comprehensive study of the HM metabolome, which includes the intermediate and end products of metabolism, can shed light on maternal status or phenotype [14,15]. The generation, analysis, and integration of large and complex data sets obtained in metabolomic studies go hand in hand with the following challenges: (i) the intrinsic complexity of the sample: a rich variety of jointly present, structurally heterogeneous compounds at concentrations that strongly vary covering several orders of magnitude; (ii) pre-analytical steps related to sampling, storage, and pre-processing (e.g., extraction, clean-up); and (iii) the diversity of platforms currently available including nuclear magnetic resonance (NMR), as well as gas chromatography (GC), liquid chromatography (LC), and capillary electrophoresis (CE) coupled to mass spectrometry (MS). The analysis of the HM metabolome has been approached employing a variety of extraction and analytical techniques to respond to a spectrum of clinically relevant questions. Several studies have compared HM metabolome with formula milk [13,16,17,18,19,20] or with milk from other mammalian species including monkey [21], donkey [17], and cow [18], whereas others have made efforts in defining the metabolome of preterm milk [13,16,22,23,24,25,26] and the evaluation of the HM metabolome during the course of lactation [15,23,27,28,29,30]. Furthermore, the influence of maternal diet [14,15,31], phenotype [14,32], obesity [30], or atopy status [33], as well as geographical location [33,34], time of the day [29,35], chemotherapy [36], or preeclampsia during pregnancy [31] on the HM metabolome have been reported.
Recent review articles that address the HM metabolome or lipidome [12,37,38,39,40,41] are available. For information on the compounds and compound families typically found in HM and their function the reader is referred to [37,38,39,40]. Readers with a particular interest in HM lipidomics are referred to a recent compilation study [41]. Technical aspects of HM analysis when performing metabolomics studies in HM have been recently described [12]. This review article gathers recent literature available on metabolomic analysis of HM, particularly focusing on untargeted approaches as indicated in Figure 1, to provide an up-to-date overview of the key factors that may influence HM metabolome coverage. Based on the information provided within the available literature, recommendations to guide study design and analytical method development of untargeted HM metabolomics assays were developed.

2. Considerations Regarding the Study Design

HM is a biofluid characterized by a dynamically varying composition according to several factors including lactation time, time of the day, throughout each feed, maternal status, and the environmental exposure. Although compositional variations have been mainly studied regarding the protein content of HM [42], changes of other compound classes such as fat or vitamins have been also reported [43,44]. Considering the intrinsic variability of HM, the complexity of obtaining representative HM samples is not negligible. Sources of variation related to sample manipulation and compositional variation can be minimized using standard operational procedures (SOPs). SOPs are fundamental to maintain quality assurance (QA) and quality control (QC) process and facilitate repeatable and reproducible research within and across laboratories. However, biologically meaningful results across studies will only be obtained if several key factors during the sample collection process are successfully controlled. This is of special importance in untargeted approaches, where the interpretation of results is especially challenging, and confounding factors introduced by a non-exhaustive sampling protocol can be wrongly attributed to differences between subjects of a studied population. Conversely, biologically meaningful information can be missed or remain unnoticed due to unwanted bias introduced during sample collection.
The information regarding study design provided in HM metabolomics studies varies considerably [13,14,15,16,17,18,19,20,21,22,23,24,26,27,28,29,30,31,32,34,35,36,45,46,47] as shown in Figure 2. Repeatedly reported factors have been grouped into three categories and are discussed in detail in the following sections: (1) maternal-infant-related factors (blue bars), (2) time-related factors (green bars), and (3) HM collection-related factors (orange bars). It should be noted that, although the importance of each factor might vary with the scientific question of each study, the authors encourage (i) the use of SOPs employed during sample collection to assure homogenous and representative sampling and (ii) the reporting of all documented factors in order to enhance comparability between results of metabolomic studies on HM. In case of HM, samples are typically collected, handled and sometimes temporary stored and transported by the mothers and not, such as it is the case for other biofluids (e.g., plasma or serum), by health professionals. During study design it is therefore very important to assure that mothers receive detailed instructions and/or training for the correct handling of collected samples. In addition, one should keep in mind that sampling protocols should neither interfere with infant feeding nor negatively impact on the mother-baby bonding. Hence, the collection of transitional and mature milk is usually preferred over colostrum, especially in studies involving mothers of preterm infants, where colostrum is usually kept exclusively for the infant’s supply.

2.1. Maternal-Infant-Related Factors

In HM metabolomic studies, gestational age is frequently reported (see Figure 2), although the impact of this factor on the HM metabolome has not been fully characterized. Studies focused on preterm milk showed that, analogously to full-term milk, its composition is dynamic throughout the first month of lactation [13,16,22]. However, after 5–7 weeks, metabolite composition of HM from mothers of preterm infants resembled that collected from mothers of full-term infants [23]. On the other hand, Marincola et al. [13] observed that HM from mothers of early preterm infants (26 weeks of gestation) differentiated from milk samples from term infants. However, the low number of samples involved in the study (n = 20 and n = 3 mothers of preterm and term infants, respectively) hindered the assessment of the statistical significance of the impact of gestational age on the milk metabolite composition. Sundekilde et al. [23] carried out a longitudinal study on milk from mothers of preterm and full-term infants covering similar lactation periods (3–14 weeks and 3–26 weeks after birth, respectively) and showed that some metabolites were present at significantly different levels in full-term milk compared to preterm milk. On the contrary, Longini et al. [16] did not observe significant differences between preterm and full-term milk within the first week after delivery, only being able to discriminate milk samples from early preterm infants (<29 weeks of gestation). It is worth noting that the effect of gestational age on the HM metabolome has been mainly studied employing NMR platforms [13,16,22,23], in which metabolite coverage is limited (see Figure 5) and some metabolite classes (e.g., lipids) are barely accessible. For this reason, and in order to further evaluate the impact of gestational age on the milk metabolome, we warrant more comprehensive metabolomic studies employing complementary analytical platforms.
Other potentially relevant, miscellaneous information about the studied population of mother-infant pairs such as infant gender, parity, and birth mode, have been frequently reported in metabolomic studies (see Figure 2), although these characteristics often remain in the background since the studies focus on other aspects. The influence of these factors on the metabolite composition of HM has not been elucidated yet, and this might be addressed in forthcoming studies. An additional factor that is not typically reported in HM metabolomics studies is maternal secretor status. Significant differences in the oligosaccharides profile of milk between so-called secretors (Se+), which are those mothers that provide a functional FUT2 gene, and non-secretors (Se−) have been reported [48]. Secretor status is mainly established based on the presence (Se+) or absence (Se−) of 2’-fucosyllactose, with a prevalence rate of approximately 80% of secretors over non-secretors [14,22,23,26,32,49]. Maternal secretor status is therefore usually determined a posteriori during data processing and analysis. Oligosaccharides are polar compounds that are present at concentrations in the mM range that will likely be preserved during sample extraction procedures employed for metabolomics studies. As their presence/absence might potentially affect clustering of milk based on maternal secretor status [24], to provide this information, when available, might be of interest.

2.2. Time-Related Factors

HM undergoes significant changes over time, having established three differentiated lactation stages: colostrum, transitional milk and mature milk. As can be seen in Figure 2, lactation time is reported in the vast majority of metabolomic studies. In particular, several studies have focused on the HM metabolome throughout lactation [15,23,27,28,29,30], all of them concluding that significant differences in the metabolic profile over time exist. Therefore, it seems reasonable to report this factor. On the other hand, although it has been demonstrated that diurnal variation affects HM fat content [50], its effect on the overall metabolite composition is a controversial issue which has not yet been adequately addressed in the available literature. Whereas no significant changes in some lipids and small polar metabolites have been observed [29,35], differences in some micronutrients (e.g., vitamins) could be evidenced [44]. The use of a pool of a 24-h expression of HM should compensate for changes due to diurnal variation, thus, obtaining more representative samples [13,24,25] in longitudinal studies. However, the drawback is that this practice is incompatible with breastfeeding of the infant, which, in turn, might raise severe ethical concerns. A feasible compromise for ameliorating diurnal variations is the use of pooled morning and evening samples [29,35].
Regarding time of collection with respect to baby’s feed, the influence of this variable has not been studied to date, but given the differences found between fore- and hindmilk [51], it seems reasonable to assume that this factor might be potentially relevant.

2.3. HM Collection-Related Factors

Any uncontrolled variable within an experiment can result in a potential source of bias. In this sense, although less attention has been payed to other factors related to the expression and storage of HM (see Figure 2, orange bars), they may be relevant to the outcomes of metabolomic studies. As can be seen, the type of sample container is indicated in 50% of the studies, whereas other specifications regarding HM expression are included scarcely. The latter factor deserves some special attention, since differences in the milk fat content between foremilk (initial milk of a feed) and hindmilk (last milk of a feed) have been reported [51]. Therefore, full expression of breast(s) is desirable in order to obtain a representative HM aliquot [52]. The influence of all other factors, to date, remains unstudied.

2.4. Pasteurization and Storage

HM banks rely on stringent protocols in which pasteurization, indispensable for minimizing the potential to transmit infectious agents, as well as freezing and long-term storage procedures are established. The pasteurization process affects some of the nutritional and biological properties of HM [53,54,55]. In this review, three studies that use milk from HM banks are included [16,20,23], but only one specifies whether or not HM has undergone pasteurization [23]. Variability of the metabolite profile of HM caused by pasteurization has not been comprehensively explored to date. Future studies focused on the systematic exploration of the effect of thermal treatment on HM are warranted.
HM is usually stored frozen employing −20 °C and −80 °C for short- and long-term storage, respectively. However, duration of storage and the effect of repeated freeze-thaw cycles are identified as additional factors with potential impact on HM composition that are missing in most published studies. In lipidomic studies, the integrity of HM samples is preserved by subjecting HM to inactivation of endogenous enzymes such as lipases in order to minimize lipolysis and lipogenesis. In this sense, immediate storage at −80 °C is advisable [41]. Particularly for metabolite composition analysis, storage at −80 °C is widespread [13,15,18,19,24,25,26,28,31,45,46], sometimes with a prior short-term storage at −20 °C [14,21,22,27,29,33,34]. Wu et al. [29] investigated the effect of storage conditions by keeping samples for different times at −20 °C and then transferring them to −80 °C versus storing samples directly at −80 °C. Variations in duration of storage at −20 °C versus −80 °C showed no detectable effect on the metabolites considered (e.g., lactose and other carbohydrates, choline and its derivatives, and a variety of amino acids) by visual inspection of sample clusters in principal component analysis scores plots. However, analysis of variance evidenced differences in butyrate, caprate, and acetate contents. However, time of storage considered in this study was limited to two weeks, which is not representative for conditions employed in clinical studies or standard home routines. It is therefore clear that further studies are required in this regard.
On the other hand, HM employed for research studies is typically stored in small aliquots. When working with raw milk, this procedure might introduce bias due to phase separation prior to the preparation of the aliquots. Hence, the homogenization of HM with a disruptor, resulting in a stable emulsion with reduced size of milk fat globules [56], as employed prior to the quantitation of macronutrients with HM analyzers, might be advisable.

3. Metabolite Extraction from HM

For metabolite extraction from HM, an array of methods has been reported. An overview of the employed approaches is shown in Figure 3. The selection of the extraction method is conditioned by the study objective and the subsequent analysis method. As in other untargeted metabolomics workflows, for HM metabolomics, the selected sample preparation approach should enable a high degree of metabolome coverage while making the sample matrix compatible with the analytical platform. Other considerations might include the available amount of sample volume and the use of one sample extraction procedure for subsequent analysis by multiple, complementary analytical platforms [13,27,28].
Liquid-liquid extraction (LLE) is the classical extraction method employed in metabolomics and lipidomics. This method, developed by Folch et al. [57] in 1957, uses a chloroform-methanol mixture (2:1, v/v), which results in two differentiate phases: an upper phase containing polar metabolites and a lower phase containing nonpolar metabolites. Subsequently, in 1959 Bligh and Dyer [58] developed a modified method using a miscible chloroform-methanol-water mixture and later separated into two phases by adding chloroform or water. Both approaches enable the separation of polar and nonpolar metabolites, thus, allowing the analysis of a wide range of metabolites and making them compatible with several analytical platforms. While the use of Bligh and Dyer LLE is widely extended for HM metabolomics studies (see Table 1) [13,16,17,18,19,24,25,29,32], only Andreas et al. [28] used a modified Folch extraction protocol for processing HM samples.
Methyl tert-butyl ether (MTBE) in combination with methanol has recently been proposed for single-phase extraction [27]. MTBE is a nontoxic and noncarcinogenic solvent and it is therefore considered a safe and environmentally friendly alternative to harmful solvents employed in traditional LLE methods, such as chloroform, which is a suspected human carcinogen. In this extraction method, a unique phase containing both, polar and nonpolar metabolites is obtained with a protein pellet at the bottom (see Figure 3). Thus, the simultaneous analysis of lipidome and metabolome in a very small amount of biological sample is achievable. This method has been successfully employed to determine polar metabolites and fatty acids (FAs) in HM by GC-MS [27,28], as well as lipids and polar metabolites by LC-MS [15,27,28], thus, increasing the metabolome coverage by the combined use of complementary analytical platforms.
Ultrafiltration makes use of centrifugal molecular weight cutoff filters. Different molecular weight cut-off filters are commercially available for this purpose and repeated centrifugation steps might be employed to remove proteins and lipids (see Table 1). Unlike single-phase extraction, ultrafiltration allows to separate polar metabolites from the HM without dilution [14,21,22,29], however this method does not have the capacity to study the global metabolome of HM. At present, this extraction method has only been used in combination with NMR analyses [14,21,22,29].
Precipitation with organic solvents separates the polar and nonpolar metabolites of the proteins that settle at the bottom of the tube which can then be easily removed by centrifugation. This simple method has been employed for the analysis of polar metabolites by GC-MS after derivatization [36] as well as for the analysis of polar and nonpolar metabolites by LC-MS without further pre-processing [18]. Furthermore, this approach has been implemented in more sophisticated workflows as recently shown by Hewelt-Belka et al. [35]. Here, the authors combined LLE and a protein precipitation and solid-phase extraction (SPE) procedure to prepare HM samples, thereby, enabling the detection of high- and low-abundant lipid species (e.g., glycerolipids and phospholipids) in one LC-MS run.

4. Analytical Platforms Employed in HM Metabolomics

As reflected in Table 1, the use of all analytical platforms that are commonly employed in untargeted metabolomics studies, such as LC-MS, GC-MS, NMR, and, to a lesser extent, CE-MS, has been reported for performing HM metabolomics. 1H-NMR is the most frequently used technique [13,14,16,20,21,22,26,28,29,31,32,33] for both, the analysis of polar and hydrophobic metabolites in HM. 13C-NMR has been reported for the detection of triacylglycerols [20]. NMR is a highly reproducible technique that allows a straightforward library matching after spectral alignment, while at the same time supporting structural elucidation of detected metabolites. However, it presents lower sensitivity, and hence, the achievable metabolome coverage is low in comparison to other analytical platforms.
All other analytical platforms rely on the use of MS detection. LC-MS provides high sensitivity and is characterized by a huge versatility due to the availability of (i) a large selection of chromatographic columns with a variety of stationary phases, that in combination with appropriate mobile phases achieve compound separation based on different retention mechanisms and (ii) an array of different instrumental configurations (i.e. different ion sources and mass analyzers). For example, reversed phase (C8, C18) LC-quadrupole time of flight MS (LC-QTOF-MS) [15,18,27,28,35,46] and LC-Orbitrap-MS [24,25] have been reported for the detection of both, polar and lipidic metabolites in HM; and hydrophilic interaction LC (HILIC)-QTOF-MS for polar metabolite detection [35]. On the other hand, for the successful screening of HM oligosaccharides, a porous graphitic carbon column installed on a LC-triple quadrupole (TQD)-MS instrument was used [24].
GC-MS is the most suitable platform for measuring volatile compounds, while other non-volatile compounds must be derivatized prior to analysis. For HM analysis, methoximation followed by silylation or methylation are commonly employed (see Table 1). The most frequently used column is the DB-5ms column for both polar and FA detection [13,17,18,19], in some cases with an integrated 10 m pre-column (deactivated fused silica) [27,28,36].
Regarding CE-MS, only one study has been reported for polar metabolite detection in HM [28]. CE provides a series of advantages over other techniques, mainly due to the small sample volumes employed and the efficient separation of polar compounds that is difficult to be achieved by LC columns. However, issues with poor reproducibility, matrix effects and sensitivity may be hindering a widely extended use of this technique for the analysis of complex biological samples such as HM.
Due to the diversified composition of HM, no single analytic technique can resolve the entire HM metabolome. Only multiplatform approaches enable a comprehensive characterization providing a high metabolome coverage including polar and nonpolar metabolites present in HM. In this sense, two studies performing a multiplatform approach were found in the literature combining LC-MS and GC-MS [18,27], and only one study that performed analysis using four different techniques (LC-MS, GC-MS, NMR, and CE-MS) was reported [28].
The use of high-end analytical platforms requires the implementation of QA and QC processes to improve data quality, repeatability, and reproducibility, especially in untargeted metabolomics. For practical guidelines on the use of QC measures in untargeted, MS-based assays the reader is referred to [59]. Pooled QC samples are prepared by mixing small aliquots of the study samples, and therefore, they are considered representative in terms of matrix composition and concentration ranges of the metabolites present in the study samples. QC samples are analyzed repeatedly throughout the analytical sequence alongside the study samples. The signal of each feature detected in QC samples can be used to model and correct systematic changes in the instrument response during the analytical sequence. Additionally, the obtained data can be used to perform intra-study reproducibility assessments and to correct for systematic variation across batches. In HM metabolomics, Smilowitz et al. [14], Andreas et al. [28], and Gay et al. [33] used QC samples for NMR studies, while Villaseñor et al. [27], Mung et al. [46], Hewelt-Belka et al. [35], and Alexandre-Gouabau et al. [24,25] used pooled HM samples for QC purposes in LC-MS-based assays. Considering the highly complex sample matrix of HM, the authors strongly recommend the implementation of QC measures, including the analysis of QC samples, to increase reproducibility and facilitate the joint analysis of data from different studies.

5. The HM Metabolome: Compound Annotation and Coverage

As in other areas of metabolomic research, compound identification is still a major bottleneck in data analysis and interpretation. The Metabolomics Standards Initiative’s (MSI) defines four levels of metabolite identification, which include: identified metabolites (level 1); putatively annotated compounds (level 2); putatively annotated compound classes (level 3); and unknown compounds (level 4) [60]. Due to the limited availability of pure analytical standards required to reach level 1, biological databanks and spectral databases are the most important resources for metabolite annotation (levels 2 and 3). A large number of databases are available today, providing different levels of information and complementary data on chemical structures, physicochemical properties, biological functions, and pathway mapping of metabolites [61]. The metabolomics community classifies these resources in several categories: (i) chemical databases; (ii) spectral libraries; (iii) pathway databases; (iv) knowledge databases; and (v) references repositories [62].
Regarding HM metabolomics, the most frequently used databases and libraries are: Human Metabolome Database (HMDB) [63], Metabolite and Chemical Entity Database (METLIN) [64], National Institute of Science and Technology (NIST) library, Fiehn RTL Library [65], LipidMAPS Structure Database (LMSD) [66], Milk Metabolome Database (MCDB) [67,68], Kyoto Encyclopedia of Genes and Genomes (KEGG) [69], MycompoundID with the evidence-based metabolome library (EML) [70], Chenomx NMR Suite Profiles and other online university databases, such as CEU-mass mediator [71,72].
Metabolite assignment in NMR spectra has been performed based on literature data and commercial resonance databases, such as Chenomx NMR Suite Profiles. Metabolite annotation was contrasted with in-house libraries containing pure compound spectra. Some of the proposed assignments were confirmed by two-dimensional NMR spectra, such as Correlation Spectroscopy (COSY) [13,29,31,32], Homonuclear Correlation Spectroscopy (TOCSY) [13,31,32,34], Diffusion-Ordered Spectroscopy (DOSY) [32], Heteronuclear Single Quantum Coherence Spectroscopy (HSQC) [32,34], and Heteronuclear Multiple Bond Correlation (HMBC) [32].
In LC-MS and CE-MS-based studies of the HM metabolome, tentative metabolite annotation has been carried out by matching of accurate masses, isotopic profiles, and/or fragmentation patterns to candidate metabolites in online databases such as KEGG, METLIN, LipidMAPS, and HMDB [18,24,25,27,28,35]. In-house built databases generated by the analysis of commercial standards are also commonly employed [24,25]. In GC-MS, retention index (RI) corrections are made by analyzing a fatty acid methyl ester (FAME) mixture standard solution and assigning a match score between the experimental FAME mixture and theoretical RI values based on the values contained in the Fiehn RTL library. Furthermore, metabolites were complementarily annotated by comparing their mass fragmentation patterns with those available in Fiehn RTL and NIST libraries [13,17,18,19,27,28,36].
A comprehensive list of annotated and/or identified metabolites in HM from untargeted metabolomics studies [14,15,17,18,19,21,22,23,24,25,26,27,28,29,31,32,33,34,35,36] is reported in Table S1. This table contains information about the metabolites reported in each reference, such as their molecular formula, IDs (LipidMAPS and/or HMDB IDs), the extraction procedure performed, the analytical platform used, and the detected metabolite class. Readers can select metabolites dynamically by filtering data according to the latter information. A total of 1187, 111, and 128 metabolites were reported using LC-MS, GC-MS, and NMR, respectively (see Figure 4). As shown in the Venn diagram, LC-MS and GC-MS allowed the detection of 36 common metabolites (mainly carbohydrates and FAs); a total of 29 metabolites overlapped between LC-MS and NMR (principally oligosaccharides); and 21 metabolites (predominantly amino acids and organic acids) were commonly reported in GC-MS and NMR based studies. Only 13 metabolites were reported by all three platforms, i.e., creatine, tyrosine, arabinose, galactose, glucose, lactose, maltose, capric acid/caprate, caprylic acid/ caprylate, citric acid/citrate, pyruvic acid/pyruvate, hippuric acid/hippurate, and myo-inositol. These metabolites were assigned to different classes including amino acids, carbohydrates, FAs, and organic acids.
Based on the available data from the literature, the distribution of metabolite classes present in HM according to each technique was assessed. As can be seen in Figure 5, the difference in detected metabolite classes as observed by LC-MS in comparison to GC-MS and NMR is evident. Using GC-MS and NMR, carbohydrates are the most reported metabolites in HM, followed by amino acids, organic acids, organooxygen compounds, and organoheterocyclic compounds, with all these metabolite classes being certainly less abundant in LC-MS studies. In the case of NMR, organonitrogen compounds have also been reported, as well as nucleosides and nucleotides on a smaller scale. In the case of lipid classes, fatty acyls have been identified by LC-MS and GC-MS with similar incidence and in lesser extent by NMR. It is indubitable that lipid classes are more comprehensively studied by LC-MS assays, where glycerophospholipids, glycerolipids, and fatty acyls are detected at relatively high abundances, followed by sphingolipids, sterol lipids, and, to a lesser extent, prenol lipids.
Table 2 shows a list of metabolites reported in > 80% of studies employing either LC-MS, GC-MS, or NMR-based assays. This table is intended to aid method development of future untargeted metabolomics workflows tailored to the study of the HM metabolome, as it shows a shortlist of metabolites that should be detected by each platform regardless of the instrumental settings employed. It should be noted that due to the high versatility of LC-MS, there is a greater variation in metabolites recorded and in return, the list of consistently reported metabolites in HM across studies is shorter than for NMR and GC-MS, where differences in experimental conditions and variations between the employed detection parameters and instruments are smaller. Again, this table represents the high orthogonality between the detected metabolites using NMR and LC-MS. While the use of LC-MS is clearly of advantage for the measurement of different lipids, NMR provides information on amino acids and small organic acids. Metabolome coverage provided by GC-MS falls in-between the other two platforms, consistently providing information on lipids, sugars, amino acids, and organic acids.

6. Conclusions and Future Perspectives

In less than a decade, 26 research papers have been published trying to shed light on the complex and dynamic composition of HM and the feasibility of different options for sample extraction and metabolite detection has been demonstrated. Due to the many factors that influence HM composition, a thorough study design including SOPs for milk extraction, collection, and storage is indispensable for obtaining biologically meaningful results. Multi-platform approaches are encouraged for providing adequate metabolome coverage, as the diversity of compounds contained in HM will not be properly reflected using one single assay. In line with metabolomics workflows tailored to other sample types, the reproducibility of HM metabolomics studies will benefit from the implementation of QA/QC procedures. Automated metabolite annotation and identification with pure chemical standards is warranted and the authors encourage the use of publicly accessible platforms for enabling the exchange of raw data for comparison between studies.

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-1989/10/2/43/s1, Table S1: List of metabolites annotated and/or identified in HM metabolomic studies.

Author Contributions

Conceptualization, G.Q., J.K., M.G., and M.V.; methodology, A.P.-L., I.T.-D., J.D.P.-R. and V.R.-G.; software, G.Q., I.T.-D. and J.K., and validation, G.Q. and J.K.; formal analysis, I.T.-D. and V.R.-G.; investigation, I.T.-D. and V.R.-G.; resources, J.K.; data curation, I.T.-D. and V.R.-G.; writing—original draft preparation, A.P.-L., I.T.-D., J.D.P.-R. and V.R.-G.; writing—review and editing, all authors; visualization, I.T.-D. and V.R.-G.; supervision, J.K. and M.V; project administration, G.Q. and J.K.; funding acquisition, J.K. All authors have read and agreed to the published version of the manuscript.

Funding

A.P.-L., I.T.-D. and J.K. received salary support by the Instituto de Salud Carlos III (Ministry of Economy and Competitiveness, Spain), grant numbers CD19/00176, CM18/00165, and CP16/00034, respectively. This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 818110.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization Breastfeeding. Available online: https://www.who.int/topics/breastfeeding/en/ (accessed on 2 July 2019).
  2. Geddes, D.; Perrella, S. Breastfeeding and human lactation. Nutrients 2019, 11, 802–806. [Google Scholar] [CrossRef] [Green Version]
  3. Owen, C.G.; Martin, R.M.; Whincup, P.H.; Davey Smith, G.; Cook, D.G. Does breastfeeding influence risk of type 2 diabetes in later life? A quantitative analysis of published evidence. Am. J. Clin. Nutr. 2006, 84, 1043–1054. [Google Scholar] [CrossRef] [Green Version]
  4. Der, G.; Batty, G.D.; Deary, I.J. Effect of breast feeding on intelligence in children: Prospective study, sibling pairs analysis, and meta-analysis. Br. Med. J. 2006, 333, 945–948. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Rozé, J.C.; Darmaun, D.; Boquien, C.Y.; Flamant, C.; Picaud, J.C.; Savagner, C.; Claris, O.; Lapillonne, A.; Mitanchez, D.; Branger, B.; et al. The apparent breastfeeding paradox in very preterm infants: Relationship between breast feeding, early weight gain and neurodevelopment based on results from two cohorts, EPIPAGE and LIFT. BMJ Open 2012, 2, 1–9. [Google Scholar] [CrossRef] [PubMed]
  6. Horta, B.L.; Loret De Mola, C.; Victora, C.G. Breastfeeding and intelligence: A systematic review and meta-analysis. Acta Paediatr. Int. J. Paediatr. 2015, 104, 14–19. [Google Scholar] [CrossRef] [PubMed]
  7. Lönnerdal, B. Bioactive proteins in human milk: Mechanisms of action. J. Pediatr. 2010, 156, S26–S30. [Google Scholar] [CrossRef]
  8. Musilova, S.; Rada, V.; Vlkova, E.; Bunesova, V. Beneficial effects of human milk oligosaccharides on gut microbiota. Benef. Microbes 2014, 5, 273–283. [Google Scholar] [CrossRef]
  9. Ballard, O.; Morrow, A.L. Human milk composition: Nutrients and bioactive factors. Pediatr. Clin. N. Am. 2013, 60, 49–74. [Google Scholar] [CrossRef] [Green Version]
  10. Alsaweed, M.; Hartmann, P.E.; Geddes, D.T.; Kakulas, F. MicroRNAs in Breastmilk and the Lactating Breast: Potential Immunoprotectors and Developmental Regulators for the Infant and the Mother. Int. J. Environ. Res. Public Health 2015, 12, 13981–14020. [Google Scholar] [CrossRef] [Green Version]
  11. Van den Berg, M.; Kypke, K.; Kotz, A.; Tritscher, A.; Lee, S.Y.; Magulova, K.; Fiedler, H.; Malisch, R. WHO/UNEP global surveys of PCDDs, PCDFs, PCBs and DDTs in human milk and benefit–risk evaluation of breastfeeding. Arch. Toxicol. 2017, 91, 83–96. [Google Scholar] [CrossRef] [Green Version]
  12. Garwolińska, D.; Namieśnik, J.; Kot-Wasik, A.; Hewelt-Belka, W. State of the art in sample preparation for human breast milk metabolomics—Merits and limitations. TrAC Trends Anal. Chem. 2019, 114, 1–10. [Google Scholar] [CrossRef]
  13. Marincola, F.C.; Noto, A.; Caboni, P.; Reali, A.; Barberini, L.; Lussu, M.; Murgia, F.; Santoru, M.L.; Atzori, L.; Fanos, V. A metabolomic study of preterm human and formula milk by high resolution NMR and GC/MS analysis: Preliminary results. J. Matern.-Fetal Neonatal Med. 2012, 25, 62–67. [Google Scholar] [CrossRef] [PubMed]
  14. Smilowitz, J.T.; Sullivan, A.O.Õ.; Barile, D.; German, J.B.; Lo, B. The human milk metabolome reveals diverse oligosaccharide profiles. J. Nutr. 2013, 143, 1709–1718. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Li, K.; Jiang, J.; Xiao, H.; Wu, K.; Qi, C.; Sund, J.; Li, D. Changes in metabolites profile of breast milk over lactation stages and their relationship with dietary intake in Chinese: HPLC-QTOFMS based metabolomic analysis. Food Funct. 2018, 9, 5189–5197. [Google Scholar] [CrossRef] [PubMed]
  16. Longini, M.; Tataranno, M.L.; Proietti, F.; Tortoriello, M.; Belvisi, E.; Vivi, A.; Tassini, M.; Perrone, S.; Buonocore, G. A metabolomic study of preterm and term human and formula milk by proton MRS analysis: Preliminary results. J. Matern.-Fetal Neonatal Med. 2014, 7058, 27–33. [Google Scholar] [CrossRef]
  17. Murgia, A.; Scano, P.; Contu, M.; Ibba, I.; Altea, M.; Demuru, M.; Porcu, A.; Caboni, P. Characterization of donkey milk and metabolite profile comparison with human milk and formula milk. LWT 2016, 74, 427–433. [Google Scholar] [CrossRef]
  18. Qian, L.; Zhao, A.; Zhang, Y.; Chen, T.; Zeisel, S.H.; Jia, W.; Cai, W. Metabolomic approaches to explore chemical diversity of human breast-milk, formula milk and bovine milk. Int. J. Mol. Sci. 2016, 17, 2128–2143. [Google Scholar] [CrossRef] [Green Version]
  19. Scano, P.; Murgia, A.; Demuru, M.; Consonni, R.; Caboni, P. Metabolite profiles of formula milk compared to breast milk. Food Res. Int. 2016, 87, 76–82. [Google Scholar] [CrossRef]
  20. Lopes, T.I.B.; Cañedo, M.C.; Oliveira, F.M.P.; Ancantara, G.B. Toward precision nutrition: Commercial infant formulas and human milk compared for stereospecific distribution of fatty acids using metabolomics. Omics J. Integr. Biol. 2018, 22, 484–492. [Google Scholar] [CrossRef]
  21. O’Sullivan, A.; He, X.; McNiven, E.M.S.; Hinde, K.; Haggarty, N.W.; Lönnerdal, B.; Slupsky, C.M. Metabolomic phenotyping validates the infant rhesus monkey as a model of human infant metabolism. J. Pediatr. Gastroenterol. Nutr. 2013, 56, 355–363. [Google Scholar] [CrossRef]
  22. Spevacek, A.R.; Smilowitz, J.T.; Chin, E.L.; Underwood, M.A.; German, J.B.; Slupsky, C.M. Infant maturity at birth reveals minor differences in the maternal milk metabolome in the first month of lactation. J. Nutr. 2015, 145, 1698–1708. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Sundekilde, U.K.; Downey, E.; Mahony, J.A.O.; Shea, C.O.; Ryan, C.A.; Kelly, A.L.; Bertram, H.C. The effect of gestational and lactational age on the human milk metabolome. Nutrients 2016, 8, 304–318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Alexandre-Gouabau, M.-C.; Moyon, T.; David-Sochard, A.; Fenaille, F.; Cholet, S.; Royer, A.-L.; Guitton, Y.; Billard, H.; Darmaun, D.; Rozé, J.-C.; et al. Comprehensive preterm breast milk metabotype associated with optimal infant early growth pattern. Nutrients 2019, 11, 528–553. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Alexandre-Gouabau, M.C.; Moyon, T.; Cariou, V.; Antignac, J.P.; Qannari, E.M.; Croyal, M.; Soumah, M.; Guitton, Y.; David-Sochard, A.; Billard, H.; et al. Breast milk lipidome is associated with early growth trajectory in preterm infants. Nutrients 2018, 10, 164–192. [Google Scholar] [CrossRef] [Green Version]
  26. Dessì, A.; Briana, D.; Corbu, S.; Gavrili, S.; Marincola, F.C.; Georgantzi, S.; Pintus, R.; Fanos, V.; Malamitsi-Puchner, A. Metabolomics of breast milk: The importance of phenotypes. Metabolites 2018, 8, 79–88. [Google Scholar] [CrossRef] [Green Version]
  27. Villaseñor, A.; Garcia-Perez, I.; Garcia, A.; Posma, J.M.; Fernández-Lópes, M.; Nicholas, A.J.; Modi, N.; Holmes, E.; Barbas, C. Breast milk metabolome characterization in a single-phase extraction, multiplatform analytical approach. Anal. Chem. 2014, 86, 8245–8252. [Google Scholar] [CrossRef]
  28. Andreas, N.J.; Hyde, M.J.; Gomez-romero, M.; Lopez-Gonzalvez, M.A.; Villaseñor, A.; Wijeyesekera, A.; Barbas, C.; Modi, N.; Holmes, E.; Garcia-Perez, I. Multiplatform characterization of dynamic changes in breast milk during lactation. Electrophoresis 2015, 36, 2269–2285. [Google Scholar] [CrossRef] [Green Version]
  29. Wu, J.; Domellöf, M.; Zivkovic, A.M.; Larsson, G.; Öhman, A.; Nording, M.L. NMR-based metabolite profiling of human milk: A pilot study of methods for investigating compositional changes during lactation. Biochem. Biophys. Res. Commun. 2016, 469, 626–632. [Google Scholar] [CrossRef]
  30. Isganaitis, E.; Venditti, S.; Matthews, T.J.; Lerin, C.; Demerath, E.W.; Fields, D.A. Maternal obesity and the human milk metabolome: Associations with infant body composition and postnatal weight gain. Am. J. Clin. Nutr. 2019, 110, 111–120. [Google Scholar] [CrossRef]
  31. Dangat, K.; Upadhyay, D.; Kilari, A.; Sharma, U.; Kemse, N.; Mehendale, S.; Lalwani, S.; Wagh, G.; Joshi, S.; Jagannathan, N.R. Altered breast milk components in preeclampsia; An in-vitro proton NMR spectroscopy study. Clin. Chim. Acta 2016, 463, 75–83. [Google Scholar] [CrossRef]
  32. Praticò, G.; Capuani, G.; Tomassini, A.; Baldassarre, E.; Delfini, M.; Miccheli, A. Exploring human breast milk composition by NMR-based metabolomics. Nat. Prod. Res. 2013, 28, 95–101. [Google Scholar] [CrossRef] [PubMed]
  33. Gay, M.C.L.; Koleva, P.T.; Slupsky, C.M.; Toit, E.; Eggesbo, M.; Johnson, C.C.; Wegienka, G.; Shimojo, N. Worldwide variation in human milk metabolome: Indicators of breast physiology and maternal lifestyle? Nutrients 2018, 10, 1151–1162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Gómez-Gallego, C.; Morales, J.M.; Monleón, D.; du Toit, E.; Kumar, H.; Linderborg, K.M.; Zhang, Y.; Yang, B.; Isolauri, E.; Salminen, S.; et al. Human breast milk NMR metabolomic profile across specific geographical locations and its association with the milk microbiota. Nutrients 2018, 10, 1355–1375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Hewelt-Belka, W.; Garwolińska, D.; Belka, M.; Bączek, T.; Namieśnik, J.; Kot-Wasik, A. A new dilution-enrichment sample preparation strategy for expanded metabolome monitoring of human breast milk that overcomes the simultaneous presence of low- and high-abundance lipid species. Food Chem. 2019, 288, 154–161. [Google Scholar] [CrossRef]
  36. Urbaniak, C.; Mcmillan, A.; Angelini, M.; Gloor, G.B.; Sumarah, M.; Burton, J.P.; Reid, G. Effect of chemotherapy on the microbiota and metabolome of human milk, a case report. Microbiome 2014, 2, 24–35. [Google Scholar] [CrossRef] [Green Version]
  37. Demmelmair, H.; Koletzko, B. Variation of metabolite and hormone contents in human milk. Clin. Perinatol. 2017, 44, 151–164. [Google Scholar] [CrossRef]
  38. Marincola, F.C.; Dessì, A.; Corbu, S.; Reali, A.; Fanos, V. Clinical impact of human breast milk metabolomics. Clin. Chim. Acta 2015, 451, 103–106. [Google Scholar] [CrossRef]
  39. Slupsky, C.M. Metabolomics in human milk research. Nestle Nutr. Inst. Workshop Ser. 2019, 90, 179–190. [Google Scholar]
  40. Bardanzellu, F.; Fanos, V.; Reali, A. “Omics” in human colostrum and mature milk: Looking to old data with new eyes. Nutrients 2017, 9, 843–867. [Google Scholar] [CrossRef] [Green Version]
  41. George, A.D.; Gay, M.C.L.; Trengove, R.D.; Geddes, D.T. Human milk lipidomics: Current techniques and methodologies. Nutrients 2018, 10, 1169–1179. [Google Scholar] [CrossRef]
  42. Lönnerdal, B.; Erdmann, P.; Thakkar, S.K.; Sauser, J.; Destaillats, F. Longitudinal evolution of true protein, amino acids and bioactive proteins in breast milk: A developmental perspective. J. Nutr. Biochem. 2017, 41, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Gidrewicz, D.A.; Fenton, T.R. A systematic review and meta-analysis of the nutrient content of preterm and term breast milk. BMC Pediatr. 2014, 14, 216–230. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Hampel, D.; Shahab-Ferdows, S.; Islam, M.M.; Peerson, J.M.; Allen, L.H. Vitamin concentrations in human milk vary with time within feed, circadian rhythm, and single-dose supplementation. J. Nutr. 2017, 147, 603–611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Mung, D.; Li, L. Development of chemical isotope labeling LC-MS for milk metabolomics: Comprehensive and quantitative profiling of the amine/phenol submetabolome. Anal. Chem. 2017, 89, 4435–4443. [Google Scholar] [CrossRef] [PubMed]
  46. Mung, D.; Li, L. Applying quantitative metabolomics based on chemical isotope labeling LC-MS for detecting potential milk adulterant in human milk. Anal. Chim. Acta 2018, 1001, 78–85. [Google Scholar] [CrossRef] [PubMed]
  47. Garcia, C.; Duan, R.D.; Brévaut-Malaty, V.; Gire, C.; Millet, V.; Simeoni, U.; Bernard, M.; Armand, M. Bioactive compounds in human milk and intestinal health and maturity in preterm newborn: An overview. Cell. Mol. Biol. 2013, 59, 108–131. [Google Scholar]
  48. Garwolińska, D.; Namieśnik, J.; Kot-Wasik, A.; Hewelt-Belka, W. Chemistry of human breast milk—A comprehensive review of the composition and role of milk metabolites in child development. J. Agric. Food Chem. 2018, 66, 11881–11896. [Google Scholar] [CrossRef]
  49. Totten, S.M.; Zivkovic, A.M.; Wu, S.; Ngyuen, U.; Freeman, S.L.; Ruhaak, L.R.; Darboe, M.K.; German, J.B.; Prentice, A.M.; Lebrilla, C.B. Comprehensive profiles of human milk oligosaccharides yield highly sensitive and specific markers for determining secretor status in lactating mothers. J. Proteome Res. 2012, 11, 6124–6133. [Google Scholar] [CrossRef]
  50. Lubetzky, R.; Littner, Y.; Mimouni, F.B.; Dollberg, S.; Mandel, D.; Lubetzky, R.; Lubetzky, R.; Littner, Y.; Mimouni, F.B.; Dollberg, S.; et al. Circadian variations in fat content of expressed breast milk from mothers of preterm infants. J. Am. Coll. Nutr. 2006, 25, 151–154. [Google Scholar] [CrossRef]
  51. Saarela, T.; Kokkonen, J.; Koivisto, M. Macronutrient and energy contents of human milk fractions during the first six months of lactation. Acta Paediatr. 2005, 94, 1176–1181. [Google Scholar] [CrossRef]
  52. Jensen, R.G.; Lammi-Keefe, C.J.; Koletzko, B. Representative sampling of human milk and the extraction of fat for analysis of environmental lipophilic contaminants. Toxicol. Environ. Chem. 1997, 62, 229–247. [Google Scholar] [CrossRef]
  53. Bertino, E.; Peila, C.; Cresi, F.; Maggiora, E.; Sottemano, S.; Gazzolo, D.; Arslanoglu, S.; Coscia, A. Donor human milk: Effects of storage and heat treatment on oxidative stress markers. Front. Pediatr. 2018, 6, 1–5. [Google Scholar] [CrossRef] [PubMed]
  54. Peila, C.; Moro, G.E.; Bertino, E.; Cavallarin, L.; Giribaldi, M.; Giuliani, F.; Cresi, F.; Coscia, A. The effect of holder pasteurization on nutrients and biologically-active components in donor human milk: A review. Nutrients 2016, 8, 447. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. García-Lara, N.R.; Vieco, D.E.; De La Cruz-Bértolo, J.; Lora-Pablos, D.; Velasco, N.U.; Pallás-Alonso, C.R. Effect of holder pasteurization and frozen storage on macronutrients and energy content of breast milk. J. Pediatr. Gastroenterol. Nutr. 2013, 57, 377–382. [Google Scholar] [CrossRef]
  56. Billard, H.; Simon, L.; Desnots, E.; Sochard, A.; Boscher, C.; Riaublanc, A.; Alexandre-Gouabau, M.C.; Boquien, C.Y. Calibration adjustment of the mid-infrared analyzer for an accurate determination of the macronutrient composition of human milk. J. Hum. Lact. 2016, 32, NP19–NP27. [Google Scholar] [CrossRef]
  57. Folch, J.; Lees, M.; Sloane Stantley, G.H. A simple method for the isolation and purification of total lipides from animal tissues. J. Biol. Chem. 1957, 266, 497–509. [Google Scholar]
  58. Bligh, E.G.; Dyer, W.J. A rapid method of total lipid extraction and purification. Can. J. Biochem. Physiol. 1959, 37, 911–917. [Google Scholar] [CrossRef] [Green Version]
  59. Broadhurst, D.; Goodacre, R.; Reinke, S.N.; Kuligowski, J.; Wilson, I.D.; Lewis, M.R.; Dunn, W.B. Guidelines and considerations for the use of system suitability and quality control samples in mass spectrometry assays applied in untargeted clinical metabolomic studies. Metabolomics 2018, 14, 72–89. [Google Scholar] [CrossRef] [Green Version]
  60. Sumner, L.W.; Samuel, T.; Noble, R.; Gmbh, S.D.; Barrett, D.; Beale, M.H.; Hardy, N. Proposed minimum reporting standards for chemical analysis Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). Metabolomics 2007, 3, 211–221. [Google Scholar] [CrossRef] [Green Version]
  61. Vinaixa, M.; Schymanski, E.L.; Neumann, S.; Navarro, M.; Salek, R.M.; Yanes, O. Mass spectral databases for LC/MS- and GC/MS-based metabolomics: State of the field and future prospects. TrAC Trends Anal. Chem. 2016, 78, 23–35. [Google Scholar] [CrossRef] [Green Version]
  62. Fiehn, O.; Barupal, D.K.; Kind, T. Extending biochemical databases by metabolomic surveys. J. Biol. Chem. 2011, 286, 23637–23643. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Wishart, D.S.; Feunang, Y.D.; Marcu, A.; Guo, A.C.; Liang, K.; Vázquez-Fresno, R.; Sajed, T.; Johnson, D.; Li, C.; Karu, N.; et al. HMDB 4.0: The human metabolome database for 2018. Nucleic Acids Res. 2018, 46, D608–D617. [Google Scholar] [CrossRef] [PubMed]
  64. Smith, C.A.; O’Maille, G.; Want, E.J.; Qin, C.; Trauger, S.A.; Brandon, T.R.; Custodio, D.E.; Abagyan, R.; Siuzdak, G. METLIN: A metabolite mass spectral database. Ther. Drug Monit. 2005, 27, 747–751. [Google Scholar] [CrossRef] [PubMed]
  65. Kind, T.; Wohlgemuth, G.; Lee, D.Y.; Lu, Y.; Palazoglu, M.; Shahbaz, S.; Fiehn, O. FiehnLib—Mass spectral and retention index libraries for metabolomics based on quadrupole and time-of-flight gas chromatography/mass spectrometry. Anal. Chem. 2009, 81, 10038–10048. [Google Scholar] [CrossRef] [Green Version]
  66. Cardiff University; Babraham Institute; University of California, S.D. LIPID MAPS Lipidomics Gateway. Available online: http://www.lipidmaps.org/ (accessed on 8 November 2019).
  67. Foroutan, A.; Guo, A.C.; Vazquez-fresno, R.; Lipfert, M.; Zhang, L.; Zheng, J.; Badran, H.; Budinski, Z.; Mandal, R.; Ametaj, B.N.; et al. Chemical composition of commercial cow’s milk. J. Agric. Food Chem. 2019, 67, 4897–4914. [Google Scholar] [CrossRef]
  68. Milk Composition Database. Available online: http://www.mcdb.ca/ (accessed on 5 November 2019).
  69. KEGG PATHWAY Database. Available online: https://www.kegg.jp/kegg/pathway.html (accessed on 8 November 2019).
  70. Li, L.; Li, R.; Zhou, J.; Zuniga, A.; Stanislaus, A.E.; Wu, Y.; Huan, T.; Zheng, J.; Shi, Y.; Wishart, D.S.; et al. MyCompoundID: Using an evidence-based metabolome library for metabolite identification. Anal. Chem. 2013, 85, 3401–3408. [Google Scholar] [CrossRef]
  71. Gil de la Fuente, A.; Godzien, J.; Saugar, S.; Garcia-Carmona, R.; Badran, H.; Wishart, D.S.; Barbas, C.; Otero, A. CEU Mass Mediator 3.0: A Metabolite Annotation Tool. J. Proteome Res. 2019, 18, 797–802. [Google Scholar] [CrossRef]
  72. CEU Mass Mediator. Available online: http://ceumass.eps.uspceu.es/mediator/ (accessed on 5 November 2019).
  73. Oliveros, J.C. Venny. An Interactive Tool for Comparing Lists with Venn’s Diagrams. Available online: http://bioinfogp.cnb.csic.es/tools/venny/index.html (accessed on 4 November 2019).
Metabolites 10 00043 g001 550
Figure 1. Flow diagram of literature selection and review process. Search “criterion 1”: term (“human milk” OR “breast milk”), AND “metabolom*”, AND “infant”; only articles. Search “criterion 2”: term (“human milk” OR “breast milk”), AND “metabolom*”, AND (“GC” OR “LC” OR “NMR” OR “CE”); only articles. Web of Science database was employed for literature search.
Figure 1. Flow diagram of literature selection and review process. Search “criterion 1”: term (“human milk” OR “breast milk”), AND “metabolom*”, AND “infant”; only articles. Search “criterion 2”: term (“human milk” OR “breast milk”), AND “metabolom*”, AND (“GC” OR “LC” OR “NMR” OR “CE”); only articles. Web of Science database was employed for literature search.
Metabolites 10 00043 g001
Metabolites 10 00043 g002 550
Figure 2. Reporting frequency of factors relevant to the human milk (HM) sampling process: Maternal-infant-related factors (blue bars), time-related factors (green bars), and HM collection-related factors (orange bars). Note: BMI = body mass index.
Figure 2. Reporting frequency of factors relevant to the human milk (HM) sampling process: Maternal-infant-related factors (blue bars), time-related factors (green bars), and HM collection-related factors (orange bars). Note: BMI = body mass index.
Metabolites 10 00043 g002
Metabolites 10 00043 g003 550
Figure 3. Sample preparation approaches employed in human milk (HM) metabolomics.
Figure 3. Sample preparation approaches employed in human milk (HM) metabolomics.
Metabolites 10 00043 g003
Metabolites 10 00043 g004 550
Figure 4. Venn diagram of metabolites reported in human milk (HM) according to the technique in [73]. Note: GC-MS, gas chromatography—mass spectrometry; LC-MS, liquid chromatography—mass spectrometry; NMR, nuclear magnetic resonance.
Figure 4. Venn diagram of metabolites reported in human milk (HM) according to the technique in [73]. Note: GC-MS, gas chromatography—mass spectrometry; LC-MS, liquid chromatography—mass spectrometry; NMR, nuclear magnetic resonance.
Metabolites 10 00043 g004
Metabolites 10 00043 g005 550
Figure 5. Distribution of metabolite classes annotated and/or identified in HM according to technique. Note: GC-MS, gas chromatography—mass spectrometry; LC-MS, liquid chromatography—mass spectrometry; NMR, nuclear magnetic resonance.
Figure 5. Distribution of metabolite classes annotated and/or identified in HM according to technique. Note: GC-MS, gas chromatography—mass spectrometry; LC-MS, liquid chromatography—mass spectrometry; NMR, nuclear magnetic resonance.
Metabolites 10 00043 g005
Table
Table 1. Sample preparation steps and platforms employed in untargeted analysis of HM metabolome.
Table 1. Sample preparation steps and platforms employed in untargeted analysis of HM metabolome.
Sample Preparation (1st. step)Sample Preparation (2nd. step)Compound ClassPlatformColumn/CapillaryReferences
Bligh & Dyer extractionDeuterated solvent addition to aqueous phasePolar metabolites1H-NMR-[13,16,29,32]
Derivatization of aqueous phase: methoximation and silylationPolar metabolites and FAsGC-MS DB-5ms[17,18,19]
Derivatization of organic phase: methylation FAsGC-MS DB-5ms[13]
Direct injection of aqueous phasePolar metabolitesLC-QTOF-MS (+)HILIC[35]
Redissolution of aqueous phase in H2O:ACN (95:5)Polar metabolitesLC-Orbitrap-MS (+, −)C18[24]
Redissolution of organic phase in
(ACN:IPA:H2O (65:30:5)		Lipidic metabolitesLC-Orbitrap-MS (+,−)C18[25]
Folch extractionDeuterated solvent addition to aqueous and organic phasesHydrophobic and polar metabolites1H-NMR-[28]
Redissolution of aqueous phase in formic acid and centrifugationPolar metabolites (amino acids)CE-TOF-MS (+) 60 m × 50 µm I.D.
Redissolution of organic phase in IPA:H2O:ACN (2:1:1) and centrifugationLipidic metabolitesUPLC-QTOF-MS (+,−) C18
Single phase extractionDerivatization: methoximation and silylationPolar metabolites and FAsGC-MS DB-5ms[27,28]
Direct injectionLipidic (and polar) metabolitesLC-QTOF-MS (+,−) C8[27,28]
UPLC-QTOF-MS (+)C18[15]
Fat extraction with n-hexane/IPADeuterated solvent additionTGs13C-NMR; 1H-NMR-[20]
Filtration 3 kDa cutoff spin filterDeuterated solvent additionPolar metabolites1H-NMR-[14,21,22,29,33]
Protein precipitationDerivatization: methoximation and silylationPolar metabolitesGC-MS DB-5ms[36]
Hybrid SPE-Phospholipid extraction and redissolution in diluted organic phase of Bligh & Dyer extractionLipidic metabolitesLC-QTOF-MS (+)C8[35]
Fat removal with CH2Cl2 and dansylation of aqueous phasePolar metabolites (amine/phenol submetabolome)Chemical isotope labelling LC-QTOF-MS (+) C18[45,46]
Direct injectionPolar metabolites and FAsUPLC-QTOF-MS (+,−) C18[18]
Fat removal by centrifugationTwo additional centrifugations and deuterated solvent additionPolar metabolites1H-NMR-[34]
Filtration 10 kDa cutoff spin filter and deuterated solvent additionPolar metabolites1H-NMR-[23,26]
HomogenizationDeuterated solvent additionPolar metabolites1H-NMR-[31]
H2O-dilutionNaBH4-reduction and PGC cartridgeOligosaccharidesUPLC-TQD-MS (+)Hypercarb®[24]
CE, capillary electrophoresis; FAs, fatty acids; GC, gas chromatography; HILIC, hydrophilic interaction liquid chromatography; IPA, 2-propanol; I.D., inner diameter; LC, liquid chromatography; MS, mass spectrometry; 13C-NMR, carbon-13 nuclear magnetic resonance; 1H-NMR, proton nuclear magnetic resonance; PGC, porous graphitic carbon; QTOF, quadrupole time of flight; TGs, triacylglycerols; TQD, triple quadrupole; UPLC, ultraperformance liquid chromatography; +, positive ionization mode; -, negative ionization mode.
Table
Table 2. Most frequently reported metabolites (>80% of studies) according to technique.
Table 2. Most frequently reported metabolites (>80% of studies) according to technique.
Metabolite classLC-MSGC-MSNMR
Fatty acylsLinoleic acid (C18:2)Oleic acid (C18:1)-
Oleic acid (C18:1)Palmitic acid (C16:0)
Palmitoleic acid (C16:1)Stearic acid (C18:0)
GlycerolipidsDG (36:1)--
GlycerophospholipidsLysoPC (16:0)--
Carbohydrates and carbohydrate conjugates-FructoseLactose
Fucose
Ribose
Organic acids and derivatives-Malic acid
Urea		Acetate
Citrate
Lactate
Organo nitrogen compounds--Choline
Amino acids, peptides, and analogues-Alanine
Glutamate
Glycine
Pyroglutamic acid
Serine
Valine		Alanine
Creatine
Glutamate
Glutamine
Isoleucine
Leucine
Tyrosine
Valine
GC-MS, gas chromatography—mass spectrometry; LC-MS, liquid chromatography—mass spectrometry; NMR, nuclear magnetic resonance; DG, diacylglycerol; PC, phosphatidylcholine.

Share and Cite

MDPI and ACS Style

Ten-Doménech, I.; Ramos-Garcia, V.; Piñeiro-Ramos, J.D.; Gormaz, M.; Parra-Llorca, A.; Vento, M.; Kuligowski, J.; Quintás, G. Current Practice in Untargeted Human Milk Metabolomics. Metabolites 2020, 10, 43. https://doi.org/10.3390/metabo10020043

AMA Style

Ten-Doménech I, Ramos-Garcia V, Piñeiro-Ramos JD, Gormaz M, Parra-Llorca A, Vento M, Kuligowski J, Quintás G. Current Practice in Untargeted Human Milk Metabolomics. Metabolites. 2020; 10(2):43. https://doi.org/10.3390/metabo10020043

Chicago/Turabian Style

Ten-Doménech, Isabel, Victoria Ramos-Garcia, José David Piñeiro-Ramos, María Gormaz, Anna Parra-Llorca, Máximo Vento, Julia Kuligowski, and Guillermo Quintás. 2020. "Current Practice in Untargeted Human Milk Metabolomics" Metabolites 10, no. 2: 43. https://doi.org/10.3390/metabo10020043

APA Style

Ten-Doménech, I., Ramos-Garcia, V., Piñeiro-Ramos, J. D., Gormaz, M., Parra-Llorca, A., Vento, M., Kuligowski, J., & Quintás, G. (2020). Current Practice in Untargeted Human Milk Metabolomics. Metabolites, 10(2), 43. https://doi.org/10.3390/metabo10020043

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

Article Metrics

Back to TopTop