Next Article in Journal
Anti-Melanogenic Effects of Fractioned Cynanchum atratum by Regulation of cAMP/MITF Pathway in a UVB-Stimulated Mice Model
Previous Article in Journal
Structural Transition of the Nucleosome during Transcription Elongation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Is There a Difference in the Proteomic Profile of Stimulated and Unstimulated Saliva Samples from Pregnant Women with/without Obesity and Periodontitis?

by
Gerson Aparecido Foratori-Junior
1,*,
Talita Mendes Oliveira Ventura
2,
Larissa Tercilia Grizzo
2,
Bruno Gualtieri Jesuino
1,
Ana Virgínia Santana Sampaio Castilho
1,
Marília Afonso Rabelo Buzalaf
2,* and
Silvia Helena de Carvalho Sales-Peres
1,*
1
Department of Pediatric Dentistry, Orthodontics and Public Health, Bauru School of Dentistry, University of São Paulo, Bauru 17012-901, Brazil
2
Department of Biological Sciences, Bauru School of Dentistry, University of São Paulo, Bauru 17012-901, Brazil
*
Authors to whom correspondence should be addressed.
Cells 2023, 12(10), 1389; https://doi.org/10.3390/cells12101389
Submission received: 2 March 2023 / Revised: 10 May 2023 / Accepted: 12 May 2023 / Published: 14 May 2023

Abstract

:
This study aimed to compare the proteomic profile of stimulated and unstimulated saliva samples from pregnant women with/without obesity and periodontitis. Pregnant women were allocated into four groups: with obesity and periodontitis (OP); with obesity but without periodontitis (OWP); with normal BMI but with periodontitis (NP); with normal BMI and without periodontitis (NWP). Stimulated saliva (SS) and unstimulated saliva (US) samples were collected, and salivary proteins were extracted and individually processed by proteomic analysis (nLC-ESI-MS/MS). Proteins involved with the immune response process, antioxidant activity, and retina homeostasis were decreased or absent in SS samples from all groups (i.e., Antileukoproteinase, Lysozyme C, Alpha-2-macroglobulin-like protein 1, Heat shock proteins—70 kDa 1-like, 1A, 1B, 6, Heat shock-related 70 kDa protein 2, Putative Heat shock 70 kDa protein 7, Heat shock cognate 71 kDa). Additionally, proteins related to the carbohydrate metabolic process and glycolytic and glucose metabolic process were absent in SS, mainly from OP and OWP (i.e., Frutose-bisphosphate aldose A, Glusoce-6-phosphate isomerase, Pyruvate kinase). Saliva stimulation decreased important proteins involved with immune response and inflammation process in all groups. Unstimulated salivary samples seem to be the best choice for the proteomic approach in pregnant women.

1. Introduction

Saliva is a biological fluid composed of more than 99% water and less than 1% protein, electrolytes, and other low molecular weight components [1]. In addition to providing biomarkers for the diagnosis of local and systemic diseases, saliva plays a key role in protecting oral tissues. Saliva collection is a non-invasive and pain-free method; therefore, several studies have sought to better understand the diagnostic potential of saliva through the identification of biomarkers, whether proteins or metabolites. Proteomic approaches based on mass spectrometry have been applied in several fields of biomedical research to identify new biomarkers [2,3,4].
Obesity is one of the major public health problems in low-, middle-, and high-income countries. The 2017 global nutrition report showed that two billion adults are overweight/obese worldwide [5]. Obesity is known to start a broad inflammatory response through the production of cytokines (interleukins, IL-1, IL-6, IL-8, and Tumor Necrosis Factor alpha, TNF-α), adipokines (leptin, adiponectin, resistin, and inhibitors of plasminogen activator-1), and other bioactive substances, such as Reactive Oxygen Species (ROS), by adipose tissue [6,7,8]. Thus, there is a greater chance of individuals with overweight/obesity to present inflammation in the oral cavity, even in the presence of a small amount of biofilm. Previous studies showed the mechanism by which obesity is linked to an increased inflammatory response to sub and supragingival biofilm, associating obesity with the high prevalence of periodontitis [9].
Evidence shows that both obesity and periodontitis are also associated with low salivary flow rates [10,11]. In the same way, although there is evidence that in the first trimester of pregnancy there is an increase in salivary flow, pregnant women during the second and third trimesters showed lower salivary flow rates when compared to non-pregnant women [12]. In addition, during the third trimester, pregnant women are more prone to periodontal inflammation due to high levels of gestational hormones [13]. In previous studies by our team, we performed proteomics and metabolomics analyses of unstimulated saliva samples from pregnant women with/without obesity and periodontitis during the third trimester of pregnancy [14,15]. One of the challenges was the difficulty in collecting saliva from pregnant women due to frequent nausea and, mainly, the reduced flow in some patients.
Collecting stimulated instead of unstimulated saliva would facilitate clinical practice, both due to the reduced collection time and the ease of obtaining an adequate amount for analysis. However, it is important to assess whether there are no losses related to protein identification in performing proteomic analysis of stimulated saliva in pregnant women with/without obesity and periodontitis. Therefore, this study aimed to compare the proteomic profile of stimulated and unstimulated saliva samples from pregnant women with/without obesity and periodontitis. The null hypothesis was that there were no differences in protein identification in the proteomic analysis of stimulated and unstimulated saliva samples from pregnant women with/without obesity and periodontitis.

2. Materials and Methods

This observational, cross-sectional, and analytical study adhered to the STROBE (Strengthening the Reporting of Observational Studies in Epidemiology) standards [16].

2.1. Ethical Statement

This study was approved by the Internal Research Ethics Committee of the Bauru School of Dentistry, University of São Paulo (CAAE 06624519.3.0000.5417), being in accordance with the protocol established by the Declaration of Helsinki (published in 1975 and revised in 2013). Following the submission of a formal written consent form, individuals were included.

2.2. Sample Selection

Sampling recruitment method, as well as inclusion and exclusion criteria, were described in detail elsewhere [14]. Initially, 126 pregnant women (27th–39th gestational week), 18–40 years old, were selected from the Primary Health Care Units of Bauru, São Paulo, Brazil, between December 2020 and March 2021, but 49 were excluded due to the following reasons: stage I of periodontitis (n = 20); overweight (intermediary BMI, 25–30 kg/m2, between eutrophics and obesity; n = 10); arterial hypertension during pregnancy (n = 5); Gestational Diabetes Mellitus (n = 4); smoker (n = 3); SARS-CoV-2 infection (n = 3); multiple tooth loss (n = 2); twin pregnancy (n = 1); and underweight (BMI < 18.5 kg/m2; n = 1). These exclusion criteria were adopted to ensure a homogenous sample, avoiding bias in the interpretation of the proteomic analysis. Pregnant women were allocated into four groups according to the presence of obesity and periodontitis: with obesity and periodontitis (OWP = 11), with obesity but without periodontitis (OWP = 27), with normal BMI but with periodontitis (NP = 10), and with normal BMI and without periodontitis (NWP = 29). According to previous in vivo individual salivary proteomic analysis by mass spectrometry [1,14,17], the sample was randomized, and 10 individuals from each group were selected for the final sample (OP = 10; OWP = 10; NP = 10; NWP = 10). As mentioned, in vivo individual salivary proteome mass spectrometry analyses were used to determine the sample size [1,14,17].

2.3. Grouping Variables

Grouping method was described elsewhere [14]. Groups OP and OWP were composed of individuals with BMI equal to or higher than 30.00 kg/m2, while groups NP and NWP were composed of individuals with BMI between 18.50–24.99 kg/m2 [14,15,18,19,20,21,22,23,24].
Periodontal evaluation was performed by one calibrated dentist (kappa = 0.95), who registered the probing pocket depth (PPD) and clinical attachment level (CAL)/attachment loss (AL). Periodontitis classifications and categorizations according to severity described by Tonetti and collaborators (2018) [25] were adopted for sampling. To ensure a more homogenous sample for proteomic analysis, participants classified as having stages I and IV of periodontitis were not considered. To avoid unnecessarily exposing pregnant women to X-rays, only clinical characteristics were used to determine the severity of periodontitis.

2.4. Saliva Collection

Saliva collection method was described elsewhere [14]. In summary, in consideration of the circadian rhythm, saliva collection was conducted in the morning (09:00–11:00) and was placed before the periodontal examination. After rinsing their mouths for 1 min with 5 mL of deionized water, the patients were asked to passively drool for 10 min into a sterilized plastic falcon tube (50 mL) immersed in ice (unstimulated whole-mouth saliva—US) [1,14,17]. For the stimulated saliva (SS) collection, patients were instructed to chew a sterilized parafilm gum for 11 min. The saliva produced in the first minute was discarded, and the whole amount produced in the following 10 min was put into a sterilized plastic falcon tube (50 mL) immersed in ice [24]. Saliva volume and flow rate were recorded.
Saliva was centrifuged at 4500× g for 15 min at 4 °C immediately following each collection to clear out all debris. Each sample’s supernatant was taken and stored at −80 °C in a freezer until the time of the proteome analysis.

2.5. Sample Preparation for Proteomic Analysis

Sample preparation for proteomic analysis was previously described in detail [1,14,17]. In summary, all samples (1000 µL each) were analyzed individually, and the sample preparation was divided into seven steps: (1) Extraction: proteins were extracted using a solution (1000 µL for each sample) containing 6 M urea, 2 M thiourea in 50 mM NH4HCO3, pH 7.8. Samples were vortexed three times for 10 min at 4 °C, sonicated for 5 min, and centrifuged at 20,817× g at 4 °C for 10 min. (2) Concentration: samples were concentrated to a volume of around 150 L in Amicon tubes (Amicon Ultra-15 Centrifugal Filter Units—Merck Millipore®, Tullagreen, County Cork, Ireland). (3) Reduction and alkylation: proteins were reduced with 5 mM dithiothreitol for 40 min at 37 °C and alkylated with 10 mM iodoacetamide for 30 min in the dark. (4) Digestion: samples were digested for 14 h at 37 °C by addition of 2% trypsin (w/w) (Thermo Scientific Pierce Trypsin Protease, Rockford, IL, USA). Digestion was stopped by the addition of 10 μL of 5% Trifluoroacetic acid. (5) Desalination and purification: samples were desalted and purified using C18 Spin columns (Thermo Scientific, Rockford, IL, USA). (6) Quantification: an aliquot of 1 μL was taken from each sample for total protein quantification by the Bradford method (Bio-Rad, Hercules, CA, USA). (7) Resuspension: samples were then resuspended in a solution containing 3% acetonitrile and 0.1% formic acid and subjected to Mass Spectrometry (nanoLC-ESI-MS/MS) (Waters, Manchester, New Hampshire, UK).

2.6. Proteomic Analysis by nLC-ESI-MS/MS

Proteomic analysis method was also previously described in detail [1,14,17]. The analysis of peptides was performed in a nanoACQUITY UPLC system (Waters, Manchester, New Hampshire, UK) coupled with a Xevo Q-TOF G2 mass spectrometer (Waters, Manchester, New Hampshire, UK). As aforementioned, for the proteomic analysis of saliva, the samples were analyzed individually. The individual analysis in the mass spectrometer is considered a highly reliable method since the system performs the technical triplicate of each sample. The nanoACQUITY UPLC system is equipped with a Trap Column, responsible for chromatography (Trap Column: 100Å, 5 µm, 180 µm × 200 mm) previously equilibrated with 99.9% phase A (0.1% formic acid in water) at a flow of 5 µL/min and an HSS T3 M-Class type column (analytical column; Acquity UPLC HSS T3 M-Class column 75 μm × 150 mm; 1.8 μm) (Waters, Manchester, New Hampshire, United Kingdom), previously equilibrated with 93% mobile phase A and mobile phase B (0.1% formic acid in ACN). Peptides were separated by a linear gradient of 7–85% mobile phase B for 70 min with 0.35 μL/min flow rate; the column temperature was maintained at 45 °C. In order to collect data using the MSE method at high energy (19–45 V), which enables data capture of both precursor and fragment ions in a single injection, the equipment was operated in positive ionic nanoelectrospray mode. Data acquisition scan range was 50–2000 Da. The lockspray was run with a solution of [Glu1] fibrinopeptide (1 pmol/L) at a flow rate of 0.5 L/min in order to assure precision and reproducibility.
ProteinLynx GlobalServer (PLGS) version 3.0.3 software (Waters Co., Manchester, UK) was used to process and search the continuous LC-MSE data. PLGS has a model of data and uses Bayesian Markov Chain Monte Carlo methods to explore the posterior probability density of the model parameters, accumulating the relevant statistics as the method proceeds. Thus, proteins were identified using the software’s ion counting algorithm, and a search was performed on the Homo sapiens database (reviewed only, UniProtKB/Swiss-Prot) downloaded in May 2022 by UniProtKB (http://www.uniprot.org/ accessed on 25 May 2022). UniProt was used to analyze each protein by its access number. Repeated and reverse proteins, as well as fragments, were excluded. All proteins identified with a confidence level greater than 95% were included in the quantitative analysis. The difference in expression between stimulated and unstimulated saliva samples for each group was expressed as p < 0.05 for the down-regulated proteins and 1-p > 0.95 for the up-regulated proteins [1,14,17]. The difference in expression between SS and US for each group was analyzed by t test (p < 0.05). The relevant comparations were performed (1-SS versus US for OP; 2-SS versus US for OWP; 3-SS versus US for NP; 4-SS versus US for NWP). Proteins were only considered uniquely expressed in the groups when they appeared in at least two of the three replicates of those individuals in each group.

2.7. Statistical Analysis and Bioinformatics

Although the sample size has been based on previous studies, the sample size was also calculated with G*Power version 3.1.9.6 (Heinrich-Heine-Universität Düsseldorf, Bavaria, Germany), utilizing information from our prior experiment [14], considering α = 0.05 and 1 − β = 0.8. Effect size was estimated to be 2.18. The estimated number of samples was 4/group. We included 10 volunteers in each group, and all samples were individually analyzed.
In light of the clinical parameters, the Kolmogorov–Smirnov test was used to determine the normality of the variables. Analysis of variance (ANOVA) with Scheffé was then used for quantitative variables with a normal distribution, while Kruskal–Wallis with Dunn was used for quantitative variables without a normal distribution.
In proteomic analysis, the difference in expression between stimulated and unstimulated saliva samples for each group was expressed as p < 0.05 for the down-regulated proteins and 1-p > 0.95 for the up-regulated proteins [1,14,17]. The difference in expression between SS and US for each group was analyzed by t test (p < 0.05). The relevant comparations were performed (1-SS versus US for OP; 2-SS versus US for OWP; 3-SS versus US for NP; 4-SS versus US for NWP).
ClueGo® plugins of the Cytoscape® 3.9.1 Software (Institute of Systems Biology, Seattle, WA, USA) were used to analyze the protein categories based on gene ontology (GO) annotation of the broad biological process, molecular function, immune system, and cell component. The functional distribution of proteins identified with differential expression (up- and down-regulated) in the comparison between SS and US for each group was carried out. Terms of significance (κ = 0.04) and distribution were according to the percentage of the number of associated genes. The mass spectrometric proteomic data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the data set identifier PXD040373.
For the protein interaction networks analysis, STRING® database (https://string-db.org/cgi/network.pl accessed on 15 February 2022) was accessed. To understand the biological processes involved only in the analysis of unstimulated saliva, unique proteins found in US samples of each group were analyzed in the STRING® database, establishing the interaction between those unique proteins present in unstimulated saliva of each group during pregnancy.

3. Results

The mean age (±SD) of the sample was 27.63 ± 4.57 years old (OP = 32.10 ± 5.02; OWP = 25.00 ± 5.09; NP = 26.00 ± 4.18; NWP = 28.40 ± 4.00). The pre-pregnancy BMI mean (±SD) was 33.6 (±3.16), 34.7 (±3.80), 21.3 (±2.59), and 23.3 (±1.71) kg/m2 for OP, OWP, NP, and NWP, respectively, while the pregnancy BMI mean during the third trimester was 36.0 (±4.40), 37.3 (±3.65), 24.5 (±2.51), and 27.3 (±3.02) kg/m2, respectively. Groups were homogeneous regarding salivary flow and the total mean amount of protein recovered (µg) both for SS and USA, showing no difference in these parameters (Table 1).
For the quantitative proteomic analysis, in the comparison between SS and US samples from the OP group, the total number of proteins identified was 108 and 192, respectively, among which 100 proteins were common to both types of saliva (Figure 1A). Eight proteins were identified exclusively in SS, while 92 proteins were uniquely identified in US (Table S1). Regarding the differentially expressed proteins in the OP group, 37 and 47 proteins were increased and decreased, respectively, in SS (Table S1). Among the main up-regulated proteins in SS were Alpha-1-antitrypsin (increased 5-fold); Neutrophil defensin 3 (increased 4-fold); Matrix metalloproteinase-9; Histatin-3 (increased 3-fold); 2 isoforms of Basic salivary proline-rich protein; 2 isoforms of POTE ankyrin domain; Carbonic anhydrase 6; Actin, alpha cardiac muscle 1; and 2 isoforms of Hemoglobin (increased 2-fold). Among the main down-regulated proteins in SS were LINE-1 type transposase domain-containing protein 1 (decreased 13-fold); Haptoglobin-related protein (decreased 10-fold); 3 isoforms of Hemoglobin (decreased 8-fold); Immunoglobulin kappa light chain; Haptoglobin (decreased 5-fold); Deleted in malignant brain tumors 1 protein; Immunoglobulin heavy constant gamma 2; Glucose-6-phosphate isomerase (decreased 4-fold); Lactoperoxidase; Protein S100-A9; Transaldolase; Alpha-2-macroglobulin; Immunoglobulin heavy constant gamma 4 (decreased 3-fold); Protein S100-A8; and Thioredoxin (decreased almost 3-fold) (Table 2).
Figure 2 shows the functional analysis for the comparison between SS and US in the OP group. Among them, we would like to highlight the categories with the highest percentages of genes in the biological process and immune system. For the biological process, the categories were defense response to a bacterium (21.01%), humoral immune response (19.33%), retina homeostasis (14.29%), and antioxidant activity (10.08%). For the immune system, the categories were complement activation (50%) and antibacterial humoral response (26.67%) (Figure 2).
In the comparison between SS and US from the OWP group, the total number of proteins identified was 94 and 185, respectively, among which 91 proteins were common to both types of saliva (Figure 1B). Three proteins were identified exclusively in SS, while ninety-four proteins were uniquely identified in US (Table S2). Regarding the differentially expressed proteins in the OWP group, 14 and 55 proteins were increased and decreased, respectively, in SS (Table S2). Among the main up-regulated proteins in SS were Statherin (increased 18-fold); 2 isoforms of Immunoglobulin (increased 3-fold); Salivary acidic proline-rich phosphoprotein ½ (increased 3-fold); and Histatin-3 (increased almost 3-fold). Among the main down-regulated proteins in SS were Transaldolase; Profilin-1 (decreased 10-fold); Myeloblastin (decreased 7-fold); 6-phosphogluconate dehydrogenase, decarboxylating (decreased 6-fold); Pyruvate kinase PKM (decreased 5-fold); 2 isoforms of Neutrophil defensin; 2 isoforms of Basic salivary proline-rich protein (decreased 4-fold); 7 isoforms of Immunoglobulin; Glyceraldehyde-3-phosphate dehydrogenase; 2 isoforms of Protein S100; Lysozyme C (decreased 3-fold) (Table 3).
Figure 3 shows the functional analysis for the comparison between SS and US in the OWP group. Among them, we would like to highlight the categories with the highest percentages of genes in the biological process and immune system. For the biological process, the categories were humoral immune response (29.79%), defense response to bacterium (27.66%), and antimicrobial humoral response (14.89%). For the immune system, the categories were humoral immune response mediated by circulating immunoglobulin (39.02%), antimicrobial humoral response (34.15%), and antimicrobial humoral immune response mediated by antimicrobial peptide (12.20%) (Figure 3).
In the comparison between SS and US from the NP group, the total number of proteins identified was 93 and 171, respectively, among which 91 proteins were common to both types of saliva (Figure 1C). Two proteins were identified exclusively in SS, while eighty proteins were uniquely identified in US (Table S3). Regarding the differentially expressed proteins in the NP group, 36 and 42 proteins were increased and decreased, respectively, in SS (Table S3). Among the main up-regulated proteins in SS were Haptoglobin (increased 9-fold); 2 isoforms of Basic salivary proline-rich protein; Histatin-3; Alpha-1-antitrypsin; Apolipoprotein A-I (increased 4-fold); Hemopexin; Cystatin-B (increased 3-fold); and BPI fold-containing family B member 2 (increased almost 3-fold). Among the main down-regulated proteins in SS were Lysozyme C (decreased 15-fold); Protein S100-A9 (decreased 7-fold); Beta-actin-like protein 2; Proline-rich protein 4; Mucin-7 (decreased 5-fold); 2 isoforms of Immunoglobulin; Hemoglobin subunit delta (decreased 3-fold) (Table 4).
Figure 4 shows the functional analysis for the comparison between SS and US in the NP group. Among them, we would like to highlight the categories with the highest percentages of genes in the biological process and immune system. For the biological process, the categories were retina homeostasis (26.92%) and defense response to bacterium (26.92%). For the immune system, the categories were antimicrobial humoral response (44%) and humoral immune response mediated by circulating immunoglobulin (36%) (Figure 4).
In the comparison between SS and US from the NWP group, the total number of proteins identified was 83 and 164, respectively, among which 81 proteins were common to both types of saliva (Figure 1D). Two proteins were identified exclusively in SS, while eighty-three proteins were uniquely identified in US (Table S4). Regarding the differentially expressed proteins in the NWP group, 18 and 48 proteins were increased and decreased, respectively, in SS (Table S4). Among the main up-regulated proteins in SS were Submaxillary gland androgen-regulated protein 3B; Beta-2-microglobulin (increased almost 3-fold); Protein S100-A9; Mucin-7; Immunoglobulin lambda-1 light chain (increased 2-fold). Among the main down-regulated proteins in SS were Putative lipocalin 1-like protein 1 (decreased 10-fold); Lipocalin-1 (decreased 9-fold); Carbonic anhydrase 6 (decreased almost 4-fold); 3 isoforms of POTE ankyrin domain; and 3 isoforms of Immunoglobulin (Table 5).
Figure 5 shows the functional analysis for the comparison between SS and US in the NWP group. Among them, we would like to highlight the categories with the highest percentages of genes in the biological process and immune system. For the biological process, the categories were defense response to bacterium (25.88%), retina homeostasis (24.71%), and humoral immune response (23.53%). For the immune system, the categories were positive regulation of B cell activation (59.26%) and antimicrobial humoral response (33.33%) (Figure 5).
Figure 6 shows the interaction networks among unique proteins identified in US from OP (Figure 6A), OWP (Figure 6B), NP (Figure 6C), and NWP (Figure 6D). In the OP group, dark blue nodes were related to protein refolding; green nodes were related to antimicrobial humoral response; pink, light blue, and red nodes were related to leukocyte activation involved in immune response, immune response, and immune system process, respectively; yellow nodes were related to carbohydrate metabolic process (Figure 6A). In the OWP group, red nodes were related to protein refold; dark blue and yellow nodes were related to glycolytic and glucose metabolic process, respectively; and green and pink nodes were related to leukocyte activation involved in immune response and immune response, respectively (Figure 6B). In the NP group, red nodes were related to protein refolding; green nodes were related to neutrophil degranulation; yellow nodes were related to the oxidation-reduction process; and dark blue nodes were related to immune response (Figure 6C). In the NWP group, red nodes were related to protein refolding; dark and light blue nodes were related to glycolytic and oxidation-reduction processes, respectively; dark green, pink, orange, and yellow nodes were related to antimicrobial humoral response, leukocyte activation involved in immune response, immune response, and neutrophil degranulation, respectively; and light green nodes were related to platelet degranulation (Figure 6D).

4. Discussion

The reduction in salivary flow in patients with obesity [10], periodontitis [11], and during the third trimester of pregnancy [12] was previously reported. Considering that collecting stimulated instead of unstimulated saliva would facilitate clinical practice, both due to the reduced collection time and the ease of obtaining an adequate amount for analysis, we sought to underly the proteomic profile of stimulated and unstimulated saliva samples from pregnant women with/without obesity and periodontitis. Our findings highlighted significant alterations in the proteomic profile of saliva when adopting different collection methods. In general, saliva stimulation decreased important proteins involved with immune response and inflammation process in all groups. Unstimulated salivary samples seem to be the best choice for the proteomic approach in pregnant women. Therefore, the null hypothesis of this study was rejected.
The proinflammatory cytokines derived from adipocytes and macrophages accumulated in the adipose tissue in obese patients can negatively affect the function of the salivary glands due to low-grade chronic inflammation in the gland [26], explaining the association between obesity and low salivary flow. Similarly, the plausible mechanisms explaining the association between low salivary flow rate and periodontitis might be due to also the increased proinflammatory cytokine levels, such as IL-2, IL-17, and TNF, triggering glandular damage and hyposalivation. Furthermore, these inflammatory cytokines accelerate the production of ROS, which, in turn, might cause structural changes in salivary gland tissue as those that occur during physiological aging [26].
As reported by Golatowski et al. (2013), the determination of the variability in whole saliva proteome is a prerequisite for the development of saliva as a diagnostic and/or prognostic human biomarker fluid. In that study, authors compared volume, protein concentrations, and proteome profile among saliva samples collected by drooling, Salivette®, and paraffin gum [27]. Our findings are in line with that study since an expected difference was found regarding the volume and flow rate among SS ad US samples, but no difference was found in protein concentrations. Golatowski et al. highlighted that passive drooling, paraffin gum, and Salivette® each allow similar coverage of the whole saliva proteome, but the specific proteins observed depended on the collection approach [27]. In this study, specific proteins related to the immune process and inflammatory response were found mainly in unstimulated saliva samples for all groups. It is important to highlight that our main objective here was to compare stimulated and unstimulated saliva samples for each group, with no intergroup comparison for quantitative proteomic analysis, and these results were discussed in detail below.
Among the proteins elevated or exclusively present in SS samples of this study, only Alpha-amylase 1C and Pancreatic alpha-amylase were similar in all groups (Tables S1–S4). Among the proteins that were decreased in SS or exclusively present in US samples, 29 proteins were similar in all groups, such as the following: 6-phosphogluconate dehydrogenase, decarboxylating; Acyl-CoA-binding protein; Alpha-2-macroglobulin-like protein 1; Antileukoproteinase; Cysteine-rich secretory protein 3; Endoplasmic reticulum chaperone BiP; Fructose-bisphosphate aldolase A; Glutathione S-transferase P; 6 isoforms of Heat shock protein (Heat shock 70 kDa protein 1-like, Heat shock 70 kDa protein 1A, 1B, and 6, Heat shock cognate 71 kDa protein, and Putative heat shock 70 kDa protein 7); 8 isoforms of Immunoglobulin (Immunoglobulin alpha-2 heavy chain, Immunoglobulin gamma-1 heavy chain, Immunoglobulin heavy constant alpha 1 and 2, Immunoglobulin heavy constant gamma 1, 2, 3, and 4); Lactoperoxidase; Lysozyme C; Polymeric immunoglobulin receptor; Protein LEG1 homolog; SH3 domain-binding glutamic acid-rich-like protein 3; Thioredoxin; and Transaldolase.
These proteins are, in general, involved with the immune response process, antioxidant activity, and retina homeostasis (i.e., Antileukoproteinase; Lysozyme C; Alpha-2-macroglobulin-like protein 1; Heat shock proteins—70 kDa 1-like, 1A, 1B, 6; Heat shock-related 70 kDa protein 2; Putative Heat shock 70 kDa protein 7; Heat shock cognate 71 kDa), showing that saliva stimulation decreased important proteins involved with immune response and inflammation process in all groups. Here, special attention should be given to Antileukoproteinase and Lysozyme C. Antileukoproteinase modulates the inflammatory and immune responses after bacterial infection and down-regulates responses to bacterial lipopolysaccharide (LPS) (UNIPROT). Similarly, Lysozyme C increases the activity of immunoagents through the monocyte-macrophage system by acting primarily as a bacteriolytic agent in tissues and bodily fluids. Our findings are in accordance with a recent study that compared the proteomic profile among stimulated and unstimulated saliva samples from head and neck cancer patients treated by radiotherapy [28]. In that study, authors reported that proteins involved with apoptosis, antibacterial, and acid-resistance were decreased in stimulated saliva in comparison to unstimulated saliva, indicating that unstimulated salivary flow seems to be the best alternative to search for biomarkers.
Other proteins that were also reduced in SS samples or exclusively present in US in all groups included 6-phosphogluconate dehydrogenase, decarboxylating; Fructose-bisphosphate aldolase A; and Transaldolase (Tables S1–S4). Thus, if only stimulated saliva samples were analyzed, there would be a lack of understanding related to the carbohydrate metabolic process and glycolytic and glucose metabolic processes in which these proteins are involved. In addition to its important role in glycolysis and gluconeogenesis, Fructose-bisphosphate aldolase A has recently been found to have non-glycolytic functions, such as binding to host cell receptors to promote invasion or activating plasminogen to possibly modify host hemostasis and improve survival [29]. Transaldolase, in turn, is important for the balance of metabolites in the pentose-phosphate pathway, and its involvement in oxidative stress and apoptosis, in multiple sclerosis, and in cancer has been discussed [30]. Additionally, 6-phosphogluconate dehydrogenase, decarboxylating catalyzes the oxidative decarboxylation of 6-phosphogluconate to ribulose 5-phosphate and CO2, with concomitant reduction of NADP to NADPH. In addition to these proteins, other proteins related to the metabolic process were also reduced or absent in samples of stimulated saliva in the groups of pregnant women with obesity (OP and OWP), such as the following: Transketolase; Pyruvate kinase; Glucose-6-phosphate isomerase; Triosephosphate isomerase; L-lactate dehydrogenase A-like 6A; L-lactate dehydrogenase C chain; and Fatty acid-binding protein (Tables S1 and S2). Therefore, in the aforementioned context, unstimulated saliva collection seems also to be a better protocol for proteomic analysis in pregnant women.
In this study, there also were some proteins decreased or absent in SS samples and in common only in periodontitis cases (OP and NP groups), for instance, ATP-dependent RNA helicase DDX55; 12 isoforms of Immunoglobulin (Immunoglobulin heavy constant mu, Immunoglobulin heavy variable 3-23, 3-30, 3-30-3, 3-30-5, 3-33, 3-53, 3-66, 3-74, Immunoglobulin kappa variable 3-11, 3D-11, Immunoglobulin mu heavy chain); Mucin-2; Plastin-1; and Submaxillary gland androgen-regulated protein 3B (Tables S1 and S3). We would especially like to highlight Plastin-1, Mucin-2, and Submaxillary gland androgen-regulated protein 3B. Plastin-1 was an example of a protein exclusively identified in US samples from periodontitis groups (OP and NP). As we mentioned in our recent study [14], to the best of our knowledge, there has never been any proof linking Plastin-1 with periodontitis. Recently, it was discovered that Plastin-1 and Plastin-3, whose mutations cause X-linked osteoporosis, are extremely similar. Furthermore, it was discovered that Plastin-1 regulates intracellular Ca2+ to encourage osteoblast development [31]. In our earlier study [14], we proposed that the existence of Plastin-1 in pregnant women with periodontitis may represent a compensation mechanism for the high levels of inflammation and bone loss, functioning to maintain bone homeostasis. Mucin-2 and Submaxillary gland androgen-regulated protein 3B are associated with periodontal diseases since the first provides a protective, lubricating barrier against particles and infectious agents at mucosal surfaces (UNIPROT), while the latter seems to bind to lipopolysaccharide of P. gingivalis, acting in promoting angiogenesis and establishing microvasculature [32]. Therefore, these proteins are potential biomarkers of periodontitis in pregnant women that are elevated or exclusively present in US, corroborating our previously discussed findings that point out that unstimulated saliva collection should be a preferable protocol for proteomic analysis.
It is also important to discuss our results focusing on pregnant women affected by obesity and periodontitis due to the combination of the two inflammatory conditions being of great interest nowadays. In this study, some proteins were elevated in the stimulated saliva of pregnant women from the OP group. Among them, Alpha-1-antitrypsin, Neutrophil defensin 1 and 3, Lactotransferrin, MMP9, and Histatin-3 were involved in the immune effector process. Nevertheless, a higher number of proteins related to not only immune and inflammatory processes but also to the metabolic process were elevated or exclusively present in unstimulated saliva samples from OP. Among proteins involved in the immune process, S100A8, S100A9, BPI fold-containing family A member 2, BPI fold containing family B member 2, Haptoglobin, Deleted in malignant brain tumors 1 protein, Lactoperoxidase, and Complement C3 might be considered antimicrobial agents strongly related to bacterium defense that were elevated in US samples from OP (Table S1); while Heat shock 70 kDa protein 1B, 2, and 6; Heat shock cognate 71 kDa protein; Alpha-1-acid glycoprotein 1; Beta-2-microglobulin; Lysozyme C; and Antileukoproteinase were examples of proteins exclusively expressed in US samples (Figure 6; Table S1). We call special attention to the importance of Protein S100-A8, Heat shock 70 kDa protein 2 and 6, Heat shock 71 kDa protein, and Haptoglobin, which were elevated or exclusively present in US samples, since these proteins were found to be potential biomarkers of the association between obesity and periodontitis during pregnancy in our previous study [14].
Proteins S100-A8 and S100-A9 present proinflammatory, antimicrobial, oxidant-scavenging, and apoptosis-inducing activities. Their proinflammatory activity includes recruitment of leukocytes, promotion of cytokine and chemokine production, and regulation of leukocyte adhesion and migration [33,34,35]. Regarding Heat shock proteins (HSPs), they are associated with the body’s metabolic and catabolic processes in addition to immune response and oxidative stress. HSP70 prevents tissue degradation by acting as a protective factor and inhibiting apoptosis [36]. Previous research indicated that HSPs were positively expressed in the periodontal pockets’ basal layer, indicating that there was an increase in the infiltration of mononuclear inflammatory cells beneath the basal layer [36,37,38]. As a result, periodontal bacteria induce the expression of HSPs in the periodontal cells, which in turn induces the production of proinflammatory cytokines by macrophages and other inflammatory cells, a mechanism that contributes to the destruction of tissue in periodontitis cases [36,37,38]. Haptoglobin has anti-inflammatory and antioxidative properties, acting as a bacteriostatic agent and, indirectly, as an antioxidant. Higher levels of Haptoglobin were associated with host defense in periodontitis cases [39,40]. Thus, these potential biomarkers of obesity plus periodontitis during pregnancy were also dependent on the type of salivary collection. Here, the unstimulated saliva samples also seem to be the best choice for the proteomic approach.
This study has some limitations. Future evaluations should be performed through different gestational trimesters, with various follow-ups after delivery, to better understand the changes in salivary proteomic profile over time, both in stimulated and unstimulated saliva samples. Ideally, future longitudinal studies with larger samples should analyze specific proteins expressed in saliva considering the different stages of periodontal diseases (including gingivitis as well) to ensure a better biological understanding regarding the progression of the diseases. Despite the limitations, to the best of our understanding, this is the first study to map the proteomic profile and to understand the differences among stimulated and unstimulated saliva samples from pregnant women with/without obesity and periodontitis, using an individual label-free quantitative shotgun proteomic analysis. The results of this study contribute to researchers around the world as it demonstrates the best saliva collection protocol that should be adopted for proteomic analysis. Furthermore, our results suggest further investigation of salivary biomarkers. The identification of biomarkers can contribute to the early diagnosis of diseases, favoring public health policies in disease prevention and health promotion since this would reduce the costs that health managers have with the sectors of greater technological density that are necessary for the rehabilitation of diseases in more advanced stages. The development of these health policies aims to achieve universal health coverage for all and to ensure people receive equitable health care and dental services.

5. Conclusions

In conclusion, there were significant differences in the proteomic profile of SS and US samples from pregnant women with/without obesity and periodontitis. Saliva stimulation decreased (or made absent) important proteins involved with immune response and inflammatory process in all groups. Unstimulated salivary samples seem to be the best choice for proteomic analysis in pregnant women. Our findings contribute in an unprecedented way to understanding the differences in the salivary proteomic profile of different flows in pregnant women with/without obesity and periodontitis, suggesting that future studies adopting the proteomic approach in the same target population should focus on the unstimulated salivary analysis.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cells12101389/s1, Table S1: Expression differences and unique proteins identified in stimulated saliva (SS) and unstimulated saliva (US) of pregnant women with obesity and periodontitis (OP); Table S2: Expression differences and unique proteins identified in stimulated saliva (SS) and unstimulated saliva (US) of pregnant women with obesity but without periodontitis (OWP); Table S3: Expression differences and unique proteins identified in stimulated saliva (SS) and unstimulated saliva (US) of pregnant women with normal BMI but with periodontitis (NP); Table S4: Expression differences and unique proteins identified in stimulated saliva (SS) and unstimulated saliva (US) of pregnant women with normal BMI and without periodontitis (NWP).

Author Contributions

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

Funding

This research was funded by the São Paulo Research Foundation (FAPESP), grant number 2018/20626-5, and Coordination of Superior Level Staff Improvement (CAPES), financial code 001. The APC was funded by São Paulo Research Foundation (FAPESP).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee from Bauru School of Dentistry, University of São Paulo (CAAE 06624519.3.0000.5417, protocol code 3.284.822; approved on 17 April 2019).

Informed Consent Statement

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

Data Availability Statement

Data presented in this study are available in the article and Supplementary Materials and made available upon request to the corresponding authors. The mass spectrometric proteomic data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the data set identifier PXD040373.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ventura, T.M.S.; Ribeiro, N.R.; Dionizio, A.S.; Sabino, I.T.; Buzalaf, M.A.R. Standardization of a protocol for shotgun proteomic analysis of saliva. J. Appl. Oral Sci. 2018, 26, e20170561. [Google Scholar] [CrossRef]
  2. Wang, X.; Kaczor-Urbanowicz, K.E.; Wong, D.T. Salivary biomarkers in cancer detection. Med. Oncol. 2017, 34, 7. [Google Scholar] [CrossRef] [PubMed]
  3. Tsuchida, S. Proteome analysis of molecular events in oral pathogenesis and virus: A review with a particular focus on periodontitis. Int. J. Mol. Sci. 2020, 21, 5184. [Google Scholar] [CrossRef] [PubMed]
  4. Sachdeva, R.; Kumar, N.; Brache, V.; Friedland, B.A.; Plagianos, M.; Zhang, S.; Kizima, L.; Cochon, L.; Tabar, A.S.T.; Blanc, A.; et al. New approaches for developing biomarkers of hormonal contraceptive use. Sci. Rep. 2023, 13, 245. [Google Scholar] [CrossRef] [PubMed]
  5. Endalifer, M.L.; Diress, G. Epidemiology, predisposing factors, biomarkers, and prevention mechanism of obesity: A systematic review. J. Obes. 2020, 2020, 6134362. [Google Scholar] [CrossRef]
  6. Furugen, R.; Hayashida, H.; Kitamura, M.; Saito, T. Relationship between adipokines and periodontitis. Jpn. Dent. Sci. Rev. 2010, 46, 159–164. [Google Scholar] [CrossRef]
  7. Huang, L.-Y.; Chiu, C.-J.; Hsing, C.-H.; Hsu, Y.-H. Interferon Family Cytokines in Obesity and Insulin Sensitivity. Cells 2022, 11, 4041. [Google Scholar] [CrossRef]
  8. Ostrowska, L.; Smarkusz-Zarzecka, J.; Gornowicz, A.; Lendzion, K.; Zyśk, B.; Pogodziński, D. Analysis of selected salivary adipokines and cytokines in patients with obesity—A pilot study. Int. J. Mol. Sci. 2023, 24, 4145. [Google Scholar] [CrossRef]
  9. Keller, A.; Rohde, J.F.; Raymond, K.; Heitmann, B.L. Association between periodontal disease and overweight and obesity: A systematic review. J. Periodontol. 2015, 86, 766–776. [Google Scholar] [CrossRef]
  10. Hatipoglu, O.; Maras, E.; Hatipoglu, F.P.; Saygin, A.G. Salivary flow rate, pH, and buffer capacity in the individuals with obesity and overweight; A meta-analysis. Niger. J. Clin. Pract. 2022, 25, 1126–1142. [Google Scholar] [CrossRef]
  11. Vallabhan, C.G.; Sivarajan, S.; Shivkumar, A.D.; Narayanan, V.; Vijayakumar, S.; Indhuja, R.S. Assessment of salivary flow rate in patients with chronic periodontitis. J. Pharm. Bioallied. Sci. 2020, 12, S308–S312. [Google Scholar] [CrossRef]
  12. Migliario, M.; Bindi, M.; Surico, D.; De Pedrini, A.; Minsenti, S.; Pezzotti, F.; Mele, B.; Foglio Bonda, P.L. Changes in salivary flow rate and pH in pregnancy. Eur. Rev. Med. Pharmacol. Sci. 2021, 25, 1804–1810. [Google Scholar] [CrossRef]
  13. Starzyńska, A.; Wychowański, P.; Nowak, M.; Sobocki, B.K.; Jereczek-Fossa, B.A.; Słupecka-Ziemilska, M. Association between maternal periodontitis and development of systematic diseases in offspring. Int. J. Mol. Sci. 2022, 23, 2473. [Google Scholar] [CrossRef]
  14. Foratori-Junior, G.A.; Ventura, T.M.O.; Grizzo, L.T.; Carpenter, G.H.; Buzalaf, M.A.R.; Sales-Peres, S.H.C. Label-free quantitative proteomic analysis reveals inflammatory pattern associated with obesity and periodontitis in pregnant women. Metabolites 2022, 12, 1091. [Google Scholar] [CrossRef] [PubMed]
  15. Foratori-Junior, G.A.; Guennec, A.L.; Fidalgo, T.K.D.S.; Cleaver, L.; Buzalaf, M.A.R.; Carpenter, G.H.; Sales-Peres, S.H.C. Metabolomic profiles associated with obesity and periodontitis during pregnancy: Cross-sectional study with proton nuclear magnetic resonance (1H-NMR)-based analysis. Metabolites 2022, 12, 1029. [Google Scholar] [CrossRef]
  16. Vandenbroucke, J.P.; von Elm, E.; Altman, D.G.; Gøtzsche, P.C.; Mulrow, C.D.; Pocock, S.J.; Poole, C.; Schlesselman, J.J.; Egger, M. Strengthening the Reporting of Observational Studies in Epidemiology (STROBE): Explanation and elaboration. Int. J. Surg. 2014, 12, 1500–1524. [Google Scholar] [CrossRef]
  17. Ventura, T.M.O.; Ribeiro, N.R.; Taira, E.A.; de Lima Leite, A.; Dionizio, A.; Rubira, C.M.F.; Santos, P.S.S.; Buzalaf, M.A.R. Radiotherapy Changes the Salivary Proteome in Head and Neck Cancer Patients: Evaluation before, during, and after Treatment. Clin. Oral Investig. 2022, 26, 225–258. [Google Scholar] [CrossRef] [PubMed]
  18. Fusco, N.D.S.; Foratori-Junior, G.A.; Missio, A.L.T.; Jesuino, B.G.; Sales-Peres, S.H.C. Systemic and oral conditions of pregnant women with excessive weight assisted in a private health system. Int. Dent. J. 2019, 69, 472–479. [Google Scholar] [CrossRef] [PubMed]
  19. Caracho, R.A.; Foratori-Junior, G.A.; Fusco, N.D.S.; Jesuino, B.G.; Missio, A.L.T.; Sales-Peres, S.H.C. Systemic conditions and oral health-related quality of life of pregnant women of normal weight and who are overweight. Int. Dent. J. 2020, 70, 287–295. [Google Scholar] [CrossRef]
  20. Foratori-Junior, G.A.; da Silva, B.M.; da Silva Pinto, A.C.; Honório, H.M.; Groppo, F.C.; de Carvalho Sales-Peres, S.H. Systemic and periodontal conditions of overweight/obese patients during pregnancy and after delivery: A prospective cohort. Clin. Oral Investig. 2020, 24, 157–165. [Google Scholar] [CrossRef]
  21. Foratori-Junior, G.A.; Jesuino, B.G.; Caracho, R.A.; Orenha, E.S.; Groppo, F.C.; Sales-Peres, S.H.C. Association between excessive maternal weight, periodontitis during the third trimester of pregnancy, and infants’ health at birth. J. Appl. Oral Sci. 2020, 28, e20190351. [Google Scholar] [CrossRef] [PubMed]
  22. Jesuino, B.G.; Foratori-Junior, G.A.; Missio, A.L.T.; Mascoli, L.S.; Sales-Peres, S.H.C. Periodontal status of women with excessive gestational weight gain and the association with their newborns’ health. Int. Dent. J. 2020, 70, 396–404. [Google Scholar] [CrossRef]
  23. Foratori-Junior, G.A.; Missio, A.L.T.; Orenha, E.S.; de Carvalho Sales-Peres, S.H. Systemic condition, periodontal status, and quality of life in obese women during pregnancy and after delivery. Int. Dent. J. 2021, 71, 420–428. [Google Scholar] [CrossRef] [PubMed]
  24. Foratori-Junior, G.A.; Mosquim, V.; Buzalaf, M.A.R.; Sales-Peres, S.H.C. Salivary cytokines levels, maternal periodontitis and infants’ weight at birth: A cohort study in pregnant women with obesity. Placenta 2021, 115, 151–157. [Google Scholar] [CrossRef] [PubMed]
  25. Tonetti, M.S.; Greenwell, H.; Kornman, K.S. Staging and grading of periodontitis: Framework and proposal of a new classification and case definition. J. Clin. Periodontol. 2018, 45, S149–S161. [Google Scholar] [CrossRef]
  26. Roa, I.; Del Sol, M. Obesity, salivary glands and oral pathology. Colomb. Med. 2018, 49, 280–287. [Google Scholar] [CrossRef]
  27. Golatowski, C.; Salazar, M.G.; Dhople, V.M.; Hammer, E.; Kocher, T.; Jehmlich, N.; Völker, U. Comparative evaluation of saliva collection methods for proteome analysis. Clin. Chim. Acta 2013, 419, 42–46. [Google Scholar] [CrossRef] [PubMed]
  28. Ventura, T.M.O.; Santos, P.S.S.; Ribeiro, N.R.; de Lima Leite, A.; Taira, E.A.; Dionizio, A.; Rubira, C.M.F.; Buzalaf, M.A.R. Is there difference in the comparative and quantitative salivary proteome between stimulated and unstimulated saliva in head and neck cancer patients treated by radiotherapy? Oral Oncol. 2021, 118, 105315. [Google Scholar] [CrossRef] [PubMed]
  29. Pirovich, D.B.; Da’dara, A.A.; Skelly, P.J. Multifunctional Fructose 1,6-Bisphosphate Aldolase as a therapeutic target. Front. Mol. Biosci. 2021, 8, 719678. [Google Scholar] [CrossRef]
  30. Samland, A.K.; Sprenger, G.A. Transaldolase: From biochemistry to human disease. Int. J. Biochem. Cell. Biol. 2009, 41, 1482–1494. [Google Scholar] [CrossRef]
  31. Wang, L.; Lan, Y.; Du, Y.; Xiang, X.; Tian, W.; Yang, B.; Li, T.; Zhai, Q. Plastin 1 Promotes osteoblast differentiation by regulating intracellular Ca2. Acta Biochim. Biophys. Sin. 2020, 52, 563–569. [Google Scholar] [CrossRef]
  32. Choi, S.; Baik, J.E.; Jeon, J.H.; Cho, K.; Seo, D.G.; Kum, K.Y.; Yun, C.H.; Han, S.H. Identification of Porphyromonas Gingivalis lipopolysaccharide-binding proteins in human saliva. Mol. Immunol. 2011, 48, 2207–2213. [Google Scholar] [CrossRef] [PubMed]
  33. Dey, A.K.; Kumar, B.; Singh, A.K.; Ranjan, P.; Thiruvengadam, R.; Desiraju, B.K.; Kshetrapal, P.; Wadhwa, N.; Bhatnagar, S.; Rashid, F.; et al. Salivary proteome signatures in the early and middle stages of human pregnancy with term birth outcome. Sci. Rep. 2020, 10, 8022. [Google Scholar] [CrossRef]
  34. Kim, H.; Kim, S.; Jeon, S.; Kim, S.; Cho, H.; Choi, Y. Diagnostic and prognostic ability of salivary MMP-9 and S100A8 for periodontitis. J. Clin. Periodontol. 2020, 47, 1191–1200. [Google Scholar] [CrossRef] [PubMed]
  35. Revenstorff, J.; Ludwig, N.; Hilger, A.; Mersmann, S.; Lehmann, M.; Grenzheuser, J.C.; Kardell, M.; Bone, J.; Kötting, N.M.; Marx, N.C.; et al. Role of S100A8/A9 in Platelet-Neutrophil Complex Formation during Acute Inflammation. Cells 2022, 11, 3944. [Google Scholar] [CrossRef] [PubMed]
  36. Motahari, P.; Pourzare Mehrbani, S.; Jabbarvand, H. Evaluation of salivary level of Heat Shock Protein 70 in patients with chronic periodontitis. J. Dent. 2021, 22, 175–179. [Google Scholar] [CrossRef]
  37. Garlet, G.P. Destructive and protective roles of cytokines in periodontitis: A re-appraisal from host defense and tissue destruction viewpoints. J. Dent. Res. 2010, 89, 1349–1363. [Google Scholar] [CrossRef]
  38. Wolf, M.; Marciniak, J.; Lossdörfer, S.; Kirschneck, C.; Brauner, I.; Götz, W.; Jäger, A. Role of HSP70 protein in human periodontal ligament cell function and physiology. Ann. Anat. 2019, 221, 76–83. [Google Scholar] [CrossRef]
  39. Haigh, B.J.; Stewart, K.W.; Whelan, J.R.; Barnett, M.P.; Smolenski, G.A.; Wheeler, T.T. Alterations in the salivary proteome associated with periodontitis. J. Clin. Periodontol. 2010, 37, 241–247. [Google Scholar] [CrossRef]
  40. Ebersole, J.L.; Machen, R.L.; Steffen, M.J.; Willmann, D.E. Systemic acute-phase reactants, C-reactive protein and haptoglobin, in adult periodontitis. Clin. Exp. Immunol. 1997, 107, 347–352. [Google Scholar] [CrossRef]
Figure 1. Venn diagrams showing the proteins identified in common between SS and US samples, as well as the number of proteins identified exclusively in SS and US for OP (A), OWP (B), NP (C), and NWP (D).
Figure 1. Venn diagrams showing the proteins identified in common between SS and US samples, as well as the number of proteins identified exclusively in SS and US for OP (A), OWP (B), NP (C), and NWP (D).
Cells 12 01389 g001
Figure 2. Functional analysis of the distribution of proteins identified with differential expression between SS and US in OP group.
Figure 2. Functional analysis of the distribution of proteins identified with differential expression between SS and US in OP group.
Cells 12 01389 g002
Figure 3. Functional analysis of the distribution of proteins identified with differential expression between SS and US in OWP group.
Figure 3. Functional analysis of the distribution of proteins identified with differential expression between SS and US in OWP group.
Cells 12 01389 g003
Figure 4. Functional analysis of the distribution of proteins identified with differential expression between SS and US in NP group.
Figure 4. Functional analysis of the distribution of proteins identified with differential expression between SS and US in NP group.
Cells 12 01389 g004
Figure 5. Functional analysis of the distribution of proteins identified with differential expression between SS and US in NWP group.
Figure 5. Functional analysis of the distribution of proteins identified with differential expression between SS and US in NWP group.
Cells 12 01389 g005
Figure 6. Interaction networks among unique proteins identified in US from OP (A), OWP (B), NP (C), and NWP (D).
Figure 6. Interaction networks among unique proteins identified in US from OP (A), OWP (B), NP (C), and NWP (D).
Cells 12 01389 g006
Table 1. Salivary and periodontal parameters of the sample.
Table 1. Salivary and periodontal parameters of the sample.
OP (n = 10)
Mean ± SD
Median
[1st–3rd Quartiles]
OWP (n = 10)
Mean ± SD
Median
[1st–3rd Quartiles]
NP (n = 10)
Mean ± SD
Median
[1st–3rd Quartiles]
NWP (n = 10)
Mean ± SD
Median
[1st–3rd Quartiles]
p
SSF (mL/min)1.64 ± 0.151.73 ± 0.141.67 ± 0.141.76 ± 0.120.236 *
Total protein SSF (ug/ptn)35.04 ± 14.1530.95 ± 20.2327.08 ± 10.3028.65 ± 13.900.668 *
USF (mL/min)0.56 ± 0.090.55 ± 0.080.57 ± 0.080.56 ± 0.060.426 *
Total protein USF (ug/ptn)32.36 ± 15.5936.22 ± 20.0934.21 ± 9.8630.44 ± 15.730.668 *
PPD (mm)2.65
[2.45–2.81]
A
2.04
[1.96–2.10]
B
2.56
[2.42–2.79]
A
2.02
[2.00–2.16]
B
<0.001
CAL (mm)2.66
[2.45–2.82]
A
2.07
[2.02–2.10]
B
2.58
[2.44–2.79]
A
2.03
[2.00–2.16]
B
<0.001
Periodontitis—n (%)
Stage II
Stage III

5 (50%)
5 (50%)
-
9 (90%)
1 (10%)
-<0.001
OP, obesity and periodontitis; OWP, obesity without periodontitis; NP, normal BMI and periodontitis; NWP, normal BMI without periodontitis; SD, standard deviation; p, significance level; SSF, stimulated salivary flow; USF, unstimulated salivary flow; PPD, probing pocket depth; CAL, clinical attachment level; * ANOVA; Kruskal–Wallis (post hoc test: Dunn). Different letters indicate a statistically significant difference (p < 0.05). Bold values indicate significance level lower than 5%.
Table 2. Main expression differences (up- or down-regulated proteins by more than 2-fold) identified in stimulated saliva (SS) and unstimulated saliva (US) of pregnant women with obesity and periodontitis (OP).
Table 2. Main expression differences (up- or down-regulated proteins by more than 2-fold) identified in stimulated saliva (SS) and unstimulated saliva (US) of pregnant women with obesity and periodontitis (OP).
Accession NumberProtein NameGeneScoreFold ChangeLog (e)SDpED
P01009Alpha-1-antitrypsinSERPINA1695.871.770.06<0.01
P59666Neutrophil defensin 3DEFA39034.851.580.18<0.01
P14780Matrix metalloproteinase-9MMP9583.861.350.12<0.01
P15516Histatin-3HTN346163.781.330.09<0.01
P02812Basic salivary proline-rich protein 2PRB214172.7210.05<0.01
P04280Basic salivary proline-rich protein 1PRB114172.590.950.07<0.01
Q6S8J3POTE ankyrin domain family member EPOTEE10872.390.870.07<0.01
A5A3E0POTE ankyrin domain family member FPOTEF10872.340.850.05<0.01
P23280Carbonic anhydrase 6CA615982.290.830.09<0.01
P68032Actin, alpha cardiac muscle 1ACTC122032.120.750.04<0.01
P02042Hemoglobin subunit deltaHBD1382.100.740.06<0.01
P69905Hemoglobin subunit alphaHBA1; HBA251462.050.720.07<0.01
Q8N4F0BPI fold-containing family B member 2BPIFB22320.49−0.720.17<0.01
P01834Immunoglobulin kappa constantIGKC51200.49−0.720.07<0.01
P01860Immunoglobulin heavy constant gamma 3IGHG39810.48−0.730.13<0.01
P01833Polymeric immunoglobulin receptorPIGR52220.47−0.760.02<0.01
P01024Complement C3C31290.45−0.80.12<0.01
P02787SerotransferrinTF14920.43−0.850.02<0.01
P0DOX6Immunoglobulin mu heavy chainIGM6740.42−0.870.09<0.01
P01591Immunoglobulin J chainJCHAIN50850.42−0.870.04<0.01
P01871Immunoglobulin heavy constant muIGHM6920.41−0.890.11<0.01
P68871Hemoglobin subunit betaHBB4140.40−0.910.03<0.01
P14618Pyruvate kinase PKMPKM1190.40−0.910.09<0.01
P30613Pyruvate kinase PKLRPKLR800.39−0.950.270.02
P0DOX2Immunoglobulin alpha-2 heavy chainIGA263360.38−0.960.02<0.01
P10599ThioredoxinTXN9340.35−1.040.15<0.01
P05109Protein S100-A8S100A815310.35−1.060.03<0.01
P01861Immunoglobulin heavy constant gamma 4IGHG47010.33−1.110.17<0.01
P01023Alpha-2-macroglobulinA2M1860.31−1.170.05<0.01
P37837TransaldolaseTALDO1890.30−1.190.1<0.01
P06702Protein S100-A9S100A95550.28−1.290.03<0.01
P22079LactoperoxidaseLPO5460.27−1.30.12<0.01
P06744Glucose-6-phosphate isomeraseGPI1340.23−1.460.11<0.01
P01859Immunoglobulin heavy constant gamma 2IGHG26800.23−1.470.08<0.01
Q9UGM3Deleted in malignant brain tumors 1 proteinDMBT13650.23−1.490.09<0.01
P00738HaptoglobinHP15260.19−1.650.07<0.01
P0DOX7Immunoglobulin kappa light chainIGK24350.17−1.770.04<0.01
P69892Hemoglobin subunit gamma-2HBG24080.12−2.150.03<0.01
P69891Hemoglobin subunit gamma-1HBG14080.11−2.170.03<0.01
P02100Hemoglobin subunit epsilonHBE14080.11−2.180.03<0.01
P00739Haptoglobin-related proteinHPR2800.10−2.320.09<0.01
Q5T7N2LINE-1 type transposase domain-containing protein 1L1TD1350.08−2.590.030.01
Note: Log (e) (“e” is a constant = 2.71); SD, standard deviation; p, statistical significance (adjusted by False Discovery Rate - FDR = 4); ED, Expression differences; ↑ = up-regulated in SS (1-p > 0.95); ↓ = down-regulated in SS (p < 0.05). Bold values indicate significance level lower than 5%.
Table 3. Main expression differences (up- or down-regulated proteins by more than 2-fold) identified in stimulated saliva (SS) and unstimulated saliva (US) of pregnant women with obesity but without periodontitis (OWP).
Table 3. Main expression differences (up- or down-regulated proteins by more than 2-fold) identified in stimulated saliva (SS) and unstimulated saliva (US) of pregnant women with obesity but without periodontitis (OWP).
Accession NumberProtein NameGeneScoreFold ChangeLog (e)SDpED
P02808StatherinSTATH23,89618.732.930.14<0.01
P01871Immunoglobulin heavy constant muIGHM3793.421.230.06<0.01
P0DOX6Immunoglobulin mu heavy chainIGM3733.251.180.09<0.01
P02810Salivary acidic proline-rich phosphoprotein 1/2PRH1; PRH211113.221.170.02<0.01
P15516Histatin-3HTN311652.891.060.1<0.01
P01860Immunoglobulin heavy constant gamma 3IGHG31080.49−0.710.08<0.01
P01591Immunoglobulin J chainJCHAIN10730.47−0.760.06<0.01
P01876Immunoglobulin heavy constant alpha 1IGHA117340.46−0.770.02<0.01
P02788LactotransferrinLTF10440.43−0.850.07<0.01
P69905Hemoglobin subunit alphaHBA1; HBA2800.40−0.910.07<0.01
P12273Prolactin-inducible proteinPIP35710.38−0.980.02<0.01
P10599ThioredoxinTXN4610.37−0.990.240.01
P05109Protein S100-A8S100A858100.34−1.080.04<0.01
P61626Lysozyme CLYZ36690.31−1.160.04<0.01
P04406Glyceraldehyde-3-phosphate dehydrogenaseGAPDH2500.31−1.180.1<0.01
P0DOY2Immunoglobulin lambda constant 2IGLC22010.29−1.240.04<0.01
P59665Neutrophil defensin 1DEFA1; DEFA1B24050.26−1.340.2<0.01
P0DOX2Immunoglobulin alpha-2 heavy chainIGA213860.26−1.350.01<0.01
P04280Basic salivary proline-rich protein 1PRB1850.25−1.40.02<0.01
P02812Basic salivary proline-rich protein 2PRB2850.25−1.40.03<0.01
P0CG04Immunoglobulin lambda constant 1IGLC11150.24−1.430.04<0.01
P01877Immunoglobulin heavy constant alpha 2IGHA214370.23−1.450.01<0.01
P59666Neutrophil defensin 3DEFA324050.23−1.450.19<0.01
P06702Protein S100-A9S100A915840.20−1.630.06<0.01
P14618Pyruvate kinase PKMPKM1180.18−1.740.11<0.01
P522096-phosphogluconate dehydrogenase,
decarboxylating
PGD1590.16−1.850.15<0.01
P24158MyeloblastinPRTN33790.13−2.050.08<0.01
P07737Profilin-1PFN14490.10−2.320.06<0.01
P37837TransaldolaseTALDO12030.10−2.350.07<0.01
Note: Log (e) (“e” is a constant = 2.71); SD, standard deviation; p, statistical significance (adjusted by False Discovery Rate - FDR = 4); ED, Expression differences; ↑ = up-regulated in SS (1-p > 0.95); ↓ = down-regulated in SS (p < 0.05). Bold values indicate significance level lower than 5%.
Table 4. Main expression differences (up- or down-regulated proteins by more than 2-fold) identified in stimulated saliva (SS) and unstimulated saliva (US) of pregnant women with normal BMI but with periodontitis (NP).
Table 4. Main expression differences (up- or down-regulated proteins by more than 2-fold) identified in stimulated saliva (SS) and unstimulated saliva (US) of pregnant women with normal BMI but with periodontitis (NP).
Accession Number Protein NameGeneScoreFold ChangeLog (e)SDpED
P00738HaptoglobinHP1149.682.270.06<0.01
P04280Basic salivary proline-rich protein 1PRB13665.261.660.25<0.01
P15516Histatin-3HTN318254.901.590.4<0.01
P01009Alpha-1-antitrypsinSERPINA11804.711.550.08<0.01
P02812Basic salivary proline-rich protein 2PRB23664.481.50.05<0.01
P02647Apolipoprotein A-IAPOA144324.101.410.05<0.01
P02790HemopexinHPX1533.251.180.05<0.01
P04080Cystatin-BCSTB12393.161.150.08<0.01
Q8N4F0BPI fold-containing family B member 2BPIFB24222.771.020.06<0.01
P02810Salivary acidic proline-rich phosphoprotein 1/2PRH1; PRH223082.510.920.16<0.01
P69892Hemoglobin subunit gamma-2HBG22712.390.870.130.01
P69891Hemoglobin subunit gamma-1HBG12712.390.870.12<0.01
P02100Hemoglobin subunit epsilonHBE12712.390.870.12<0.01
P02808StatherinSTATH22862.270.820.22<0.01
P02679Fibrinogen gamma chainFGG1182.230.80.13<0.01
P59666Neutrophil defensin 3DEFA310622.120.750.06<0.01
P20742Pregnancy zone proteinPZP2122.100.740.09<0.01
P59665Neutrophil defensin 1DEFA1; DEFA1B10622.100.740.06<0.01
P0DTE8Alpha-amylase 1CAMY1C39942.080.730.01<0.01
P22079LactoperoxidaseLPO820.49−0.720.11<0.01
P0CG38POTE ankyrin domain family member IPOTEI2570.47−0.750.07<0.01
P02788LactotransferrinLTF580.47−0.750.260.01
P0CG39POTE ankyrin domain family member JPOTEJ1760.47−0.760.1<0.01
P68133Actin, alpha skeletal muscleACTA19040.45−0.80.07<0.01
Q8NHQ9ATP-dependent RNA helicase DDX55DDX551690.44−0.820.260.02
P63267Actin, gamma-enteric smooth muscleACTG29040.43−0.850.07<0.01
Q5VSP4Putative lipocalin 1-like protein 1LCN1P19670.43−0.850.08<0.01
Q6S8J3POTE ankyrin domain family member EPOTEE3660.39−0.940.05<0.01
A5A3E0POTE ankyrin domain family member FPOTEF3660.39−0.950.06<0.01
P68032Actin, alpha cardiac muscle 1ACTC19040.37−10.04<0.01
Q9BYX7Putative beta-actin-like protein 3POTEKP1090.36−1.020.05<0.01
P62736Actin, aortic smooth muscleACTA29040.35−1.050.05<0.01
P63261Actin, cytoplasmic 2ACTG113840.34−1.090.02<0.01
P60709Actin, cytoplasmic 1ACTB13840.33−1.10.02<0.01
P02042Hemoglobin subunit deltaHBD4990.32−1.130.05<0.01
A0M8Q6Immunoglobulin lambda constant 7IGLC73200.23−1.480.08<0.01
Q8TAX7Mucin-7MUC78420.20−1.630.02<0.01
Q16378Proline-rich protein 4PRR49840.18−1.690.080.01
Q562R1Beta-actin-like protein 2ACTBL23620.17−1.750.02<0.01
P01859Immunoglobulin heavy constant gamma 2IGHG2850.17−1.80.06<0.01
P06702Protein S100-A9S100A93440.13−2.020.04<0.01
P61626Lysozyme CLYZ32230.06−2.770.03<0.01
Note: Log (e) (“e” is a constant = 2.71); SD, standard deviation; p, statistical significance (adjusted by False Discovery Rate - FDR = 4); ED, Expression differences; ↑ = up-regulated in SS (1-p > 0.95); ↓ = down-regulated in SS (p < 0.05). Bold values indicate significance level lower than 5%.
Table 5. Main expression differences (up- or down-regulated proteins by more than 2-fold) identified in stimulated saliva (SS) and unstimulated saliva (US) of pregnant women with normal BMI and without periodontitis (NWP).
Table 5. Main expression differences (up- or down-regulated proteins by more than 2-fold) identified in stimulated saliva (SS) and unstimulated saliva (US) of pregnant women with normal BMI and without periodontitis (NWP).
Accession NumberProtein NameGeneScoreFold ChangeLog (e)SDpED
P02814Submaxillary gland androgen-regulated protein 3BSMR3B51,4082.941.080.02<0.01
P61769Beta-2-microglobulinB2M1322.941.080.1<0.01
P06702Protein S100-A9S100A91812.530.930.08<0.01
Q8TAX7Mucin-7MUC714972.270.820.04<0.01
P0DOX8Immunoglobulin lambda-1 light chainIGL132562.080.730.06<0.01
Q9BYX7Putative beta-actin-like protein 3POTEKP2540.50−0.70.09<0.01
P01036Cystatin-SCST416,4800.50−0.70.02<0.01
A0M8Q6Immunoglobulin lambda constant 7IGLC721520.48−0.740.23<0.01
P07737Profilin-1PFN18590.46−0.770.16<0.01
P61626Lysozyme CLYZ37430.45−0.80.16<0.01
P02810Salivary acidic proline-rich phosphoprotein 1/2PRH1; PRH239310.43−0.850.02<0.01
Q562R1Beta-actin-like protein 2ACTBL26260.43−0.850.09<0.01
A5A3E0POTE ankyrin domain family member FPOTEF3370.42−0.860.07<0.01
P37837TransaldolaseTALDO11990.42−0.870.23<0.01
Q96DA0Zymogen granule protein 16 homolog BZG16B22,2150.42−0.870.03<0.01
Q6S8J3POTE ankyrin domain family member EPOTEE3370.41−0.880.05<0.01
P0CG38POTE ankyrin domain family member IPOTEI830.36−1.020.09<0.01
P01861Immunoglobulin heavy constant gamma 4IGHG41220.32−1.130.29<0.01
P01859Immunoglobulin heavy constant gamma 2IGHG2300.31−1.170.29<0.01
P0CG39POTE ankyrin domain family member JPOTEJ830.30−1.220.17<0.01
P23280Carbonic anhydrase 6CA611770.25−1.380.03<0.01
P31025Lipocalin-1LCN158070.10−2.30.03<0.01
Q5VSP4Putative lipocalin 1-like protein 1LCN1P145920.10−2.330.03<0.01
Note: Log (e) (“e” is a constant = 2.71); SD, standard deviation; p, statistical significance (adjusted by False Discovery Rate - FDR = 4); ED, Expression differences; ↑ = up-regulated in SS (1-p > 0.95); ↓ = down-regulated in SS (p < 0.05). Bold values indicate significance level lower than 5%.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Foratori-Junior, G.A.; Ventura, T.M.O.; Grizzo, L.T.; Jesuino, B.G.; Castilho, A.V.S.S.; Buzalaf, M.A.R.; Sales-Peres, S.H.d.C. Is There a Difference in the Proteomic Profile of Stimulated and Unstimulated Saliva Samples from Pregnant Women with/without Obesity and Periodontitis? Cells 2023, 12, 1389. https://doi.org/10.3390/cells12101389

AMA Style

Foratori-Junior GA, Ventura TMO, Grizzo LT, Jesuino BG, Castilho AVSS, Buzalaf MAR, Sales-Peres SHdC. Is There a Difference in the Proteomic Profile of Stimulated and Unstimulated Saliva Samples from Pregnant Women with/without Obesity and Periodontitis? Cells. 2023; 12(10):1389. https://doi.org/10.3390/cells12101389

Chicago/Turabian Style

Foratori-Junior, Gerson Aparecido, Talita Mendes Oliveira Ventura, Larissa Tercilia Grizzo, Bruno Gualtieri Jesuino, Ana Virgínia Santana Sampaio Castilho, Marília Afonso Rabelo Buzalaf, and Silvia Helena de Carvalho Sales-Peres. 2023. "Is There a Difference in the Proteomic Profile of Stimulated and Unstimulated Saliva Samples from Pregnant Women with/without Obesity and Periodontitis?" Cells 12, no. 10: 1389. https://doi.org/10.3390/cells12101389

APA Style

Foratori-Junior, G. A., Ventura, T. M. O., Grizzo, L. T., Jesuino, B. G., Castilho, A. V. S. S., Buzalaf, M. A. R., & Sales-Peres, S. H. d. C. (2023). Is There a Difference in the Proteomic Profile of Stimulated and Unstimulated Saliva Samples from Pregnant Women with/without Obesity and Periodontitis? Cells, 12(10), 1389. https://doi.org/10.3390/cells12101389

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