Next Article in Journal
Physical Activity and Inhibitory Control: The Mediating Role of Sleep Quality and Sleep Efficiency
Previous Article in Journal
Web-Based Assessment of the Phenomenology of Autobiographical Memories in Young and Older Adults
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Brain Sciences
Volume 11
Issue 5
10.3390/brainsci11050661
brainsci-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Proteomic Profile of Saliva in Parkinson’s Disease Patients: A Proof of Concept Study

by
Monika Figura
1,*,
Ewa Sitkiewicz
2,
Bianka Świderska
2,
Łukasz Milanowski
1,
Stanisław Szlufik
1,
Dariusz Koziorowski
1 and
Andrzej Friedman
1
1
Department of Neurology, Faculty of Health Sciences, Medical University of Warsaw, 03-242 Warsaw, Poland
2
Mass Spectrometry Laboratory, Institute of Biochemistry and Biophysics Polish Academy of Sciences, 02-106 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Brain Sci. 2021, 11(5), 661; https://doi.org/10.3390/brainsci11050661
Submission received: 10 April 2021 / Revised: 11 May 2021 / Accepted: 13 May 2021 / Published: 18 May 2021
(This article belongs to the Section Neurodegenerative Diseases)
Browse Figures
Versions Notes

Abstract

:
Parkinson’s disease (PD) is a progressive neurodegenerative disorder. It affects many organs. Lewy bodies—a histopathological “hallmark” of PD—are detected in about 75% of PD submandibular gland samples. We hypothesize that saliva can be a source of biomarkers of PD. The aim of the study was to evaluate and compare the salivary proteome of PD patients and healthy controls (HC). Salivary samples from 39 subjects (24 PD patients, mean age 61.6 ± 8.2; 15 HC, mean age 60.9 ± 6.7) were collected. Saliva was collected using RNA-Pro-Sal kits. Label-free LC-MS/MS mass spectrometry was performed to characterize the proteome of the saliva. IPA analysis of upstream inhibitors was performed. A total of 530 proteins and peptides were identified. We observed lower concentrations of S100-A16, ARP2/3, and VPS4B in PD group when compared to HC. We conclude that the salivary proteome composition of PD patients is different than that of healthy controls. We observed a lower concentration of proteins involved in inflammatory processes, exosome formation, and adipose tissue formation. The variability of expression of proteins between the two groups needs to be considered.

1. Introduction

Despite being the second most frequent neurodegenerative disorder, Parkinson’s disease (PD) is still diagnosed in its late motor stage. The diagnosis becomes possible only when the pathology has spread throughout the nervous system, causing damage to over 50% of the substantia nigra cells [1]. The only route to early intervention in PD treatment is through the development of biomarkers. The current availability of biomarkers facilitating early PD diagnosis is limited. Diagnosis is established with limited accuracy on the basis of clinical symptoms of bradykinesia with concomitant muscular rigidity or rest tremor or postural instability [2]. Diagnostic tools such as a DAT scan, transcranial ultrasonography of the substantia nigra, olfactory test, or autonomic assessment may be of help, but due to high costs and limited accessibility they cannot be used for routine screening of all suspected cases.
Looking for biomarkers of PD-related degeneration seems to be the crucial next step in its early diagnosis. Current efforts such as the BIOFIND initiative or Systemic Synuclein Sampling Study are directed towards the identification of tissue that may serve as a biomarker source [3,4]. Organs such as salivary glands, gastroenteric plexuses, adrenal glands, or the urinary system seem to be involved in PD pathology long before motor symptoms develop [5,6,7,8]. They are also more feasible for sample collection than the CNS nervous tissue. Staining for alpha-synuclein deposits (Lewy bodies and neurites) performed in samples from peripheral tissues is currently being investigated as a biomarker of PD. Protein misfolding cyclic amplification (PMCA) and real-time quaking-induced conversion (RT-QuIC) techniques seem to be promising, new, ultrasensitive tools for the detection of very small amounts of alpha-synuclein prone to forming aggregates [9].
Saliva production is abnormal in PD, with higher total protein concentrations [10]. In addition to drooling, which is commonly associated with PD, xerostomia is also frequent among PD patients, being reported in 60.8% of cases [11]. Detailed findings on salivary protein changes in PD were summarized in a paper by our group [12]. We concluded that higher concentrations of oligomeric alpha-synuclein seem to be the most promising salivary biomarker of PD [13,14]. Oligomeric alpha-synuclein seems to be a key player in PD-related neurodegeneration. It is involved in mitochondrial dysfunction and inefficient removal of misfolded proteins by proteasomes, synaptic dysfunction, and neuroinflammation, among others [15,16,17]. Its higher expression in saliva supports the idea that the neurodegenerative process in PD is generalized and may be reflected in saliva composition.
Human saliva has several functions. It contains many proteins and peptides as its components. Its organic components include amylases, cystatins, hormones, lysozyme, lipase, lactoferrins, mucins, peroxidase, and growth factors, among other proteins [18,19].
Unstimulated saliva production is secreted mostly by the submandibular glands. The submandibular glands’ saliva differs from the parotid glands—it is more viscous and mucin-rich [19]. The submandibular glands of PD patients also have the highest prevalence of Lewy-type synucleinopathy, ranging from 74% to 100%, depending on the tissue collection method [20,21,22,23]. This could provide a rationale for the preference of unstimulated saliva over the stimulated, watery saliva produced by parotid glands in the search for PD biomarkers.
LC-MS/MS mass spectrometry is a useful tool for biomarker candidate searches as it can identify large numbers of proteins across a large dynamic range and results in hypothesis-free analysis driven by the obtained data. In this study, we hypothesized that the involvement of salivary glands in synucleinopathy, as well as aberrations in protein secretion, may be reflected by changes in the saliva proteome of PD patients. Such changes in the composition of saliva could then be developed into a biomarker. The proteomic approach to saliva was initially introduced in studies of dental and oral diseases and is now applied in generalized conditions such as oncological disorders, addictions, or pediatric disorders [24,25,26,27]. The aim of this study is to compare the proteome of PD saliva vs. healthy controls (HC), as well as to further analyze specific pathways leading to different expression of proteins.

2. Materials and Methods

2.1. Study Group Characteristics and Material Collection

A total of 50 subjects were recruited for the study. Eleven out of fifty patients (9 controls and 2 in the PD group) were unable to provide enough saliva during the collection procedure.
The remaining 39 subjects recruited for the study were divided into two groups: 24 PD patients and 15 HC. All patients in the PD group were diagnosed in accordance with UK Brain Bank criteria for Parkinson’s disease. Detailed characteristics of the study groups are presented in Table 1. The mean duration of the disease in the PD group was 8.2 +/− 4.3 years. In total, 7 patients in the PD group were treated with a dopamine agonist (5-ropinirole, 5 piribedil, 1-rotigotine) in addition to levodopa treatment.
The study was approved by the Ethics Committee of the Medical University of Warsaw (KB/239/2015, with amendment KB/25/A/2016), and therefore was performed in accordance with the ethical standards laid down in the 1964 Declaration of Helsinki and its later amendments. All participants signed informed consent forms prior to their inclusion in the study. Participants were recruited between April 2018 and February 2019. Patients were recruited in the Department of Neurology, Faculty of Health Sciences.
Exclusion criteria included cigarette smoking, previous injection of botulinum toxin to salivary glands, treatment with anticholinergic medications, clinical diagnosis of neurodegenerative diseases other than PD, and known malignancies. Among PD patients, only levodopa and oral dopamine agonist treatments were allowed. HC were subjects with no clinical evidence of neurodegeneration, matched by age and sex to the PD group. The control group was screened and confirmed negative for symptoms of bradykinesia, tremor, imbalance, or rigidity, as well as REM behavior disorder history. The control group was recruited among invited individuals (site personnel and non-blood-related members of patients’ families or caregivers). No patients who had a periodontist, active inflammation, or oral cancer history were included.
Subjects were asked to refrain from drinking, eating, or using oral hygiene procedures for at least 2 h before the procedure and to rinse their mouths with tap water 30–60 min prior to collection. Saliva samples were collected in the morning hours, using RNA-Pro-Sal kits. The collection procedure with this device was described in detail by Chiang et al. [28]. Samples were immediately frozen at −80 °C after collection and later processed in the Mass Spectrometry Laboratory at the Institute of Biochemistry and Biophysics PAS.

2.2. Analytical Methods

2.2.1. Sample Preparation for Mass Spectrometry Analysis

Mass spectrometry experiments were performed at the Mass Spectrometry Laboratory at the Institute of Biochemistry and Biophysics PAS. Freshly prepared urea buffer (9 M urea, 200 mM Tris-HCl pH 8.0 in MS-grade water) was aliquoted into 85 μL volumes and fully dried in a SpeedVac. The urea was redissolved in 50 μL of saliva to obtain 85 μL of protein sample in urea buffer. The protein concentration was measured with a BCA Protein Assay Kit (Thermo). Protein eluates were processed using single-pot, solid-phase-enhanced sample preparation (SP3) with some modifications (ultrasensitive proteome analysis using paramagnetic bead technology). Then, 40 μg of protein from each sample was transferred to a new tube and filled up to 100 μL with urea buffer. Cysteine bridges were reduced by 1 h incubation with 20 mM tris(2-carboxyethyl) phosphine (TCEP) at 37 °C followed by 30 min incubation at room temperature with 40 mM iodoacetamide (IAA). A magnetic bead mix was prepared by combining equal parts of Sera-Mag Carboxyl hydrophilic and hydrophobic particles (09-981-121 and 09-981-123, GE Healthcare). The bead mix was washed three times with MS-grade water and resuspended in a working concentration of 10 μg/μL. Then, 250 μg of the prepared bead mix, along with 5 μL of 10% formic acid and 800 μL of acetonitrile, were added to each sample. Proteins bound to beads were washed with 75% ethanol, isopropanol, and acetonitrile, followed by overnight digestion with 2 μg of trypsin/Lys-C mix (Promega). After digestion, peptides were washed with acetonitrile and eluted from the beads by subsequent incubation with MS-grade water and 2% DMSO with sonication during each step. Pulled aliquots were dried in a SpeedVac and resuspended in 40 μL 2% acetonitrile and 0.1% formic acid.

2.2.2. Mass Spectrometry

Here, 2 ug of each saliva sample was analyzed using a nanoAcquity UPLC (Waters) directly coupled to a QExactive mass spectrometer (Thermo Scientific, Bremen, Germany). Peptides were trapped on a C18 precolumn (180 µm × 20 mm, Waters) with 0.1% FA in water as a mobile phase and transferred to a nanoAcquity BEH C18 column (75 µm × 250 mm, 1.7 µm, Waters) using ACN gradient (0–35% ACN in 160 min) in the presence of 0.1% FA at a flow rate of 250 nL/min. Measurements were performed in data-dependent mode with top 12 precursors selected for MS2. Full MS scans covering the range of 300–1650 m/z were acquired at a resolution of 70,000, with a maximum injection time of 60 ms and an AGC target value of 1e6. MS2 scans were acquired at a resolution of 17,500 and an AGC target value of 5e5. Dynamic exclusion was set to 30 s.

2.2.3. Data Analysis

Obtained data were pre-processed with Mascot Distiller software (Matrixscience) and protein identification was performed using a Mascot Server 2.5 (Matrixscience) against the Homo sapiens protein sequences (20,490 sequences) deposited in the Swiss-Prot database (201,903, 559,634 sequences; 201,129,965 residues). The parameters were set as follows: enzyme—Trypsin; missed cleavages—2; fixed modifications—carbamidomethyl (C); variable modifications—oxidation (M); instrument—HCD. To reduce mass errors, peptide and fragment mass tolerance settings were established separately for each file after an off-line mass recalibration [29]. The assessment of confidence was based on a target–decoy database search strategy, as described by Elias et al., which provided q-value estimates for each peptide spectrum match [30,31]. All queries with q-values > 0.01, subset proteins, and proteins identified with one peptide were discarded from further analysis. The mass recalibration, FDR computations, and data filtering were done with Mscan software, developed in-house [32].

2.2.4. Quantitative MS Data Processing

The lists of identified peptides were merged into one common list and overlayed onto 2-D heatmaps generated from LC-MS spectra and the volumes were obtained from the assigned peaks (a more detailed description of data extraction procedures can be found in [33]). The abundance of each peptide was determined as the height of a 2-D fit to the monoisotopic peak of the tagged isotopic envelope. Quantitative values were then exported into text files for statistical analysis with Diffprot software for non-parametric statistical analysis of differential proteomics data (Malinowska et al., 2012). Diffprot is an in-house software for statistical significance assessment. In this program, the statistical validity of the regulation or expression status of a protein represented by its calculated protein ratio is based solely on the statistical analysis of the datasets from a given experiment, without assumptions on the character of the distribution of peptide ratios in a dataset (e.g., its normality). The probability of obtaining a given protein ratio by random selection from the dataset was tested by calculating protein ratios for a large number of permuted decoy datasets in which the peptide–protein assignment was scrambled. Unfortunately, calculated p-values in this dataset were non-significant for identified proteins, so it was decided to report only raw p-values. Diffprot was run with the following parameters: number of random peptide sets = 106; clustering of peptide sets—only when 90% identical; normalization by LOWESS, min-pep 4, quantification based on unique peptides. Only proteins with fold changes >1.5 were taken into consideration during further analysis.
The ROC plots were obtained by plotting all sensitivity values (true positive rate— TPR) on the y-axis against their equivalent (1-specificity) values (false positive rate—FPR) for all available thresholds on the x-axis using ROCit 2.1.1 software (authors: R.A Khan, T. Brandenburger). The area under the curve (AUC) was calculated to provide a summary of overall potential marker effectiveness. An optimal point was chosen as described by Youden [34].
Proteins with at least 1.5-fold change between groups were analyzed using Ingenuity Pathway Analysis (IPA; QIAGEN, Germany) software to identify relevant biological pathways and upstream regulators.

3. Results

A total of 1328 peptides corresponding to 530 proteins were identified. We observed a −10.47-fold change in the concentration of protein S100-A16 in the PD group vs. healthy control. We also observed changes in concentrations of proteins from the annexin family (annexin A2 (−4.4-fold change) and annexin A8 (−3.84)) in PD vs. control. The resistin concentration was 4.04 times lower in PD than in control. The proteins with the highest and lowest fold changes for MSP/HC peak areas (>1.5/<−1.5) are presented in Supplementary Table S1. Protein S100-A16, actin-related protein 2/3 complex subunit 1A (ARPC1A), and vacuolar protein sorting-associated protein 4B (VPS4B) had p-values < 0.05 in Diffprot software analysis [29]. Unfortunately, we observed heterogenous expression of salivary proteins among samples, with the mean percentage of proteins detected in each sample at around 32%. In effect, the data did not achieve enough statistical power to give significant results after FDR. Therefore, in Supplementary Table S1, we show raw p-values for proteins with a fold change >1.5 or <−1.5 in PD versus HC saliva samples. Figure 1 represents the results of AUC ROC calculated for 2 proteins selected based on high fold changes and low p-values—S100A16 and ARPC1A. The values were AUC = 0.7, specificity = 0.67, sensitivity = 0.91 for S100A16; AUC = 0.62, specificity = 1, sensitivity = 0.4 for ARPC1A. The third interesting protein with a good p-value, VPS4B, was characterized by lower AUC ROC of 0.54, with 100% sensitivity and specificity of 0.4.
To visualize the proteomic findings in samples, we provide a volcano plot reporting p-values against fold changes (Figure 2).
Next, we investigated networks of proteins with fold changes above 1.5 using IPA to further explore molecular processes possibly leading to observed changes. We performed upstream regulator analysis to predict molecules that may have been causing observed protein expression changes. Table 2 contains a list of upstream regulators identified with the highest and lowest predicted z-scores.

4. Discussion

To the best of our knowledge, this is one of the very first studies to assess the salivary proteomes of PD patients and healthy subjects. The results of our analysis, as well as previous studies on PD saliva, indicate that the salivary proteome of PD patients might differ from that of healthy controls. A study by Kumari et al. took a metabolomic approach to the assessment of saliva from PD patients. The authors described increased concentrations of N-acetylglutamate, acetoin, acetate, alanine, phenylalanine, tyrosine, histidine, glycine, acetoacetate, taurine, TMAO, GABA, fucose, propionate, isoleucine, and valine in PD patients versus healthy controls (HC). They also indicated that subgroup analysis revealed differences in metabolite concentrations in saliva among different stages of the disease [35]. Masters et al. reported preliminary results from a study involving 3 PD patients and one control. They reported upregulation of S100-A9 and S100-A8 proteins in the saliva of PD patients versus the control [36]. In our larger cohort, we observed lower concentrations of the S100A16 protein in the PD group vs. control. It also had the highest AUC of all the identified proteins. Higher detection of the S100B family of proteins is associated with numerous neurodegenerative processes, including Alzheimer’s disease, PD, and amyotrophic lateral sclerosis [37]. It was also proposed as a biomarker of neural injury and clinical outcome in stroke and trauma [38,39]. S100B protein overexpression has an established role in PD. It was associated with age of onset of PD and neuronal and glial damage [40,41]. To our best knowledge, no studies on salivary S100A16 protein measurements have been performed. S100A16 participates in calcium ion binding and adipocyte differentiation. In clinical settings, its role has mainly been investigated in oncogenesis. Specifically, S100A16 expression was upregulated in tumors of the lungs, thyroid gland, pancreas, bladder, and ovaries [42]. Active neoplastic process was, however, an exclusion criterion in our cohort. The S100A16 protein’s higher expression was also described as promoting adipocyte differentiation in animal models [43,44]. We can only speculate that lower salivary levels of this protein can reflect malnutrition and lower adipose tissue levels in the PD group [45,46].
ARP2/3 participates in polymerization of actin upon stimulation via nucleation-promoting factor. The Arp2/3 complex takes part in the formation of branched actin networks in the cytoplasm. This process provides the force for cell motility. The Arp2/3 complex is also involved in homologous recombination repair after DNA damage [47]. To our best knowledge, changes of expression of this protein were not associated previously with neurodegeneration. We observed lower expression of ARP2/3 in PD vs. HC. Animal studies prove that accumulation of alpha-synuclein can induce DNA single-strand and double-strand breaks [48,49], which are not limited to the central nervous system. We hypothesize that lower expression of the Arp2/3 complex may reflect an inefficient repair mechanism of DNA damage in the PD cohort.
VPS4B is involved in the endosomal multivesicular body pathway. It is required for the exosomal release of SDCBP, CD63, and syndecan [50]. In our cohort, lower expression of VPS4B in saliva was observed in the PD group. Simultaneous inhibition of VSP4A and VSP4B reduces exosome secretion [50]. Exosomes play a crucial role in alpha-synuclein propagation [51,52]. A pilot study by Rani et al. indicated that PD patients have increased secretion of exosomes from neuronal endings in salivary glands [53]. A lower concentration of VPS4B could perhaps reflect its higher uptake in the exosome formation process in salivary glands. In our study, although VPS4B had a very high sensitivity of 100%, this classifier’s effectiveness is insufficient.
Although the proteomic data on PD are limited, there have been several genomic studies on this disease [54,55,56]. The GWAS study by Nalls et al. revealed 90 PD-related genetic risk loci [55]. In our proteomic study, we did not identify any proteins encoded by these genes. However, the study by Nalls et al. showed the VPS13C loci as a potential PD risk factor. The VPS13C protein, similarly to the VPS4B protein discovered in our study, is involved in the endosomal multivesicular bodies pathway. It also increases PINK1/Parkin-dependent mitophagy. Additionally, Lesage et al. identified VPS13C as a gene responsible for autosomal recessive PD [57].
Interestingly, some of the proteins generally associated with inflammation had lower concentrations in PD group than in control group. This proves that oral inflammation related to worse oral hygiene in PD is not the main factor determining differences in saliva composition between groups [58]. Regarding proteins with high fold change (>4) but no significant difference in concentration between cohorts, we identified proteins from the annexin family and resistin. Inhibition of NF-κ B pathways observed in analysis of upstream regulators may coincide with lower concentrations of salivary resistin. NF-κ B plays a pivotal role in PD-related neuroinflammation. Numerous substances working through inhibition of NF-κ B pathways were investigated regarding their neuroprotective role in PD models [59,60,61,62]. Its expression is also increased in the substantia nigra of PD patients [63]. Our results indicate that this process is not reflected by saliva composition.
The high number of identified upstream regulators with high z-scores may reflect the complex nature of the disease, which is partially reflected in the salivary composition. The molecules identified as possibly driving differential expression of proteins between PD and HC groups were mainly cytokines, transcription regulators, chemicals, and kinase inhibitors. We observed an activation of an upstream regulator previously associated with PD pathogenesis: PD98059. PD98059 is an ERK1/2 signaling inhibitor with implications for antidyskinetic effects in PD [64,65].
Inhibition of prolactin-related proteins may be in line with previously reported effects of levodopa treatment on prolactin secretion [66,67,68].
We did not detect alpha-synuclein in our samples, which may be due to method limitations. Recent findings did not confirm that salivary alpha-synuclein can differentiate between HC and PD patients [4].
A limitation of our study was the high variability of salivary protein expression, with only around 30% of the total number of proteins identified in each sample. The saliva proteome may be significantly affected by such general aspects as patients’ oral cavity status, type of diet, nutrition, and comorbidities. Treatment received for PD (levodopa, domperidone, botulinum toxin) may also influence saliva production and composition [69,70]. These variables led to vast changes in the results observed among our subjects. We took measures to unify collection procedures and exclude some possibly interfering factors. Our goal was, however, to include heterogenous patients with PD without artificially narrowing study groups and to reflect the real-life proteomic profile of saliva. The findings of this study should be validated independently with Western blotting or enzyme-linked immunosorbent assay to increase its accuracy.
It remains of interest whether changes in the salivary proteome may appear ahead of motor symptoms of PD. In particular, markers of inflammation could serve as its early indicators. We can only hypothesize that this may be the case, since gastroenteric symptoms appear in PD years before motor manifestation. Salivary glands and gastroenteric plexuses are positive for Lewy pathology, with similar high frequency [21]. This “liquid biopsy” could then result in early causative treatment.
This study’s focus on saliva is both its limitation and advantage. Our study reveals several interesting directions for future studies of the salivary proteome in PD. It suggests that inflammatory processes and imbalances in cellular pathways associated with PD pathology may be reflected in the composition of saliva. The differences observed between the salivary proteomic profiles of PD patients and HCs add to the current knowledge about the state of the widespread pathology underlying the disease. However, very little is known about salivary changes in neurodegeneration; therefore, some findings are difficult to interpret. Further studies will have to be conducted to evaluate the levels of selected salivary proteins at different stages of PD. A more detailed assessment of received therapies, dental status, and motor subtype of PD may all also play a role in the proteomic profile of saliva. We believe that a proteomic and bioinformatic approach highlights the potential of salivary diagnostics in understanding the pathophysiology of PD and opens the way for future biomarker research.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/brainsci11050661/s1; Table S1: The protein clusters with fold changes >1.5 or <−1.5 in Parkinson patients (PD) versus healthy control (HC) saliva samples.

Author Contributions

Conceptualization, M.F., A.F., and D.K.; methodology, M.F., E.S., and B.Ś.; software, E.S. and B.Ś.; validation, M.F., E.S., and B.Ś.; formal analysis, E.S. and B.Ś.; investigation, M.F., Ł.M., and S.S.; resources, M.F., S.S., E.S., and B.Ś.; data curation, M.F., E.S., and B.Ś.; writing—original draft preparation, M.F., E.S., B.Ś., and Ł.M.; writing—review and editing, M.F., A.F., E.S., B.Ś., and D.K.; visualization, M.F. and B.Ś.; supervision, A.F., D.K., E.S., and M.F.; funding acquisition, M.F., A.F., and D.K. All authors have read and agreed to the published version of the manuscript.

Funding

The study was funded by a grant provided by the Medical University of Warsaw no. NZP/PM1/17.

Institutional Review Board Statement

The study was approved by the Ethics Committee of the Medical University of Warsaw (KB/239/2015, with amendment KB/25/A/2016), and therefore was performed in accordance with the ethical standards laid down in the 1964 Declaration of Helsinki and its later amendments.

Informed Consent Statement

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

Data Availability Statement

The data supporting this study are available in PRIDE using accession number PXD023489.

Acknowledgments

We would like to acknowledge Michał Kistowski for his help in statistical analysis of the data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Emamzadeh, F.N.; Surguchov, A. Parkinson’s disease: Biomarkers, treatment, and risk factors. Front. Neurosci. 2018, 12, 612. [Google Scholar] [CrossRef]
  2. Rizzo, G.; Copetti, M.; Arcuti, S.; Martino, D.; Fontana, A.; Logroscino, G. Accuracy of clinical diagnosis of Parkinson disease: A systematic review and meta-analysis. Neurology 2016, 86, 566–576. [Google Scholar] [CrossRef] [PubMed]
  3. Goldman, J.G.; Andrews, H.; Amara, A.; Naito, A.; Alcalay, R.N.; Shaw, L.M.; Taylor, P.; Xie, T.; Tuite, P.; Henchcliffe, C. Cerebrospinal fluid, plasma, and saliva in the BioFIND study: Relationships among biomarkers and Parkinson’s disease features. Mov. Disord. 2018, 33, 282–288. [Google Scholar] [CrossRef] [Green Version]
  4. Chahine, L.M.; Beach, T.G.; Seedorff, N.; Caspell-Garcia, C.; Coffey, C.S.; Brumm, M.; Adler, C.H.; Serrano, G.E.; Linder, C.; Mosovsky, S. Feasibility and Safety of Multicenter Tissue and Biofluid Sampling for α-Synuclein in Parkinson’s Disease: The Systemic Synuclein Sampling Study (S4). J. Parkinson’s Dis. 2018, 8, 517–527. [Google Scholar] [CrossRef] [Green Version]
  5. Gelpi, E.; Navarro-Otano, J.; Tolosa, E.; Gaig, C.; Compta, Y.; Rey, M.J.; Martí, M.J.; Hernández, I.; Valldeoriola, F.; Reñé, R. Multiple organ involvement by alpha-synuclein pathology in Lewy body disorders. Mov. Disord. 2014, 29, 1010–1018. [Google Scholar] [CrossRef] [PubMed]
  6. Del Tredici, K.; Hawkes, C.H.; Ghebremedhin, E.; Braak, H. Lewy pathology in the submandibular gland of individuals with incidental Lewy body disease and sporadic Parkinson’s disease. Acta Neuropathol. 2010, 119, 703–713. [Google Scholar] [CrossRef]
  7. Hawkes, C.H.; Del Tredici, K.; Braak, H. Parkinson’s disease: A dual-hit hypothesis. Neuropathol. Appl. Neurobiol. 2007, 33, 599–614. [Google Scholar] [CrossRef] [PubMed]
  8. Savica, R.; Boeve, B.F.; Mielke, M.M. When do α-synucleinopathies start? An epidemiological timeline: A review. JAMA Neurol. 2018, 75, 503–509. [Google Scholar] [CrossRef] [PubMed]
  9. Paciotti, S.; Bellomo, G.; Gatticchi, L.; Parnetti, L. Are we ready for detecting α-synuclein prone to aggregation in patients? The case of “protein-misfolding cyclic amplification” and “real-time quaking-induced conversion” as diagnostic tools. Front. Neurol. 2018, 9, 415. [Google Scholar] [CrossRef] [Green Version]
  10. Masters, J.M.; Noyce, A.J.; Warner, T.T.; Giovannoni, G.; Proctor, G.B. Elevated salivary protein in Parkinson’s disease and salivary DJ-1 as a potential marker of disease severity. Parkinsonism Relat. Disord. 2015, 21, 1251–1255. [Google Scholar] [CrossRef] [Green Version]
  11. Cersosimo, M.G.; Raina, G.B.; Calandra, C.R.; Pellene, A.; Gutiérrez, C.; Micheli, F.E.; Benarroch, E.E. Dry mouth: An overlooked autonomic symptom of Parkinson’s disease. J. Parkinson’s Dis. 2011, 1, 169–173. [Google Scholar] [CrossRef] [PubMed]
  12. Figura, M.; Friedman, A. In search of Parkinson’s disease biomarkers—is the answer in our mouths? A systematic review of the literature on salivary biomarkers of Parkinson’s disease. Neurol. I Neurochir. Pol. 2020, 54, 14–20. [Google Scholar] [CrossRef] [PubMed]
  13. Vivacqua, G.; Latorre, A.; Suppa, A.; Nardi, M.; Pietracupa, S.; Mancinelli, R.; Fabbrini, G.; Colosimo, C.; Gaudio, E.; Berardelli, A. Abnormal salivary total and oligomeric alpha-synuclein in Parkinson’s disease. PLoS ONE 2016, 11, e0151156. [Google Scholar]
  14. Kang, W.; Chen, W.; Yang, Q.; Zhang, L.; Zhang, L.; Wang, X.; Dong, F.; Zhao, Y.; Chen, S.; Quinn, T.J. Salivary total α-synuclein, oligomeric α-synuclein and SNCA variants in Parkinson’s disease patients. Sci. Rep. 2016, 6, 28143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Plotegher, N.; Gratton, E.; Bubacco, L. Number and Brightness analysis of alpha-synuclein oligomerization and the associated mitochondrial morphology alterations in live cells. Biochim. Biophys. Acta 2014, 1840, 2014–2024. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Sharon, R.; Bar-Joseph, I.; Frosch, M.P.; Walsh, D.M.; Hamilton, J.A.; Selkoe, D.J. The formation of highly soluble oligomers of alpha-synuclein is regulated by fatty acids and enhanced in Parkinson’s disease. Neuron 2003, 37, 583–595. [Google Scholar] [CrossRef] [Green Version]
  17. Bengoa-Vergniory, N.; Roberts, R.F.; Wade-Martins, R.; Alegre-Abarrategui, J. Alpha-synuclein oligomers: A new hope. Acta Neuropathol 2017, 134, 819–838. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Khurshid, Z.; Moin, S.F.; Khan, R.S.; Agwan, M.A.S.; Alwadaani, A.H.; Zafar, M.S. Human salivary protein extraction from RNAPro· SAL™, Pure· SAL™, and passive drooling method. Eur. J. Dent. 2017, 11, 385–389. [Google Scholar]
  19. Iorgulescu, G. Saliva between normal and pathological. Important factors in determining systemic and oral health. J. Med. Life 2009, 2, 303. [Google Scholar]
  20. Beach, T.G.; Adler, C.H.; Serrano, G.; Sue, L.I.; Walker, D.; Dugger, B.N.; Shill, H.A.; Driver-Dunckley, E.; Caviness, J.N.; Intorcia, A. Prevalence of submandibular gland synucleinopathy in Parkinson’s disease, dementia with Lewy bodies and other Lewy body disorders. J. Parkinson’s Dis. 2016, 6, 153–163. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Beach, T.G.; Adler, C.H.; Dugger, B.N.; Serrano, G.; Hidalgo, J.; Henry-Watson, J.; Shill, H.A.; Sue, L.I.; Sabbagh, M.N.; Akiyama, H. Submandibular gland biopsy for the diagnosis of Parkinson disease. J. Neuropathol. Exp. Neurol. 2013, 72, 130–136. [Google Scholar] [CrossRef] [Green Version]
  22. Adler, C.H.; Dugger, B.N.; Hentz, J.G.; Hinni, M.L.; Lott, D.G.; Driver-Dunckley, E.; Mehta, S.; Serrano, G.; Sue, L.I.; Duffy, A. Peripheral Synucleinopathy in Early Parkinson’s Disease: Submandibular Gland Needle Biopsy Findings. Mov. Disord. 2016, 31, 250–256. [Google Scholar] [CrossRef] [Green Version]
  23. Adler, C.H.; Dugger, B.N.; Hinni, M.L.; Lott, D.G.; Driver-Dunckley, E.; Hidalgo, J.; Henry-Watson, J.; Serrano, G.; Sue, L.I.; Nagel, T. Submandibular gland needle biopsy for the diagnosis of Parkinson disease. Neurology 2014, 82, 858–864. [Google Scholar] [CrossRef] [Green Version]
  24. Castagnola, M.; Scarano, E.; Passali, G.C.; Messana, I.; Cabras, T.; Iavarone, F.; Di Cintio, G.; Fiorita, A.; De Corso, E.; Paludetti, G. Salivary biomarkers and proteomics: Future diagnostic and clinical utilities. Acta Otorhinolaryngol. Ital. 2017, 37, 94–101. [Google Scholar] [CrossRef]
  25. Pappa, E.; Kousvelari, E.; Vastardis, H. Saliva in the “Omics” era: A promising tool in paediatrics. Oral Dis 2019, 25, 16–25. [Google Scholar] [CrossRef] [Green Version]
  26. Batista, T.B.D.; Chaiben, C.L.; Penteado, C.A.S.; Nascimento, J.M.C.; Ventura, T.M.O.; Dionizio, A.; Rosa, E.A.R.; Buzalaf, M.A.R.; Azevedo-Alanis, L.R. Salivary proteome characterization of alcohol and tobacco dependents. Drug Alcohol. Depend. 2019, 204, 107510. [Google Scholar] [CrossRef]
  27. Tierney, C.; Bazou, D.; Lê, G.; Dowling, P.; O’Gorman, P. Saliva-omics in plasma cell disorders- Proof of concept and potential as a non-invasive tool for monitoring disease burden. J. Proteom. 2021, 231, 104015. [Google Scholar] [CrossRef]
  28. Chiang, S.H.; Thomas, G.A.; Liao, W.; Grogan, T.; Buck, R.L.; Fuentes, L.; Yakob, M.; Laughlin, M.J.; Schafer, C.; Nazmul-Hossain, A. RNAPro• SAL: A device for rapid and standardized collection of saliva RNA and proteins. Biotechniques 2015, 58, 69–76. [Google Scholar] [CrossRef]
  29. Malinowska, A.; Kistowski, M.; Bakun, M.; Rubel, T.; Tkaczyk, M.; Mierzejewska, J.; Dadlez, M. Diffprot—software for non-parametric statistical analysis of differential proteomics data. J. Proteom. 2012, 75, 4062–4073. [Google Scholar] [CrossRef]
  30. Elias, J.E.; Gygi, S.P. Target-decoy search strategy for increased confidence in large-scale protein identifications by mass spectrometry. Nat. Methods 2007, 4, 207–214. [Google Scholar] [CrossRef]
  31. Elias, J.E.; Gygi, S.P. Target-decoy search strategy for mass spectrometry-based proteomics. In Proteome Bioinformatics; Springer: Berlin/Heidelberg, Germany, 2010; pp. 55–71. [Google Scholar]
  32. Rubel, T. Available online: http://proteom.ibb.waw.pl/mscan/ (accessed on 14 May 2021).
  33. Bakun, M.; Karczmarski, J.; Poznanski, J.; Rubel, T.; Rozga, M.; Malinowska, A.; Sands, D.; Hennig, E.; Oledzki, J.; Ostrowski, J. An integrated LC-ESI-MS platform for quantitation of serum peptide ladders. Application for colon carcinoma study. Proteom. Clin. Appl. 2009, 3, 932–946. [Google Scholar] [CrossRef]
  34. Youden, W.J. Index for rating diagnostic tests. Cancer 1950, 3, 32–35. [Google Scholar] [CrossRef]
  35. Kumari, S.; Goyal, V.; Kumaran, S.S.; Dwivedi, S.; Srivastava, A.; Jagannathan, N. Quantitative metabolomics of saliva using proton nmr spectroscopy in patients with parkinson’s disease and healthy controls. Neurol. Sci. 2020, 41, 1–10. [Google Scholar] [CrossRef]
  36. Masters, J.; Noyce, A.; Lynham, S.; Warner, T.; Giovannoni, G.; Proctor, G. A comprehensive characterisation of the salivary proteome of patients with Parkinson’s disease. In Proceedings of the 20th International Congress of Parkinson’s Disease and Movement Disorders, Berlin, Germany, 19–23 June 2016. [Google Scholar]
  37. Serrano, A.; Apolloni, S.; Rossi, S.; Lattante, S.; Sabatelli, M.; Peric, M.; Andjus, P.; Michetti, F.; Carrì, M.T.; Cozzolino, M. The S100A4 transcriptional inhibitor niclosamide reduces pro-inflammatory and migratory phenotypes of microglia: Implications for amyotrophic lateral sclerosis. Cells 2019, 8, 1261. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Qiu, X.; Ping, S.; Kyle, M.; Longo, J.; Chin, L.; Zhao, L.-R. S100 Calcium-Binding Protein A9 Knockout Contributes to Neuroprotection and Functional Improvement after Traumatic Brain Injury. J. Neurotrauma 2019, 37, 950–965. [Google Scholar] [CrossRef]
  39. Kaca-Oryńska, M.; Tomasiuk, R.; Friedman, A. Neuron-specific enolase and S 100B protein as predictors of outcome in ischaemic stroke. Neurol. I Neurochir. Pol. 2010, 44, 459–463. [Google Scholar] [CrossRef]
  40. Fardell, C.; Zettergren, A.; Ran, C.; Carmine Belin, A.; Ekman, A.; Sydow, O.; Bäckman, L.; Holmberg, B.; Dizdar, N.; Söderkvist, P. S100B polymorphisms are associated with age of onset of Parkinson’s disease. BMC Med. Genet. 2018, 19. [Google Scholar] [CrossRef] [Green Version]
  41. Papuć, E.; Rejdak, K. Increased Cerebrospinal Fluid S100B and NSE Reflect Neuronal and Glial Damage in Parkinson’s Disease. Front. Aging Neurosci. 2020, 12, 156. [Google Scholar] [CrossRef]
  42. Marenholz, I.; Heizmann, C.W. S100A16, a ubiquitously expressed EF-hand protein which is up-regulated in tumors. Biochem. Biophys. Res. Commun. 2004, 313, 237–244. [Google Scholar] [CrossRef]
  43. Liu, Y.; Zhang, R.; Xin, J.; Sun, Y.; Li, J.; Wei, D.; Zhao, A.Z. Identification of S100A16 as a novel adipogenesis promoting factor in 3T3-L1 cells. Endocrinology 2011, 152, 903–911. [Google Scholar] [CrossRef]
  44. Zhang, R.; Kan, J.B.; Lu, S.; Tong, P.; Yang, J.; Xi, L.; Liang, X.; Su, D.; Li, D.; Liu, Y. S100A16-induced adipogenesis is associated with up-regulation of 11 β-hydroxysteroid dehydrogenase type 1 (11β-HSD1). Biosci. Rep. 2019, 39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Sheard, J.M.; Ash, S.; Mellick, G.D.; Silburn, P.A.; Kerr, G.K. Malnutrition in a sample of community-dwelling people with Parkinson’s disease. PLoS ONE 2013, 8, e53290. [Google Scholar]
  46. Sheard, J.M.; Ash, S.; Silburn, P.A.; Kerr, G.K. Prevalence of malnutrition in Parkinson’s disease: A systematic review. Nutr. Rev. 2011, 69, 520–532. [Google Scholar] [CrossRef] [PubMed]
  47. Schrank, B.R.; Aparicio, T.; Li, Y.; Chang, W.; Chait, B.T.; Gundersen, G.G.; Gottesman, M.E.; Gautier, J. Nuclear ARP2/3 drives DNA break clustering for homology-directed repair. Nature 2018, 559, 61–66. [Google Scholar] [CrossRef] [PubMed]
  48. Vasquez, V.; Mitra, J.; Hegde, P.M.; Pandey, A.; Sengupta, S.; Mitra, S.; Rao, K.; Hegde, M.L. Chromatin-bound oxidized α-synuclein causes strand breaks in neuronal genomes in in vitro models of Parkinson’s disease. J. Alzheimer’s Dis. 2017, 60, S133–S150. [Google Scholar] [CrossRef]
  49. Gonzalez-Hunt, C.P.; Sanders, L.H. DNA damage and repair in Parkinson’s disease: Recent advances and new opportunities. J. Neurosci. Res. 2020, 99, 180–189. [Google Scholar] [CrossRef]
  50. Baietti, M.F.; Zhang, Z.; Mortier, E.; Melchior, A.; Degeest, G.; Geeraerts, A.; Ivarsson, Y.; Depoortere, F.; Coomans, C.; Vermeiren, E. Syndecan–syntenin–ALIX regulates the biogenesis of exosomes. Nat. Cell Biol. 2012, 14, 677–685. [Google Scholar] [CrossRef]
  51. Grey, M.; Dunning, C.J.; Gaspar, R.; Grey, C.; Brundin, P.; Sparr, E.; Linse, S. Acceleration of α-synuclein aggregation by exosomes. J. Biol. Chem. 2015, 290, 2969–2982. [Google Scholar] [CrossRef] [Green Version]
  52. Yu, H.; Sun, T.; An, J.; Wen, L.; Liu, F.; Bu, Z.; Cui, Y.; Feng, J. Potential Roles of Exosomes in Parkinson’s Disease: From Pathogenesis, Diagnosis, and Treatment to Prognosis. Front. Cell Dev. Biol. 2020, 8. [Google Scholar] [CrossRef] [Green Version]
  53. Rani, K.; Mukherjee, R.; Singh, E.; Kumar, S.; Sharma, V.; Vishwakarma, P.; Bharti, P.S.; Nikolajeff, F.; Dinda, A.K.; Goyal, V. Neuronal exosomes in saliva of Parkinson’s disease patients: A pilot study. Parkinsonism Relat. Disord. 2019, 67, 21–23. [Google Scholar] [CrossRef]
  54. Fung, H.-C.; Scholz, S.; Matarin, M.; Simón-Sánchez, J.; Hernandez, D.; Britton, A.; Gibbs, J.R.; Langefeld, C.; Stiegert, M.L.; Schymick, J. Genome-wide genotyping in Parkinson’s disease and neurologically normal controls: First stage analysis and public release of data. Lancet Neurol. 2006, 5, 911–916. [Google Scholar] [CrossRef]
  55. Nalls, M.A.; Blauwendraat, C.; Vallerga, C.L.; Heilbron, K.; Bandres-Ciga, S.; Chang, D.; Tan, M.; Kia, D.A.; Noyce, A.J.; Xue, A. Identification of novel risk loci, causal insights, and heritable risk for Parkinson’s disease: A meta-analysis of genome-wide association studies. Lancet Neurol. 2019, 18, 1091–1102. [Google Scholar] [CrossRef]
  56. Nalls, M.A.; Pankratz, N.; Lill, C.M.; Do, C.B.; Hernandez, D.G.; Saad, M.; DeStefano, A.L.; Kara, E.; Bras, J.; Sharma, M. Large-scale meta-analysis of genome-wide association data identifies six new risk loci for Parkinson’s disease. Nat. Genet. 2014, 46, 989–993. [Google Scholar] [CrossRef] [PubMed]
  57. Lesage, S.; Drouet, V.; Majounie, E.; Deramecourt, V.; Jacoupy, M.; Nicolas, A.; Cormier-Dequaire, F.; Hassoun, S.M.; Pujol, C.; Ciura, S. Loss of VPS13C function in autosomal-recessive Parkinsonism causes mitochondrial dysfunction and increases PINK1/Parkin-dependent mitophagy. Am. J. Hum. Genet. 2016, 98, 500–513. [Google Scholar] [CrossRef] [Green Version]
  58. Van Stiphout, M.A.; Marinus, J.; van Hilten, J.J.; Lobbezoo, F.; de Baat, C. Oral Health of Parkinson’s Disease Patients: A Case-Control Study. Parkinson’s Dis. 2018, 2018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Youdim, M.; Grünblatt, E.; Mandel, S. The Pivotal Role of Iron in NF-κB Activation and Nigrostriatal Dopaminergic Neurodegeneration: Prospects for Neuroprotection in Parkinson’s Disease with Iron Chelators. Ann. N. Y. Acad. Sci. 1999, 890, 7–25. [Google Scholar] [CrossRef] [PubMed]
  60. Dehmer, T.; Heneka, M.T.; Sastre, M.; Dichgans, J.; Schulz, J.B. Protection by pioglitazone in the MPTP model of Parkinson’s disease correlates with IκBα induction and block of NFκB and iNOS activation. J. Neurochem. 2004, 88, 494–501. [Google Scholar] [CrossRef] [PubMed]
  61. Kim, B.-W.; Koppula, S.; Kumar, H.; Park, J.-Y.; Kim, I.-W.; More, S.V.; Kim, I.-S.; Han, S.-D.; Kim, S.-K.; Yoon, S.-H. α-Asarone attenuates microglia-mediated neuroinflammation by inhibiting NF kappa B activation and mitigates MPTP-induced behavioral deficits in a mouse model of Parkinson’s disease. Neuropharmacology 2015, 97, 46–57. [Google Scholar] [CrossRef]
  62. Xu, W.; Zheng, D.; Liu, Y.; Li, J.; Yang, L.; Shang, X. Glaucocalyxin B alleviates lipopolysaccharide-induced Parkinson’s disease by inhibiting TLR/NF-κB and activating Nrf2/HO-1 pathway. Cell. Physiol. Biochem. 2017, 44, 2091–2104. [Google Scholar] [CrossRef] [Green Version]
  63. Soós, J.; Engelhardt, J.I.; Siklós, L.; Havas, L.; Majtényi, K. The expression of PARP, NF-κB and parvalbumin is increased in Parkinson disease. Neuroreport 2004, 15, 1715–1718. [Google Scholar] [CrossRef]
  64. Chen, G.; Nie, S.; Han, C.; Ma, K.; Xu, Y.; Zhang, Z.; Papa, S.; Cao, X. Antidyskinetic effects of MEK inhibitor are associated with multiple neurochemical alterations in the striatum of hemiparkinsonian rats. Front Neurosci 2017, 11, 112. [Google Scholar] [CrossRef] [Green Version]
  65. Chen, W.-F.; Wu, L.; Du, Z.-R.; Chen, L.; Xu, A.-L.; Chen, X.-H.; Teng, J.-J.; Wong, M.-S. Neuroprotective properties of icariin in MPTP-induced mouse model of Parkinson’s disease: Involvement of PI3K/Akt and MEK/ERK signaling pathways. Phytomedicine 2017, 25, 93–99. [Google Scholar] [CrossRef]
  66. Lavin, P.J.; Gawel, M.J.; Das, P.K.; Alaghband-Zadeh, J.; Rose, F.C. Effect of Levodopa on Thyroid Function and Prolactin Release: A Study in Patients With Parkinson’s Disease. Arch. Neurol. 1981, 38, 759–760. [Google Scholar] [CrossRef] [PubMed]
  67. Eisler, T.; Thorner, M.; MacLeod, R.; Kaiser, D.; Calne, D. Prolactin secretion in Parkinson disease. Neurology 1981, 31, 1356. [Google Scholar] [CrossRef] [PubMed]
  68. Lawton, N.; MacDermot, J. Abnormal regulation of prolactin release in idiopathic Parkinson’s disease. J. Neurol. Neurosurg. Psychiatry 1980, 43, 1012–1015. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Tumilasci, O.R.; Cersosimo, M.; Belforte, J.E.; Micheli, F.E.; Benarroch, E.E.; Pazo, J.H. Quantitative study of salivary secretion in Parkinson’s disease. Mov. Disord. 2006, 21, 660–667. [Google Scholar] [CrossRef]
  70. Tiigimäe-Saar, J.; Tamme, T.; Rosenthal, M.; Kadastik-Eerme, L.; Taba, P. Saliva changes in Parkinson’s disease patients after injection of Botulinum neurotoxin type A. Neurol. Sci. 2018, 39, 871–877. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Brainsci 11 00661 g001 550
Figure 1. Area under the curve (AUC) receiver operating characteristic curve (ROC) analysis for S100A16 and ARPC1A proteins in HC versus PD patients.
Figure 1. Area under the curve (AUC) receiver operating characteristic curve (ROC) analysis for S100A16 and ARPC1A proteins in HC versus PD patients.
Brainsci 11 00661 g001
Brainsci 11 00661 g002 550
Figure 2. A volcano plot from the mass spectrometry data is shown. The x-axis shows the fold change values and the y axis shows the −log10 p-values showing statistical significance. Q96FQ6–S100A16 protein; Q92747–ARPC1A protein; O75351–VPS4B protein. The horizontal dashed line shows where p = 0.05 (−log10(0.05)~ = 1.3) is, while the vertical dashed line shows where the fold change is 1.5 (log2 (1.5)~ = 0.58). The absolute 1.5-fold change and p-value of 0.05 are used as the threshold cutoff.
Figure 2. A volcano plot from the mass spectrometry data is shown. The x-axis shows the fold change values and the y axis shows the −log10 p-values showing statistical significance. Q96FQ6–S100A16 protein; Q92747–ARPC1A protein; O75351–VPS4B protein. The horizontal dashed line shows where p = 0.05 (−log10(0.05)~ = 1.3) is, while the vertical dashed line shows where the fold change is 1.5 (log2 (1.5)~ = 0.58). The absolute 1.5-fold change and p-value of 0.05 are used as the threshold cutoff.
Brainsci 11 00661 g002
Table
Table 1. Clinical characteristic of study groups. LED: levodopa equivalent dose; UPDRS: Unified Parkinson’s Disease Rating Scale; ns: non- significant; na: non applicable.
Table 1. Clinical characteristic of study groups. LED: levodopa equivalent dose; UPDRS: Unified Parkinson’s Disease Rating Scale; ns: non- significant; na: non applicable.
PD GroupHealthy Controlp-Value
Number of patients2415
Sex M/F15/99/5ns
Age61.6 ± 8.260.9 ± 6.7ns
Levodopa dose1352 ± 763na
LED1524 ± 786na
UPDRS part III OFF34.6 ± 12.4na
Table
Table 2. Upstream regulators identified using information from the Ingenuity Pathway Analysis (IPA) software (Qiagen, Germany).
Table 2. Upstream regulators identified using information from the Ingenuity Pathway Analysis (IPA) software (Qiagen, Germany).
Upstream RegulatorMolecule TypePredicted Activation StateActivation z-Scorep-Value of OverlapTarget Molecules in Dataset
PD98059chemical - kinase inhibitorActivated2.7379.13 × 10−5ANXA1⇓,CAMP⇓, FN1⇓, GRN⇓,ITGAM⇓, KRT19⇓, NQO1⇓, RETN⇓, S100A4⇓, S100A8⇓,SOD1⇓
SRFtranscription regulatorActivated2.6462.72 × 10−4CAMP⇓, ELANE⇓, ITGAM⇓, LTF⇓, MYH9⇓, S100A8⇓, S100A9⇓, SLPI⇓
SB203580chemical - kinase inhibitorActivated2.5968.58 × 10−5ANXA1⇓, ANXA3⇓, CAMP⇓, FN1⇓, HMGB1⇑, ITGAM⇓, NQO1⇓, RETN⇓, SLPI⇓, TGM1⇓
RICTORotherActivated2.2361.11 × 10−2PSMA1⇓, PSMA5⇓, PSMB1⇓, PSMB3⇓, RPSA⇓
GATA3transcription regulatorActivated2.2161.22 × 10−3FN1⇓, PPL⇓, S100A8⇓, S100A9⇓, SLPI⇓, TGM1⇓
mifepristonechemical drugActivated2.1211.75 × 10−4ANXA1⇓, CAST⇓,GSTM1⇓, HSPD1⇓, ITGAM⇓, JUP⇓, LTF⇓, PSMA1⇓
tert-butyl-hydroquinonechemical reagentInhibited−2.1756.90 × 10−5GSS⇓, ITGAM⇓, ME1⇓, NQO1⇓, PSMA5⇓
arsenic trioxidechemical drugInhibited−2.1761.08 × 10−3FN1⇓, ITGAM⇓, ME1⇓, NQO1⇓, PDIA4⇓, S100A8⇓, VCP⇓
IL17AcytokineInhibited−2.1968.72 × 10−3CAMP⇓, MPO⇓, S100A12⇓, S100A8⇓, S100A9⇓
SP1transcription regulatorInhibited−2.2001.91 × 10−5CAMP⇓, CES1⇓, ELANE⇓, FN1⇓, GSS⇓, ITGAM⇓, KRT16⇓, KRT19⇓, KRT4⇓, PADI4⇑, SOD1⇓, TGM1⇓
IL2cytokineInhibited−2.2133.35 × 10−2ANXA1⇓,FCGR3A/FCGR3B⇑,GRN⇓, HSPD1⇓, PSMB1⇓, S100A4⇓, S100A8⇓
butyric acidchemical - endogenous mammalianInhibited−2.2272.81 × 10−8ALDH1A1⇓, ANXA1⇓, ANXA5⇓, CAMP⇓, CEACAM5⇓, ELANE⇓, GRN⇓, HMGB1⇑, ITGAM⇓, KRT13⇓, MVP⇓, PCMT1⇑, TGM1⇓, TP53I3⇓
curcuminchemical drugInhibited−2.2367.96 × 10−3CAMP⇓, HMGB1⇑, ITGAM⇓, NQO1⇓, RETN⇓, SOD1⇓
OSMcytokineInhibited−2.2725.08 × 10−7ANXA1⇓, ANXA3⇓, CAMP⇓, CDA⇓, FN1⇓, GCA⇓, KLK13⇑, KRT16⇓, KRT19⇓, LRRFIP1⇓, S100A12⇓, S100A8⇓, S100A9⇓, SLPI⇓
HSF1transcription regulatorInhibited−2.3606.56 × 10−4CCT2⇓, CCT3⇑, FKBP4⇓, HMGB1⇑, HSPD1⇓, TCP1⇓
1,2-dithiol-3-thionechemical reagentInhibited−2.5036.27 × 10−8CCT3⇑, GSTM1⇓, NQO1⇓, PDIA4⇓, PSMA1⇓, PSMA5⇓, PSMB1⇓, PSMB3⇓, SOD1⇓, VCP⇓
EGFgrowth factorInhibited−2.5161.16 × 10−5FN1⇓, ITGAM⇓, KRT16⇓, KRT19⇓, KRT5⇓, LTF⇓, PFKM⇓, PSMB1⇓, PSMB3⇓, S100A4⇓, S100A9⇓, TGM1⇓
tetradecanoylphorbol acetatechemical drugInhibited−2.6081.43 × 10−5ANXA1⇓, AZU1⇓, CES1⇓, ELANE⇓, GSTM1⇓, ITGAM⇓, LRRFIP1⇓, MPO⇓, MYH9⇓, NQO1⇓, RETN⇓, S100A14⇓, S100A8⇓, S100A9⇓, SLPI⇓, SOD1⇓, STATH⇑, TGM1⇓
NFkB (complex)complexInhibited−2.6213.90 × 10−3ALDH7A1⇓, CAMP⇓, CAPNS1⇓, FN1⇓, HSPA9⇓, ITGAM⇓, KRT19⇓, PKP1⇓, SLPI⇓
KLF4transcription regulatorInhibited−2.7042.55 × 10−7ALDH1A1⇓, DSP⇓, FN1⇓, ITGAM⇓, KRT13⇓, KRT19⇓, PFKP⇓, PPL⇓, S100A14⇓, SLPI⇓, TGM1⇓
CEBPAtranscription regulatorInhibited−2.7551.11 × 10−5ANXA1⇓, CAMP⇓, ELANE⇓, ITGAM⇓, LTF⇓, MPO⇓, PPL⇓, RETN⇓, S100A8⇓, S100A9⇓, SOD1⇓
Lipopolysaccharidechemical drugInhibited−2.9519.57 × 10−9ANXA1⇓, ANXA3⇓, ANXA5⇓, AZU1⇓, CAMP⇓, CLIC3⇓, ELANE⇓, FN1⇓, GCA⇓, HMGB1⇑, ITGAM⇓, ITIH4⇓, KRT13⇓, KRT4⇓, LGALS3⇓, LRRFIP1⇓, LTF⇓, MPO⇓, MYH9⇓, NQO1⇓, ORM2⇓, PDIA4⇓, PFKP⇓, RETN⇓, S100A12⇓, S100A8⇓, S100A9⇓, SLPI⇓, SOD1⇓, TRIM29⇓
PRLcytokineInhibited−3.0711.07 × 10−6ANXA3⇓, ANXA5⇓, FN1⇓, GSTM1⇓, HSPD1⇓, KRT19⇓, KRT5⇓, PDIA4⇓, RPSA⇓, SOD1⇓
NFE2L2transcription regulatorInhibited−3.2632.58 × 10−8CCT3⇑, FN1⇓, GSS⇓, GSTM1⇓, HSPA9⇓, ME1⇓, NQO1⇓, PDIA4⇓, PSMA1⇓, PSMA5⇓, PSMB1⇓, PSMB3⇓, SOD1⇓, VCP⇓
beta-estradiolchemical - endogenous mammalianInhibited−3.4286.34 × 10−8ADK ⇓,ALDH7A1⇓, ANXA1⇓, ANXA3⇓, CAST⇓,CCT2⇓, DSP⇓, FN1⇓, HNRNPD⇓, HSPA9⇓, HSPD1⇓, ITIH4⇓, KRT13⇓, KRT16⇓, KRT19⇓, KRT4⇓, KRT5⇓, LGALS3⇓, LTF⇓, MPO⇓, PADI4⇑, PDIA4⇓, PPL⇓, PSMA1⇓, PSMB1⇓, S100A4⇓, S100A9⇓, SLPI⇓, TP53I3⇓, TRIM29⇓
Knowledge base to explain observed proteomic changes between the patients and the control group. The table contains regulators with the highest and lowest Z-score (cut-off at > 2-activated or < −2-inhibited) identified from the IPA software (Qiagen, Germany) with p-values overlapping < 0.05. Targets in the experimental results are labeled as ⇓ for downregulated or ⇑ for upregulated proteins between the patient and the control group. The colors of the target genes indicate the literature prediction of the influence of the upstream regulator on the given protein: red—activated; blue—inhibited; black—affected with no proven direction of the relationship.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Figura, M.; Sitkiewicz, E.; Świderska, B.; Milanowski, Ł.; Szlufik, S.; Koziorowski, D.; Friedman, A. Proteomic Profile of Saliva in Parkinson’s Disease Patients: A Proof of Concept Study. Brain Sci. 2021, 11, 661. https://doi.org/10.3390/brainsci11050661

AMA Style

Figura M, Sitkiewicz E, Świderska B, Milanowski Ł, Szlufik S, Koziorowski D, Friedman A. Proteomic Profile of Saliva in Parkinson’s Disease Patients: A Proof of Concept Study. Brain Sciences. 2021; 11(5):661. https://doi.org/10.3390/brainsci11050661

Chicago/Turabian Style

Figura, Monika, Ewa Sitkiewicz, Bianka Świderska, Łukasz Milanowski, Stanisław Szlufik, Dariusz Koziorowski, and Andrzej Friedman. 2021. "Proteomic Profile of Saliva in Parkinson’s Disease Patients: A Proof of Concept Study" Brain Sciences 11, no. 5: 661. https://doi.org/10.3390/brainsci11050661

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