The studies were performed on a set of clinical samples including 20 sialoliths, which were divided into 3 groups based on spectroscopic studies: 6 calcified stones, 4 lipid stones and 10 mixed stones (Table S1). Each stone was separately crushed into the powder using mortar. The obtained material was used to extract the proteins. The diagram of the entire process, starting from a stone sample and ending with data analysis, is shown in Fig. 2. For the qualitative analysis and, simultaneously, for the creation of the spectral library required for the quantitative SWATH analysis, the pooled clinical samples were processed following protocol 1 (see materials and methods), and the IDA spectra were recorded. In addition, a three-stage enrichment of the library was carried out—the first step by using alkaline pH chromatographic separation of trypsin-digested peptides in FASP methodology. The second type of enrichment involved the separation of proteins extracted from stones on a gel into 12 fractions in the SageELF apparatus, which performs automatic elution of each fraction. Finally, they were individually digested according to protocol 2 clinical samples. Spectra in the IDA format were recorded for each fraction from the chromatographic and gel separation and individual clinical samples. For each clinical sample, the possibility of digesting the powdered material directly on the membrane (in a microcon, 10 kDa) was checked according to the FASP methodology - protocol 2. The aim was to obtain as much protein extract as possible for digestion and the possibility of direct digestion on the sialolith pellet. Due to its structure, the powder did not block the membrane, and it was possible to perform all standard FASP steps for each sample freely. The two protocols were compared quantitatively, and the results were similar. Fragmentation spectra were analyzed using PeaksStudio (qualitative analysis) and Protein Pilot (spectral library for SWATH-MS analysis).
Trypsin-digested peptides were separated using HPLC chromatography at alkaline pH, which increased from about 200 IDs to 694 and made the most outstanding contribution to the spectral library. In gel separation combined with automatic elution in the SageELF system (see materials and methods) combined with further digestion (FASP) of the thirteen fractions obtained increased the ID number in a library to 794 human proteins.
Quantitative analysis
The quantitative analysis was carried out per the SWATH-MS approach, and the first step was constructing a spectral library using spectra from IDA experiments. The combined fraction of protein extracts from all sialoliths obtained according to protocol 1 was used to create the library.
The results of MS analysis of clinical samples performed in DIA mode were analysed in PeakView 2.2 software using constructed spectral library. The quantitative analysis allowed us to identify up-regulated and down-regulated proteins. For each analysed clinical sample, sets of proteins which were statistically significant (q < 0.05) and the log2(FC) ≥ 0.45 for up-regulated proteins and log2(FC)≤-0.45 for down-regulated proteins were chosen. In the case of calcified stones, 109 proteins were quantitatively analysed (Table S3), among which up- and down-regulated proteins detected in more than 50% of samples from this group are presented in the form of a heat map (Fig. 5A). For all samples from the LIP group, 84 statistically significant proteins were quantified in total samples from the LIP group. Most of the proteins were down-regulated for this group's samples (Table S4). The next step was selecting the set of statistically significant up-regulated and down-regulated proteins detected in more than 50% of samples from the LIP group. Heatmap for these 22 proteins presents the level of regulation of proteins (Fig. 5B).
For all samples from the MIX group, 114 statistically significant proteins were quantified. In this case, the type of regulation of proteins is more varied in this group (Table S5). The set of statistically significant up-regulated and down-regulated proteins detected in more than 50% of samples from the MIX group composed of 26 proteins is presented as a heatmap for this set presents the level of regulation of proteins (Fig. 5C).
DISCUSSION
Although salivary stone disease is a common pathological state, it can cause serious consequences, such as severe pain, discomfort, nerve injuries and the necessity of re-surgical operation to remove sialoliths. Despite many theories, there is no confirmed cause of pathogenesis leading to salivary stone formation. Thanks to the use of qualitative, quantitative and bioinformatic analysis of MS data, we were able to examine protein compositions of sialoliths collected from patients taking into account their classification to calcified, lipid and mixed groups based on spectroscopic methods. Using advanced software for proteomic analysis, we had the chance to identify human and bacterial proteins.
Bacterial proteins.
MS data registered in IDA mode were processed using PeaksSTUDIO software, and thanks to the Multi-round search option, raw spectra were searched first against the human and next against the bacterial database. Detected bacterial proteins belonged to 44 bacteria species. The most numerous bacteria group for which the proteins were identified was Actinomyces, which belongs to gram-positive opportunistic pathogens common in the oral cavity, especially in the gums, causing various oral infections [75]. Many Actinomyces species can cause actinomycosis associated with swelling and the formation of abscesses [76]. Actinomyces viscosus is a pathobiont which can colonize the oral cavity of even 70% of adults [77]. Its presence is connected with periodontal disease – A. viscosus was isolated from root surface caries and dental calculus [78]. A. naeslundii can also cause periodontal disease and is one of the first bacteria which can colonize the oral cavity and then cover the surface of teeth [79], [80]. A. radicidentis was also identified in infected root canals of teeth [81]. It is also part of the biofilm found on oral surfaces [82]. Capnocytophaga sp. oral is a gram-negative opportunistic pathogen involved in pathogenesis leading to periodontal disease [83]. It can often be isolated from periodontal pockets and abscesses [84]. The following identified bacteria was Eikenella corrodens, a Gram-negative bacteria common in the oral cavity but can act as an opportunistic pathogen [85]. Its presence can cause the formation of abscesses, including the area of submandibular glands [86]. Another large group of detected bacteria is Fusobacterium. According to the current statement, gram-negative Fusobacterium species should be permanently treated as pathogens – they cause several human diseases, including periodontal disease [75], [87]. The species responsible for that pathological state are, for example, F. nucleatum and F. polymorphum. Moreover, they can form aggregates with other bacteria in the oral cavity [88], [89]. Usually, the Gram-negative Hamephilus species are commensal organisms, including the mouth. They are also part of the salivary microbiome [90], [91]. However, detected H. haemolyticus and H. parainfluenzae are opportunistic pathogens whose presence leads to the formation of abscesses [92]. 4 species from Neisseria group were detected in sialoliths: N. bacilliformis, N. macacae, N. sicca and N. sp. oral. These Gram-negative bacteria can be found on mucosal surfaces.
In most cases, they are commensals but living in the oral cavity, and they can act as opportunistic pathogens, causing some diseases and infections [93], [94]. Porphyromonas sp. oral is gram-negative bacteria commonly present in the oral cavity and found in the salivary microbiome [95], [96]. As pathobiont in the case of disturbed homeostasis, Porphyromonas can cause different diseases, for example, periodontitis [97], [98]. Pseudopropionibacterium propionicum, a Gram-positive bacteria and opportunistic pathogen, form biofilm on external root surfaces of teeth, suggesting P. propionicum as the cause of endodontic pathological states. Its presence is the leading cause of actinomycosis [99]–[101]. Gram-positive Rothia dentocariosa is a standard part of the oral and respiratory tract microbiome [102]. This bacteria was associated with periodontal disease, and according to one of the hypotheses, this pathological state caused by R. dentocariosa can lead to infections of other tissues [103]. They identified proteins from Gram-negative Tannerella forsythia and Treponema denticola in clinical sialolith samples and Porphyromonas gingivalis (not detected in this research but shown in other) from the Red Complex. They are the primary virulent pathogens which cause chronic periodontitis [104]. Other species of bacteria were also detected based on the analyzed proteins, usually belonging to the natural microbiome of the oral cavity. However, they transformed into pathogens under stressful conditions and were allowed to cause variable diseases. They were primarily Gram-negative species: Aggregatibacter aphrophilus, Desulfobulbus oralis, Fretibacterium sp., Kingella potus, Ottowia sp. oral, Selenomonas noxia and Selenomonas sp. oral. There were also Gram-positive species: Peptostreptococcus stomatis and Streptococcus mitis. As we can see, bacteria are commonly present in the oral cavity (gums, saliva, surface of teeth), and often, they are causes of periodontal diseases. It can also suggest their potential influence on pathogenesis leading to the development of salivary stone disease.
Some of the identified bacteria species, for which the proteins were extracted and detected from clinical samples of salivary stones, were also shown in previous research describing the presence of bacterial pathogens in sialoliths and their potential influence on pathogenesis and mineralisation leading to the formation of deposits in salivary glands. There were detected precisely the same species: Actinomyces viscosus, Eikenella corrodens, Fusobacterium nucleatum, Haemophilus parainfluenzae and Streptococcus mitis. What is more, identified bacteria, if they were not the same species, they were from the same groups: Actinomyces, Capnocytophaga, Eikenella, Haemophylus, Kingella, Neisseria, Peptostreptococcus, Porphyromonas, Rothia and Streptococcus [22], [45], [49]–[51]. All of the identified bacteria groups are common in the oral cavity. In most cases, they are opportunistic pathogens which can cause various infections and diseases when the environmental homeostasis is disturbed. There were also selected bacteria species that occur uniquely in different types of sialoliths and are standard for all 3 salivary stones (Table 3). These sets are numerous, taking into account limited samples in each group. A similar analysis should be conducted using more clinical samples to confirm the uniqueness and repeatability of these species. Besides, the protein extraction protocol should be improved to obtain a more significant number of proteins, especially bacterial proteins.
Human proteins
Qualitative analysis of identified human proteins in the studied groups of sialoliths, supported by enrichment and functional analysis, showed that the most often repeated proteins between groups with the highest frequency are- neutrophil defensin 3 (DEFA3), protein S100A9 and S100-A8, cathepsin G (CTSG), myeloperoxidase (MPO), lactotransferrin (LTF), eosinophil cationic protein (RNASE3) and lysozyme C (LYZ) - according to the GO enrichment mainly involved in the defence response to the bacterium, regulated exocytosis and neutrophil degranulation. Most of the proteins identified and shown in Fig. 3 are multifunctional, suggesting a complex process influenced by many factors. Neutrophil defensin 3 has antibiotic, fungicide and antiviral activities. The group of neutrophil defensins can kill microorganisms by permeabilizing their plasma membrane [105] and are present in the granules of neutrophils and the epithelia of mucosal surfaces such as the oral cavity [106]. The decreased level of neutrophil defensin 3 was noted in the case of dental caries. It is connected with antimicrobial activity because, in the pathogenesis of periodontitis, antimicrobial proteins such as DEFA3 can cooperate with other inflammatory proteins and regulate distinct inflammatory pathways [107], [108]. Neutrophil defensin 3 is also associated with dyslipidemia causing lipids imbalance and influencing the lipid and calcium balance in sialolithiasis [109]. On the other hand, proteins S100A8 and S100A9 are important calcium- and zinc-binding proteins and form a complex called calprotectin [110] which has antimicrobial properties thanks to the metal sequestration process conducted in the presence of calcium by chelation [111]–[113]. A low level of S100-A9 is the cause of neutrophil extracellular trap formation (NETs) in the presence of bacteria [114]. The formation of these structures directly influences calcification, leading to sialolith generation. On the other hand, increased level of proteins S100-A8 and S100-A9 is treated as a marker of periodontitis [115] indicating that an imbalance of these 2 proteins can cause oral pathological states. Cathepsin G (CTSG) is a protein which also plays antibacterial activity. Moreover, this property can also allow fighting against biofilms, which can play an important role in sialolith formation [116]–[118]. Cathepsin G belongs to the neutrophil serine proteases family. This protein was first identified as a degradative enzyme that acts at inflammatory sites in 2 ways - intracellularly (degradation of pathogens) and extracellularly (breakdown of extracellular matrix components) [119]. CTSG is also localized in NETs – important for forming sialoliths - because of its high affinity for chromatin [120]. This protein causes changes in the level of calcium ions during exposure to endothelial cells [121]. On the other hand, the cathepsin family was also detected as a mediator of the metabolism of lipids [122], so an imbalance of lipids and calcium can also be caused by the level disturbance of cathepsin G. Myeloperoxidase is a peroxidase enzyme expressed mostly in neutrophil granulocytes. Their antimicrobial activity is connected with the secretion of hypohalous acids [123], [124]. Increased level of MPO was noted in saliva in patients with periodontitis [125]. Focusing on the balance of calcium and lipid, first, myeloperoxidase has a high affinity for calcium [126] and second, during the inflammation process in the presence of pathogenic bacteria, MPO is responsible for initiating the initiation of lipid peroxidation, changing their features [127]. Myeloperoxidase can also be found in NETs [128]. Lactotransferrin (LTF) is a multifunctional protein of the transferrin family and has antimicrobial activity. This feature depends on the extracellular cation concentration [129] possibly connected with the calcium ions' level. LTF is mostly present in secretory fluids, such as saliva [130]. Lactoferrin in saliva decreases bacterial growth, biofilm development and inflammatory processes [131], [132]. LTS is also a biomarker of salivary gland pathological states [133]. This protein also takes part in the creation of NETs [128]. Eosinophil cationic protein (RNASE3) is a heparin-binding ribonuclease with cytotoxic abilities [134]. This feature is used against many pathogens – lipid bilayers of pathogenic microorganisms are destabilized [135]. The concentration of eosinophil cationic protein is increased in plasma and other body fluids, also in saliva, during inflammation [136]. There are many studies about the role of eosinophil cationic protein as a biomarker of asthma [137]. Lysozyme (LYZ) is another antimicrobial enzymatic protein. It causes hydrolyzation of peptidoglycan - the crucial component of the cell walls of Gram-positive bacteria. Thanks to that, the bacterial cells are lysed [138]. Moreover, lysozyme and calcium cause an imbalance of calcium concentration [139]. Besides, LYZ is also part of neutrophil extracellular traps [128].
Quantitative analysis.
The SWATH-MS approach was used for the quantitative analysis of salivary stones. The crucial stage of this approach was constructing a well-developed spectral library. For this purpose, a pooled sample of peptides obtained after FASP digestion of proteins extracted from salivary stones was fractionated with high-performance liquid chromatography in basic pH and SageELF methodology. Thanks to that, the number of proteins included in the spectral library was maximised. After adding the spectra of 60 HPLC fractions and 12 SageELF fractions into the final version of the spectral library, there are almost 3,5 times more proteins than after processing only IDA spectra of clinical sialolith samples. There is no optimal control sample for salivary stones to compare the level of protein regulation. During previous studies, maxillary bones and teeth were used as control [58], but the same way was not possible during our research. However, during this project stage, we decided to follow the approach successfully used during our published preliminary studies – control was a pooled sample of sialolith peptides fractions. The exact amount of peptides was taken from each clinical salivary stones sample, and after mixing and final clean-up with use, the StageTips procedure control sample was ready for analysis. To check the applicability of pooled samples as a control, principal component analysis (PCA) was conducted (Figure S2). As expected, the pooled sample is placed in the middle of the graph. Relative quantitative analysis was performed considering the original classification of sialolith samples.
The comparison of the selected sets of the most common statistically significant proteins quantified for each group (it was assumed that these proteins must be present in more than 50% of samples in each group) is presented in Figure 6 as a Venn diagram. Thanks to that, it was possible to detect the number of unique proteins for different groups.
Calcified sialoliths.
A STRING-generated network presenting experimentally verified protein-protein interactions was prepared in Fig. 7A. The network presents the type of regulation of proteins, their frequency among the samples in the group and their uniqueness by marking the proteins with yellow circles for 33 quantified proteins, all unique for the CAL group. Almost all of them were up-regulated in each of the samples from the group. Only one protein - immunoglobulin J chain (IGJ) – was down-regulated in one sample (sample 12). This protein is crucial to forming a complex of immunoglobulins M and A, which then is secreted into the mucosa, where this complex plays an initial role in the immune response [140]. Variable levels of IGJ can cause this process disorder. The biggest group of unique up-regulated proteins is the keratin family (KRT1, KRT2, KRT4, KRT5, KRT13). Keratins are structural fibrous proteins which are responsible for structural balance. They play a crucial role in the repair of epidermal barriers. During this process, their level is increased, causing changes in the lipids metabolism, and this state is possible in the case of calcified sialoliths [141]. A following protein family is a group of proteins S100, which representatives were mentioned earlier. The median values of fold change for these proteins are as follows: S100A8 = 4.30, S100A12 = 5.25, S100P = 11.76, S100A9 = 1.76 (protein S100A9 is not shown on the STRING-generated network, because it was not quantified in more than 50% of samples in CAL group – S100A9 was quantified in 3 calcified samples). These proteins can bind calcium, so their increased level can cause the domination of calcium in the structure of sialoliths. The essential thing is that the S100 family is typical in NETs [128]. Two common in salivary glands and saliva proteins from BPI fold-containing family (BPIFA2, BPIFB1), responsible for an antimicrobial activity causing changes in the cellular response to lipopolysaccharide by binding to LPS, are also up-regulated, pointing to the presence of bacteria [142], [143]. Mucin 8 (MUC8) was proposed earlier as one candidate for the biomarker of sialolithiasis [59]. MUC8 was not detected during this research, but mucin 7 (MUC7), typical for oral cavity protein, was quantified. MUC7 was up-regulated, similar to the case of MUC8 during inflammatory processes [60]. Another unique CAL group protein was stomatin (STOM). The function of this protein is regulating ion channel activity in membranes, so a higher level of STOM can influence the concentration of calcium ions [144]. Changes in calcium concentration can also be associated with the up-regulation of transketolase (TKT). This enzymatic protein can form complexes with calcium [145]. The last 2 up-regulated unique proteins – haptoglobin (HP) and fibrinogen (FGB) – both have antimicrobial activity, which is connected with the presence of different bacteria [146], [147]. Non-unique but persistent protein is also up-regulated keratin 19 (KRT19), belonging to the mentioned earlier keratin family. Another protein common for the calcified group is down-regulated neutrophil elastase (ELANE), a serine proteinase secreted by neutrophils during inflammation. This protein has a high affinity to DNA. It can be found in neutrophil extracellular traps. Down-regulation of ELANE can suggest that the influence of NETs on biocalcification is reduced [148]. The high frequency of up-regulated immunoglobulin mu heavy chain (P0DOX6) clearly shows the activity of the immune system against the present bacteria.
In order to enrich the analysis and obtain additional information on unique and classified as significant based on the SWATH analysis, enrichment analysis was performed with the use of g:Profiler bioinformatic tool in the following annotation category: Gene Ontology (Biological Process, Cellular Component, Molecular Function), KEGG, Reactome and type of regulation: up-regulation, down-regulation or varied level of regulation of proteins among the samples in the group (Fig. 8). The uniqueness of proteins for the group was marked with yellow names of proteins. Analysing Biological Process GO terms, most of them are connected with the body’s defence against bacteria, which proteins were identified earlier. About 50% of analysing proteins are involved in these processes and up-regulated for all calcified samples. Only 4 proteins: neutrophil elastase (ELANE), statherin (STATH), eosinophil cationic protein (RNASE3) and azurocidin (AZU1) are down-regulated for all of the samples. Other 4 proteins: immunoglobulin J chain (IGJ), myeloperoxidase (MPO), cathepsin G (CTSG), and immunoglobulin lambda constant 2 (P0DOY2), have varied values of fold change among the samples in the calcified group. These changes are directly connected with pathological states caused by pathogenic bacteria.
Focussing on the Cellular Component GO terms, all of them are associated with creating extracellular components, referring to the hypothesis about the role of neutrophil extracellular traps in the biocalcification leading to the sialoliths formation. In this case, only keratin 4 (KRT4) is not involved.
Most of the proteins are up-regulated for all of the samples. 7 proteins: neutrophil elastase (ELANE), statherin (STATH), submaxillary gland androgen-regulated protein 3B (SMR3B), eosinophil cationic protein (RNASE3), azurocidin (AZU1), haemoglobin subunit beta (HBB) and neutrophil gelatinase-associated lipocalin (LCN2), which are involved in the formation of cellular components, are down-regulated. Other 6 proteins: carbonic anhydrase 1 (CA1), immunoglobulin J chain (IGJ), myeloperoxidase (MPO), cathepsin G (CTSG), immunoglobulin lambda constant 2 (P0DOY2) and alpha-amylase 1A (P0DUB6) have a different level of regulation among the samples in the group. For Biological Process and Cellular Component GO terms, the values of -log(q-value) are high, pointing to the high statistical significance of the identified terms. However, it can be caused by numerous groups of analysing proteins. Detected Molecular Function GO terms are mostly connected to antimicrobial activities, but fewer proteins were linked to these terms compared to previous results. One of the most important conclusions is detecting the Neutrophil extracellular trap formation KEGG pathway based on quantified proteins. It is evidence of the crucial role of NETs in salivary stones formation, but the group of proteins involved in this process is relatively small. Several proteins also influence the Salivary secretion KEGG pathway, so modifying this process may have an essential role in pathogenesis. What is more, identified Neutrophil degranulation Reactome pathway can be indicative of the ongoing process of NETs formation because, during the degranulation of neutrophils, their chromatin is released. Most of the proteins associated with the Neutrophil degranulation Reactome pathway are up-regulated, so we suppose this process intensifies.
Lipid sialoliths.
For the chosen set of the most common, statistically significant proteins quantified in more than 50% of lipid salivary stones samples, the STRING-generated network presenting experimentally verified protein-protein interactions was prepared as Fig. 7B. In the group of unique proteins for the lipid type of salivary stones, there are 4 quantified proteins. These proteins were marked with yellow circles on the STRING-generated network. The group of LIP unique proteins includes transthyretin (TTR), lactotransferrin (LTF), matrix Gla protein (MGP) and submaxillary gland androgen-regulated protein 3A (SMR3A). Transthyretin and submaxillary gland androgen-regulated protein 3A were down-regulated in all lipid samples. Lactotransferrin and matrix Gla protein were up-regulated in only one LIP20 sample but down-regulated in the others. This sample in the PCA analysis deviated the most from the group and was in the range of mixed stones. Transthyretin (TTR) can form oligomers, which are responsible for increasing the concentration of calcium ions in cells, so down-regulation of TTR can be caused a lower level of calcium in lipid salivary stones [149]. There is a similar case with matrix Gla protein (MGP) – this protein has a high affinity to calcium ions, so decreasing the level of MGP can influence calcium balance [150].
Another example is lactotransferrin (LTF), which was mentioned before. As we know, this antimicrobial protein and its activity depend on the level of extracellular cations. It was assumed that in lipid sialoliths, the concentration of calcium ions is lower, and the down-regulation of this protein in most samples causes antimicrobial activity inhibition. Submaxillary gland androgen-regulated protein 3A (SMR3A) is secreted only by submaxillary glands into saliva, so its decreased level in all samples can be evidence of activity and functions of salivary glands are disturbed [151]. Non-unique but down-regulated for a significant part of all of the lipid samples is the eosinophil cationic protein (RNASE3), mentioned earlier, responsible for defence activity against the pathogens. Azurocidin 1 (AZU1), also known as cationic antimicrobial protein CAP37 or heparin-binding protein, is an important multifunctional inflammatory mediator, and its activity is mainly directed against Gram-negative bacteria [152]. Azurocidin 1 has inhibitory activity during periodontitis, so its variable level can cause different pathological states [153]. The activating function of azurocidin 1 concerning the macrophages is connected to the immobilisation of calcium ions, which is decreased in this group of sialoliths. However, in consequence, the concentration of Ca2+ in cells of salivary gland tissue can be higher, and that imbalance causes different dysfunctions [154]. This protein is also part of NETs [128]. Focusing on the different types of regulation of immunoglobulin alpha-2 heavy chain (P0DOX2), the body’s response to pathological bacteria is variable. Statherin (STATH) was mentioned earlier as a protein responsible for stabilising saliva supersaturated with calcium salts by inhibiting the precipitation of calcium phosphate salts. A variable level of STATH can indicate calcium imbalance in lipid stones.
Also, for this group g:Profiler bioinformatic tool enrichment analysis was used, taking into account the following annotation category: Gene Ontology (Biological Process, Cellular Component, Molecular Function), KEGG, Reactome and type of regulation: up-regulation, down-regulation or varied level of regulation of proteins among the samples in the group (Fig. 9). The uniqueness of proteins for the group was marked with yellow names of proteins. Focusing on Biological Process GO terms, about half of the proteins are associated with the body’s defence against bacteria, so this part of selected proteins is smaller than in the case of the CAL group. Most of them have a variable level of regulation among the sample in the group. Only 2 proteins, eosinophil cationic protein (RNASE3) and myeloperoxidase (MPO), are down-regulated for all of the samples in the group. This imbalance directly shows the influence of bacteria on the pathological state leading to the formation of sialoliths. However, comparing this analysis with results for calcified stones, it is suggested that the body’s response against the pathogens is at a lower level. Several proteins were connected with the body’s activity against fungus, so there is the possibility that oral fungal infections can influence the development of salivary stone disease, but similar terms were not connected with proteins selected for calcified stones. Some of the proteins are responsible for the regulation of endopeptidase activity. Most often, these proteins have a variable level of regulation, 3 proteins: submaxillary gland androgen-regulated protein 3B (SMR3B), leukocyte elastase inhibitor (SERPINB1), submaxillary gland androgen regulated protein 3A (SMR3A) are down-regulated among the whole group. Endopeptidase activity is associated with lipids, so disturbing levels of these proteins can result from the imbalance of lipids [155]. These functions were not detected in the case of calcified sialoliths. Detected Cellular Component GO terms show the connection between the formation of salivary stones and neutrophil extracellular matrix. In this process, all of the proteins from the LIP group are involved, and they usually have a variable type of regulation, but 5 proteins: eosinophil cationic protein (RNASE3), transthyretin (TTR), submaxillary gland androgen-regulated protein 3B (SMR3B), myeloperoxidase (MPO) and leukocyte elastase inhibitor (SERPINB1) are down-regulated for all of the samples. Analysing Molecular Function GO terms, several proteins are again involved in peptidases and endopeptidases activity. 3 of these proteins are down-regulated, the same as those detected for Biological Process GO. Again, the Neutrophil extracellular trap formation KEGG pathway and Neutrophil degranulation Reactome pathway were detected, confirming the pivotal role of NETs in biocalcification leading to the salivary stones formation. The values of -log(q-value) are lower, but it is probably the result of the less numerous group of analysing proteins.
Mixed sialoliths.
For the chosen set of the most common, statistically significant proteins quantified in more than 50% of mixed salivary stones samples, the STRING-generated network presenting experimentally verified protein-protein interactions is presented in Fig. 7C. In the group of unique proteins for the mixed type of salivary stones, there is only 1 quantified protein – fibrinogen alpha chain (FGA). This protein is marked with a yellow circle on the STRING-generated network. FGA was down-regulated for all of the samples in the group. Fibrinogen is allowed to protect the neutrophils against the cytotoxic effects caused by, for example, the presence of bacteria and, consequently, the formation of neutrophil extracellular traps is delayed. The reduced level of fibrinogen alpha chain can cause the abnormal activity of neutrophils, leading to the formation of NETs [156], [157].
Non-unique but down-regulated in many samples was eosinophil cationic protein (RNASE3) and immunoglobulin gamma-1 heavy chain (P0DOX5), pointing to the disturbance in the immune response. Other non-unique frequent proteins have a different level of regulation among the MIX samples. These include, for example, haemoglobin subunit beta (HBB), which can harm oral mucosa because of the presence of the reactive heme group [158]. Cathepsin G (CTSG) and leukocyte elastase inhibitor (SERPINB1) are engaged in the immune response of the body against pathogens [159]. Zinc-alph-2-glycoprotin (AZGP1) is responsible for the degradation of lipids [160], but on the other hand, carbonic anhydrase (CA1) is a protein present in saliva, and its primary function is controlling the process of calcification [161].
With use g:Profiler bioinformatic tool enrichment analysis was performed taking into account the following annotation category: Gene Ontology (Biological Process, Cellular Component, Molecular Function), KEGG, Reactome and type of regulation: up-regulation, down-regulation or varied level of regulation of proteins among the samples in the group (Fig. 10). The uniqueness of proteins for the group was marked with yellow names of proteins.
Analysing Biological Process GO terms, about half of the proteins are associated with the body’s defence against bacteria, so this part of selected proteins is smaller than in the case of the CAL group but roughly equal to the part of the LIP group. Most of them have a variable level of regulation among the sample in the group. Only 4 proteins, eosinophil cationic protein (RNASE3), azurocidin (AZU1), fibrinogen alpha chain (FGA) and neutrophil elastase (ELANE), are down-regulated for all of the samples in the group. Again, as in the case of lipid sialoliths, we can suppose that the body’s response against the pathogens is lower than calcified salivary stones. Terms describing the body’s activity against fungus were not detected in this case. Some proteins are responsible for the regulation of endopeptidase activity. However, this group of proteins is less numerous, and the values of -log(q-value) are lower than in the enrichment analysis of the LIP group. The Biological Processes connected with the regulation of endopeptidase activity were not identified for MIX sialoliths. Cellular Component GO terms indicate the role of neutrophil extracellular matrix in the sialolithiasis, and analysing the values of -log(q-value), these terms have higher significance than in the case of the lipid group but lower than the calcified group. For these terms, most of the proteins from the MIX group are involved. They usually have a variable type of regulation, but 5 proteins: eosinophil cationic protein (RNASE3), azurocidin (AZU1), fibrinogen alpha chain (FGA), submaxillary gland androgen-regulated protein 3B (SMR3B) and neutrophil elastase (ELANE) are down-regulated for all of the samples. About Molecular Function GO terms, fewer proteins are involved in peptidases and endopeptidase activity compared to the lipid sialoliths. The significance is also on the lower level. Neutrophil extracellular trap formation and Salivary secretion KEGG pathways and Neutrophil degranulation Reactome pathway were detected again. Proteins connected with these 2 terms have mainly variable levels of regulation.
Comparison of CAL, LIP and MIX sialoliths
The unique protein for the calcified and lipid group was immunoglobulin lambda-like polypeptide 5 (IGLL5), engaged in the immune system response. In all of the calcified samples for which this protein was quantified (5 of 6), its level was up-regulated, and the level of its regulation in lipid samples was varied – up-regulated in 1 sample and down-regulated in 2 samples. The unique proteins for calcified and mixed sialolith groups included 8 proteins. Carbonic anhydrase 1 (CA1), a protein responsible for calcification, has a variable type of regulation for both CAL and MIX groups, pointing to the imbalance of calcium. Up-regulated immunoglobulin mu heavy chain (P0DOX6) and variable immunoglobulin lambda constant 2 (P0DOY2) are evidence of the immune system's activity caused by bacteria. Albumin (ALB) can bind calcium ions, so the up-regulation of this protein in the CAL group and variable level in the MIX group is completely understandable [162]. The same types of regulation as for albumin were detected for mentioned earlier keratin, type I cytoskeletal 19 (KRT19). This protein-fixing structural balance can change the metabolism of lipids. Serotransferrin (TF) with the same type of regulation belongs to the transferrin family typical in saliva. Hence, its activity depends on the concentration of ions and calcium ions, which is crucial in calcified sialoliths, where TF is up-regulated. Alpha-amylase 1A (P0DUB6), common in saliva, is a calcium-binding protein, so its variable level in salivary stones indicates calcium imbalance [163]. Neutrophil gelatinase-associated lipocalin (LCN2) takes part in the regulation of immune response in the presence of bacteria, and as a member of the lipocalin family, this protein is responsible for the transport of hydrophobic molecules, such as lipids. Down-regulation of LCN2 can cause lipids deficiency in calcified salivary stones. Besides, LCN2 is also part of NETs [128].
On the other hand, the set of unique proteins for lipid and mixed sialolith groups included 4 proteins. Variable levels of protein S100-A9 (S100A9) and zinc-alpha-2-glycoprotein (AZGP1) can influence the balance of calcium and lipids, respectively. Besides, S100A9 modifies the formation of NETs. Haemoglobin subunit alpha (HBA2), such as another described earlier subunit, haemoglobin subunit beta (HBB), can cause damage to the oral mucosa, accelerating inflammation. Leukocyte elastase inhibitor (SERPINB1) and its down-regulation in the LIP group and variable level in the MIX group show that the body’s response to the presence of bacteria is clearly at a lower level in the case of lipid sialoliths. The unique proteins for these 2 pairs of sialolith groups are present in Table 2, considering the regulation type.
The Standard set of proteins for each type of salivary stone included 13 proteins: azurocidin (AZU1), cystatin-SN (CST1), cystatin-S (CST4), cathepsin G (CTSG), neutrophil elastase (ELANE), haemoglobin subunit beta (HBB), histone H3,1 (HIST1H3J), myeloperoxidase (MPO), eosinophil cationic protein (RNASE3), submaxillary gland androgen-regulated protein 3B (SMR3B), statherin (STATH), immunoglobulin alpha-2 heavy chain (P0DOX2), immunoglobulin gamma-1 heavy chain (P0DOX5). The level of regulation of each protein for all samples is presented as a heatmap (Fig. 11).
For the standard statistically significant proteins quantified in more than 50% of samples in each group, the STRING-generated network presenting experimentally verified protein-protein interactions was prepared (Fig. 12). The network presents the type of regulation of proteins in each group and their frequency among the samples. Most of the proteins were mentioned above. The activity of azurocidin (AZU1) concerning the role of macrophages during the body's immune response depends on calcium ions inside the cells. Down-regulation of AZU1 in the case of calcium and mixed sialoliths can be associated with immobilising a more significant amount of calcium in stones, thus reducing intracellular Ca2+. Without a proper amount of calcium ions, azurocidin cannot fulfil its functions, even in the presence of pathogens. The role of azurocidin in the creation of NETs is also reduced. A disturbed calcium balance in lipid stones leads to a variable level of AZU1 in these sialoliths. Myeloperoxidase (MPO) has a high affinity for calcium, so different concentrations of MPO in CAL and MIX groups point to the imbalance of this mineral compound in salivary stones and the presence of bacteria. Down-regulation of MPO in the case of lipid stones shows the lower concentration of calcium in these sialoliths, and immobilisation of lipids, so their peroxidation is reduced, reduced role of this protein during NETosis and the smaller number of identified bacteria in lipid sialoliths. Neutrophils secrete neutrophil elastase (ELANE) during inflammation, so down-regulation of this protein in calcified and mixed groups suggests a lower level of immune response mediated by this protein and reduced part of ELANE in the neutrophil extracellular trap. Cathepsin G (CTSG), as a mediator of metabolisms of calcium and lipids, quantified in all of the groups on the different levels, can cause calcium-lipid imbalance and influence the formation of NETs. At the same time, its presence clearly shows the influence of bacteria on the calcification process. Haemoglobin subunit beta (HBB) is down-regulated for calcified and lipid sialoliths, suggesting that the risk of harming oral mucosa is reduced. The role of statherin (STATH) is the protection of saliva from the precipitation of calcium phosphate salts in high concentrations. Reduced levels of STATH in calcified stones can be associated with lower calcium salts in saliva concentration caused by the immobilization of calcium in sialoliths. Submaxillary gland androgen-regulated protein 3B (SMR3B) is a poorly described protein, but there is a prediction that SMR3B can inhibit the activity of peptidases and endopeptidases, which, as it was mentioned above, are connected with the presence of lipid [164]. Variable types of regulation of this protein in LIP and MIX stones indicate some imbalance of this mineral compound in these stones when the lipid concentration is higher in the LIP group, inhibiting peptidases and endopeptidases can be difficult. This way, the lowest level of submaxillary gland androgen-regulated protein 3B in lipid sialoliths is explained. Cystatin-SN (CST1) and cystatin-S (CST4) are similar proteins common in saliva. Their primary function is the inhibition of human cathepsins, which is why the level of cathepsin G (CTSG) is variable [53]. Cystatins S and SN can bind calcium, so their up-regulation in CAL sialoliths show a high concentration of this mineral compound. Different level in LIP and MIX salivary stones probably causes an imbalance of calcium in these groups. Histones are nuclear proteins that tightly pack DNA into chromatin [165]. During NETosis, chromatin is decondensed, and the level of decondensation can be regulated by histone H3,1 (HIST1H3J) [166]. Variable types of regulation of HIST1H3J in all groups may influence chromatin decondensation, leading to the creation of neutrophil extracellular matrix and calcification. Eosinophil cationic protein (RNASE3), immunoglobulin alpha-2 heavy chain (P0DOX2) and immunoglobulin gamma-1 heavy chain (P0DOX5) are mainly responsible for the body’s defence against pathogens, so down-regulation of eosinophil cationic protein and immunoglobulin gamma-1 heavy chain and variable level of immunoglobulin alpha-2 heavy chain shows, that the immune response in the presence of bacteria is disturbed. That effect can have a strong influence on the formation of sialoliths. 6 proteins have a common type of regulation comparing pairs CAL-MIX (azurocidin, myeloperoxidase, neutrophil elastase) and LIP-MIX (statherin, cystatin-SN, cystatin-S), proving that mixed sialoliths have features of both calcified and lipid stones. The following 6 proteins are down-regulated (eosinophil cationic protein, submaxillary gland androgen-regulated protein 3B, immunoglobulin gamma-1 heavy chain) or have variable levels (cathepsin G, immunoglobulin alpha-2 heavy chain, histone H3,1) for all of the groups. Haemoglobin subunit beta is regulated differently in the case of mixed salivary stones.
With the use of g:Profiler bioinformatic tool enrichment analysis was performed taking into account the following annotation categories: Gene Ontology (Biological Process, Cellular Component, Molecular Function), KEGG, Reactome and type of regulation: up-regulation, down-regulation or varied level of regulation of proteins among the samples in the groups (Fig. 13).
Biological Process GO terms detected for the set of common 13 proteins for all types of sialoliths are mainly associated with the body’s immune defence against pathogenic bacteria. 7 quantified proteins evidence it: azurocidin (AZU1), cystatin-S (CST4), cystatin-SN (CST1), submaxillary gland androgen-regulated protein 3B (SMR3B), myeloperoxidase (MPO), neutrophil elastase (ELANE) and cathepsin G (CTSG). In Cellular Component GO enrichment, detected terms are connected to extracellular components, pointing to the creation of neutrophil extracellular traps during calcification leading to the formation of sialoliths. There are 6 proteins responsible for the regulation of this process: statherin (STATH), azurocidin (AZU1), eosinophil cationic protein (RNASE3), myeloperoxidase (MPO), neutrophil elastase (ELANE) and cathepsin G (CTSG). Focusing on the Molecular Functional GO enrichment, summarizing all of the proteins are associated with modulation of peptidases and endopeptidases, which activity can be disturbed by the imbalance of lipids. Also, for this set of proteins, the Neutrophil extracellular trap formation KEGG pathway was detected, highlighting the role of NETosis in forming stones. In this pathway, there are engaged 5 proteins: azurocidin (AZU1), myeloperoxidase (MPO), neutrophil elastase (ELANE), cathepsin G (CTSG) and histone H3,1 (HIST1H3J). The influence on the Salivary secretion KEGG pathway has statherin (STATH), Cystatin-SN (CST1) and cystatin-S (CST4), and this can also have negative consequences. The Reactome pathways concern mainly immune response, but the most important for NETosis is Neutrophil degranulation Reactome pathway engaging 6 proteins: azurocidin (AZU1), eosinophil cationic protein (RNASE3), myeloperoxidase (MPO), neutrophil elastase (ELANE), cathepsin G (CTSG) and histone H3,1 (HIST1H3J).