Next Article in Journal
A Multilevel Bayesian Approach to Improve Effect Size Estimation in Regression Modeling of Metabolomics Data Utilizing Imputation with Uncertainty
Previous Article in Journal
Target Analysis of Volatile Organic Compounds in Exhaled Breath for Lung Cancer Discrimination from Other Pulmonary Diseases and Healthy Persons
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolic Profiles of Whole, Parotid and Submandibular/Sublingual Saliva

1
Centro Universitario Odontoiatria, University of Parma, Via Gramsci 14, 43126 Parma, Italy
2
Department of Medicine and Surgery, University of Parma, Via Volturno 39, 43125 Parma, Italy
3
Transfusion Medicine Unit, Azienda USL-IRCCS di Reggio Emilia, Viale Umberto I 50, 43123 Reggio Emilia, Italy
*
Author to whom correspondence should be addressed.
Metabolites 2020, 10(8), 318; https://doi.org/10.3390/metabo10080318
Submission received: 15 July 2020 / Revised: 31 July 2020 / Accepted: 5 August 2020 / Published: 6 August 2020

Abstract

:
The detection of salivary molecules associated with pathological and physiological alterations has encouraged the search of novel and non-invasive diagnostic biomarkers for oral health evaluation. While genomic, transcriptomic, and proteomic profiles of human saliva have been reported, its metabolic composition is a topic of research: metabolites in submandibular/sublingual saliva have never been analyzed systematically. In this study, samples of whole, parotid, and submandibular/sublingual saliva from 20 healthy donors, without dental or periodontal diseases, were examined by nuclear magnetic resonance. We identified metabolites which are differently distributed within the three saliva subtypes (54 in whole, 49 in parotid, and 36 in submandibular/sublingual saliva). Principal component analysis revealed a distinct cluster for whole saliva and a partial overlap for parotid and submandibular/sublingual metabolites. We found exclusive metabolites for each subtype: 2-hydroxy-3-methylvalerate, 3-methyl-glutarate, 3-phenylpropionate, 4-hydroxyphenylacetate, 4-hydroxyphenyllactate, galactose, and isocaproate in whole saliva; caprylate and glycolate in submandibular/sublingual saliva; arginine in parotid saliva. Salivary metabolites were classified into standard and non-proteinogenic amino acids and amines; simple carbohydrates; organic acids; bacterial-derived metabolites. The identification of a salivary gland-specific metabolic composition in healthy people provides the basis to invigorate the search for salivary biomarkers associated with oral and systemic diseases.

Graphical Abstract

1. Introduction

Human whole saliva (WS) is a mixture of fluids produced by parotid (20%), submandibular (65–70%), sublingual (7% to 8%), minor (<10%) salivary glands, and by gingival sulcus (crevicular fluid) [1].
Chemical-physical properties and volume of WS can grossly vary among people, as well as in the same person, according to endogenous and exogenous factors (e.g., age, gender, circadian rhythm, psychological state, nutrition, diseases, drugs, and environmental exposures). Moreover, qualitative variations of saliva, mostly related to the presence and concentration of specific categories of molecules, have been reported [2,3,4].
Variation in salivary flow (e.g., by stimulation with citric acid) is associated with changes in most metabolites’ concentrations. The concentration of acetate in unstimulated saliva is markedly higher than in the stimulated one. By contrast, lactate has more elevated levels in stimulated than in unstimulated saliva. Glucose has a similar concentration in the two types of fluids [5].
It is important to highlight that the inter-individual variability of salivary metabolic profiles seems to be higher than the intra-individual one. Such observation has led to hypothesize that under standardized conditions, an individual metabolic phenotype is relatively stable [2].
Despite the relatively well-described composition of WS, parotid saliva (PS) and submandibular/sublingual saliva (SM/SL), in terms of nucleic acids and proteins [6,7], the metabolites composition of these fluids is still subject of research.
The origin of the metabolites in WS is quite diversified: some molecules are produced by human metabolic processes, others by oral microorganisms, and several are of exogenous origin [2,8,9,10]. Being unlikely a contamination by exogenous and microbial molecules, the majority of the metabolites in PS [5] and SM/SL are presumably host-derived. As for the metabolites of human origin in WS, they either originate within the salivary glands or are released from alive or desquamated oral mucosal cells [8].
To the best of our knowledge, SM/SL metabolites have never been systematically described: only few examples of metabolic profiles have been published so far [11]. In particular, Yamada-Nosaka and co-workers [12] recorded broad and not well resolved proton Nuclear Magnetic Resonance (NMR) spectra of SM/SL, most likely due to the presence of the viscous mucous component.
Currently, the increasing number of studies on salivary metabolites [11] point to interpret the metabolite content of each type of saliva, in light of physiological and pathological changes that characterize each salivary gland.
Indeed, in-depth analysis of the human salivary metabolome may significantly boost: (1) the research of salivary biomarkers for oral and systemic diseases [13,14,15,16] and (2) the interpretation of the metabolic alterations occurring in physiological conditions (e.g., effects of physical exercise, weight changes, activation of specific metabolic pathways) [17,18].
Here, we report and compare the metabolic composition of human unstimulated PS, SM/SL, and WS derived from a cohort of young and healthy volunteers in physiological conditions. Quantitative and qualitative differences between salivary types are discussed, with special emphasis on resident microflora contribution.

2. Results

For each participant (n = 20), a sample of unstimulated PS, SM/SL, and WS was collected separately, in this exact order, as described in the Material and Methods section (Figure 1).
The salivary flow rates (0.15 ± 0.16 mL/min for PS; 0.20 ± 0.09 mL/min for SM/SL; 0.49 ± 0.28 mL/min for WS, expressed as mean value ± SD) show a wide inter-individual variability, mainly for PS. Our flow rate values of unstimulated saliva are in accordance with flow rates reported in the literature: 0.1, 0.1, and 0.6 mL/min for PS, SM/SL, and WS, respectively [19,20]. Using a sialometry test, we preliminary measured WS flow rate obtaining a good correlation (r = 0.86) with the value measured at the end of the whole sampling procedure. This fact indicates that the flow rate is not significantly altered, even after the prolonged sample collection required for some participants.

2.1. Metabolite Content of Saliva Subtypes: Emerging Differences

The 1D 1H-NMR spectra of WS, SM/SL, and PS samples (Figure 2) highlight different and characteristic metabolites patterns.
Since the total signal area under each NMR spectrum is proportional to the total metabolite content, the comparison of the values derived from the whole cohort reveals a similar total metabolite content for WS and PS and a sensibly smaller value for SM/SL, with a ratio of 1:1:0.3 (WS:PS:SM/SL).
Principal component analysis (PCA) of WS, PS, and SM/SL was applied to binned NMR spectra of all samples (20 samples for each of the 3 salivary types). In Figure 3, WS appears as a well separated cluster of scores, while PS and SM/SL display a partial cluster overlap that might suggest some similarity between their metabolic profiles. As expected, all clusters are characterized by a spread of scores, very likely due to the contribution of inter-individual variability. The principal data variance is given by PC1, and the separation between WS and the other two salivary subtypes is accounted for by lactate, propionate, maltose, 2-aminoadipate, and taurine.

2.2. Salivary Metabolites

We identified 66 metabolites with average concentrations higher than 5 μM. Heatmap analysis of all salivary metabolites profiles highlights the heterogeneity of WS, PS, and SM/SL composition and shows that 54 are in WS, 49 in PS, and 36 in SM/SL (Figure 4). Overall, 32 metabolites (48%) are common to the three saliva subtypes. Notably, it has been possible to single out a number of metabolites uniquely present in each salivary subtype: 2-hydroxy-3-methylvalerate, 3-methyl-glutarate, 3-phenylpropionate, 4-hydroxyphenylacetate, 4-hydroxyphenyllactate, galactose, and isocaproate in WS; arginine in PS; caprylate and glycolate in SM/SL.

2.3. Classes of Metabolites

To describe quantitatively each type of saliva, we selected an ensemble of metabolites, focusing on molecules with relevant concentrations and/or differential expression in the three salivary types. The median value and the range of concentration for each of the selected metabolites are reported in Table 1. Table 1—Section A and Figure 4A indicate the presence of the majority of the standard amino acids, non-proteinogenic amino acids, e.g., pyroglutamate and taurine, and biogenic amines, e.g., cadaverine, creatine, homoserine, and putrescine. When the highest median value (Table 1, in bold) is found in WS, most frequently, we observe that the related metabolite concentrations are significantly higher than in SM/SL. Only for lysine, phenylalanine, proline, and putrescine, their WS concentrations are significantly higher also than the ones found for PS. In the case of alanine, creatine, glutamine, and taurine, instead, the PS concentrations result significantly higher than in WS and SM/SL.
In Table 1—Section B, the relevant simple carbohydrates are listed. PS presents high levels of glucose and maltose that are significantly higher than in WS and SM/SL (Figure 4B).
Finally, Table 1—Section C summarizes the main organic acids. Once more, the highest median value belongs more frequently to WS. Exceptions are lactate and citrate that are significantly higher in PS than in WS and SM/SL (Figure 4C).
Figure 4D and Table 1—Section A and C report the distribution of metabolites that, according to the “Human Metabolome Database” (www.hmbd.ca), are mainly referable to bacterial metabolism. Noteworthy, WS presents high abundance of short chain fatty acids (SCFAs, i.e., formate, acetate, propionate, and butyrate), products of amino acid degradation such as putrescine and 5-aminopentanoate, and metabolic products of aromatic amino acids fermentation such as 4-hydroxy-phenyllactate and 3-phenylpropionate [21]. Yet, some of these metabolites are detected only in gland saliva, possibly reflecting either a microbial contamination or a host gland contribution. In fact, cadaverine and homoserine, detected in PS and SM/SL, are absent in WS; glycolate is detected only in SM/SL.

2.4. Salivary Cell Count

Prokaryotic cell counts are significantly higher in WS than in PS and SM/SL, reflecting the prokaryotic metabolite proportion observed in the three salivary types (Figure 5). The median value of WS cell distribution (Figure 5) is consistent with the value reported by Sender et al. [22].
In all types of saliva, eukaryotic cell counts range between 0.1 and 0.5 × 106 cells/mL, and they are approximately three orders of magnitude lower than the prokaryotic ones.

3. Discussion

Mapping the human oral metabolome, with emphasis on the metabolic composition of saliva subtypes, is expected to provide hints to clarify the physiologic and pathologic processes of the salivary glands and oral cavity. In this study, we were able to identify and quantify a considerably high number of metabolites: 54, 49, and 36 in WS, PS, and SM/SL, respectively.
Free salivary amino acids are known to be primarily produced by endogenous and exogenous proteases of salivary glands, exfoliating cells, and oral microflora [23]. Proteolytic amino acid-degrading bacteria dissect proteins and peptides into amino acids and convert them in short chain fatty acids [24], contributing, together with the saccharolytic bacteria, to the organic acid content of saliva.
Thus, we hypothesize that the high concentration of amino acids and organic acids found in WS (Table 1—Sections A and C) may reflect bacterial metabolic pathways.
Indeed, glutamine, glycine, and proline, the most represented residues in salivary Proline-Rich Proteins (PRPs) [25], are among the amino acids found at the highest concentrations (Table 1—Section A) being produced by proteolytic processes. Pyroglutamic acid, the common N-terminal of acidic PRP, is present as well at significant concentrations in saliva samples (Table 1—Section A). Moreover, salivary mucins, a heterogeneous group of glycoproteins synthesized and secreted by the submandibular, sublingual, and minor salivary glands, contribute to proline salivary concentration [26]. Arginine, ornithine, and lysine, which originate from proteins and peptides lysis, are metabolized by oral cavity bacteria, and contribute to the salivary content of putrescine, by decarboxylation of ornithine, an intermediate in the degradation of arginine, and cadaverine, formed by lysine decarboxylation [27]. Our data (Table 1—Section A) point to the presence of cadaverine only in gland salivary types.
We found a significant concentration of arginine only in PS (Table 1—Section A) in agreement with Van Wuyckhuyse and co-workers [28] that found a concentration of free arginine and lysine in PS of caries-free adults significantly higher than in caries-susceptible individuals. Noteworthy, arginine deiminase system is a relevant source of alkali generation by means of ammonia production. Accordingly, this enzymatic system of both saliva and dental plaque results more active in caries-free people when compared to caries-active individuals, likely contributing to the neutralization of plaque acids and to caries resistance [29].
Taurine, a beta-sulfonic amino acid, is probably the most abundant free amino acid in mammalian tissues. As suggested by Revenga-Parra and co-workers [30], the determination of its concentration in various body fluids seems relevant for the early diagnosis of Alzheimer’s disease, growth retardation, diabetes mellitus, epilepsy, sepsis, and some types of cancers. In rat SM glands, taurine is suggested to act as a regulator of the saliva ionic strength [31] and, in human saliva, appears to be correlated to physical stress [18]. Our study shows that such metabolite is present in all three types of saliva with a marked prevalence in PS (Table 1—Section A), thus, suggesting a role in salivary glands function.
Glucose, a blood component, passes in saliva through the salivary gland apparatus in proportion to its blood concentration. A statistically significant positive correlation has been found between fasting salivary glucose and fasting blood glucose [32]. Our study shows that glucose level in PS not only is considerably higher than in WS, as reported by Wang and co-workers [33], but is also higher than the concentration found in SM/SL, suggesting that parotid gland is the primary route of entry (Table 1—Section B). Therefore, we can assume that, in fasting conditions, such as in our cohort, the WS glucose level is severely influenced by the oral microorganism’s glucose metabolism as well as by fasting blood glucose concentration.
Salivary glycoproteins are a suite of macromolecules that, while contributing with specific functional roles to the oral cavity defense, constitute endogenous nutrients for the resident oral microflora, thus, being also responsible for the microbial plaque growth [34]. Because glycoproteins degradation is the result of the combined action of various microbial glucosidases [24], we conclude that the prevalent presence of monosaccharides such as fucose, N-acetylglucosamine, and galactose in WS samples (Table 1—Section B) should primarily be associated with the microbial saccharolytic activity on the oligosaccharide chains linked to glycoproteins.
Interestingly, the unexpected presence of salivary maltose (Table 1—Section B), a significant source of carbohydrate to oral bacteria, might be attributed to the digestive action of α-amylase, an enzyme produced by serous cells of parotid glands, with minor contributions from other glands’ enzymes [35]. On the basis of the results reported here, we can hypothesize that the concentrations of maltose in glandular saliva reflect the concentration of gland α-amylase, while its presence in WS is drastically reduced because of microbial utilization.
Saccharolytic microflora converts sugars to lactic, formic, acetic, succinic, and other organic acids through the glycolytic pathway [24,36]. On the other hand, in subgingival sites, asaccharolytic and/or proteolytic bacteria metabolize nitrogenous compounds derived from gingival crevicular fluid, creating an environment rich in SCFAs and ammonia [24,36].
Organic acids, in our saliva samples, are preferentially present in anionic form and frequently display their highest concentrations in WS (Table 1—Section C).
Consistent with the data reported by Gardner and co-workers [8], we find that acetate is the most abundant metabolite in WS samples, being present at relevant concentrations also in PS and SM/SL saliva. Interestingly, formate, acetate, propionate, and butyrate metabolites are present at a sensibly lower concentration in glandular saliva as compared to WS, a fact that might be due to their reduced bacterial contamination (Figure 5).
Lactate and citrate are more concentrated in PS than in WS and SM/SL, suggesting that the parotid gland is a relevant route of entry of these metabolites into the oral cavity [8]. However, WS lactate concentration may reflect also the contribution of microorganisms and oral mucosa cells [37]. The median lactate concentration that we measured in the WS of individuals with low dental plaque score is in agreement with the value reported by Gardner and colleagues [8].
Eventually, it is worth mentioning that we found a significant concentration of three metabolites of bacterial origin (glycolate, cadaverine, and homoserine) in PS and SM/SL saliva, but not in WS. We interpret the presence of those metabolites as an indication of microbial contamination of saliva. Particularly in the case of the glandular saliva, even if the fluids are collected at the close proximity of the excretory ducts, it is not possible to exclude a bacterial contamination of the outlet of the terminal portion of the salivary ducts, also in the absence of clinical signs of glandular infections. Because the outlets of the Stensen and Wharton ducts are very close to those portions of the dental arches where plaque, calculus, and periodontal diseases are more frequent (vestibular area of maxillary molars and lingual area of mandibular incisors) [38], even following rigorous procedures for saliva collection, it is possible that the area of sponge application is contaminated by dental and periodontal bacterial species. The hypothesis of some bacterial contamination is sustained by the high number of prokaryotic cells measured in all types of saliva (Figure 5). The absence of those metabolites in WS, on the other hand, may be the result of a pronounced dilution effect and/or additional degradation processes.
Overall, we have been able to identify and evaluate the concentration of a relevant number of metabolites in human saliva and to highlight qualitative and quantitative differences between WS, PS, and SM/SL saliva. Particularly, for the first time, we provided a metabolic profile of SM/SL.
We believe that mapping the human salivary metabolome is central for understanding most of the physiologic and pathologic oral metabolic pathways, including those related to the host–microbiome relationships.

4. Materials and Methods

The present study was approved by the Ethical Committee of the “Area Vasta Emilia Nord” (AVEN) (protocol number: 808/2018/SPER/UNIPR METASAL3). Written consent was obtained from all volunteers who participated in this study, according to the Declaration of Helsinki.

4.1. Subjects of the Study

A cohort of twenty healthy volunteers (10 males, 10 females), aged 19–25 years, qualified for saliva collection after oral clinical examination, interview for data acquisition on general medical history, and salivary flow rate assessment by sialometry test (modified Saxon Test) [39]. None of them revealed a full-mouth plaque score (FMPS) and/or full-mouth bleeding score (FMBS) higher than 25%.

4.2. Saliva Collection

For each participant, a sample of PS, SM/SL, and WS was collected separately and in the absence of stimulation. Participants were asked to refrain from eating, smoking, and performing intense physical activity for at least 12 h before salivary sampling and to drink only water. Furthermore, it was requested not to carry out oral hygiene (tooth brushing and flossing) in the 45 min before saliva collection. The procedure took place between 8:00 a.m. and 10:00 a.m. to minimize the influence of the circadian rhythm on salivary composition. Immediately before collection, patients rinsed their mouth with water for 1 min. For PS and SM/SL collection, the outlets of the Stensen and Wharton ducts were isolated and gently cleaned with a sterile gauze. A sterile sponge, capable of absorbing saliva flow, was positioned on the outlet of the ducts. Periodically, the sponge was squeezed and a syringe was used to collect the saliva within a vial. WS was collected by the passive drooling method (Figure 1).
During collection, salivary samples were transferred to a tube containing NaN3 (0.5% final concentration) and kept on ice until a volume of 5.4 mL, of each salivary type, was obtained and then frozen at −80 °C.

4.3. Sample Preparation and 1H-NMR Spectra Collection and Analysis

Each frozen saliva sample was thawed at room temperature and centrifuged at 15,000× g for 10 min at 4 °C to remove eukaryotic and prokaryotic cells, cellular debris, and mucins, according to Gardner et al. [40]. The supernatants were protein-depleted by ultra-filtration, using Amicon Ultra-4 Centrifugal filters (3000 MWCO, Merck Millipore) at 4000× g for 120 min at 10 °C, and lyophilized.
For 1H-NMR measurements, each of the lyophilized samples was suspended in potassium phosphate buffer (50 mM, pH 7.4) and 3-trimethylsilyl propanoic acid (TSP) was added as the chemical shift reference (0.00 ppm) and quantitative internal standard.
High-resolution one-dimensional (1D) 1H-NMR spectra acquisition and processing were carried out according to Pertinhez et al. [41]. Metabolites identification and quantification were carried out using Chenomx NMR Suite 8.3 software (Chenomx Inc., Edmonton, AL, Canada).
Heatmap analysis was carried out on targeted metabolites, with concentrations higher than 5 μM at least for one saliva subtype. Heatmaps were generated using MetaboAnalystR (https://www.metaboanalyst.ca) [42], with normalization referenced to TSP and autoscaling.

4.4. Cell Counting

Eukaryotic cells (oral epithelial cells and leucocytes) and prokaryotic cells were counted according to Gardner et al. [40], to estimate their possible contribution to the metabolic profile of WS, SM/SL, and PS.

4.5. Statistical Analysis

To compare the metabolite composition of each saliva subtype, the upper-tailed Mann–Whitney test (Origin 2019 software) was applied. P < 0.05 was considered statistically significant. The saliva subtype with the highest median concentration of each metabolite is shown in bold in Table 1.
Unsupervised multivariate analysis. To produce an overview of the overall variability, NMR spectra datasets were analyzed by principal component analysis (PCA), using the PCA module of MestreNova 11.0 software (Mestrelab Research, Santiago de Compostela, Spain).

Author Contributions

M.M.—Contributed to conception, design, data interpretation, drafted and critically revised the manuscript. E.Q.—Contributed to data acquisition and interpretation, drafted the manuscript. R.A.—Contributed to sample collection. M.E.P.—Contributed to sample collection. B.G.—Contributed to data acquisition and interpretation. M.V.V.—Contributed to sample collection. G.S.—Contributed to sample collection. E.C.—Contributed to data acquisition and interpretation and statistical analysis. E.F.—Contributed to data interpretation, statistical analyses, drafted and critically revised the manuscript. T.C.—Contributed to data acquisition and interpretation. A.S.—Contributed to drafted critically revised the manuscript. T.A.P.—Contributed to conception, design, data acquisition and interpretation, drafted and critically revised the manuscript. All authors gave their final approval and agree to be accountable for all aspects of the work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

We gratefully acknowledge Paolo Govoni, University of Parma, for the use of microscope.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De Almeida, P.D.V.; Grégio, A.M.; Machado, M.A.; de Lima, A.A.; Azevedo, L.R. Saliva composition and functions: A comprehensive review. J. Contemp. Dent. Pract. 2008, 9, 72–80. [Google Scholar]
  2. Wallner-Liebmann, S.; Tenori, L.; Mazzoleni, A.; Dieber-Rothenederm, M.; Konrad, M.; Hofmann, P.; Luchinat, C.; Turano, P.; Zatloukal, K. Individual Human Metabolic Phenotype Analyzed by 1H NMR of Saliva Samples. J. Proteom. Res. 2016, 15, 1787–1793. [Google Scholar] [CrossRef] [PubMed]
  3. Setti, G.; Pezzi, M.E.; Viani, M.V.; Pertinhez, T.A.; Cassi, D.; Magnoni, C.; Bellini, P.; Musolino, A.; Vescovi, P.; Meleti, M. Salivary MicroRNA for Diagnosis of Cancer and Systemic Diseases: A Systematic Review. Int. J. Mol. Sci. 2020, 21, 907. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Slavish, D.C.; Graham-Engeland, J.E.; Smyth, J.M.; Engeland, C.G. Salivary markers of inflammation in response to acute stress. Brain Behav. Immun. 2015, 44, 253–269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Figueira, J.; Gouveia-Figueira, S.; Öhman, C.; Lif Holgerson, P.; Nording, M.L.; Öhman, A. Metabolite quantification by NMR and LC-MS/MS reveals differences between unstimulated; stimulated; and pure parotid saliva. J. Pharm. Biomed. Anal. 2017, 140, 295–300. [Google Scholar] [CrossRef] [PubMed]
  6. Fábián, T.K.; Fejérdy, P.; Csermely, P. Salivary genomics, transcriptomics and proteomics: The emerging concept of the oral ecosystem and their use in the early diagnosis of cancer and other diseases. Curr. Genom. 2008, 9, 11–21. [Google Scholar] [CrossRef] [PubMed]
  7. Schipper, R.G.; Silletti, E.; Vingerhoeds, M.H. Saliva as research material: Biochemical, physicochemical and practical aspects. Arch. Oral Biol. 2007, 52, 1114–1135. [Google Scholar] [CrossRef] [PubMed]
  8. Gardner, A.; Parkes, H.G.; So, P.W.; Carpenter, G.H. Determining bacterial and host contributions to the human salivary metabolome. J. Oral Microb. 2019, 11, 1617014. [Google Scholar] [CrossRef] [Green Version]
  9. Pereira, J.L.; Duarte, D.; Carneiro, T.J.; Ferreira, S.; Cunha, B.; Solares, D.; Cosa, A.L.; Gil, A.M. Saliva NMR metabolomics: Analytical issues in pediatric oral health research. Oral Dis. 2019, 25, 1545–1554. [Google Scholar] [CrossRef]
  10. Dame, Z.T.; Aziat, F.; Mandal, R.; Krishnamurthy, R.; Bouatra, S.; Borzouie, S.; Guo, A.C.; Sajed, T.; Deng, L.; Lin, H.; et al. The human saliva metabolome. Metabolomics 2015, 11, 1864–1883. [Google Scholar] [CrossRef]
  11. Gardner, A.; Carpenter, G.; So, P.W. Salivary Metabolomics: From Diagnostic Biomarker Discovery to Investigating Biological Function. Metabolites 2020, 10, 47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Yamada-Nosaka, A.; Fukutomi, S.; Uemura, S.; Hashida, T.; Fujishita, M.; Kobayashi, Y.; Kyogoku, Y. Preliminary nuclear magnetic resonance studies on human saliva. Arch. Oral Biol. 1991, 36, 697–701. [Google Scholar] [CrossRef]
  13. Lohavanichbutr, P.; Zhang, Y.; Wang, P.; Gu, H.; Nagana Gowda, G.A.; Djukovic, D.; Buas, M.F.; Raftery, D.; Che, C. Salivary metabolite profiling distinguishes patients with oral cavity squamous cell carcinoma from normal controls. PLoS ONE 2018, 13, e0204249. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Mikkonen, J.J.W.; Singh, S.P.; Akhi, R.; Salo, T.; Lappalainen, R.; González-Arriagada, W.A.; Lopes, M.A.; Kullaa, A.M.; Myllymaa, S. Potential role of nuclear magnetic resonance spectroscopy to identify salivary metabolite alterations in patients with head and neck cancer. Oncol. Lett. 2018, 16, 6795–6800. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Kumari, S.; Goyal, V.; Kumaran, S.S.; Dwivedi, S.N.; Srivastava, A.; Jagannathan, N.R. Quantitative metabolomics of saliva using proton NMR spectroscopy in patients with Parkinson’s disease and healthy controls. Neurol. Sci. 2020, 41, 1201–1210. [Google Scholar] [CrossRef] [PubMed]
  16. Meleti, M.; Cassi, D.; Vescovi, P.; Setti, G.; Pertinhez, T.A.; Pezzi, M.E. Salivary biomarkers for diagnosis of systemic diseases and malignant tumors. A systematic review. Med. Oral Patol. Oral Cir. Bucal 2020, 25, e299–e310. [Google Scholar] [CrossRef]
  17. Pitti, E.; Petrella, G.; Di Marino, S.; Summa, V.; Perrone, M.; D’Ottavio, S.; Bernardini, A.; Cicero, D.O. Salivary metabolome and soccer match: Challenges for understanding exercise induced changes. Metabolites 2019, 9, 141. [Google Scholar] [CrossRef] [Green Version]
  18. Cicero, D.O.; Di Marino, S.; Dinallo, V.; Pero, M.; Summa, V.; Desideri, A.; Bernardini, A.; Perondi, S.; D’Ottavio, S. A small sided game session affects salivary metabolite levels in young soccer players. Biomed. Spectrosc. Imaging 2016, 5, 55–70. [Google Scholar] [CrossRef]
  19. Percival, R.S.; Challacombe, S.j.; Marshi, P.D. Flow rates of resting whole and stimulated parotid saliva in relation to age and gender. J. Dent. Res. 1994, 73, 1416–1420. [Google Scholar] [CrossRef]
  20. Ship, J.A.; Nolan, N.E.; Puckett, S.A. Longitudinal analysis of parotid and submandibular salivary flow rates in healthy, different-aged adults. J. Gerontol. A Biol. Sci. Med. Sci. 1995, 50, M285–M289. [Google Scholar] [CrossRef] [Green Version]
  21. Smith, E.A.; Macfarlane, G.T. Formation of phenolic and indolic compounds by anaerobic bacteria in the human large intestine. Microb. Ecol. 1997, 33, 180–188. [Google Scholar] [CrossRef] [PubMed]
  22. Sender, R.; Fuchs, S.; Milo, R. Revised estimates for the number of human and bacteria cells in the body. PLoS Biol. 2016, 14, e1002533. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Thomadaki, K.; Helmerhorst, E.J.; Tian, N.; Sun, X.; Siqueira, W.L.; Walt, D.R.; Oppenheim, F.G. Whole-saliva proteolysis and its impact on salivary diagnostics. J. Dent. Res. 2011, 90, 1325–1330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Takahashi, N. Oral microbiome metabolism: From "who are they?" to "what are they doing?". J. Dent. Res. 2015, 94, 1628–1637. [Google Scholar] [CrossRef] [PubMed]
  25. Levine, M. Susceptibility to dental caries and the salivary proline-rich proteins. Int. J. Dent. 2011, 953412. [Google Scholar] [CrossRef]
  26. Zalewska, A.; Zwierz, K.; Zółkowski, K.; Gindzieński, A. Structure and biosynthesis of human salivary mucins. Acta Biochim. Pol. 2000, 47, 1067–1079. [Google Scholar] [CrossRef] [Green Version]
  27. Cooke, M.; Leeves, N.; White, C. Time profile of putrescine; cadaverine; indole and skatole in human saliva. Arch. Oral Biol. 2003, 48, 323–327. [Google Scholar] [CrossRef]
  28. Van Wuyckhuyse, B.C.; Perinpanayagam, H.E.; Bevacqua, D.; Raubertas, R.F.; Billings, R.J.; Bowen, W.H.; Tabak, L.A. Association of free arginine and lysine concentrations in human parotid saliva with caries experience. J. Dent. Res. 1995, 74, 686–690. [Google Scholar] [CrossRef]
  29. Reyes, E.; Martin, J.; Moncada, G.; Neira, M.; Palma, P.; Gordan, V.; Oyarzo, J.F.; Yevenes, I. Caries-free subjects have high levels of urease and arginine deiminase activity. J. Appl. Oral Sci. 2014, 22, 235–240. [Google Scholar] [CrossRef]
  30. Revenga-Parra, M.; Martínez-Periñán, E.; Moreno, B.; Pariente, F.; Lorenzo, E. Rapid taurine and lactate biomarkers determination with disposable electrochemical detectors. Electrochim. Acta 2017, 240, 506–513. [Google Scholar] [CrossRef]
  31. Mozaffari, M.S.; Borke, J.L. Taurine in submandibular gland of the rat: Effect of muscarinic drugs. J. Histochem. Cytochem. 2002, 50, 527–532. [Google Scholar] [CrossRef] [PubMed]
  32. Gupta, S.; Nayak, M.T.; Sunitha, J.D.; Dawar, G.; Sinha, N.; Rallan, N.S. Correlation of salivary glucose level with blood glucose level in diabetes mellitus. J. Oral Maxillofac. Pathol. 2017, 21, 334–339. [Google Scholar] [PubMed]
  33. Wang, B.; Du, J.; Zhu, Z.; Ma, Z.; Wang, S.; Shan, Z. Evaluation of parotid salivary glucose level for clinical diagnosis and monitoring type 2 diabetes mellitus patients. Biomed. Res. Int. 2017, 2569707. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Kejriwal, S.; Bhandary, R.; Thomas, B.; Kumari, S. Estimation of levels of salivary mucin, amylase and total protein in gingivitis and chronic periodontitis patients. J. Clin. Diagn. Res. 2014, 8, ZC56–ZC60. [Google Scholar] [CrossRef] [PubMed]
  35. Nikitkova, A.E.; Haase, E.M.; Scannapieco, F.A. Taking the starch out of oral biofilm formation: Molecular basis and functional significance of salivary α-amylase binding to oral streptococci. Appl. Environ. Microbiol. 2013, 79, 416–423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Takahashi, N. Microbial ecosystem in the oral cavity: Metabolic diversity in an ecological niche and its relationship with oral diseases. Int. Congr. Ser. 2005, 1284, 103–112. [Google Scholar] [CrossRef]
  37. Nagler, R.M.; Lischinsky, S.; Diamond, E.; Klein, I.; Reznick, A.Z. New insights into salivary lactate dehydrogenase of human subjects. J. Lab. Clin. Med. 2001, 137, 363–369. [Google Scholar] [CrossRef]
  38. Albandar, J.M.; Kingman, A. Gingival recession, gingival bleeding, and dental calculus in adults 30 years of age and older in the United States, 1988-1994. J. Periodontol. 1999, 70, 30–43. [Google Scholar] [CrossRef]
  39. Gomes, P.D.; Juodzbalys, G.; Fernandes, M.H.; Guobis, Z. Diagnostic approaches to sjögren’s syndrome: A literature review and own clinical experience. J. Oral Maxillofac. Res. 2012, 3, e3. [Google Scholar] [CrossRef]
  40. Gardner, A.; Parkes, H.G.; Carpenter, G.H.; So, P.W. Developing and standardizing a protocol for quantitative proton nuclear magnetic resonance (1H NMR) spectroscopy of saliva. J. Proteom. Res. 2018, 17, 1521–1531. [Google Scholar] [CrossRef] [Green Version]
  41. Pertinhez, T.A.; Casali, E.; Lindner, L.; Spisni, A.; Baricchi, R.; Berni, P. Biochemical assessment of red blood cells during storage by 1H nuclear magnetic resonance spectroscopy. Identification of a biomarker of their level of protection against oxidative stress. Blood Transfus. 2014, 12, 548–556. [Google Scholar] [PubMed]
  42. Chong, J.; Wishart, D.S.; Xia, J. Using MetaboAnalyst 4.0 for comprehensive and integrative metabolomics data analysis. Curr. Protoc. Bioinform. 2019, 68, e86. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Experimental set-up used for parotid (A), submandibular/sublingual (B), and whole saliva (C) collection.
Figure 1. Experimental set-up used for parotid (A), submandibular/sublingual (B), and whole saliva (C) collection.
Metabolites 10 00318 g001
Figure 2. 1H-NMR spectra of unstimulated WS, PS, and SM/SL saliva from the same participant, acquired at 25 °C. The left region of the spectra shows the vertical scale increased by a factor of 10.
Figure 2. 1H-NMR spectra of unstimulated WS, PS, and SM/SL saliva from the same participant, acquired at 25 °C. The left region of the spectra shows the vertical scale increased by a factor of 10.
Metabolites 10 00318 g002
Figure 3. PCA score plot of 1H-NMR spectra of whole saliva (blue), parotid saliva (red), and submandibular/sublingual saliva (green) samples. The contribution of the three principal components of the total variance are PC1 = 62.6%, PC2 = 6.8%, and PC3 = 5.1%.
Figure 3. PCA score plot of 1H-NMR spectra of whole saliva (blue), parotid saliva (red), and submandibular/sublingual saliva (green) samples. The contribution of the three principal components of the total variance are PC1 = 62.6%, PC2 = 6.8%, and PC3 = 5.1%.
Metabolites 10 00318 g003
Figure 4. Heatmap analysis of metabolites grouped according to the categories: (A) amino acids, (B) carbohydrates, (C) organic acids, and (D) selected prokaryotic metabolites. PS (red), SMS (green), and WS (blue) samples.
Figure 4. Heatmap analysis of metabolites grouped according to the categories: (A) amino acids, (B) carbohydrates, (C) organic acids, and (D) selected prokaryotic metabolites. PS (red), SMS (green), and WS (blue) samples.
Metabolites 10 00318 g004
Figure 5. Eukaryotic (A) and prokaryotic (B) cell count distributions obtained with the three salivary types of all the subjects of the study. The significance level of the independent t-tests is set at 0.05 (* p < 0.05).
Figure 5. Eukaryotic (A) and prokaryotic (B) cell count distributions obtained with the three salivary types of all the subjects of the study. The significance level of the independent t-tests is set at 0.05 (* p < 0.05).
Metabolites 10 00318 g005
Table 1. Concentrations of selected metabolites in whole, parotid, and submandibular/sublingual saliva.
Table 1. Concentrations of selected metabolites in whole, parotid, and submandibular/sublingual saliva.
METABOLITE a WS (μM) PS (μM) SM/SL (μM)
SECTION A: STANDARD AND NON–PROTEINOGENIC AMINO ACIDS
Alanine27.4 (9.4–212.6)69.1 (41.5–219.1) ° § 20.7 (8.1–74.9)
ArginineND 24.2 (5.9–56.3)ND
Aspartate19.4 (9.8–76.7)25.7 (7.9–64.0) §7.3 (1.8–18.7)
Cadaverine bND3.2 (0.7–17.6)2.5 (0–16.6)
Creatine10.8 (6.7–49.4)48.9 (23.7–116.7) ° § 16.6 (6.4–48.1)
Glutamate108.1 (30.6–250.1) §57.5 (27.7–224.0)19.4 (8.6– 71.2)
Glutamine28.4 (6.6–134.6)63.3 (9.1–151.4) ° § 21.6 (2.3–70.6)
Glycine82.5 (12.1–306.6) §66.6 (3.6–192.3)14.4 (4.5–76.5)
Histidine20.1 (4.1–49.9) §17.6 (7.8–70.7)5.2 (2.9–21.1)
Homoserine bND36.4 (0–158.4)32.4 (3.2–57.9)
Isoleucine5.0 (0.9–44.3) §3.8 (2.0–13.4)1.2 (0.7–4.1)
Leucine13.2 (3.5–56.9) §12.0 (6.3–30.7)4.8 (2.1–9.5)
Lysine66.7 (17.6–164.8) * § 19.4 (10.1–97.5) 6.2 (2.2–37.9)
Phenylalanine16.8 (6.4–48.6) § *10.9 (5.0–41.1) 4.6 (2.0–12.7)
Proline64.1 (24.8–446.8) * § 41.1 (10.8–156.0) 7.2 (2.8–74.4)
Putrescineb38.6 (8.5–96.4) * § 5.1 (0.6–27.3) 0.9 (0.5–17.7)
Pyroglutamate12.9 (2.9–70.5) § 9.4 (0–32.3)7.6 (3.3–14.7)
Taurine46.2 (2.8–132.0)121.4 (0–342.2) ° §60.2 (3.3–146.9)
Threonine4.6 (2.4–31.4)ND7.3 (2.9–22.3)
Tyrosine34.9 (10.5–93.5) § 28.2 (14.38–90.9)10.1 (3.0–30.8)
Valine9.3 (2.9–59.2)12.4 (3.0–38.0) § 4.7 (1.1–13.1)
SECTION B: SIMPLE CARBOHYDRATES
Fucose34.8 (11.5–275.8) * § 5.6 (0.9–57.9) 4.8 (2.4–16.3)
Galactose18.9 (6.3–173.9)ND ND
Glucose11.8 (6.8–137.8)204.6 (81.8–697.8) ° § 46.8 (7.6–211.8)
Maltose1.3 (0.2–52.8)296.9 (103.5–1587.9) ° § 76.5 (16.4–420.4)
N-acetylglucosamine26.3 (2.2–141.5) §ND10.5 (1.4–40.8)
SECTION C: ORGANIC ACIDS
2-Aminoadipate186.0 (77.2–530.3) * § 25.0 (2.4–117.3) 11.6 (2.0–92.8)
2-Hydroxybutyrate13.6 (1.4–53.1) * § 7.6 (4.3–12.8) 2.5 (0.9–6.5)
3-Phenylpropionate b10.0 (2.0–38.1)NDND
4-hydroxyphenylacetate b8.6 (1.9–19.0)ND ND
4-hydroxyphenyllactate b4.4 (0.7–23.9)NDND
5-Aminopentanoate b100.5 (0–386.3) *12.1 (0–82.2)ND
Acetate b2277.9 (734.1–4322.8) * § 470.7 (81.8–3145.0) 237.4 (54.0–1370.4)
Butyrate b20.9 (3.2–77.3) * § 5.1 (0.8–14.8) 4.0 (1.1–18.6)
Citrate12.7 (0.7–33.0)35.3 (16.0–125.8) ° §20.4 (5.9–146.9)
Formate b37.7 (8.7–234.0) * § 15.1 (6.8–106.7) 17.2 (6.5–97.9)
Glycolate bNDND7.8 (1.8–55.4)
Isobutyrate b18.6 (3.76–47.6) * § 1.3 (0.2–9.9) 1.7 (0.3–5.2)
Lactate123.1 (23.7–517.5)714.8 (408.8–1683.9) ° § 162.7 (84.0–444.5)
Phenylacetate b15.9 (4.2–46.2) * § 1.3 (0–25.7) 1.6 (0.3–3.3)
Propionate b261.8 (64.7–627.9) * § 31.2 (4.4–200.8) 17.5 (4.3–151.7)
Pyruvate18.8 (4.0–52.8) §12.1 (1.4–232.0)5.1 (2.1–14.2)
Succinate b16.1 (9.9–39.2) * § 12.1 (3.8–23.3)5.1 (2.2–11.8)
a In each section, metabolites are presented in alphanumerical order. The values reported for each metabolite are median concentrations in WS, PS and SM/SL. Numbers in brackets indicate the concentration range. Numbers in bold are the highest median values obtained in the three saliva types; For each metabolite, a Mann–Whitney test has been applied for comparing the concentration pool producing the highest median value with that of the other saliva samples, when detected. The significance level has been set at 0.05; The symbols °, *, and § are associated with the highest median values and indicate that the related metabolite concentration tends to be significantly higher than in WS, PS, or SMS, respectively; b Metabolite of bacterial origin, according to the “Human Metabolome Database” (www.hmbd.ca); ND is the abbreviation for not detected.

Share and Cite

MDPI and ACS Style

Meleti, M.; Quartieri, E.; Antonelli, R.; Pezzi, M.E.; Ghezzi, B.; Viani, M.V.; Setti, G.; Casali, E.; Ferrari, E.; Ciociola, T.; et al. Metabolic Profiles of Whole, Parotid and Submandibular/Sublingual Saliva. Metabolites 2020, 10, 318. https://doi.org/10.3390/metabo10080318

AMA Style

Meleti M, Quartieri E, Antonelli R, Pezzi ME, Ghezzi B, Viani MV, Setti G, Casali E, Ferrari E, Ciociola T, et al. Metabolic Profiles of Whole, Parotid and Submandibular/Sublingual Saliva. Metabolites. 2020; 10(8):318. https://doi.org/10.3390/metabo10080318

Chicago/Turabian Style

Meleti, Marco, Eleonora Quartieri, Rita Antonelli, Margherita E. Pezzi, Benedetta Ghezzi, Maria Vittoria Viani, Giacomo Setti, Emanuela Casali, Elena Ferrari, Tecla Ciociola, and et al. 2020. "Metabolic Profiles of Whole, Parotid and Submandibular/Sublingual Saliva" Metabolites 10, no. 8: 318. https://doi.org/10.3390/metabo10080318

APA Style

Meleti, M., Quartieri, E., Antonelli, R., Pezzi, M. E., Ghezzi, B., Viani, M. V., Setti, G., Casali, E., Ferrari, E., Ciociola, T., Spisni, A., & Pertinhez, T. A. (2020). Metabolic Profiles of Whole, Parotid and Submandibular/Sublingual Saliva. Metabolites, 10(8), 318. https://doi.org/10.3390/metabo10080318

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