ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 02 May 2025

Sec. Oral Microbes and Host

Volume 15 - 2025 | https://doi.org/10.3389/fcimb.2025.1576906

Characterisation of the periodontal proteome in gingival crevicular fluid and saliva using SWATH-MS

  • 1. Oral Sciences Research Group, Special Needs Unit, Department of Surgery and Medical-Surgical Specialties, School of Medicine and Dentistry, Universidade de Santiago de Compostela, Health Research Institute of Santiago (IDIS), Santiago de Compostela, Spain

  • 2. Oral Pathology and Rehabilitation Research Unit (UNIPRO), University Institute of Health Sciences (IUCS-CESPU), Gandra, Portugal

  • 3. Proteomic Unit, Health Research Institute of Santiago de Compostela (IDIS), Santiago de Compostela, Spain

Article metrics

View details

4

Citations

2,5k

Views

682

Downloads

Abstract

Introduction:

Proteomic techniques are useful to analyse the periodontal proteome in gingival crevicular fluid (GCF) and saliva. However, few investigations have assessed and compared the GCF and salivary proteomes. Therefore, this research aims to analyse the proteome structure and compare protein expression in these fluids between individuals with periodontal health and those with periodontitis.

Methods:

GCF and saliva were collected from 44 periodontally healthy subjects and 41 with periodontitis (stages III-IV). Samples were analysed using sequential window acquisition of all theoretical mass spectra (SWATH-MS), and proteins were identified employing the UniProt database. The periodontal proteome structure was assessed using principal component analysis (PCA). Differential protein expression was defined as an adjusted p-value <0.05 combined with a fold-change ≥2 (upregulated) or ≤0.5 (downregulated).

Results:

250 abundant proteins were quantified in GCF and 377 in saliva (238 in common). The proteome structure was different in periodontitis compared to periodontal health in both oral fluids. In GCF, 63 (25.2%) proteins were differentially expressed, with 38 upregulated and 25 downregulated in periodontitis. The most overexpressed proteins were haemoglobin subunits (Hbs) beta (fold-change of 5.06) and alpha (4.35), carbonic anhydrase 1 (4.28), and protein S100-P (4.27). Among the underexpressed proteins, 14 were keratins, with type II cytoskeletal 6B being the most downregulated (0.10), together with glyceraldehyde-3-phosphate dehydrogenase (0.12) and zymogen granule protein 16 homolog B (0.13).

In saliva, 59 (15.7%) proteins were differentially expressed, with 55 upregulated and four downregulated in periodontitis. Twenty-nine proteins showed a fold-change ≥4, highlighting beta-2-microglobulin (44.14), keratin, type I cytoskeletal 13 (36.23), neutrophil defensin 1 (25.08), proteins S100-A9 (12.30), A8 (10.61), A12 (4.76), and P (4.72), annexin A1 (9.34), lysozyme C (4.98), immunoglobulin heavy constant alpha 1 (4.45), resistin (4.37), and Hbs beta (4.20) and alpha (4.06). The most downregulated protein was lipocalin-1 (0.35). Fourteen proteins were differentially expressed in GCF and saliva, where seven were keratins being underexpressed in GCF but overexpressed in saliva.

Conclusion:

Periodontitis alters the periodontal proteome structure and the expression of numerous abundant proteins in GCF and saliva. However, proteins expressed vary qualitatively and quantitatively, indicating different expression patterns between oral fluids.

Introduction

Periodontitis is a chronic inflammatory disease that affects the gingiva, the alveolar bone, and the periodontal ligament surrounding the teeth (Kinane et al., 2017). Periodontitis not only has an impact on oral health but is also related to numerous systemic diseases, being especially relevant in their association with cardiovascular diseases or diabetes (Monsarrat et al., 2016). Today, periodontitis remains one of the most prevalent diseases worldwide, affecting almost 45% of the world’s population (Eke et al., 2020).

Several proteins have been shown to play an essential role in the initiation, progression, and severity of periodontitis (Bostanci and Belibasakis, 2018). However, other relevant proteins have not yet been identified, and their associations with the disease remain undiscovered (Bostanci and Belibasakis, 2018). Analyses using high-throughput proteomic techniques enable the periodontal proteome to be studied in greater detail, contributing to the determination of extensive and comprehensive protein profiles associated with periodontal diseases (Bostanci and Bao, 2017; Bostanci and Belibasakis, 2018).

Considering the range of proteomic techniques available, sequential window acquisition of all theoretical mass spectra (SWATH-MS) stands out as a particularly effective approach (Anjo et al., 2017). SWATH-MS is an untargeted, data-independent acquisition (DIA), label-free quantification technique designed to eliminate the variability introduced by chemical modifications, expand the dynamic range of quantification, and optimise the detection of data present in a sample (Frederick and Ciborowski, 2016). This method was developed to combine the proteome coverage potential of data-dependent acquisition (DDA) approaches with the quantitative consistency of targeted DIA proteomics (e.g., selected reaction monitoring – SRM) (Ludwig et al., 2018). These features make SWATH-MS a powerful strategy for discovering novel biomarkers (Anjo et al., 2017).

The most commonly used sample types for studying the periodontal proteome are gingival crevicular fluid (GCF) and saliva. Both can be easily collected in a non-invasive manner, making them particularly suitable for potential translation into clinical diagnostic tools (Guzman et al., 2014). The analysis of these fluids using various proteomic approaches has enabled the identification and quantification of numerous proteins with significantly different expression profiles depending on periodontal status (Ahmad et al., 2025). However, few studies have simultaneously analysed the GCF and salivary proteome collected from the same individuals (Grant et al., 2022; Tang et al., 2019). Additionally, SWATH-MS has previously proven useful in identifying novel GCF and salivary biomarkers with strong diagnostic potential for periodontitis (Blanco-Pintos et al., 2024, 2025a). Nevertheless, to date, SWATH-MS has not been employed to comprehensively study and compare the proteomic profile and protein expression between periodontal health and periodontitis.

Accordingly, this cross-sectional study aimed to (1) analyse the periodontal proteome structure and (2) compare differential protein expression in GCF and saliva between individuals with periodontal health and those with periodontitis using SWATH-MS.

Material and methods

The complete analysis protocol applied in the present study is detailed in Figure 1.

Figure 1

Selection of study groups

A convenience sample of 90 eligible participants, comprising 45 periodontally healthy subjects (control group) and 45 individuals affected by untreated periodontitis (periodontitis group) was recruited between 2018 and 2021 from 250 consecutive patients in the general population who were referred for assessment of their oral health status to the School of Medicine and Dentistry (Universidade de Santiago de Compostela, Spain) and the Instituto Superior de Ciências da Saúde-Norte, Cooperativa de Ensino Superior, Politécnico e Universitário (CESPU, Gandra, Paredes, Portugal).

Patients were selected if they fulfilled the following inclusion criteria: 1) aged 20 to 75; 2) had at least 10 natural teeth; 3) no previous periodontal treatment; 4) had no medical history of diabetes mellitus, hepatic or renal disease, or other severe medical conditions or transmittable illnesses; 5) had not taken systemic antimicrobials in the previous six months; 6) had no intake of antiinflammatory medication in the previous four months; 7) did not routinely use oral antiseptics; 8) had no history of alcohol or drug abuse; 9) were not pregnant or breastfeeding; and 10) used no orthodontic appliances.

Two experienced and previously calibrated dentists performed the periodontal diagnostics, with intra- and interexaminer intraclass correlation coefficients (ICC) greater than 0.950. Bleeding on probing (BOP) and the bacterial plaque level (BPL) were recorded for the entire mouth on a binary scale (presence/absence) at six sites per tooth. The probing pocket depth (PPD) and clinical attachment loss (CAL) were recorded using a PCP-UNC 15 probe at six sites per tooth. PPD was measured as the distance from the gingival margin to the base of the gingival sulcus. A positive value for gingival recession (REC) was assessed for the apical displacement of the gingival margin in relation to the cementoenamel junction. Considering an average biological width (BW) of 2 mm (Mulla et al., 2023), CAL was defined for three conditions: (1) PPD <4 mm and REC of 0 mm, CAL was equal to 0 mm; (2) PPD ≥4 mm and REC of 0 mm, CAL was equal to PPD subtracting BW; and (3) if REC >0 mm, CAL was PPD plus REC.

Standardised teeth radiographs were acquired to assess the alveolar bone status. The diagnosis of periodontitis was based on the clinical and radiographic information thus obtained. The control group included subjects who were periodontally healthy according to the criteria in the Classification of Periodontal Diseases and Conditions (Caton et al., 2018; Tonetti et al., 2018), i.e., BOP <10%, no location with a PPD ≥4 mm, and no radiographic evidence of alveolar bone loss. This classification was also used for the periodontitis group, which was composed of patients diagnosed with stage III or IV and grades B or C generalised periodontitis.

A questionnaire was used to evaluate the smoking habits of the participants, collecting information on their smoking status (non-smoker, former smoker, current smoker), the duration of their smoking status (as a former or current smoker), and the number of cigarettes consumed per day. Non-smokers were considered those who had never smoked, current smokers were those who had been smoking for at least one year, and former smokers were those who stopped doing so at least one year before, irrespective of how long they had been smoking.

The research was conducted following the principles of the Declaration of Helsinki (revised in 2013) on human experimentation studies (World Medical Association, 2013). The protocol was approved by the Galician Clinical Research Ethics Committee (registration numbers 2018/295 and 2021/417) and the Instituto Superior de Ciências da Saúde-Norte, CESPU (registration number 35/CEIUCS/2019). All the participants provided written informed consent to their involvement in the study.

Collection of the GCF and saliva samples and preparation for the mass spectrometry analysis

The GCF and saliva samples were collected one or two weeks after the initial examination. The GCF samples from the control and periodontal subjects were collected and pooled from five non-adjacent proximal sites. The samples were taken from teeth in quadrants 1 and 3 in the periodontally healthy controls. In patients with periodontitis, samples were collected from the sites with the deepest PPD values. Two paper strips were inserted into the gingival sulcus or periodontal pocket for 30 seconds. In cases of visible contamination with blood, samples were discarded, and new sites were sampled. The strips were collected and inserted into labelled tubes containing 300 mL of 0.01 M phosphate-buffered saline (PBS) with a pH of 7.2.

Immediately before the collection of the GCF samples, one ml of unstimulated whole saliva was collected from each participant using the spitting method at least one hour after the last meal and after the last brushing (Navazesh and Christensen, 1982). Subsequently, the GCF and saliva samples were stored at -80°C until further analysis.

The total protein concentrations of the samples were quantified using the reducing-agent and detergent-compatible (RC DC™) Protein Assay Kit (BIO-RAD, Hercules, CA, USA), following the manufacturer’s protocol. Approximately 50 μg of protein from each sample was precipitated using the methanol/chloroform (MeOH/CHCl3) method to remove potential contaminants such as deoxyribonucleic acid (DNA), ribonucleic acid (RNA), and lipids. Purified GCF and salivary proteins were loaded onto a 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) gel to concentrate them into a band for in-gel digestion. Electrophoresis was interrupted once the front had penetrated 3 mm into the gel (Bonzon-Kulichenko et al., 2011; Perez-Hernandez et al., 2013). The protein band was detected with SYPRO Ruby fluorescent staining (Lonza, Basel, Switzerland) and cut out. We then proceeded to manual in-gel tryptic digestion, as previously explained (Shevchenko et al., 1996), but with slight modifications. The gel pieces were reduced with 10 mM dithiothreitol (Sigma-Aldrich, St. Louis, Missouri, USA) in 50 mM ammonium bicarbonate (Sigma-Aldrich) and alkylated with 55 mM iodoacetamide (Sigma-Aldrich) in 50 mM ammonium bicarbonate.

The gel pieces were then rinsed with 50 mM ammonium bicarbonate in 50% MeOH (HPLC grade, Scharlau, Barcelona, Spain), dehydrated by adding acetonitrile (ACN) (HPLC grade, Scharlau), and dried in a SpeedVac (Savant™). Modified porcine trypsin (Promega, Madison, Wisconsin, USA) was added to the dried gel pieces at a final concentration of 20 ng/μL in 20 mM ammonium bicarbonate and incubated at 37°C for 16 hours. Peptides were extracted three times by way of a 20-minute incubation in 40 μl of 60% ACN in 0.5% acetic acid (HCOOH). The resulting peptide extracts were pooled, concentrated in a SpeedVac, and stored at -20°C until use.

Protein quantification by SWATH-MS

Construction of MS/MS spectral libraries

To obtain a robust representation of the peptides and MS/MS present in the GCF and saliva study groups, four pools of around 10 samples from each group (periodontal health and periodontitis) were prepared after protein digestion using equal amounts of peptides from each sample. In more detail, after peptide resuspension in water plus 0.1% of formic acid with sonication per 10 minutes (the resuspension was made to obtain a µg/µl peptides solution), from each of the resuspended samples/peptides, we used three µl for create the corresponding pool to obtain the spectral library (Blanco-Pintos et al., 2024, 2025a, 2025b; Bugallo-Casal et al., 2025; Casado-Fernández et al., 2024; López-Valverde et al., 2024; Mondelo-Macía et al., 2024; Nóvoa et al., 2025; Rodrigues et al., 2024; Torres Iglesias et al., 2024; Vázquez-Mera et al., 2024). Four µl (around four µg) of each pool were separated using reversed-phase (RP) LC to identify the proteins and create the SWATH library.

The gradient was created with a micro-LC system (Eksigent Technologies nanoLC 400, SCIEX, Redwood City, California, USA) coupled with a TripleTOF® (time-of-flight) 6600 high-speed mass spectrometer (SCIEX) via a microflow ionisation source. A silica-based RP analytical column (C18CL 150×0.30 mm) with a particle size of 3 µm and a pore size of 120 Å (Eksigent, SCIEX) was employed. The samples were pre-washed on a YMC-TRIART C18 pre-column (YMC Technologies, Teknokroma, Barcelona, Spain), which had the same particle and pore sizes, and was connected in-line with the analytical column. The charge pump provided a 0.1% formic acid solution in water at a rate of 10 μl/min. The micropump generated a flow rate of 5 μl/min and operated under gradient elution conditions, employing 0.1% formic acid in water as mobile phase A and 0.1% formic acid in ACN as mobile phase B. Peptides were separated using a 40-minute gradient from 2% to 90% of mobile phase B. The injection volume was 4 μl (approximately 4 μg of protein).

The TripleTOF 6600 (SCIEX) was used for data acquisition, employing a DDA workflow. The source and interface conditions comprised an ion-source voltage (ISVF) of 5500V, a curtain gas (CUR) of 25, a collision energy (CE) of 10, and an ion-source gas 1 (GS1) of 25. The instrument was operated using the software Analyst TF 1.7.1 (SCIEX). Switching criteria were set for ions with a mass-to-charge ratio (m/z) ranging from 350 to 1400, a charge state of 2-5, a mass tolerance of 250 parts per million (ppm), and an abundance threshold of more than 200 counts per second (cps). Selected ions meeting these conditions were excluded for 15 seconds. The instrument was automatically calibrated every four hours using commercial PepCalMix tryptic peptides (SCIEX) as the external calibrants.

After the MS/MS analysis, the data files were processed using ProteinPilotTM 5.0.1 (SCIEX), which employs the ParagonTM algorithm for database search and Progroup™ for data grouping. Files were searched for peptide/protein identification against the UniProt human-specific database (UniProt Consortium, 2014) (https://www.uniprot.org/, UniProt release 2022_02, containing 20387 human proteins), specifying iodoacetamide at cysteine alkylation as a variable modification and methionine oxidation as a fixed modification. All known contaminants (keratins epidermal, trypsin…) were manually excluded. The false discovery rate (FDR) was set at ≤1% for both peptides and proteins (Aggarwal and Yadav, 2016). The MS/MS spectra of the identified peptides were used to generate a spectral library, which subsequently enabled the extraction of SWATH peaks using MicroApp MS/MSALL 2.0 for PeakView 2.2 (SCIEX). Only peptides with a confidence score above 99% (obtained by searching against the UniProt database) were included in the spectral library.

Relative quantification

The GCF and saliva samples were analysed individually, using a DIA method. These were then examined using the equipment and gradients described above for spectral library construction, although a SWATH-MS acquisition strategy was employed instead of a DDA approach. The method involves repeating an acquisition cycle. In particular, 100 TOF MS/MS scans (400 to 1500 m/z, high sensitivity mode, and 50 millisecond acquisition time) were performed using sequential, overlapping precursor isolation windows of varying widths (1 m/z overlap) that cover the entire mass range between 400 and 1250 m/z, obtained with a previous TOF MS scan (400–1500 m/z and 50 millisecond acquisition time) for each cycle. The total cycle time was 6.3 seconds. A SWATH variable window spreadsheet (SCIEX) was used to optimise the width of the 100 variable windows for each set of samples according to the ion density found in the DDA runs.

The extraction of specific data from the trace-ion chromatograms and fragments obtained using the SWATH-MS methodology was performed with the PeakView 2.2 software (SCIEX), utilising MicroApp MS/MSALL 2.0 and the previously created spectral library. Only peak area values integrated based on signal intensity from 10 peptides per protein and 7 fragments per peptide were selected. Shared and modified peptides were excluded from the processing. Five-minute windows and widths of 30 ppm were used to extract the ion chromatograms. SWATH-MS quantification was performed for all the proteins in the ion library identified by ProteinPilot™ with an FDR ≤1%. The retention times (RTs) of the selected peptides for each protein were realigned in each run according to the peptide-indexed RTs obtained in each sample throughout the chromatogram. The areas of each protein in each sample were acquired by combining the peak areas of the 10 peptides and the seven corresponding fragment ions.

The integrated peak areas (processed PeakView® mrkvw files) were exported directly to the MarkerView™ software (SCIEX) for the relative quantitative examination, which is an approach that has been used previously for SWATH-MS data analyses (Agra-Bermejo et al., 2020; Álvarez et al., 2020; Camino et al., 2020, 2022; Chantada-Vázquez et al., 2020; García-Vence et al., 2020; Novelle et al., 2021; Ortea et al., 2018). This export generated three files containing quantitative information on individual ions, the total intensity of the different ions for a particular peptide, and the total intensity of the different peptides for a specific protein. MarkerView™ utilises processing algorithms that accurately identify chromatographic and spectral peaks directly within the raw SWATH data. The software’s alignment of the data compensates for slight variations in mass and RT values, ensuring precise comparisons of identical compounds in different samples. To homogenise the data obtained, a most likely ratio (MLR) normalisation was performed (Colangelo et al., 2015).

Statistical analysis

For the calculation of the sample size, we considered the design of independent groups (periodontal health vs periodontitis), the possibility of using a nonparametric two-tailed test of mean differences between two independent groups, an effect size of 0.65, an alpha error of 0.05 and a statistical power of 0.80, a minimum sample size of 40 subjects in each group was required. Due to the possible loss of patients for various reasons, each study group initially consisted of 45 participants. Sample size calculation was performed with the G*Power program (version 3.1.9.4) (Faul et al., 2007).

Five of the 90 subjects selected were excluded for non-compliance with the project phases. Therefore, the clinical groups were finally composed of 44 periodontally healthy subjects and 41 patients with periodontitis. One GCF sample from those with periodontitis had to be excluded due to a processing error.

The R software (version 4.2.2) (R Core Team, 2022) was used to perform the statistical data analysis of the GCF and saliva samples from the study groups (control and periodontitis).

Comparison of clinical characteristics between control subjects and patients with periodontitis

The Shapiro-Wilk test was used to assess the distribution of quantitative variables. For variables with a normal distribution, the Student’s t-test was applied; otherwise, the Mann-Whitney U test was used. Fisher’s exact test and the Chi-squared test were used to assess the association between qualitative variables and the clinical conditions. A significance level of p<0.05 was established.

Periodontal proteome structure

Principal component analysis (PCA) was employed to visualise the clustering of the GCF and saliva samples in each analysis group. The FactoMineR software (version 2.8) (Lê et al., 2008) was used to perform the PCA, and then the amount of variability explained within each component was assessed using the factoextra package (version 1.0.7) (Kassambara and Mundt, 2020). The vegan software (version 2.6.4) (Oksanen et al., 2022) was used to perform two evaluations: 1) an analysis of the similarities (ANOSIM) between groups; and 2) an analysis of variance, using a similarity matrix via the adonis2 function of the non-parametric test permutational multivariate analysis of variance (PERMANOVA) to that end (Anderson, 2001). A significance level of p<0.05 was established.

Differential protein expression

Concerning all the proteins quantified in GCF and saliva, we evaluated the mean expression difference between the clinical groups using the non-parametric Mann-Whitney U test. The p-values obtained were adjusted using the mutoss package (version 0.1-13) (Team et al., 2023) with the Benjamini-Hochberg correction (Benjamini and Hochberg, 1995). The differential expression was statistically significant if the adjusted p-value was <0.05 (-log10 adjusted p-value=1.30).

Subsequently, the effect size of each protein in all the phases was calculated using the effsize package (version 0.8.1) (Torchiano, 2020) and Cohen’s d (Cohen, 1988). The expression ratio of each protein between conditions was calculated, and a fold-change ≥2 (log2 fold-change ≥1) or ≤0.5 (log2 fold-change ≤-1) was considered significant.

Volcano plots were created from the adjusted p-values and the differential expression of the log2 fold-change. Proteins with an adjusted p-value <0.05 and a log2 fold-change ≥1 (upregulated) or ≤-1 (downregulated) were differentially expressed.

Results

Clinical characteristics

Considering the clinical parameters, full-mouth BPL, BOP, PPD, and CAL levels were significantly higher in patients with periodontitis than in control subjects. Similarly, BOP, PPD and CAL levels from sampled GCF sites were significantly higher in periodontal patients (Table 1).

Table 1

Clinical characteristicsStudy groups
Control group (n=44)Periodontitis group (n=41)*p-value
Age, years31.86 ± 10.2653.80 ± 10.09<0.001
No. male (%)17 (38.64)21 (51.22)NS
No. female (%)27 (61.36)20 (48.78)
Smoking dimensions
No. patients non-smokers (%)27 (61.36)13 (31.71)0.004
No. patients smokers (%)15 (34.09)17 (41.46)
No. patients ex-smokers (%)2 (4.55)11 (26.83)
No. cigarettes/day11.26 ± 7.4916.71 ± 10.29NS
No. patients <10 cigarettes/day (%)7 (15.91)7 (17.07)NS
No. patients 10–20 cigarettes/day (%)8 (18.18)15 (36.59)
No. patients >20 cigarettes/day (%)2 (4.55)6 (14.63)
No. years of smoking11.75 ± 6.6728.39 ± 10.95<0.001
No. patients <5 years smoking (%)2 (4.55)1 (2.44)NS
No. patients ≥5 years smoking (%)15 (34.09)27 (65.85)
No. years without smoking7.00 ± 2.837.27 ± 5.83NS
No. patients <5 years without smoking (%)0 (0.00)4 (9.76)NS
No. patients ≥5 years without smoking (%)2 (4.55)7 (17.07)
Clinical parameters
No. of teeth27.66 ± 1.1221.41 ± 4.91<0.001
No. sites evaluated165.95 ± 6.72131.12 ± 29.04<0.001
BPL, %14.23 ± 11.0351.10 ± 27.15<0.001
No. patients BPL ≤30% (%)39 (88.64)8 (19.51)<0.001
No. patients BPL >30% (%)5 (11.36)33 (80.49)
BOP, %3.98 ± 2.6550.95 ± 19.18<0.001
No. patients BOP <10% (%)44 (100)0 (0.00)<0.001
No. patients BOP 10-30% (%)0 (0.00)6 (14.63)
No. patients BOP >30% (%)0 (0.00)35 (85.37)
PPD, mm1.75 ± 0.293.59 ± 0.58<0.001
% sites PPD ≤3 mm100 ± 0.0055.48 ± 17.49<0.001
% sites PPD 4–5 mm0.00 ± 0.0033.21 ± 10.62<0.001
% sites PPD ≥6 mm0.00 ± 0.0011.31 ± 9.21<0.001
CAL, mm0.01 ± 0.033.34 ± 1.80<0.001
% sites CAL 0 mm99.77 ± 1.0031.02 ± 20.74<0.001
% sites CAL 1–2 mm0.23 ± 1.0012.64 ± 7.99<0.001
% sites CAL 3–4 mm0.00 ± 0.0022.03 ± 8.28<0.001
% sites CAL ≥5 mm0.00 ± 0.0034.31 ± 25.62<0.001
Sampled sites
BOP, %0.85 ± 3.1987.25 ± 16.48<0.001
PPD, mm1.91 ± 0.406.09 ± 0.75<0.001
CAL, mm0.00 ± 0.006.08 ± 1.75<0.001
Periodontitis stage and grade
No. patients stage III (%)NA29 (70.73)NA
No. patients stage IV (%)NA12 (29.27)
No. patients grade B (%)NA17 (41.46)NA
No. patients grade C (%)NA24 (58.54)

Age, gender, smoking habit, and clinical characteristics of GCF and saliva sampled patients associated with periodontal status in the control group and periodontitis.

No., number; BPL, bacterial plaque level; BOP, bleeding on probing; PPD, probing pocket depth; mm, millimetre; CAL, clinical attachment loss; NA, not applicable; NS, not significant.

*One GCF sample from the periodontitis group had to be excluded due to a processing error resulting in 40 patients.

Mann-Whitney U test (non-normal distribution).

CAL values different from 0 for traumatic reasons, not of periodontal origin.

Values are presented as means ± standard deviations, and in some cases, the number of subjects or the percentage of sites is also provided. A significant p-value of <0.05 was established.

Proteins identified and quantified in GCF and saliva samples by SWATH-MS

The in silico spectral library contained 37985 transitions, 4556 peptides, and 427 proteins. Of these, 250 GCF and 377 salivary proteins were quantified by SWATH-MS using the previously defined criteria. Interestingly, 238 proteins were quantified in both oral fluids, representing 95.2% of the GCF proteins and 63.1% of those in saliva (Appendix S1).

Periodontal proteome structure in GCF and saliva

The periodontal proteome structure of the GCF and salivary proteins was evaluated using a PCA. The first two principal components or dimensions accounted for 37.3% of the total variance in GCF and 23.0% in saliva (Appendix S2).

The PCA plot of the GCF and saliva samples revealed distinct clustering between samples with good periodontal health and those with periodontitis. These findings were confirmed by ANOSIM and adonis2 tests, which showed significant differences in GCF and saliva between individuals with periodontitis and those with healthy periodontium (Figures 2, 3).

Figure 2

Figure 3

Differential expression of the GCF proteins

Considering the expression of the 250 GCF proteins quantified and the previously established criteria, 63 (25.2%) were found to be differentially expressed in periodontitis compared to periodontal health. Of these, 38 (60.3%) were shown to be upregulated, and 25 (39.7%) were downregulated (Figure 4).

Figure 4

Among the 38 overexpressed proteins, four showed a ≥4-fold higher expression in periodontitis compared to periodontal health (log2 fold-change ≥2). Of these, haemoglobin subunits (Hb) beta and alpha showed the highest fold-change (5.06 and 4.35, respectively), followed by carbonic anhydrase 1 (4.28) and protein S100-P (4.27). The remaining 34 proteins showed at least 2-fold higher expression in periodontitis, highlighting immunoglobulins (Ig) heavy constant gamma 3 (3.78), kappa light chain (2.19), and heavy variable 5-51 (2.18), plasma protease C1 inhibitor (3.61), haptoglobin (3.35), and another Hb, delta (3.26) (Table 2).

Table 2

UniProt accessionProtein nameAdjusted p-valueEffect sizeFold-changeMag Log2 fold-change
P68871Hb beta1.4938E-080.625.06≥2 (N=4)
P69905Hb alpha1.1748E-080.604.35
P00915Carbonic anhydrase 11.7237E-090.714.28
P25815Protein S100-P2.3003E-081.594.27
P01860Ig heavy constant gamma 37.795E-071.173.781  to− 2 (N=34)
P05155Plasma protease C1 inhibitor1.7237E-091.323.61
P00738Haptoglobin7.4947E-081.083.35
P02042Hb delta6.3584E-080.643.26
P02768Albumin2.2729E-070.903.14
P43652Afamin6.4177E-071.042.87
P01008Antithrombin-III2.5737E-081.272.79
P02730Band 3 anion transport protein5.2674E-050.732.77
P02749Beta-2-glycoprotein 13.809E-060.892.65
P00918Carbonic anhydrase 23.1194E-060.462.64
P02774Vitamin D-binding protein2.2729E-070.942.63
P32119Peroxiredoxin-27.4947E-080.612.62
P08697Alpha-2-antiplasmin3.4816E-081.362.61
P02766Transthyretin7.3368E-081.172.61
P02751Fibronectin3.4816E-081.332.55
P19823ITI H22.2156E-050.822.49
Q14624ITI H43.5623E-060.992.46
P02765Alpha-2-HS-glycoprotein3.1194E-061.092.45
P01042Kininogen-15.8831E-060.962.44
P01023Alpha-2-macroglobulin2.6176E-071.232.41
P10909Clusterin0.03290.662.34
P00751Complement factor B9.0221E-060.932.33
P02679Fibrinogen gamma chain3.5623E-061.062.27
P01011Alpha-1-antichymotrypsin2.5979E-071.352.25
P0DOX7Ig kappa light chain0.00010.892.19
A0A0C4DH38Ig heavy variable 5-510.01620.762.18
P13716Delta-aminolevulinic acid dehydratase1.2837E-050.772.18
P02787Serotransferrin1.2263E-050.932.16
P0C0L5Complement C4-B1.7525E-061.062.15
P01024Complement C31.527E-050.982.12
P04004Vitronectin1.1241E-051.082.10
P02790Hemopexin0.00030.842.10
P01009Alpha-1-antitrypsin6.7049E-071.102.05
P19827ITI H10.00170.682.00

GCF proteins differentially upregulated in periodontitis compared to periodontal health (N=38 proteins).

Log, logarithmic; Mag, magnitude; N, number.

Abbreviations related to protein name: Hb, haemoglobin subunit; Ig, immunoglobulin; ITI, inter-alpha-trypsin inhibitor heavy chain.

Differentially upregulated proteins were considered if the adjusted p-value was <0.05 and a fold-change ≥2 (log2 fold-change ≥1).

Magnitude Log2 fold-change: ≥2 (red)=fold-change ≥4; 1 to 2 (purple)=fold-change 2 − 4.

Interestingly, of the 25 downregulated proteins, 14 were keratins, accounting for 56.0%. Considering a fold-change ≤0.25 (log2 fold-change ≤-2), which means a ≤4-fold lower expression in periodontitis, four were keratins: type II cytoskeletal 6B (0.10), type II cytoskeletal 6A (0.13), type II cytoskeletal 5 (0.21), and type I cytoskeletal 16 (0.24). The remaining 11 downregulated keratins exhibited a fold-change of 0.25-0.44. Other protein types that showed a fold-change ≤0.25 were glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (0.12) and zymogen granule protein 16 homolog B (ZG16B) (0.13) (Table 3).

Table 3

UniProt accessionProtein nameAdjusted p-valueEffect sizeFold-changeMag Log2 fold-change
P04259Keratin, type II cytoskeletal 6B1.2006E-07-0.840.10≤-2 (N=6)
P04406GAPDH2.0474E-11-1.720.12
P02538Keratin, type II cytoskeletal 6A1.4675E-08-0.850.13
Q96DA0ZG16B2.0474E-11-1.590.13
P13647Keratin, type II cytoskeletal 51.4901E-07-0.820.21
P08779Keratin, type I cytoskeletal 161.2894E-06-0.840.24
P02533Keratin, type I cytoskeletal 141.0958E-06-0.890.25-1 −  to -2 (N=19)
Q6P5S2Protein LEG1 homolog6.7049E-07-0.870.26
P13646Keratin, type I cytoskeletal 131.3100E-05-0.720.27
P19013Keratin, type II cytoskeletal 43.3994E-05-0.720.29
P04264Keratin, type II cytoskeletal 17.4947E-08-0.770.33
P01833Polymeric Ig receptor9.7776E-07-0.810.33
P48668Keratin, type II cytoskeletal 6C0.0040-0.710.34
P07355Annexin A21.1010E-05-0.800.34
P13645Keratin, type I cytoskeletal 101.8719E-06-0.610.35
Q04695Keratin, type I cytoskeletal 170.0003-0.760.36
P08185Corticosteroid-binding globulin0.0005-0.760.37
P49913Cathelicidin antimicrobial peptide4.2083E-08-1.340.38
Q01546Keratin, type II cytoskeletal 2 oral0.0002-0.510.38
P35527Keratin, type I cytoskeletal 97.3206E-07-0.650.41
P35908Keratin, type II cytoskeletal 2 epidermal9.9088E-05-0.490.44
P20700Lamin-B10.0012-0.710.44
P02545Prelamin-A/C0.0007-0.930.45
P0886540S ribosomal protein SA2.1657E-05-0.960.45
P19961Alpha-amylase 2B0.0003-0.720.46

GCF proteins differentially downregulated in periodontitis compared to periodontal health (N=25 proteins).

Log, logarithmic; Mag, magnitude; N, number.

Abbreviations related to protein name: GAPDH, glyceraldehyde-3-phosphate dehydrogenase; Ig, immunoglobulin; LEG, liver-enriched gene; ZG16B, zymogen granule protein 16 homolog B.

Differentially downregulated proteins were considered if the adjusted p-value was <0.05 and a fold-change ≤0.50 (log2 fold-change ≤-1).

Magnitude Log2 fold-change: ≤-2 (green)=fold-change ≤0.25; -1 to -2 (blue)=fold-change 0.25 − 0.50.

Differential expression of the salivary proteins

In saliva, of the 377 proteins quantified, 59 (15.7%) were differentially expressed, with 55 (93.2%) being upregulated and only four (6.8%) downregulated (Figure 5).

Figure 5

Consequently, the most differentially expressed proteins were upregulated, with beta-2-microglobulin exhibiting the greatest fold-change (44.14). Furthermore, in contrast to GCF, nine keratins were overexpressed in saliva in periodontitis, with seven exhibiting a fold-change of ≥ 4 and two between 2 and 4. Of these, keratin, type I cytoskeletal 13 showed the highest fold-change (36.23), and the remaining eight values ranged from 2.65 to 23.31. Other noteworthy proteins include the S100 family: A9 (12.30), A8 (10.61), A12 (4.76), P (4.72), and A6 (2.67), in addition to Ig heavy constant alpha 1 (4.45), alpha-2 heavy chain (2.59), and Hb beta (4.20) and alpha (4.06). Considering the remaining overexpressed proteins, of particular interest were neutrophil defensin 1 (25.08), annexin A1 (9.34), lysozyme C (4.98), resistin (4.37), cathepsin G (3.28), myeloperoxidase (2.77), and azurocidin (2.75) (Table 4).

Table 4

UniProt accessionProtein nameAdjusted p-valueEffect sizeFold-changeMag Log2 fold-change
P61769Beta-2-microglobulin3.2279E-081.4544.14≥2 (N=29)
P13646Keratin, type I cytoskeletal 132.1223E-051.2036.23
Q9UBC9Small proline-rich protein 34.8555E-050.5825.77
P59665Neutrophil defensin 13.2279E-080.7625.08
P13647Keratin, type II cytoskeletal 56.7217E-051.1623.31
Q01546Keratin, type II cytoskeletal 2 oral5.9948E-051.1716.90
P19013Keratin, type II cytoskeletal 40.00021.2016.39
P06702Protein S100-A98.006E-090.9712.30
P08779Keratin, type I cytoskeletal 168.6934E-061.1811.55
P04264Keratin, type II cytoskeletal 11.3533E-061.6211.18
P05109Protein S100-A81.1339E-081.0810.61
P04083Annexin A10.00031.209.34
P10599Thioredoxin3.2903E-051.448.95
P28676Grancalcin6.1142E-060.848.94
Q04695Keratin, type I cytoskeletal 173.327E-051.026.10
O14818Proteasome subunit alpha type-70.00441.046.09
Q02487.2Isoform 2 of desmocollin-23.3862E-071.175.44
P07305Histone H1.02.5758E-050.815.10
P61626Lysozyme C1.4364E-071.044.98
P04792Heat shock protein beta-10.00061.024.85
P80511Protein S100-A124.7261E-061.014.76
P25815Protein S100-P7.7803E-071.294.72
P01876Ig heavy constant alpha 13.6244E-061.064.45
Q9HD89Resistin1.3885E-070.914.37
P19105Myosin regulatory light chain 12A2.0242E-051.254.37
P68871Hb beta0.00110.584.20
Q9UBG3Cornulin6.3468E-050.994.18
P02652Apolipoprotein A-II0.00120.754.10
P69905Hb alpha0.00090.554.06
Q9NUQ9CYFIP-related Rac1 interactor B3.6474E-071.253.771  to− 2 (N=26)
P04114Apolipoprotein B-1003.9948E-050.873.66
P0886540S ribosomal protein SA0.00630.933.66
P62805Histone H42.4866E-051.183.46
P08311Cathepsin G0.00090.773.28
P0DUB6Alpha-amylase 1A0.00090.893.15
Q8N1N4Keratin, type II cytoskeletal 782.5758E-051.213.10
P60842Eukaryotic initiation factor 4A-I0.00051.203.07
P02671Fibrinogen alpha chain1.2471E-051.003.01
P07355Annexin A20.00121.142.99
P02647Apolipoprotein A-I0.00060.712.89
P05164Myeloperoxidase0.00030.892.77
P20160Azurocidin0.00060.752.75
P06703Protein S100-A61.9327E-051.192.67
P08727Keratin, type I cytoskeletal 190.00270.842.65
P12429Annexin A30.00060.852.62
P0DOX2Ig alpha-2 heavy chain2.7128E-050.992.59
P68431Histone H3.10.01090.322.54
P62987Ubiquitin-60S ribosomal protein L400.00021.282.43
Q9BTM1Histone H2A.J0.00290.832.31
P49913Cathelicidin antimicrobial peptide0.00051.012.30
Q99879Histone H2B type 1-M0.00060.932.26
P43490Nicotinamide phosphoribosyltransferase0.00061.032.24
P13639Elongation factor 20.00211.142.12
P6225814-3–3 protein epsilon0.00650.912.11
P35579Myosin-91.2471E-051.232.02

Salivary proteins differentially upregulated in periodontitis compared to periodontal health (N=55 proteins).

Log, logarithmic; Mag, magnitude; N, number. Abbreviations related to protein name: CYFIP, cytoplasmic fragile X messenger ribonucleoprotein 1-interacting protein; Hb, haemoglobin subunit; Ig, immunoglobulin.

Differentially upregulated proteins were considered if the adjusted p-value was <0.05 and a fold-change ≥2 (log2 fold-change ≥1).

Magnitude Log2 fold-change: ≥2 (red)=fold-change ≥4; 1 to 2 (purple)=fold-change 2 − 4.

The fold-change of the downregulated proteins was considerably lower, with no protein having a value of ≤0.25. The most underexpressed protein was lipocalin-1 with a fold-change of 0.35 (Table 5).

Table 5

UniProt accessionProtein nameAdjusted p-valueEffect sizeFold-changeMag Log2 fold-change
P31025Lipocalin-11.9327E-05-0.810.35-1 −  to -2 (N=4)
P30086Phosphatidylethanolamine-binding protein 11.4203E-05-1.050.47
P1080960 kDa heat shock protein, mitochondrial5.9948E-05-0.620.48
Q96DA0ZG16B0.0001-0.760.49

Salivary proteins differentially downregulated in periodontitis compared to periodontal health (N=4 proteins).

Log, logarithmic; Mag, magnitude; N, number.

Abbreviations related to protein name: kDa, kilodalton; ZG16B, zymogen granule protein 16 homolog B.

Differentially downregulated proteins were considered if the adjusted p-value was <0.05 and a fold-change ≤0.50 (log2 fold-change ≤-1).

Magnitude Log2 fold-change: -1 to -2 (blue)=fold-change 0.25 − 0.50.

Comparison of differentially expressed proteins in GCF and saliva

Comparing the results obtained in GCF and saliva, 14 proteins were differentially expressed in both sample types, representing 22.2% and 23.7% of the total proteins differentially expressed, respectively.

As previously noted, Hb beta and alpha were overexpressed in GCF and saliva (GCF: 5.06 and 4.35; saliva: 4.20 and 4.06, respectively), in addition to protein S100-P (4.27; 4.72). Furthermore, ZG16B was found to be downregulated in both sample types (0.13; 0.49).

However, the remaining 10 proteins were underexpressed in GCF but overexpressed in saliva. Interestingly, seven of these proteins were keratins (types I cytoskeletal 13, 16, and 17, and types II cytoskeletal 1, 2 oral, 4, and 5) with a fold-change in GCF of 0.21-0.38, whereas in saliva it was 6.10-36.23. The remaining three proteins were annexin A2 (0.34; 2.99), cathelicidin antimicrobial peptide (0.38; 2.30), and 40S ribosomal protein SA (0.45; 3.66) (Table 6).

Table 6

UniProt accessionProtein nameGCFSaliva
Adjusted p-valueFold-changeMag.Adjusted p-valueFold-changeMag.
P68871Hb beta1.4938E-085.06UP0.00114.20UP
P69905Hb alpha1.1748E-084.350.00094.06
P25815Protein S100-P2.3003E-084.277.7803E-074.72
Q96DA0ZG16B2.0474E-110.13DW0.00010.49DW
P13647Keratin, type II cytoskeletal 51.4901E-070.216.7217E-0523.31UP
P08779Keratin, type I cytoskeletal 161.2894E-060.248.6934E-0611.55
P13646Keratin, type I cytoskeletal 131.3100E-050.272.1223E-0536.23
P19013Keratin, type II cytoskeletal 43.3994E-050.290.000216.39
P04264Keratin, type II cytoskeletal 17.4947E-080.331.3533E-0611.18
Q04695Keratin, type I cytoskeletal 170.00030.363.327E-056.10
Q01546Keratin, type II cytoskeletal 2 oral0.00020.385.9948E-0516.90
P07355Annexin A21.1010E-050.340.00122.99
P49913Cathelicidin antimicrobial peptide4.2083E-080.380.00052.30
P0886540S ribosomal protein SA2.1657E-050.450.00633.66

Proteins differentially expressed in GCF and saliva in periodontitis compared to periodontal health (N=14 proteins).

DW, downregulated; GCF, gingival crevicular fluid; Mag, magnitude of fold-change; UP, upregulated.

Abbreviations related to protein name: Hb, haemoglobin subunit; ZG16B, zymogen granule protein 16 homolog B.

Differentially proteins were considered if the adjusted p-value was <0.05 and a fold-change ≥2 (log2 fold-change ≥1; upregulated) or ≤0.5 (log2 fold-change ≤-1; downregulated).

Discussion

The analysis of GCF and saliva samples using SWATH-MS enabled the quantification of 250 abundant proteins in GCF and 377 in saliva. In both fluids, the structure of the periodontal proteome differed between individuals with periodontitis and those with a healthy periodontium. Regarding protein expression, approximately 25% of GCF proteins and 16% of salivary proteins were found to be differentially expressed. In both fluids, most of these proteins were overexpressed in periodontitis—60% vs. 40% in GCF and 93% vs. 7% in saliva, respectively. Furthermore, although 95% of GCF proteins and 63% of salivary proteins were quantified in both sample types, the most differentially expressed proteins varied, and only about 22% of them were shared between the two oral fluids.

Methodological differences in periodontal proteome analysis

To date, the periodontal proteome has been evaluated using different proteomic techniques (Rizal et al., 2020). However, this is the first study to use SWATH-MS to analyse the structure and differential protein expression of the periodontal proteome.

Sample size is a relevant limitation in previous proteomics research, with most studies including ≤20 subjects per clinical group (Bellei et al., 2022; Bostanci et al., 2018; Casarin et al., 2023; Grant et al., 2022; Guzman et al., 2018; Mertens et al., 2018; Tang et al., 2019; Yi et al., 2021). In contrast, the present study included a notably larger sample size, with more than 40 patients. Increasing the sample size is crucial for enhancing statistical robustness and improving the ability to detect true differences in protein expression levels between clinical conditions (Ahmad et al., 2025).

The periodontal proteome has been primarily studied in gingival crevicular fluid (GCF) and saliva samples (Guzman et al., 2014), as both oral fluids are valuable for investigating the oral and systemic environments (Ahmad et al., 2025). However, there is a tendency to justify the selection of one fluid over the other, since most studies have evaluated the periodontal proteome in either GCF or saliva (Bellei et al., 2022; Bostanci et al., 2018; Casarin et al., 2023; Hartenbach et al., 2020; Kim et al., 2021; Yi et al., 2021). Only two studies have analysed both sample types (Grant et al., 2022; Tang et al., 2019). These studies focused on assessing the diagnostic accuracy of proteins without comparing the structural and expression differences between the GCF and saliva proteomes.

Analysis of GCF and saliva samples using SWATH-MS yielded a spectral library of 427 proteins, with 250 and 377 proteins quantified in each fluid, respectively, based on stringent criteria that required at least 10 peptides per protein and seven fragments per peptide. Although this number of quantified proteins is lower compared to previous studies, those studies applied less stringent criteria, quantifying proteins based on as few as two unique peptides (Bostanci et al., 2018; Grant et al., 2022; Kim et al., 2021; Silbereisen et al., 2024; Yi et al., 2021).

Although we could have employed more extensive libraries or adopted a library-free DIA approach (Demichev et al., 2020; Rosenberger et al., 2014), building our own spectral library ensured accurate protein identification by focusing on the most abundant proteins. This approach lays the groundwork for further validation using less sensitive methods, which could ultimately support clinical applications.

Comparing GCF and saliva, saliva samples have a higher number of proteins (Grant et al., 2022; Tang et al., 2019). This is because saliva also comprises proteins from the salivary glands (Guzman et al., 2014). Moreover, 238 proteins were common to both oral fluids, representing 95% of those in GCF and 63% in saliva. This reflects that most GCF components are excreted from the pocket through saliva and can be identified in both fluids (Bostanci and Bao, 2017).

Grant et al. reported that about 30% of the total proteins were present in both fluids (Grant et al., 2022). This lower proportion may be attributed to the different proteomic techniques used (isobaric tag for relative and absolute quantitation -iTRAQ-) and the analysis of stimulated saliva samples (Grant et al., 2022). Conversely, we collected unstimulated saliva as the stimulation of salivary flow may produce inappropriate glycoprotein glycosylation (Buzalaf et al., 2020).

The criteria applied in previous publications to consider a protein as differentially expressed were based on a p-value <0.05, a minimum fold-change, or both (Bellei et al., 2022; Bostanci et al., 2018; Casarin et al., 2023; Hartenbach et al., 2020; Kim et al., 2021; Shin et al., 2019). However, our study was more demanding, being the first to apply the Benjamini-Hochberg correction for controlling false positives (Benjamini and Hochberg, 1995). This variability in the criteria applied, as well as the different proteomics techniques employed, complicates the comparison of the results obtained with other studies.

Using a more permissive statistical threshold may increase the number of differentially expressed proteins detected. However, in scenarios involving distinct clinical states—such as periodontal health and advanced periodontitis—it is more relevant to prioritise the most strongly differentially expressed proteins. Conversely, when comparing closely related conditions, such as periodontal health and gingivitis, applying a less stringent threshold may be necessary to detect subtle variations in protein expression.

Periodontal proteome structure in GCF and saliva

The periodontal proteome structure was different between periodontitis and periodontal health in both sample types, consistent with previous research (Shin et al., 2019; Yi et al., 2021). Our research group demonstrated that smoking status does not alter the distribution of the GCF periodontal proteome in periodontitis (Blanco-Pintos et al., 2025b). Therefore, periodontitis is a more relevant condition for modifying the structure of the periodontal proteome. However, future research should investigate the impact of other clinical variables, such as age, on the spatial distribution of the periodontal proteome.

Differential expression of GCF proteins

The functions and possible relationship to periodontitis of the most relevant GCF proteins differentially expressed are specified in Table 7.

Table 7

MagUniProt accessionProtein nameRelation with periodontitis
UPREGULATEDP68871Hb betaBleeding caused by inflammation (Guzman et al., 2014)
P69905Hb alpha
P02042Hb delta
P00915Carbonic anhydrase 1Unclear
P25815Protein S100-PUnclear
P01860Ig heavy constant gamma 3Bacterial defence line (Megha and Mohanan, 2021)
P0DOX7Ig kappa light chain
A0A0C4DH38Ig heavy variable 5-51
P05155Plasma protease C1 inhibitorAntiinflammatory activity (Beinrohr et al., 2011)
P00738HaptoglobinBind free haemoglobin. Acute phase protein (Smith and McCulloh, 2015)
DOWNREGULATEDP04259Keratin, type II cytoskeletal 6BMaintain epithelial integrity (Groeger and Meyle, 2019)
P02538Keratin, type II cytoskeletal 6A
P13647Keratin, type II cytoskeletal 5
P08779Keratin, type I cytoskeletal 16
P02533Keratin, type I cytoskeletal 14
P13646Keratin, type I cytoskeletal 13
P19013Keratin, type II cytoskeletal 4
P04264Keratin, type II cytoskeletal 1
P48668Keratin, type II cytoskeletal 6C
P13645Keratin, type I cytoskeletal 10
Q04695Keratin, type I cytoskeletal 17
Q01546Keratin, type II cytoskeletal 2 oral
P35527Keratin, type I cytoskeletal 9
P35908Keratin, type II cytoskeletal 2 epidermal
P04406GAPDHGAIT complex (Sampath et al., 2003)
Q96DA0ZG16BUnclear

Functions and relationship to periodontitis of the most relevant GCF proteins differentially expressed in periodontitis compared to periodontal health.

GAIT, gamma-interferon-activated inhibitor of translation; Mag, magnitude of fold-change.

Abbreviations related to protein name: GAPDH, glyceraldehyde-3-phosphate dehydrogenase; Hb, haemoglobin subunit; Ig, immunoglobulin; ZG16B, zymogen granule protein 16 homolog B.

Differentially upregulated proteins were considered if the adjusted p-value was <0.05 and a fold-change ≥2 (log2 fold-change ≥1) or ≤0.5 (log2 fold-change ≤-1).

Hb beta and alpha were the most overexpressed proteins, with Hb delta also being up to five times increased in periodontitis. It is reported that, despite discarding highly visible blood-contaminated GCF samples, there is a greater tendency for spontaneous bleeding in periodontitis due to inflammation of periodontal tissues (Guzman et al., 2014).

Consequently, haptoglobin was also found to be upregulated (fold-change >3), because it binds to free haemoglobin. Furthermore, haptoglobin is an acute-phase protein, and its expression is increased in infectious and inflammatory processes (Smith and McCulloh, 2015). Previous proteomic studies have also demonstrated the overexpression of these proteins in periodontitis (Bellei et al., 2022; Carneiro et al., 2014; Tang et al., 2019).

Moreover, the Ig’s heavy constant gamma 3, kappa light chain, and heavy variable 5–51 were also shown to be upregulated (fold-change >2). Igs are essential components of the immune system for recognising and neutralising pathogens by binding to their antigens (Megha and Mohanan, 2021). Therefore, in periodontitis, their expression is increased due to bacterial dysbiosis and inflammation of periodontal tissues (Bellei et al., 2022; Carneiro et al., 2014; Marinho et al., 2019; Yi et al., 2021).

Considering the remaining upregulated proteins, the plasma protease C1 inhibitor (fold-change >3.5) is particularly noteworthy. This protein functions as an essential antiinflammatory protein by regulating the complement system (Beinrohr et al., 2011). However, carbonic anhydrase 1 and protein S100-P, they were overexpressed in periodontitis (fold-change around 4.3 for both), as in previous studies (Baliban et al., 2012; Silva-Boghossian et al., 2013); however, their relationship with periodontitis is still unclear.

Conversely, the most differentially expressed proteins were downregulated, as keratin, type II cytoskeletal 6B showed a fold-change of 0.10, implying a 10-fold lower expression in periodontitis compared to periodontal health. Moreover, the majority (56%) of the 25 proteins underexpressed in periodontitis were different types of keratins. Keratins are essential structural proteins in the epithelial barrier, providing resistance to microbial invasion (Groeger and Meyle, 2019). The decrease in keratins in GCF periodontitis samples may be related to periodontal tissue degradation, compromising the integrity of the epithelium and making the tissues more vulnerable to bacterial invasion (Groeger and Meyle, 2019).

Apart from keratins, the most underexpressed proteins include GAPDH and ZG16B with also a fold-change of around 0.10. GAPDH is a key enzyme in glycolysis (Tisdale, 2002), but it is also part of the gamma interferon-activated inhibitor of translation (GAIT) complex. This complex may contribute to the resolution of the inflammatory response following cytokine activation (Sampath et al., 2003). Therefore, its underexpression may be associated with reduced antiinflammatory activity of the GAIT complex. However, the function of ZG16B and its potential relationship with periodontitis remain unclear.

Differential expression of the salivary proteins

In saliva, the most differentially expressed protein was beta-2-microglobulin, which was approximately 44 times more abundant in patients with periodontitis than in healthy subjects. This protein is a component of major histocompatibility complex (MHC) class I molecules and is released during inflammatory processes due to its crucial role in the immune system and antigen presentation (Bernier, 1980; Sreejit et al., 2014) (Table 8).

Table 8

MagUniProt accessionProtein nameRelation with periodontitis
UPREGULATEDP61769Beta-2-microglobulinPresentation of antigens and regulator of the immune system (Bernier, 1980; Sreejit et al., 2014)
P59665Neutrophil defensin 1Antibacterial activity, induces the production of proinflammatory cytokines (Ericksen et al., 2005; Wang et al., 2016)
P04083Annexin A1Control proinflammatory stimuli (Casarin et al., 2023)
P61626Lysozyme CAntibacterial activity, enhances the activity of immunoagents (Surna et al., 2009)
Q9HD89ResistinStimulate proinflammatory reaction (Tripathi et al., 2020)
P08311Cathepsin GAntibacterial activity, regulate inflammatory reaction (Gao et al., 2018)
P05164MyeloperoxidaseAntimicrobial and anti/proinflammatory activity (Aratani, 2018)
P20160AzurocidinAntimicrobial activity, induces cytokine release (Soehnlein and Lindbom, 2009)
P13646Keratin, type I cytoskeletal 13Maintain epithelial integrity (Groeger and Meyle, 2019)
P13647Keratin, type II cytoskeletal 5
Q01546Keratin, type II cytoskeletal 2 oral
P19013Keratin, type II cytoskeletal 4
P08779Keratin, type I cytoskeletal 16
P04264Keratin, type II cytoskeletal 1
Q04695Keratin, type I cytoskeletal 17
Q8N1N4Keratin, type II cytoskeletal 78
P08727Keratin, type I cytoskeletal 19
P06702Protein S100-A9Heterodimer S100A8/A9 (calprotectin) has antimicrobial and proinflammatory activity (Wang et al., 2018)
P05109Protein S100-A8
P80511Protein S100-A12Innate immune response with proinflammatory function (Meijer et al., 2012)
P25815Protein S100-PUnclear
P06703Protein S100-A6Unclear
P01876Ig heavy constant alpha 1Bacterial defence line (Megha and Mohanan, 2021)
P0DOX2Ig alpha-2 heavy chain
P68871Hb betaBleeding caused by inflammation (Guzman et al., 2014)
P69905Hb alpha
DWP31025Lipocalin-1Possible relation with innate immunity (Ong et al., 2006).

Functions and relationship to periodontitis of the most relevant salivary proteins differentially expressed in periodontitis compared to periodontal health.

DW, downregulated; Mag, magnitude of fold-change.

Abbreviations related to protein name: Hb, haemoglobin subunit; Ig, immunoglobulin.

Differentially upregulated proteins were considered if the adjusted p-value was <0.05 and a fold-change ≥2 (log2 fold-change ≥1) or ≤0.5 (log2 fold-change ≤-1).

Other proteins that stand out are neutrophil defensin 1, annexin A1, lysozyme C, resistin, cathepsin G, myeloperoxidase, and azurocidin, which were upregulated between 2.8 and 25.1-fold in periodontitis. These proteins participate in the immune system by having proinflammatory (promoting cytokine secretion) and/or antiinflammatory functions (Aratani, 2018; Casarin et al., 2023; Gao et al., 2018; Soehnlein and Lindbom, 2009; Surna et al., 2009; Tripathi et al., 2020; Wang et al., 2016). Additionally, some of them have also demonstrated antimicrobial activity (Ericksen et al., 2005) (Table 8). Thus, overexpression of these proteins can indicate the dysregulation of the host immune response in periodontitis.

Similarly, other proteins with proinflammatory and antimicrobial functions that were overexpressed were those belonging to the S100 family, A8, A9, and A12, with a fold-change >4.5 (Meijer et al., 2012; Wang et al., 2018). Moreover, S100-P and S100-A6 were also found to be upregulated, as in previous studies (Grant et al., 2022; Salazar et al., 2013), although their possible relationship with periodontitis has not yet been discovered.

Furthermore, as in GCF, two Igs (heavy constant alpha 1 and alpha-2 heavy chain), as well as Hb beta and alpha, were overexpressed in periodontitis. However, keratins showed a different expression pattern, as nine keratins were found to be upregulated in saliva.

Conversely, only four salivary proteins were found to be downregulated in periodontitis. Moreover, the difference in expression with respect to healthy subjects was lower than that shown by the overexpressed proteins. Lipocalin-1 was the most underexpressed saliva protein with a nearly three times lower expression in periodontitis. This protein is involved in taste perception but also has a possible relationship with innate immunity (Ong et al., 2006).

Periodontal proteome expression in GCF and saliva: similarities and differences

As mentioned above, the protein composition of GCF and saliva samples was similar, as 95% of abundant GCF proteins and 63% of abundant salivary proteins were quantified in both oral fluids. However, the proportion of differentially expressed proteins was higher in GCF than in saliva (25% vs. 16%).

Furthermore, although more upregulated than downregulated proteins were found in both sample types in periodontitis, the distribution was more balanced in GCF (60% vs. 40%), whereas saliva showed a strong predominance of overexpressed proteins (93% vs. 7%). Additionally, the most differentially expressed proteins in GCF were underexpressed in periodontitis, with expression levels up to 10-fold lower. In contrast, in saliva, 11 proteins showed a 10-fold or greater increase in expression in periodontitis, with the highest difference reaching a 44-fold increase.

When comparing the differentially expressed proteins in both sample types, keratins are particularly noteworthy, as previously mentioned. Of the 14 proteins differentially expressed in both oral fluids, seven were keratins that were underexpressed in GCF but overexpressed in saliva (Table 9). This contrasting pattern may result from keratin exfoliation from the periodontal pocket due to epithelial barrier breakdown, leading to increased concentrations in saliva (Bostanci and Bao, 2017).

Table 9

MagUniProt accessionProtein nameRelation with periodontitis
GCFSaliva
UPUPP68871Hb betaBleeding caused by inflammation (Guzman et al., 2014)
P69905Hb alpha
P25815Protein S100-PUnclear
DOWNREGULATEDDWQ96DA0ZG16BUnclear
UPREGULATEDP13647Keratin, type II cytoskeletal 5Maintain epithelial integrity (Groeger and Meyle, 2019)
P08779Keratin, type I cytoskeletal 16
P13646Keratin, type I cytoskeletal 13
P19013Keratin, type II cytoskeletal 4
P04264Keratin, type II cytoskeletal 1
Q04695Keratin, type I cytoskeletal 17
Q01546Keratin, type II cytoskeletal 2 oral
P07355Annexin A2Annexin A2/S100-A10 heterotetramer has proinflammatory function (Bharadwaj et al., 2013)
P49913Cathelicidin antimicrobial peptideAntimicrobial protein, integral component of the innate immune system (Murakami et al., 2004)
P0886540S ribosomal protein SAFacilitate bacteria elimination, antiinflammatory activity (Sun et al., 2021)

Functions and relationship to periodontitis of the GCF and salivary proteins differentially expressed in periodontitis compared to periodontal health in both oral fluids.

DW, downregulated; GCF, gingival crevicular fluid; Mag, magnitude of fold-change; UP, upregulated.

Abbreviations related to protein name: Hb, haemoglobin subunit; ZG16B, zymogen granule protein 16 homolog B.

Differentially upregulated proteins were considered if the adjusted p-value was <0.05 and a fold-change ≥2 (log2 fold-change ≥1) or ≤0.5 (log2 fold-change ≤-1).

A similar opposing expression pattern between GCF and saliva was observed for annexin A2, cathelicidin antimicrobial peptide, and 40S ribosomal protein SA. These proteins play roles in the immune response, exhibiting antiinflammatory, proinflammatory, and antimicrobial activities (Bharadwaj et al., 2013; Murakami et al., 2004; Sun et al., 2021).

Implications for practice and research perspectives

To better understand the periodontal proteome in GCF and saliva across periodontal health and periodontitis, we included a heterogeneous population with respect to age and smoking status. Proteomic analyses have shown that both factors influence protein expression (Blanco-Pintos et al., 2025b; Zhang et al., 2022). This diversity enhances the generalisability of the findings by helping to isolate the specific impact of periodontitis on the proteome. However, the influence of clinical variables on the periodontal proteome in GCF and saliva should be further investigated. This would allow for a more accurate evaluation and comparison of their effects on the proteome in periodontitis across both oral fluids.

Additionally, we analysed two distinct clinical conditions: periodontal health and advanced periodontitis (stages III–IV). Future research should also investigate the periodontal proteome in earlier disease stages, such as initial periodontitis (stages I–II) and gingivitis, in both GCF and saliva. This would help clarify how the proteome is progressively altered from health to disease, and how these molecular changes are reflected in both oral fluids.

From a clinical perspective, differentially expressed proteins cannot be considered biomarkers without diagnostic accuracy analyses (Gürsoy and Kantarci, 2022), such as those previously published by our research group (Blanco-Pintos et al., 2024, 2025a). Consistently, nearly all proteins demonstrating excellent diagnostic accuracy for detecting periodontitis were differentially expressed in both GCF and saliva. This includes proteins such as ZG16B and carbonic anhydrase 1, whose relationship with periodontitis remains unclear. As a result, functional studies aimed at elucidating the role of these proteins in the pathogenesis of periodontitis are of particular interest.

Moreover, the proteins with the highest diagnostic capacity differed between GCF and saliva, just as the most differentially expressed proteins varied between these fluids. This study is the first to demonstrate that while similar proteins are quantified in both GCF and saliva, their expression patterns differ qualitatively and quantitatively, underscoring the potential need for different molecular biomarkers in each oral fluid.

In previous studies from our research group, the diagnostic accuracy of GCF proteins was found to be higher than that of salivary proteins (Blanco-Pintos et al., 2024, 2025a). This difference may be attributed to the higher protein composition in saliva, which leads to greater inter-subject variability (Guzman et al., 2014). Clinically, both oral fluids are applicable, as while GCF biomarkers are probably more relevant in clinical practice alongside traditional diagnostic parameters, saliva, in addition to the clinical setting, is more suitable for large-scale population screenings and self-assessment (Guzman et al., 2014).

However, these results require validation in future studies. Commonly used validation techniques, such as enzyme-linked immunosorbent assay (ELISA), Western blot, and protein arrays, have certain limitations, including higher detection limits for some proteins compared to proteomic approaches, the requirement for highly specific antibodies, and an inability to provide a comprehensive analysis of the proteome (Jayasena et al., 2016). Therefore, it may be more appropriate to use SWATH-MS in an independent patient cohort or to employ other proteomic techniques such as SRM (Bostanci et al., 2018).

In conclusion, the periodontal proteome structure in GCF and saliva differs between individuals with periodontitis and those with periodontal health. Differential expression analysis reveals that a considerable proportion of abundant proteins in both sample types are altered in periodontitis, being higher in GCF. Moreover, the most differentially expressed proteins vary between these oral fluids both qualitatively and quantitatively. Therefore, although periodontitis induces changes in the periodontal proteome in GCF and saliva, the expression patterns differ between these fluids.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ebi.ac.uk/pride/archive/, PXD043474 (GCF samples) and PXD043491 (saliva samples).

Ethics statement

The studies involving humans were approved by Galician Clinical Research Ethics Committee (registration numbers 2018/295 and 2021/417) Instituto Superior de Ciências da Saúde-Norte, CESPU (registration number 35/CEIUCS/2019). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

BT: Conceptualization, Investigation, Methodology, Visualization, Writing – original draft. RA: Investigation, Methodology, Writing – review & editing. SB: Investigation, Visualization, Writing – original draft. SN: Investigation, Visualization, Writing – original draft. RM: Investigation, Methodology, Resources, Writing – review & editing. BS: Data curation, Formal analysis, Methodology, Writing – review & editing. BC: Data curation, Methodology, Software, Writing – review & editing. TI: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Validation, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study has been funded by the Instituto de Salud Carlos III (ISCIII) through the project PI21/00588 and co-funded by the European Union (EU).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2025.1576906/full#supplementary-material

References

  • 1

    AggarwalS.YadavA. K. (2016). False discovery rate estimation in proteomics. Methods Mol. Biol.1362, 119128. doi: 10.1007/978-1-4939-3106-4_7

  • 2

    Agra-BermejoR. M.Cacho-AntonioC.Rozados-LuisA.Couselo-SeijasM.FernandezA. L.Martinez-CereijoJ. M.et al. (2020). CD5L, macrophage apoptosis inhibitor, was identified in epicardial fat-secretome and regulated by isoproterenol from patients with heart failure. Front. Physiol.11. doi: 10.3389/fphys.2020.00620

  • 3

    AhmadP.Escalante-HerreraA.MarinL. M.SiqueiraW. L. (2025). Progression from healthy periodontium to gingivitis and periodontitis: Insights from bioinformatics-driven proteomics - A systematic review with meta-analysis. J. Periodontal. Res. 60, 829. doi: 10.1111/jre.13313

  • 4

    ÁlvarezV. J.BravoS. B.Chantada-VazquezM. P.ColónC.De CastroM. J.MoralesM.et al. (2020). Characterization of new proteomic biomarker candidates in mucopolysaccharidosis Type IVA. Int. J. Mol. Sci.22, 226. doi: 10.3390/ijms22010226

  • 5

    AndersonM. J. (2001). A new method for non-parametric multivariate analysis of variance. Austral. Ecol.26, 3246. doi: 10.1111/j.1442-9993.2001.01070.pp.x

  • 6

    AnjoS. I.SantaC.ManadasB. (2017). SWATH-MS as a tool for biomarker discovery: From basic research to clinical applications. Proteomics17, 1600278n/a. doi: 10.1002/pmic.201600278

  • 7

    ArataniY. (2018). Myeloperoxidase: Its role for host defense, inflammation, and neutrophil function. Arch. Biochem. Biophys.640, 4752. doi: 10.1016/j.abb.2018.01.004

  • 8

    BalibanR. C.SakellariD.LiZ.DiMaggioP. A.GarciaB. A.FloudasC. A. (2012). Novel protein identification methods for biomarker discovery via a proteomic analysis of periodontally healthy and diseased gingival crevicular fluid samples. J. Clin. Periodontol.39, 203212. doi: 10.1111/j.1600-051X.2011.01805.x

  • 9

    BeinrohrL.Murray-RustT. A.DyksterhuisL.ZávodszkyP.GálP.PikeR. N.et al. (2011). Serpins and the complement system. Methods Enzymol.499, 5575. doi: 10.1016/B978-0-12-386471-0.00004-3

  • 10

    BelleiE.BertoldiC.MonariE.BergaminiS. (2022). Proteomics disclose the potential of gingival crevicular fluid (GCF) as a source of biomarkers for severe periodontitis. Materials15, 2161. doi: 10.3390/ma15062161

  • 11

    BenjaminiY.HochbergY. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc B. Stat. Methodol.57, 289300. doi: 10.1111/j.2517-6161.1995.tb02031.x

  • 12

    BernierG. M. (1980). β2-microglobulin: Structure, function and significance. Vox. Sang.38, 323327. doi: 10.1111/j.1423-0410.1980.tb04500.x

  • 13

    BharadwajA.BydounM.HollowayR.WaismanD. (2013). Annexin A2 heterotetramer: structure and function. Int. J. Mol. Sci.14, 62596305. doi: 10.3390/ijms14036259

  • 14

    Blanco-PintosT.Regueira-IglesiasA.KuzI.Sánchez-BarcoA.Seijas-OteroN.Chantada-VázquezM. P.et al. (2025b). Impact of smoking habit on the subgingival proteome in patients with periodontitis. J. Periodontol.96, 217229. doi: 10.1002/JPER.24-0062

  • 15

    Blanco-PintosT.Regueira-IglesiasA.RelvasM.Alonso-SampedroM.BravoS. B.Balsa-CastroC.et al. (2025a). Diagnostic accuracy of novel protein biomarkers in saliva to detect periodontitis using untargeted ‘SWATH’ mass spectrometry. J. Clin. Periodontol.52, 199214. doi: 10.1111/jcpe.14103

  • 16

    Blanco-PintosT.Regueira-IglesiasA.RelvasM.Alonso-SampedroM.Chantada-VázquezM. P.Balsa-CastroC.et al. (2024). Using SWATH-MS to identify new molecular biomarkers in gingival crevicular fluid for detecting periodontitis and its response to treatment. J. Clin. Periodontol.51, 13421358. doi: 10.1111/jcpe.14037

  • 17

    Bonzon-KulichenkoE.Pérez-HernándezD.NúñezE.Martínez-AcedoP.NavarroP.Trevisan-HerrazM.et al. (2011). A robust method for quantitative high-throughput analysis of proteomes by 18O labeling. Mol. Cell. Proteomics10, M110.003335. doi: 10.1074/mcp.M110.003335

  • 18

    BostanciN.BaoK. (2017). Contribution of proteomics to our understanding of periodontal inflammation. Proteomics17, 1500518n/a. doi: 10.1002/pmic.201500518

  • 19

    BostanciN.BelibasakisG. N. (2018). Gingival crevicular fluid and its immune mediators in the proteomic era. Periodontol. 2000.76, 6884. doi: 10.1111/prd.12154

  • 20

    BostanciN.SelevsekN.WolskiW.GrossmannJ.BaoK.WahlanderA.et al. (2018). Targeted proteomics guided by label-free quantitative proteome analysis in saliva reveal transition signatures from health to periodontal disease. Mol. Cell. Proteomics.17, 13921409. doi: 10.1074/mcp.RA118.000718

  • 21

    Bugallo-CasalA.MuiñoE.BravoS. B.HervellaP.Arias-RivasS.Rodríguez-YáñezM.et al. (2025). NOTCH3 variant position affects the phenotype at the pluripotent stem cell level in CADASIL. Neuromol. Med.27, 18. doi: 10.1007/s12017-025-08840-6

  • 22

    BuzalafM. A. R.OrtizA. C.CarvalhoT. S.FidelesS. O. M.AraújoT. T.MoraesS. M.et al. (2020). Saliva as a diagnostic tool for dental caries, periodontal disease and cancer: is there a need for more biomarkers? Expert. Rev. Mol. Diagn.20, 543555. doi: 10.1080/14737159.2020.1743686

  • 23

    CaminoT.Lago-BaameiroN.BravoS. B.Molares-VilaA.SueiroA.CoutoI.et al. (2022). Human obese white adipose tissue sheds depot-specific extracellular vesicles and reveals candidate biomarkers for monitoring obesity and its comorbidities. Transl. Res.239, 85102. doi: 10.1016/j.trsl.2021.01.006

  • 24

    CaminoT.Lago-BaameiroN.BravoS. B.SueiroA.CoutoI.SantosF.et al. (2020). Vesicles shed by pathological murine adipocytes spread pathology: Characterization and functional role of insulin resistant/hypertrophied adiposomes. Int. J. Mol. Sci.21, 2252. doi: 10.3390/ijms21062252

  • 25

    CarneiroL. G.NouhH.SalihE. (2014). Quantitative gingival crevicular fluid proteome in health and periodontal disease using stable isotope chemistries and mass spectrometry. J. Clin. Periodontol.41, 733747. doi: 10.1111/jcpe.12262

  • 26

    Casado-FernándezL.Laso-GarcíaF.PiniellaD.Gómez-de FrutosM. C.Otero-OrtegaL.BravoS. B.et al. (2024). The proteomic signature of circulating extracellular vesicles following intracerebral hemorrhage: Novel insights into mechanisms underlying recovery. Neurobiol. Dis.201, 106665. doi: 10.1016/J.NBD.2024.106665

  • 27

    CasarinR. C. V.SalmonC. R.StolfC. S.PazH. E. S.RangelT. P.DominguesR. R.et al. (2023). Salivary annexin A1: A candidate biomarker for periodontitis. J. Clin. Periodontol.50, 942951. doi: 10.1111/jcpe.13803

  • 28

    CatonJ. G.ArmitageG.BerglundhT.ChappleI. L. C.JepsenS.KornmanK. S.et al. (2018). A new classification scheme for periodontal and peri-implant diseases and conditions – Introduction and key changes from the 1999 classification. J. Clin. Periodontol.45, S1S8. doi: 10.1111/jcpe.12935

  • 29

    Chantada-VázquezM. P.Castro LópezA.García VenceM.Vázquez-EstévezS.Acea-NebrilB.CalatayudD. G.et al. (2020). Proteomic investigation on bio-corona of Au, Ag and Fe nanoparticles for the discovery of triple negative breast cancer serum protein biomarkers. J. Proteomics.212, 103581. doi: 10.1016/j.jprot.2019.103581

  • 30

    CohenJ. (1988). Statistical power analysis for the behavioral sciences. 2.ed (New York: Erlbaum).

  • 31

    ColangeloC. M.IvosevG.ChungL.AbbottT.ShifmanM.SakaueF.et al. (2015). Development of a highly automated and multiplexed targeted proteome pipeline and assay for 112 rat brain synaptic proteins. Proteomics15, 12021214. doi: 10.1002/pmic.201400353

  • 32

    DemichevV.MessnerC. B.VernardisS. I.LilleyK. S.RalserM. (2020). DIA-NN: neural networks and interference correction enable deep proteome coverage in high throughput. Nat. Methods17, 4144. doi: 10.1038/s41592-019-0638-x

  • 33

    EkeP. I.BorgnakkeW. S.GencoR. J. (2020). Recent epidemiologic trends in periodontitis in the USA. Periodontol. 2000.82, 257267. doi: 10.1111/prd.12323

  • 34

    EricksenB.WuZ.LuW.LehrerR. I. (2005). Antibacterial activity and specificity of the six human {alpha}-defensins. Antimicrob. Agents Chemother.49, 269275. doi: 10.1128/AAC.49.1.269-275.2005

  • 35

    FaulF.ErdfelderE.LangA.BuchnerA. (2007). G*Power 3: A flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav. Res. Methods39, 175191. doi: 10.3758/bf03193146

  • 36

    FrederickK.CiborowskiP. (2016). “SWATH-MS: Data acquisition and analysis,” in Proteomic Profiling and Analytical Chemistry. Eds. CiborowskiP.SilberringJ. (Amsterdam: Elsevier), 161173. doi: 10.1016/B978-0-444-63688-1.00009-4

  • 37

    GaoS.ZhuH.ZuoX.LuoH. (2018). Cathepsin G and its role in inflammation and autoimmune diseases. Arch. Rheumatol.33, 498504. doi: 10.5606/ArchRheumatol.2018.6595

  • 38

    García-VenceM.Chantada-VázquezM. P.Cameselle-TeijeiroJ. M.BravoS. B.NúñezC. (2020). A novel nanoproteomic approach for the identification of molecular targets associated with thyroid tumors. Nanomaterials10, 2370. doi: 10.3390/nano10122370

  • 39

    GrantM. M.TaylorJ. J.JaedickeK.CreeseA.GowlandC.BurkeB.et al. (2022). Discovery, validation, and diagnostic ability of multiple protein-based biomarkers in saliva and gingival crevicular fluid to distinguish between health and periodontal diseases. J. Clin. Periodontol.49, 622632. doi: 10.1111/jcpe.13630

  • 40

    GroegerS.MeyleJ. (2019). Oral mucosal epithelial cells. Front. Immunol.10. doi: 10.3389/fimmu.2019.00208

  • 41

    GürsoyU. K.KantarciA. (2022). Molecular biomarker research in periodontology: A roadmap for translation of science to clinical assay validation. J. Clin. Periodontol.49, 556561. doi: 10.1111/jcpe.13617

  • 42

    GuzmanY. A.SakellariD.ArsenakisM.FloudasC. A. (2014). Proteomics for the discovery of biomarkers and diagnosis of periodontitis: A critical review. Expert. Rev. Proteomics.11, 3141. doi: 10.1586/14789450.2014.864953

  • 43

    GuzmanY. A.SakellariD.PapadimitriouK.FloudasC. A. (2018). High-throughput proteomic analysis of candidate biomarker changes in gingival crevicular fluid after treatment of chronic periodontitis. J. Periodontal. Res.53, 853860. doi: 10.1111/jre.12575

  • 44

    HartenbachF. A. R. R.VelasquezÉ.NogueiraF. C. S.DomontG. B.FerreiraE.ColomboA. P. V. (2020). Proteomic analysis of whole saliva in chronic periodontitis. J. Proteomics.213, 103602. doi: 10.1016/j.jprot.2019.103602

  • 45

    JayasenaT.PoljakA.BraidyN.ZhongL.RowlandsB.MuenchhoffJ.et al. (2016). Application of targeted mass spectrometry for the quantification of sirtuins in the central nervous system. Sci. Rep.6, 35391. doi: 10.1038/srep35391

  • 46

    KassambaraA.MundtF. (2020). factoextra: Extract and visualize the results of multivariate data analyses. R package. Available online at: https://CRAN.R-project.org/package=factoextra (Accessed November 14, 2024).

  • 47

    KimJ. S.ChoI.KimK.HwangY. S. (2021). Identification of galectin-10 as a biomarker for periodontitis based on proteomic analysis of gingival crevicular fluid. Mol. Med. Rep.23, 123. doi: 10.3892/mmr.2020.11762

  • 48

    KinaneD. F.StathopoulouP. G.PapapanouP. N. (2017). Periodontal diseases. Nat. Rev. Dis. Primers.3, 1703917038. doi: 10.1038/nrdp.2017.38

  • 49

    S.JosseJ.HussonF. (2008). FactoMineR: An R package for multivariate analysis. J. Stat. Software25, 118. doi: 10.18637/jss.v025.i01

  • 50

    López-ValverdeL.Vázquez-MosqueraM. E.Colón-MejerasC.BravoS. B.Barbosa-GouveiaS.ÁlvarezJ. V.et al. (2024). Characterization of the plasma proteomic profile of Fabry disease: Potential sex- and clinical phenotype-specific biomarkers. Transl. Res.269, 4763. doi: 10.1016/j.trsl.2024.02.006

  • 51

    LudwigC.GilletL.RosenbergerG.AmonS.CollinsB. C.AebersoldR. (2018). Data-independent acquisition-based SWATH-MS for quantitative proteomics: a tutorial. Mol. Syst. Biol.14 Suppl, e8126e8n/a. doi: 10.15252/msb.20178126

  • 52

    MarinhoM. C.PachecoA. B. F.CostaG. C. V.OrtizN. D.ZajdenvergL.SansoneC. (2019). Quantitative gingival crevicular fluid proteome in type 2 diabetes mellitus and chronic periodontitis. Dis.25, 588595. doi: 10.1111/odi.12996

  • 53

    MeghaK. B.MohananP. V. (2021). Role of immunoglobulin and antibodies in disease management. Int. J. Biol. Macromol.169, 2838. doi: 10.1016/j.ijbiomac.2020.12.073

  • 54

    MeijerB.GearryR. B.DayA. S. (2012). The role of S100A12 as a systemic marker of inflammation. Int. J. Inflam.2012, 907078907076. doi: 10.1155/2012/907078

  • 55

    MertensB.OrtiV.VialaretJ.GibertP.Relaño-GinésA.LehmannS.et al. (2018). Assessing a multiplex-targeted proteomics approach for the clinical diagnosis of periodontitis using saliva samples. Bioanalysis10, 3545. doi: 10.4155/bio-2017-0218

  • 56

    Mondelo-MacíaP.García-GonzálezJ.León-MateosL.AbaloA.BravoS.Chantada VazquezM. P.et al. (2024). Identification of a proteomic signature for predicting immunotherapy response in patients with metastatic non-small cell lung cancer. Mol. Cell. Proteomics.23, 100834. doi: 10.1016/J.MCPRO.2024.100834

  • 57

    MonsarratP.BlaizotA.KémounP.RavaudP.NabetC.SixouM.et al. (2016). Clinical research activity in periodontal medicine: A systematic mapping of trial registers. J. Clin. Periodontol.43, 390400. doi: 10.1111/jcpe.12534

  • 58

    MullaS. A.PatilA.MaliS.JainA.SharmaD.JaiswalH. C.et al. (2023). Exploring the biological width in dentistry: A comprehensive narrative review. Curēus15, e42080. doi: 10.7759/cureus.42080

  • 59

    MurakamiM.Lopez-GarciaB.BraffM.DorschnerR. A.GalloR. L. (2004). Postsecretory processing generates multiple cathelicidins for enhanced topical antimicrobial defense. J. Immunol.172, 30703077. doi: 10.4049/jimmunol.172.5.3070

  • 60

    NavazeshM.ChristensenC. M. (1982). A comparison of whole mouth resting and stimulated salivary measurement procedures. J. Dent. Res.61, 11581162. doi: 10.1177/00220345820610100901

  • 61

    NovelleM. G.BravoS. B.DeshonsM.IglesiasC.García-VenceM.AnnellsR.et al. (2021). Impact of liver-specific GLUT8 silencing on fructose-induced inflammation and omega oxidation. iScience24, 102071. doi: 10.1016/j.isci.2021.102071

  • 62

    NóvoaE.da Silva LimaN.Gonzalez-RellanM. J.Chantada-VazquezM. P.VerheijJ.RodriguezA.et al. (2025). Mitochondrial antiviral signaling protein enhances MASLD progression through the ERK/TNFα/NFκβ pathway. Hepatology. 81:1535–1552. doi: 10.1097/hep.0000000000000930

  • 63

    OksanenJ.SimpsonG.BlanchetF.KindtR.LegendreP.MinchinP.et al. (2022). vegan: Community ecology package. R package. Available online at: https://CRAN.R-project.org/package=vegan (Accessed November 14, 2024).

  • 64

    OngS. T.Shan HoJ. Z.HoB.DingJ. L. (2006). Iron-withholding strategy in innate immunity. Immunobiology211, 295314. doi: 10.1016/j.imbio.2006.02.004

  • 65

    OrteaI.Ruiz-SánchezI.CañeteR.Caballero-VillarrasoJ.CañeteM. D. (2018). Identification of candidate serum biomarkers of childhood-onset growth hormone deficiency using SWATH-MS and feature selection. J. Proteomics.175, 105113. doi: 10.1016/j.jprot.2018.01.003

  • 66

    Perez-HernandezD.Gutiérrez-VázquezC.JorgeI.López-MartínS.UrsaA.Sánchez-MadridF.et al. (2013). The intracellular interactome of tetraspanin-enriched microdomains reveals their function as sorting machineries toward exosomes. J. Biol. Chem.288, 1164911661. doi: 10.1074/jbc.M112.445304

  • 67

    R Core Team (2022). R: A language and environment for statistical computing (Vienna, Austria: R Foundation for Statistical Computing). Available at: https://www.R-project.org/ (Accessed November 14, 2024).

  • 68

    RizalM. I.SoerosoY.SulijayaB.AssiddiqB. F.BachtiarE. W.BachtiarB. M. (2020). Proteomics approach for biomarkers and diagnosis of periodontitis: systematic review. Heliyon6, e04022. doi: 10.1016/j.heliyon.2020.e04022

  • 69

    RodriguesJ. S.ChenloM.BravoS. B.Perez-RomeroS.Suarez-FariñaM.SobrinoT.et al. (2024). dsRNAi-mediated silencing of PIAS2beta specifically kills anaplastic carcinomas by mitotic catastrophe. Nat. Commun.15, 3736. doi: 10.1038/s41467-024-47751-1

  • 70

    RosenbergerG.KohC. C.GuoT.RöstH. L.KouvonenP.CollinsB. C.et al. (2014). A repository of assays to quantify 10,000 human proteins by SWATH-MS. Sci. Data.1, 140031. doi: 10.1038/sdata.2014.31

  • 71

    SalazarM. G.JehmlichN.MurrA.DhopleV. M.HoltfreterB.HammerE.et al. (2013). Identification of periodontitis associated changes in the proteome of whole human saliva by mass spectrometric analysis. J. Clin. Periodontol.40, 825832. doi: 10.1111/jcpe.12130

  • 72

    SampathP.MazumderB.SeshadriV.FoxP. L. (2003). Transcript-selective translational silencing by gamma interferon is directed by a novel structural element in the ceruloplasmin mRNA 3′ untranslated region. Mol. Cell. Biol.23, 15091519. doi: 10.1128/MCB.23.5.1509-1519.2003

  • 73

    ShevchenkoA.WilmM.VormO.MannM. (1996). Mass spectrometric sequencing of proteins from silver-stained polyacrylamide gels. Anal. Chem.68, 850858. doi: 10.1021/ac950914h

  • 74

    ShinM.KimY.ShinY. J.KoB. J.KimS.KimH. (2019). Deep sequencing salivary proteins for periodontitis using proteomics. Clin. Invest.23, 35713580. doi: 10.1007/s00784-018-2779-1

  • 75

    SilbereisenA.BaoK.WolskiW.NanniP.KunzL.AfacanB.et al. (2024). Probing the salivary proteome for prognostic biomarkers in response to non-surgical periodontal therapy. J. Clin. Periodontol.52, 5667. doi: 10.1111/jcpe.13990

  • 76

    Silva-BoghossianC. M.ColomboA. P. V.TanakaM.RayoC.XiaoY.SiqueiraW. L. (2013). Quantitative proteomic analysis of gingival crevicular fluid in different periodontal conditions. PloS One8, e75898. doi: 10.1371/journal.pone.0075898

  • 77

    SmithA.McCullohR. J. (2015). Hemopexin and haptoglobin: allies against heme toxicity from hemoglobin not contenders. Front. Physiol.6. doi: 10.3389/fphys.2015.00187

  • 78

    SoehnleinO.LindbomL. (2009). Neutrophil-derived azurocidin alarms the immune system. J. Leukoc. Biol.85, 344351. doi: 10.1189/jlb.0808495

  • 79

    SreejitG.AhmedA.ParveenN.JhaV.ValluriV. L.GhoshS.et al. (2014). The ESAT-6 protein of Mycobacterium tuberculosis interacts with beta-2-microglobulin (β2M) affecting antigen presentation function of macrophage. PloS Pathog.10, e1004446. doi: 10.1371/journal.ppat.1004446

  • 80

    SunQ.LiN.JiaL.GuoW.JiangH.LiuB.et al. (2021). Ribosomal protein SA-positive neutrophil elicits stronger phagocytosis and neutrophil extracellular trap formation and subdues pro-inflammatory cytokine secretion against Streptococcus suis serotype 2 infection. Front. Immunol.11. doi: 10.3389/fimmu.2020.585399

  • 81

    SurnaA.KubiliusR.SakalauskieneJ.VitkauskieneA.JonaitisJ.SaferisV.et al. (2009). Lysozyme and microbiota in relation to gingivitis and periodontitis. Med. Sci. Monit.15, CR66CR73.

  • 82

    TangH.YuanC.MaZ.ZhuC.TongP.GallagherJ. E.et al. (2019). The potentiality of salivary peptide biomarkers for screening patients with periodontal diseases by mass spectrometry. Clin. Chim. Acta495, 278286. doi: 10.1016/j.cca.2019.04.076

  • 83

    TeamM. C.BlanchardG.DickhausT.HackN.KonietschkeF.RohmeyerK.et al. (2023). mutoss: Unified multiple testing procedures (Wien: R package). Available at: https://CRAN.R-project.org/package=mutoss (Accessed November 14, 2024).

  • 84

    TisdaleE. J. (2002). Glyceraldehyde-3-phosphate dehydrogenase is phosphorylated by protein kinase ciota/lambda and plays a role in microtubule dynamics in the early secretory pathway. J. Biol. Chem.277, 3334. doi: 10.1074/jbc.M109744200

  • 85

    TonettiM. S.GreenwellH.KornmanK. S. (2018). Staging and grading of periodontitis: Framework and proposal of a new classification and case definition. J. Clin. Periodontol.45, S149S161. doi: 10.1111/jcpe.12945

  • 86

    TorchianoM. (2020). effsize: Efficient effect size computation (Wien: R package). Available at: https://CRAN.R-project.org/package=effsize (Accessed November 14, 2024).

  • 87

    Torres IglesiasG.López-MolinaM. P.BotellaL.Laso-GarcíaF.ChamorroB.Fernández-FournierM.et al. (2024). Differential protein expression in extracellular vesicles defines treatment responders and non-responders in multiple sclerosis. Int. J. Mol. Sci.25, 10761. doi: 10.3390/ijms251910761

  • 88

    TripathiD.KantS.PandeyS.EhteshamN. Z. (2020). Resistin in metabolism, inflammation, and disease. Febs. J.287, 31413149. doi: 10.1111/febs.15322

  • 89

    UniProt Consortium (2014). Activities at the universal protein resource (UniProt). Nucleic Acids Res.42, D191D198. doi: 10.1093/nar/gkt1140

  • 90

    Vázquez-MeraS.Martelo-VidalL.Miguéns-SuárezP.BravoS. B.Saavedra-NievesP.AriasP.et al. (2024). Exploring CD26–/lo subpopulations of lymphocytes in asthma phenotype and severity: A novel CD4+ T cell subset expressing archetypical granulocyte proteins. Allergy79, 30053021. doi: 10.1111/all.16327

  • 91

    WangF.QiaoL.LvX.TrivettA.YangR.OppenheimJ. J.et al. (2016). Alarmin human α defensin HNP1 activates plasmacytoid dendritic cells by triggering NF-κB and IRF1 signaling pathways. Cytokine83, 5360. doi: 10.1016/j.cyto.2016.03.015

  • 92

    WangS.SongR.WangZ.JingZ.WangS.MaJ. (2018). S100A8/A9 in inflammation. Front. Immunol.9. doi: 10.3389/fimmu.2018.01298

  • 93

    World Medical Association (2013). World medical association declaration of Helsinki: Ethical principles for medical research involving human subjects. JAMA310, 21912194. doi: 10.1001/jama.2013.281053

  • 94

    YiJ.ShenY.YangY.ShenC.LiuB.QiaoL.et al. (2021). Direct MALDI-TOF profiling of gingival crevicular fluid sediments for periodontitis diagnosis. Talanta225, 121956. doi: 10.1016/j.talanta.2020.121956

  • 95

    ZhangX.XiaoX.MuY.LiuY. R.LinX.LiQ.et al. (2022). A comparative proteomic analysis to define the influencing factors on gingival crevicular fluid using LC-MS/MS. J. Proteomics.252, 104421. doi: 10.1016/j.jprot.2021.104421

Summary

Keywords

periodontitis, gingival crevicular fluid, saliva, oral fluids, protein expression, SWATH-MS, mass spectrometry, proteomics

Citation

Blanco-Pintos T, Regueira-Iglesias A, Suárez-Rodríguez B, Seijas-Otero N, Relvas M, Bravo SB, Balsa-Castro C and Tomás I (2025) Characterisation of the periodontal proteome in gingival crevicular fluid and saliva using SWATH-MS. Front. Cell. Infect. Microbiol. 15:1576906. doi: 10.3389/fcimb.2025.1576906

Received

28 February 2025

Accepted

07 April 2025

Published

02 May 2025

Volume

15 - 2025

Edited by

Lucinda Janete Bessa, Egas Moniz Center for Interdisciplinary Research (CiiEM), Portugal

Reviewed by

Ganesan Saibaba, Hebrew University of Jerusalem, Israel

Stefania Bergamini, University of Modena and Reggio Emilia, Italy

Updates

Copyright

*Correspondence: Inmaculada Tomás, ; Alba Regueira-Iglesias,

†ORCID: Alba Regueira-Iglesias, orcid.org/0000-0002-6549-7738; Inmaculada Tomás, orcid.org/0000-0002-3317-0853

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics