# Untargeted Metabolomic Approach to Study the Impact of Aging on Salivary Metabolome in Women

Pauline Bosman, Valérie Pichon, Ana Carolina Acevedo, Laëtitia Le Pottier,

Jacques Olivier Pers, Hélène Chardin, Audrey Combès

► To cite this version:

Pauline Bosman, Valérie Pichon, Ana Carolina Acevedo, Laëtitia Le Pottier, Jacques Olivier Pers, et al.. Untargeted Metabolomic Approach to Study the Impact of Aging on Salivary Metabolome in Women. Metabolites, 2022, 12 (10), pp.986. 10.3390/metabo12100986. hal-05038095

HAL Id: hal-05038095

https://hal.univ-brest.fr/hal-05038095v1

Submitted on 17 Apr 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Untargeted Metabolomic Approach to Study the Impact of Aging on Salivary Metabolome in Women

Pauline Bosman<fnref n="1" />, Valérie Pichon <fnref n="1" />,<fnref n="2" /> ID, Ana Carolina Acevedo <fnref n="3" />,<fnref n="4" />, Laëtitia Le Pottier <fnref n="5" />, Jacques Olivier Pers <fnref n="5" />,<fnref n="6" /> ID,
Hélène Chardin <fnref n="1" />,<fnref n="4" /> and Audrey Combès <fnref n="1" />,* ID

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Bosman, P.; Pichon, V.; Acevedo, A.C.; Le Pottier, L.; Pers, J.O.; Chardin, H.; Combès, A.

Untargeted Metabolomic Approach
to Study the Impact of Aging on
Salivary Metabolome in Women.

Metabolites **2022**, *12*, 986. https://doi.org/10.3390/metabo12100986

Academic Editors: Tatiana K.
S. Fidalgo and Ana Paula Valente

Received: 23 September 2022

Accepted: 15 October 2022

Published: 18 October 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** Despite the growing interest in salivary metabolomics, few studies have investigated the impact of aging on the salivary metabolome. The alterations in metabolic pathways that occur with aging are likely to be observed in pathologies affecting older people and may interfere with the search for salivary biomarkers. It is therefore important to investigate the age-related changes occurring in the salivary metabolome. Using reversed phase liquid chromatography and hydrophilic interaction chromatography coupled to mass spectrometry used in positive and negative ionization modes, the salivary metabolic profiles of young (22 to 45 years old) and older people (55 to 92 years old) were obtained. Those profiles were compared with the use of XCMS online to highlight the under or overexpression of some metabolites with aging. A total of 60 metabolites showed differential expression with age. The identification of 26 of them was proposed by the METLIN database and, among them, 17 were validated by standard injections. Aging seemed to affect most of the main metabolic pathways (amino acid metabolism, Krebs cycle, fatty acid synthesis, and nucleic acid synthesis). Moreover, most of the metabolites that were over- or under-expressed with age in this study have already been identified as being potential biomarkers of diseases affecting older people, such as in Alzheimer's disease. Special attention should be paid in the search for biomarkers of pathologies affecting the elderly to differentiate age-related changes from disease-related changes.

**Keywords:** LC-HRMS; salivary metabolome; aging

## 1. Introduction

With increases in the aging of the population, the age structure of the world population has changed. Approximately 10% of the world's population is now over the age of 60 and this figure will almost double by 2050 [1,2]. Aging leads to various changes in the human body that are more or less deleterious and can lead to loss of function and low resistance to chronic disease [2,3].

The oral fluid, more commonly known as saliva, is a complex biological fluid essential for maintaining oral health and nutrition [4]. Saliva is a viscous, slightly acidic liquid composed mainly of water and containing organic and inorganic compounds [5,6]. Saliva is secreted by the major (submandibular, sublingual, and parotid) and minor salivary glands [7]. With age, histological changes occur in the structure of the salivary glands, leading to a decrease in the proportion of acini and increase in the proportions of adipose and fibrous tissue. These histological changes naturally result in the hypofunction of the salivary glands, and subsequent decreases in secretion and alterations of the composition

<footnote n="1">¹ Laboratoire des Sciences Analytiques, Bioanalytiques et Miniaturisation, UMR 8231 CBI CNRS, ESPCI Paris, PSL Université, 75005 Paris, France</footnote>

<footnote n="2">² Sorbonne Université, 75006 Paris, France</footnote>

<footnote n="3">³ Laboratory of Oral Histopathology, Health Sciences Faculty, University of Brasília, Brasília DF CEP 70910-900, Brazil</footnote>

<footnote n="4">⁴ Université Paris Cité, 75006 Paris, France</footnote>

<footnote n="5">⁵ LBAI, UMR 1227, Université de Brest, Inserm, 29200 Brest, France</footnote>

<footnote n="6">⁶ University Hospital of Brest, 29200 Brest, France</footnote>

<footnote>* Correspondence: audrey.combes@espci.fr</footnote>
of saliva [1,8–11]. In addition, there is a decrease in the number of olfactory and gustatory receptors in the elderly, as well as a decrease in the neuronal stimulation of the glands, which leads to decreases in the intensity of stimuli and secretion reflexes, which then reinforce the hypofunction of the glands [10,11]. Some contradictions appear in the literature concerning the decrease in salivary flow. These contradictions are largely due to variations in study design. Indeed, some studies differentiate between groups of young and older people who are too close in age, other studies exclude all individuals taking medication and only consider healthy people, and others include minimal and hormone medications and/or people with disease. Finally, there are variations in the salivary collection methods, with studies based on unstimulated, stimulated, or parotid saliva. After a literature review, Affoo et al. [12] concluded that salivary flow rates in the elderly are significantly lower than in younger people, and this decrease is mainly attributed to the activity of the sublingual and submandibular glands. The parotid gland, on the other hand, seems to be less affected by age [1,8–12].

Despite the rise of metabolomics and growing interest in the salivary metabolome for
biomarker discovery, only a few studies have been performed on the impact of aging on the
salivary metabolome. Some authors have found alterations in several important metabolic
pathways, such as glycolysis and gluconeogenesis, the pentose phosphate pathway, the
tricarboxylic acid cycle, the urea cycle, nucleotide metabolism, glutathione metabolism,
and acetylcholine metabolism [11,13]. Some alterations in amino acid metabolism may
be associated with diseases that often develop with age, such as sarcopenia, which is a
gradual loss of muscle mass and function, and osteoporosis, which is characterized by
reduced bone strength [14,15]. Therefore, to discover biomarkers for pathologies affecting
the elderly, we first need to characterize the natural aging process.

Thus, the aim of the present work was to identify alterations in the salivary metabolome related to the natural aging process. As sex differences exist, especially in the aging process, we chose to focus our study on the salivary metabolome of women. Thus, the saliva of two groups of healthy women, young (22 to 45 years old) and older (55 to 92 years old), were analyzed. The metabolic profiles obtained by LC-HRMS analysis were then compared, to highlight the metabolites over- or under-expressed in saliva; these expression differences were correlated with age-related metabolic pathways and disease.

**2. Materials and Methods**

**2.1. *Chemicals***

HPLC-grade ACN and MeOH were supplied by Carlo Erba (Val de Reuil, France) and formic acid by Sigma-Aldrich (Saint Quentin Fallavier, France). High-purity water was dispensed by a Milli-Q purification system (Millipore, Saint Quentin en Yvelines, France). Alanine, GABA, lysine, threitol, ornithine, tyrosine, phenylalanine, and fucose were purchased from Sigma Aldrich (St Louis, MO, USA). Pyroglutamic acid was purchased from EGA Chemie (Steinheim, Germany). Threonine, serine, glutamic acid, and leucine/isoleucine were purchased from Merck (Darmstadt, Germany). Succinic acid was purchased from Prolabo (Paris, France). Isovaleric acid was purchased from Fluka (St Louis, MO, USA). Hypoxanthine was purchased from Extrasynthese (Genay, France) and linoleic acid was purchased from Interchim (Montluçon, France).

**2.2. *Selection Criteria for Participants***

To participate in the study, participants self-reported their age and whether or not
they had any chronic pathologies or take any medications. This study involved a total of
50 healthy women who were divided into two groups: the control group, which included
24 women aged 22–45 years (32.8 ± 7.4), and the experimental group, including 26 women
aged 55–92 years (61.5 ± 7.8). This study design would allow us to see the potential effects
of aging on salivary metabolites. As aging is a phenomenon that can vary according to
each individual and as menopause, which may be a major determinant of metabolomic
alterations, happens between 45 and 55 years of age, a gap of 10 years was deliberately left
between the two groups.

2.3. *Collection and Storage of Salivary Samples*

Our group previously reported on the importance of standardizing sample collection and treatment to allow for proper comparison of metabolomic profiles [16]. In view of the different partners, the samples were collected either by aspiration (24 samples) or by spitting (26 samples), two methods that were previously shown to give similar results for metabolomic analysis [16]. However, the following rules were defined to maximize standardization of sample collection: (i) All the volunteers were asked to refrain from having any food or drink or brushing their teeth for at least 2 h before saliva collection; (ii) The samples were collected without prior rinsing of the mouth; (iii) The collection time was set at 5 min; and (iv) Subjects were asked not to swallow, so that all oral fluid produced during the collection time would be collected. After collection, all saliva samples were centrifuged at 4 °C, 4500 × g for 15 min and were stored at −20 °C until further analyses. Prior to LC-MS analyses, the samples were diluted by a factor 2 with ACN for HILIC and with a mixture of ACN/deionized water, 10/90, *v/v* for RP-LC.

2.4. Instrumentation and Analytical Conditions

Chromatographic separation and detection by mass spectrometry were performed according to Bosman et al. [16]. The analyses were performed using a liquid chromatograph (1260 Infinity, Agilent Technologies, Waldbronn Germany), coupled to a tandem mass spectrometer (microOTOF2, Bruker Daltonics, Bremen, Germany), and equipped with an electrospray ionization source. All data acquisitions were controlled using TOF Control software (version 3.4, Bruker Daltonics), and Hystar software (version 3.2, Bruker Daltonics) was used to interface the HPLC and MS systems. An Acclaim 120 C18 column (Thermo Scientific, Saint Quentin Fallavier, France, 150 × 2.1 mm, 3 µm) and an Acclaim HILIC 10 (Thermo Scientific, 150 × 2.1 mm, 3 µm) were used for RP-LC and HILIC separations, respectively. The temperature of columns was set at 40 °C, the mobile phase flow rate at 0.2 mL/min, and injection volume at 10 µL, for both chromatographic modes. A mobile phase composed of deionized water (A) and ACN (B) were both acidified with 0.1% of formic acid for RP-LC and HILIC methods, respectively. For RP-LC, a gradient consisting of 95% of A for 5 min, followed by an increase of 2%/min of B for 45 min, a plateau at 95% of B for 2 min, return to starting conditions, i.e., 5% of B, in 1 min, and ending with 10 min of equilibration was used. For HILIC, we used a gradient consisting of 98% of B for 5 min, followed by an increase of 0.3%/min of A for 45 min, then a plateau at 15% of A for 2 min, return to initial conditions in 1 min, and finally ending with 10 min equilibration. The MS detection was made alternatively in positive and negative mode in scan mode (m/z ranging from 50 to 500). The capillary voltage was set at 4000 V, nebulizer pressure at 20 psi, drying gas at 8 L/min, dry temperature at 200 °C, ion energy at 2.5 eV, and collision cell energy at 5 eV. A stepping in the basic mode that alternates every 0.25 s between 70 and 500 Vpp for RF collision and between 50 and 120 µs for the transfer time was used to facilitate the transfer of ions on the range of m/z. As previously demonstrated, the widest possible coverage of salivary metabolites, and thus, the most comprehensive comparison possible for metabolomic profiles was obtained when each saliva sample was analyzed using four analytical methods [16]. Two chromatographic separation modes (RP-LC and HILIC) coupled with detection by HRMS used in positive or negative ionization modes, were respectively named RP-LC+, RP-LC−, HILIC+, and HILIC−, having been thus implemented. The online software XCMS was then used to perform four comparisons, one for each analytical method used, between LC-MS profiles obtained in an untargeted way for women under 45 years and those obtained for women older than 55 years.
## 2.5. *Quality Control for Untargeted Metabolomics*

Each saliva sample was injected twice into the LC-MS system. After verifying that the two LC-MS injections of the same sample gave the same LC-MS profile, only one of the replicates was used to compare the LC-MS profiles between the two groups, using XCMS online [17]. Otherwise, a blank was inserted between each injection of a saliva sample to verify the absence of carry-over between the different samples. Finally, for quality control, we used a solution containing 90 compounds that had already been described in the literature and was confirmed by our group as being present in saliva [16], i.e., amino acids, fatty acids, and other small molecules, with *m/z* between 60.021 and 430.381 and belonging to a wide range of polarity (lop P values in the range −4.70–10.51). Indeed, in a previous study [16], we demonstrated that these compounds, which are representative of salivary metabolites in terms of molecular mass, structural family, and polarity, induced a response from the analytical system that was proportional to the analyte concentration and allowed us to obtain limits of quantification to the order of a hundred ng/mL. Thus, a solution containing the 90 selected standards was injected at the beginning, middle, and end of each sequence of injections in the LC-MS system, regardless of the method used (RP-LC or HILIC with MS detection in positive or negative mode). This allowed us to verify the absence of significant fouling (less than 25% decrease in the observed signal) of the MS detector. This control ensured the absence of bias in the measured signal, and thus, a proper comparison of the observed metabolomic profiles.

## 2.6. *Data Analysis*

The LC-HRMS profiles were converted into mzML files using Hystar software, version 3.2 (Bruker Technology, Bremen, Germany). They were then uploaded to XCMS Online (xcmsonline.scripps.edu) [17]. All LC-MS profiles were processed to extract the molecular features, using an allowance of 10 ppm on the experimental *m/z*, a minimum S/N ratio of 10, and the “center-wave” algorithm. The obi warp algorithm was used for the retention time correction and alignment, with a maximum of 5 s allowed for deviation and a step size of 0.5 *m/z*.

## 2.7. *Metabolic Pathways*

The metabolites identified using the METLIN database [18] were searched for in the Human Metabolome Database (hmdb.ca), KEGG network (genome.jp), and PubChem (pubchem.ncbi.nlm.nih.gov). Stock solutions with concentrations ranging from 0.1 to 1 mg/mL were prepared in appropriate solvents (water, ACN, MeOH, and mixtures of these solvents), depending on the solubility of the compound, and stored at 4 °C. Working solutions were prepared from these stock solutions at concentrations between 0.25 and 18 µg/mL and injected.

# 3. Results and Discussion

## 3.1. *Comparison of the LC-MS Profiles of Young and Elderly Women*

Due to the alignment algorithm, XCMS detected hundreds of features (defined by the same retention time during chromatographic separations and the same *m/z* ratio during the MS analysis) that were present in both groups. Among these features, a significant inter-group difference (up- or down-regulated, fold change minimum set at 1.5, and *p*-value inferior to 0.05) was observed using MS detection in positive mode for 26 features after RP-LC separation, as shown in Figure 1 in a cloud plot representing under- and over-expressed compounds in RP-LC-MS+ (results not shown for the three other methods), and 13 features after HILIC separation. Then, the extracted ion chromatograms (EIC) for the 39 corresponding features were generated and manually verified. Two features were just noise and did not correspond to a chromatographic peak, so they were eliminated. An additional step of manual verification of the raw mass spectra, corresponding with a given retention time, made it possible to eliminate six additional features which were eluted from the chromatographic column with the same retention time and that corresponded with an
ion in mass spectra belonging to the isotopic pattern of another ion already listed under its
monoisotopic form. These verifications were necessary and supported that software should
not be blindly used in without manual verification. The total list of 31 features (couple
of *m/z* and retention time), showing statistical differential expression with age (12 were
highlighted from the comparison of the LC-MS profile obtained after HILIC versus 19
highlighted after RP-LC separation and detection in positive mode) is given in Table S1
(Supplementary material). The same approach was applied to the LC-MS profile obtained
after separations on RP-LC and HILIC mode coupled with a MS detection in negative
mode. XCMS online highlighted 24 and 15 features after RP-LC and HILIC separation,
respectively. The manual verifications of EIC and raw mass spectra led us to eliminate four
features that were noise and did not correspond to a chromatographic peak from HILIC
analysis and three that did not correspond to monoisotopic mass from RP-LC−. The list
of the 32 features (couple of *m/z* and retention time), showing significative differential
expression with age (11 after HILIC and 21 after RP-LC separation), is shown in Table S1
(Supplementary material).

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Cloud plots showing under- and over-expressed features with age in RPLC/MS+. The red bubbles over the chromatogram represent the over-expressed features and the green bubbles under the under-expressed. The bigger the bubble is, the bigger the fold change, and the darker the bubble, the more significative the *p*-value. * corresponds to features that were removed after manual verification.</figcaption></figure>

Among these 63 features, only seven ions (same *m/z*) were found to be in common with
the RP and HILIC separations (two and four in positive and negative modes, respectively,
as seen in Figure 2, and one in common between RP-LC+ and HILIC−). This proves the
complementarity of the four analytical methods and the usefulness in implementing them
all. A closer look at these ions (they are highlighted in Table S1) demonstrated that their
differences in expression were similar (same direction of variation and same magnitude
of change), regardless of the chromatographic separation mode used. For example, the
ion with a *m/z* of 274.875 detected in the MS positive ionization mode decreased with
age by a factor of 1.7 and 1.5 depending on the chromatographic separation in HILIC and
RP-LC, respectively.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Venn diagram showing the number of features detected according to the analytical method used. RPLC+: reversed-phase liquid chromatography coupled to mass spectrometry used in positive ionization mode; RPLC -: reversed-phase liquid chromatography coupled to mass spectrometry used in negative ionization mode; HILIC+: hydrophilic interaction liquid chromatography coupled to mass spectrometry used in positive ionization mode; HILIC -: hydrophilic interaction liquid chromatography coupled to mass spectrometry used in negative ionization mode.</figcaption></figure>

Furthermore, one of the known effects of aging is decreasing the salivary flow rate, and thus, the volume of saliva produced during collection time [10,11], which could lead to a bias in the comparison of metabolic profiles. This could lead to = generalized over- or under-expression of all the compounds detected in the saliva of elderly women. However, it is interesting to note that, among the 57 features with different expression levels highlighted listed in Table S1 (Supplementary material), 35 are overexpressed and 22 are under-expressed in the group of elderly people. This clearly demonstrates the absence of such bias in this study.

3.2. *Identification of the over = or Under-Expressed Compounds*

To identify these 57 features showing differential expression with age, the METLIN database was used and identification was proposed for 26 features. For 17 of them, analytical standards were available, and they were thus injected using the same analytical method(s) that revealed their difference in expression with aging. The superposition of the EIC corresponding to the analytical standard with that of the corresponding compound (same *m/z*) in the saliva of young and elderly people confirmed their identification (same *m/z* and same retention time). Figure 3 shows an example of these EIC overlays for three compounds: succinic acid (*m/z* 117.090 or 119.051 according to the detection mode), hypoxanthine (*m/z* 137.047), and linoleic acid (*m/z* 279.240).

The 17 compounds for which identification could be confirmed using these comparisons with available analytical standards are summarized in Table 1. The nine remaining compounds are listed in Table 2 with their respective putative identification.

<table><caption>**Table 1.** Salivary metabolites for which the identification was confirmed.</caption><thead><tr><th>Measured<br>m/z</th><th>Fold Change<br>(over Expressed)</th><th>Raw Formula</th><th>Confirmed<br>Identification *</th><th>Log P</th><th>Analytical<br>Method</th><th>Retention Time<br>(min)</th></tr></thead><tbody><tr><td>88.041</td><td>1.8</td><td>C<sub>3</sub>H<sub>7</sub>NO<sub>2</sub></td><td>Alanine</td><td>-3.2</td><td>RPLC-</td><td>1.73</td></tr><tr><td>102.058</td><td>2.7</td><td>C<sub>4</sub>H<sub>9</sub>NO<sub>2</sub></td><td>γ-aminobutyric acid (GABA)</td><td>-2.3</td><td>RPLC-</td><td>1.75</td></tr><tr><td>103.056</td><td>2.4</td><td>C<sub>5</sub>H<sub>10</sub>O<sub>2</sub></td><td>Isovaleric acid</td><td>1.21</td><td>RPLC+</td><td>5.04</td></tr><tr><td>104.037</td><td>2.6</td><td>C<sub>3</sub>H<sub>7</sub>NO<sub>3</sub></td><td>Serine</td><td>-3.9</td><td>RPLC-</td><td>1.69</td></tr><tr><td>119.051</td><td>1.9</td><td rowspan="2">C<sub>4</sub>H<sub>6</sub>O<sub>4</sub></td><td rowspan="2">Succinic acid **</td><td rowspan="2">-0.4</td><td rowspan="2">RPLC+<br>HILIC-</td><td>2.66</td></tr><tr><td>117.020</td><td>3.0</td><td>6.61</td></tr><tr><td>120.083</td><td>2.3</td><td>C<sub>4</sub>H<sub>9</sub>NO<sub>3</sub></td><td>Threonine</td><td>-3.5</td><td>RPLC+</td><td>5.04</td></tr><tr><td>123.046</td><td>2.0</td><td>C<sub>4</sub>H<sub>10</sub>O<sub>4</sub></td><td>Threitol</td><td>-2.5</td><td>RPLC+</td><td>2.65</td></tr><tr><td>128.037</td><td>4.6</td><td>C<sub>5</sub>H<sub>7</sub>NO<sub>3</sub></td><td>Pyroglutamic acid</td><td>-0.89</td><td>RPLC-</td><td>2.54</td></tr><tr><td>130.089</td><td>3.1</td><td>C<sub>6</sub>H<sub>13</sub>NO<sub>2</sub></td><td>Leucine/Isoleucine</td><td>-1.6</td><td>RPLC-</td><td>3.12</td></tr><tr><td>131.083</td><td>2.3</td><td>C<sub>5</sub>H<sub>12</sub>N<sub>2</sub>O<sub>2</sub></td><td>Ornithine</td><td>-3.7</td><td>RPLC-</td><td>1.68</td></tr><tr><td>137.047</td><td>3.0</td><td>C<sub>5</sub>H<sub>4</sub>N<sub>4</sub>O</td><td>Hypoxanthine</td><td>-0.05</td><td>RPLC+</td><td>2.32</td></tr><tr><td>145.101</td><td>2.8</td><td>C<sub>6</sub>H<sub>14</sub>N<sub>2</sub>O<sub>2</sub></td><td>Lysine</td><td>-3.2</td><td>RPLC-</td><td>1.56</td></tr><tr><td>146.050</td><td>2.5</td><td>C<sub>5</sub>H<sub>9</sub>NO<sub>4</sub></td><td>Glutamic acid</td><td>-3.2</td><td>RPLC-</td><td>1.75</td></tr></tbody></table>
<table><caption>**Table 1.** *Cont.*</caption><thead><tr><th>Measured<br>m/z</th><th>Fold Change<br>(over Expressed)</th><th>Raw Formula</th><th>Confirmed<br>Identification *</th><th>Log P</th><th>Analytical<br>Method</th><th>Retention Time<br>(min)</th></tr></thead><tbody><tr><td>165.056</td><td>2.2</td><td>C<sub>6</sub>H<sub>12</sub>O<sub>5</sub></td><td>Fucose</td><td>-1.9</td><td>RPLC+</td><td>2.65</td></tr><tr><td>166.087</td><td>2.5</td><td>C<sub>9</sub>H<sub>11</sub>NO<sub>2</sub></td><td>Phenylalanine</td><td>-1.2</td><td>RPLC+</td><td>5.03</td></tr><tr><td>182.083</td><td>2.6</td><td>C<sub>9</sub>H<sub>11</sub>NO<sub>3</sub></td><td>Tyrosine</td><td>-1.5</td><td>RPLC+</td><td>2.65</td></tr><tr><td>279.236</td><td>1.9</td><td>C<sub>18</sub>H<sub>32</sub>O<sub>2</sub></td><td>Linoleic acid</td><td>6.42</td><td>RPLC-</td><td>53.75</td></tr></tbody></table>

* The identification of compounds was confirmed by comparing the *m/z* and retention time of the unknown compound with those of an analytical standard injected under the same analytical conditions (LC separation and MS detection mode). ** For these compounds, the identification was confirmed with two different analytical methods.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Extracted ion chromatograms for succinic acid (m/z 117.020), hypoxanthine (m/z 137.047), and linoleic acid (m/z 279.240) analytical standards and for the corresponding compounds in the saliva of young and elderly people: (a) succinic acid in HILIC/MS-, (b) succinic acid in RPLC/MS+, (c) hypoxanthine in RPLC/MS+, and (d) linoleic acid in RPLC/MS-.</figcaption></figure>

<table><caption>**Table 2.** Potential salivary metabolites identified with METLIN.</caption><thead><tr><th>Measured<br>m/z</th><th>Fold Change (over<br>or under Expressed)</th><th>Raw Formula</th><th>Putative Identification<br>(Error in ppm)</th><th>Analytical<br>Method</th><th>Retention Time<br>(min)</th></tr></thead><tbody><tr><td>136.076</td><td>2.1 (Over)</td><td>C<sub>8</sub>H<sub>9</sub>NO</td><td>Phenylacetamide (5)</td><td></td><td>2.65</td></tr><tr><td>138.068</td><td>2.0 (Over)</td><td>C<sub>6</sub>H<sub>13</sub>N</td><td>Cyclohexylammonium (4)</td><td></td><td>16.41</td></tr><tr><td>173.094</td><td>2.7 (Over)</td><td>C<sub>10</sub>H<sub>14</sub>O</td><td>Thymol (4)</td><td></td><td>2.03</td></tr><tr><td>195.088</td><td>1.9 (Over)</td><td>C<sub>8</sub>H<sub>10</sub>N<sub>4</sub>O<sub>2</sub></td><td>Caffein (4)</td><td>RPLC+</td><td>16.41</td></tr><tr><td>223.026</td><td>1.9 (Over)</td><td>C<sub>11</sub>H<sub>8</sub>N<sub>2</sub>O</td><td>Fuberidazole (3)</td><td></td><td>1.62</td></tr><tr><td>384.150</td><td>4.1 (Over)</td><td>C<sub>14</sub>H<sub>25</sub>NO<sub>11</sub></td><td>N-acetylactosamine or<br>B-mannose-N-acetylglucosamine (1)</td><td></td><td>1.82</td></tr><tr><td>98.984</td><td>2.5 (Over)</td><td>H<sub>3</sub>PO<sub>4</sub></td><td>Phosphoric acid (4)</td><td>HILIC+</td><td>2.27</td></tr><tr><td>256.103</td><td>2.6 (Over)</td><td>C<sub>9</sub>H<sub>15</sub>N<sub>5</sub>O<sub>4</sub></td><td>Tetrahydroneoptin (10)</td><td>RPLC-</td><td>2.27</td></tr></tbody></table>
### 3.3. *Evolution of Saliva Metabolome with Age and Potential Affection of Metabolism*

The Human Metabolome DataBase (HMDB) verified that all the compounds listed in Table 1 had already been described as present in saliva. Among the 17 compounds overexpressed in the saliva of the elderly women, more than half were amino acids or amino acid derivatives. Indeed, nine of them were free proteinogenic amino acids (alanine, serine, threonine, leucine/isoleucine, lysine, glutamic acid, phenylalanine, and tyrosine). Among these amino acids, four were amino acids (leucine/isoleucine, lysine, phenylalanine, and threonine) that could not be produced by the human organism (essential amino acids) and therefore need to be provided by the diet. Their presence in higher concentrations in the saliva of elderly women was a sign of an imbalance between protein synthesis (anabolism) and proteolysis (catabolism), in favor of proteolysis. This imbalance has been described as a characteristic of sarcopenia, which affects many elderly people [15,19]. The higher concentration of amino acids in the saliva of elderly women differs from the results obtained in Teruya et al. [13] study, in which the levels of amino acids decreased with age, but correlates partially with the results from Tanaka et al. [20], who observed higher levels of lysine. Ornithine is a non-proteinogenic amino acid that plays a crucial role in the urea cycle. This allows the excretion of the nitrogen generated by the breakdown of amino acids in protein and other nitrogen-containing molecules from the body in urine. It is also interesting to note the overexpression of pyroglutamic acid, an uncommon amino acid derivative in which the free amino group of glutamic acid, which is itself overexpressed in the saliva of elderly people, cyclizes to form a lactam. It has been shown in an animal model that pyroglutamic acid led to the release of GABA from the cerebral cortex [21], a neurotransmitter also overexpressed in the saliva of elderly women.

Hypoxanthine is a naturally occurring purine derivative (by spontaneous deamination of adenine) and an intermediate in the metabolism of adenosine and in the formation of nucleic acids by the nucleotide salvage pathway.

Concerning the two fatty acids belonging to the fatty acyl family that were overexpressed in saliva, linoleic acid is a doubly unsaturated fatty acid, also known as an omega-6 fatty acid. It is an essential fatty acid that is used in the biosynthesis of prostaglandins (via arachidonic acid) and cell membranes [22]. Isovaleric acid is a natural methyl-branched fatty acid found in a wide variety of plants and essential oils.

Fucose is the only carbohydrate that was detected as being overrepresented in the saliva of elderly women. It is a monosaccharide, hexose deoxy sugar, and it is a common component of many N- and O-linked glycans and glycolipids. As a free sugar, L-fucose is normally found at very low levels of concentration in the saliva of mammals [23]. However, high concentrations of free L-fucose in serum and urine have been described to be a biomarker for different pathologies that affect the liver (cancer, cirrhosis, alcoholic liver disease, and gastric ulcers) [24,25]. The increase in this sugar in the saliva of elderly people could therefore be linked to a decrease in the metabolic activity of the liver in connection with aging [26]. In the same pathway, threitol is the main end-product of D-xylose metabolism in humans.

Another major metabolic pathway that seems to be affected by aging is the Krebs cycle. Succinic acid, a dicarboxylic acid that is an important component of the Krebs cycle, and thus, of the formation of adenosine triphosphate (ATP), was found in higher concentrations in the saliva of elderly women. Succinic acid acts as a signaling molecule reflecting the cellular metabolic state and regulates many processes in living cells. It has also been shown that succinic acid is an oncometabolite [27] that promotes tumor development and progression by inhibiting the activity of several 2-oxoglutarate-dependent enzymes. This superfamily of enzymes participates in many biological processes, including adaptation to hypoxia, epigenetic reprogramming, or the maturation and stability of collagen fibers.

Given all the metabolites identified (putatively or confirmed), aging seems to affect most of the main metabolic pathways of the organism (proteolysis, Krebs cycle, fatty acid synthesis, and nucleic acid synthesis). However, almost all of the salivary metabolites that were found to be overexpressed with age have already been suggested to be associated
with pathologies of the oral sphere (periodontitis, dental caries, oral cancers, etc.) or with systemic diseases, such as pancreatic, liver or breast cancers, Alzheimer's disease, or Lewy body disease [26,28,29].

**4. Conclusions**

The metabolic profiles of the saliva of women younger than 45 years and older than 55 years were obtained with the use of RP-LC and HILIC, coupled to HRMS, used in positive and negative ionization modes. When these metabolic profiles were compared using the statistical tool XCMS Online, our results showed that 35 features (coupled with retention time and *m/z*) were overexpressed and 22 under-expressed in older women. Using the METLIN database, an identification was proposed for 26 metabolites and confirmed for 17 of them. Among these 17 metabolites that showed differential expression with aging, 11 amino acids and amino acid derivatives, a purine derivative, 2 carbohydrate and carbohydrate derivatives, 2 fatty acids, and a dicarboxylic acid were found. Given the identified metabolites, the metabolic pathways that seemed to be affected with age were mainly amino acid metabolism, the Krebs cycle, fatty acid metabolism, and nucleic acid synthesis. Some of the metabolites identified in this study as being over- or under-expressed with age have been put forward as potential biomarkers of pathologies of the elderly, which raises questions about whether these biomarkers are related to age or to the pathology. Our results show, therefore, the essential need to characterize the natural process of aging before the identification of biomarkers for pathologies affecting older people.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/metabo12100986/s1, Table S1: *m/z* values of features statistically over- or under-expressed with age after analysis in RP-LC or HILIC coupled to HRMS in positive or negative mode.

**Author Contributions:** Conceptualization, P.B., A.C.A., H.C. and A.C.; methodology, H.C. and A.C.; formal analysis, P.B., H.C. and A.C.; investigation, P.B., A.C.A., H.C. and A.C.; resources, P.B., J.O.P., L.L.P. and A.C.A.; writing—original draft preparation, P.B. and A.C.; writing—review and editing, P.B., A.C.A., H.C., J.O.P., L.L.P., V.P. and A.C.; supervision, A.C.A., H.C. and A.C.; funding acquisition, V.P., H.C. and A.C. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by the "Plateforme d'expertise pour le suivi de toxines à l'état de traces en milieu réel" project from DIM analytics from the Ile de France region.

**Institutional Review Board Statement:** The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee (Hospital Center Regional University Morvan De Brest, CHRU Brest, CCPP Ouest VI, BB-0033-00037) in 2015. The Brest Biological Resources Center (CRB Santé du CHRU de Brest) provided the saliva samples for this study.

**Informed Consent Statement:** Informed consent was obtained from all subjects involved in the study.

**Data Availability Statement:** All data and materials are stored in the Laboratory of Bioanalytical Sciences and Miniaturization (LSABM), UMR CBI 8231, ESPCI Paris, CNRS, PSL University, Paris, France. The data are not publicly available due to ethical concerns.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. Nassar, M.; Hiraishi, N.; Islam, M.S.; Otsuki, M.; Tagami, J. Age-Related Changes in Salivary Biomarkers. *J. Dent. Sci.* **2014**, *9*, 85-90. [CrossRef]

2. Khalil, N.A. Saliva Can Be an Indicator for Aging. A Review. *Curr. Sci. Int.* **2022**, *11*, 84-98.

3. Kudryashova, K.S.; Burka, K.; Kulaga, A.Y.; Vorobyeva, N.S.; Kennedy, B.K. Aging Biomarkers: From Functional Tests to Multi-Omics Approaches. *Proteomics* **2020**, *20*, 1900408. [CrossRef] [PubMed]

4. Humphrey, S.P.; Williamson, R.T. A Review of Saliva: Normal Composition, Flow, and Function. *J. Prosthet. Dent.* **2001**, *85*, 162-169. [CrossRef]

5. Carpenter, G.H. The Secretion, Components, and Properties of Saliva. *Annu. Rev. Food Sci. Technol.* **2013**, *4*, 267-276. [CrossRef]
6. Yoshizawa, J.M.; Schafer, C.A.; Schafer, J.J.; Farrell, J.J.; Paster, B.J.; Wong, D.T.W. Salivary Biomarkers: Toward Future Clinical and Diagnostic Utilities. *Clin. Microbiol. Rev.* **2013**, *26*, 781–791. [CrossRef]

7. Roblegg, E.; Coughran, A.; Sirjani, D. Saliva: An All-Rounder of Our Body. *Eur. J. Pharm. Biopharm.* **2019**, *142*, 133–141. [CrossRef]

8. Nagler, R.M.; Hershkovich, O. Age-Related Changes in Unstimulated Salivary Function and Composition and Its Relations to Medications and Oral Sensorial Complaints. *Aging Clin. Exp. Res.* **2005**, *17*, 358–366. [CrossRef]

9. Vissink, A.; Spijkervet, F.K.L.; Amerongen, A.V.N. Aging and Saliva: A Review of the Literature. *Spec. Care Dent.* **1996**, *16*, 95–103. [CrossRef]

10. Xu, F.; Laguna, L.; Sarkar, A. Aging-Related Changes in Quantity and Quality of Saliva: Where Do We Stand in Our Understanding? *J. Texture Stud.* **2019**, *50*, 27–35. [CrossRef]

11. Toan, N.K.; Ahn, S.-G. Aging-Related Metabolic Dysfunction in the Salivary Gland: A Review of the Literature. *IJMS* **2021**, *22*, 5835. [CrossRef] [PubMed]

12. Affoo, R.H.; Foley, N.; Garrick, R.; Siqueira, W.L.; Martin, R.E. Meta-Analysis of Salivary Flow Rates in Young and Older Adults. *J. Am. Geriatr. Soc.* **2015**, *63*, 2142–2151. [CrossRef] [PubMed]

13. Teruya, T.; Goga, H.; Yanagida, M. Human Age-Declined Saliva Metabolic Markers Determined by LC-MS. *Sci. Rep.* **2021**, *11*, 18135. [CrossRef] [PubMed]

14. Picca, A.; Calvani, R.; Cesari, M.; Landi, F.; Bernabei, R.; Coelho-Júnior, H.J.; Marzetti, E. Biomarkers of Physical Frailty and Sarcopenia: Coming up to the Place? *Int. J. Mol. Sci.* **2020**, *21*, 5635. [CrossRef]

15. Miyamoto, K.; Hirayama, A.; Sato, Y.; Ikeda, S.; Maruyama, M.; Soga, T.; Tomita, M.; Nakamura, M.; Matsumoto, M.; Yoshimura, N.; et al. A Metabolomic Profile Predictive of New Osteoporosis or Sarcopenia Development. *Metabolites* **2021**, *11*, 278. [CrossRef] [PubMed]

16. Bosman, P.; Pichon, V.; Acevedo, A.C.; Chardin, H.; Combes, A. Development of Analytical Methods to Study the Salivary Metabolome: Impact of the Sampling. *Anal. Bioanal. Chem.* **2022**, *414*, 6899–6909. [CrossRef]

17. Tautenhahn, R.; Patti, G.J.; Rinehart, D.; Siuzdak, G. XCMS Online: A Web-Based Platform to Process Untargeted Metabolomic Data. *Anal. Chem.* **2012**, *84*, 5035–5039. [CrossRef]

18. Montenegro-Burke, J.R.; Guijas, C.; Siuzdak, G. METLIN: A Tandem Mass Spectral Library of Standards. *Methods Mol. Biol.* **2020**, *2104*, 149–163. [CrossRef]

19. Papadopoulou, S.K. Sarcopenia: A Contemporary Health Problem among Older Adult Populations. *Nutrients* **2020**, *12*, 1293. [CrossRef]

20. Tanaka, S.; Machino, M.; Akita, S.; Yokote, Y.; Sakagami, H. Changes in Salivary Amino Acid Composition during Aging. *In Vivo* **2010**, *24*, 853–856.

21. Antonelli, T.; Carlà, V.; Lambertini, L.; Moroni, F.; Bianchi, C. Pyroglutamic Acid Administration Modifies the Electrocorticogram and Increases the Release of Acetylcholine and GABA from the Guinea-Pig Cerebral Cortex. *Pharmacol. Res. Commun.* **1984**, *16*, 189–197. [CrossRef]

22. Kaduce, T.L.; Spector, A.A.; Bar, R.S. Linoleic Acid Metabolism and Prostaglandin Production by Cultured Bovine Pulmonary Artery Endothelial Cells. *Arterioscler. Off. J. Am. Heart Assoc. Inc.* **1982**, *2*, 380–389. [CrossRef] [PubMed]

23. Dallmann, R.; Viola, A.U.; Tarokh, L.; Cajochen, C.; Brown, S.A. The Human Circadian Metabolome. *Proc. Natl. Acad. Sci. USA* **2012**, *109*, 2625–2629. [CrossRef] [PubMed]

24. Sakai, T.; Yamamoto, K.; Yokota, H.; Hakozaki-Usui, K.; Hino, F.; Kato, I. Rapid, Simple Enzymatic Assay of Free L-Fucose in Serum and Urine, and Its Use as a Marker for Cancer, Cirrhosis, and Gastric Ulcers. *Clin. Chem.* **1990**, *36*, 474–476. [CrossRef]

25. Yamauchi, M.; Kimura, K.; Maezawa, Y.; Ohata, M.; Mizuhara, Y.; Hirakawa, J.; Nakajima, H.; Toda, G. Urinary Level of L-Fucose as a Marker of Alcoholic Liver Disease. *Alcohol. Clin. Exp. Res.* **1993**, *17*, 268–271. [CrossRef] [PubMed]

26. Liebsch, C.; Pitchika, V.; Pink, C.; Samietz, S.; Kastenmüller, G.; Artati, A.; Suhre, K.; Adamski, J.; Nauck, M.; Völzke, H.; et al. The Saliva Metabolome in Association to Oral Health Status. *J. Dent. Res.* **2019**, *98*, 642–651. [CrossRef] [PubMed]

27. Losman, J.-A.; Koivunen, P.; Kaelin, W.G. 2-Oxoglutarate-Dependent Dioxygenases in Cancer. *Nat. Rev. Cancer* **2020**, *20*, 710–726. [CrossRef]

28. Sugimoto, M.; Wong, D.T.; Hirayama, A.; Soga, T.; Tomita, M. Capillary Electrophoresis Mass Spectrometry-Based Saliva Metabolomics Identified Oral, Breast and Pancreatic Cancer-Specific Profiles. *Metabolomics* **2009**, *6*, 78–95. [CrossRef]

29. Tsuruoka, M.; Hara, J.; Hirayama, A.; Sugimoto, M.; Soga, T.; Shankle, W.R.; Tomita, M. Capillary Electrophoresis-Mass Spectrometry-Based Metabolome Analysis of Serum and Saliva from Neurodegenerative Dementia Patients. *Electrophoresis* **2013**, *34*, 2865–2872. [CrossRef]