Skip to main content

Label-free quantitative proteomic analysis of gingival crevicular fluid to identify potential early markers for root resorption



Orthodontically-induced root resorption is an iatrogenic effect and it cannot be examined regularly due to the harmful effects of sequential doses of radiation with more frequent radiography. This study aims to compare protein abundance (PA) of pre-treatment and during orthodontic treatment for root resorption and to determine potential early markers for root resorption.


Ten subjects (n = 10) who had upper and lower fixed appliances (MBT, 3 M Unitek, 0.022″ × 0.028″) were recruited for this study. Human gingival crevicular fluid (GCF) was obtained using periopaper strips at pre-treatment (T0), 1 month (T1), 3 months (T3), and 6 months (T6) of orthodontic treatment. Periapical radiographs of the upper permanent central incisors were taken at T0 and T6 to measure the amount of root resorption. Identification of changes in PA was performed using liquid chromatography-tandem mass spectrometry. Student’s t-test was then performed to determine the significance of the differences in protein abundance before and after orthodontic treatment.


Our findings showed that all ten subjects had mild root resorption, with an average resorption length of 0.56 ± 0.30 mm. A total of 186 proteins were found to be commonly present at T0, T1, T3, and T6. There were significant changes in the abundance of 16 proteins (student’s t-test, p ≤ 0.05). The increased PA of S100A9, immunoglobulin J chain, heat shock protein 1A, immunoglobulin heavy variable 4–34 and vitronectin at T1 suggested a response to stress that involved inflammation during the early phase of orthodontic treatment. On the other hand, the increased PA of thymidine phosphorylase at T3 suggested growth promotion and, angiogenic and chemotactic activities.


The identified proteins can be potential early markers for root resorption based on the increase in their respective PA and predicted roles during the early phase of orthodontic treatment. Non-invasive detection of root resorption using protein markers as early as possible is extremely important as it can aid orthodontists in successful orthodontic treatment.

Peer Review reports


Orthodontic treatment has iatrogenic effects that can be detrimental to the stability of teeth and shorten their roots. Root resorption occurs when the root is shortened, which often results in tooth extraction (in severe cases) [1]. Orthodontic treatment applies a mechanical force [2] that is capable of creating external root resorption in the patients, hence shortening the roots [3, 4]. Gingival crevicular fluid (GCF) is an oral cavity fluid and its content may be representative of the body’s health [5,6,7]. GCF contains proteins involved in oral pathologic phenomenon, such as internal root resorption [8, 9], as well as a diverse population of cells, desquamated epithelial cells, organic ions, inorganic ions, enzymes, and bacteria from adjacent plaque. GCF has attracted researchers for its potential use as a diagnostic fluid for periodontal diseases and drug analyses due to its non-invasive sampling [8,9,10]. As such, we believe that GCF is a reliable source for root resorption biomarkers.

Radiography is currently the primary approach in determining root resorption. A radiograph is an image produced on a sensitive plate or film by X-rays, gamma rays, or similar radiation [11, 12]. Frequent radiographs may have negative effects on patients [1, 13, 14]. As radiography cannot be conducted regularly, it limits its utility as a tool for monitoring a patient’s teeth. Therefore, identification and the subsequent utilization of biomarkers for root resorption would be a useful and attractive alternative to monitoring root resorption. A non-invasive way to determine the severity of root resorption is by detecting the presence of protein biomarkers in GCF. Comparative proteomics using liquid chromatography-mass spectrometry (LC-MS) can be used to screen and compare GCF proteome profiles and identify potential root resorption biomarkers. This study is a preliminary investigation to identify potential diagnostic biomarkers for root resorption during the early phase of orthodontic treatment. Eventually, in practice, detection of the identified biomarkers will be performed using protein assay or enzyme-linked immunosorbent assay (ELISA) technique. These biomarkers will allow continuous monitoring of root resorption and better management of the condition. Our study aims to compare the protein abundance (PA) of GCF at pre-treatment (T0), after 1 month (T1), 3 months (T3), and 6 months (T6) of orthodontic treatment. This is will allow us to identify potential early root resorption biomarkers.



Subjects were patients (n = 10) who have undergone treatment at the Postgraduate Orthodontic Clinic 1, Faculty of Dentistry, Universiti Kebangsaan Malaysia, Kuala Lumpur. All wore upper and lower fixed appliances (MBT, 3 M Unitek, 0.022″ × 0.028″). The exclusion criteria were patients who were pregnant, smokers, who had previous orthodontic treatment, a non-vital tooth, a treated root or fractured root, and patients who had consumed anti-inflammatory drugs in the month preceding the study.

Ethical approval and sample collection

This study was approved by the Ethics Research Committee, Universiti Kebangsaan Malaysia (UKM PPI/111/8/JEP-2018-268), and patients were briefed and written informed consent was obtained. GCF sampling was carried out using Periopaper strips (Proflow, USA). Strips were inserted at the gingival sulcus mesial and distal side of the upper permanent central incisor at approximately 1–2 mm of gingiva depth for 60 s. This was repeated three times with intervals of 60 s in between. Each crevicular sulcus on the test tooth was dried with cotton rolls, and a saliva ejector was used to remove the remaining saliva prior to insertion of the periopaper strip. Sampling was conducted at T0, T1, T3, and T6. The periopaper strips were then placed into microcentrifuge tubes containing protease inhibitors.


Periapical radiographs of the upper permanent central incisor were taken at T0 and T6 to measure the amount of root resorption by the same radiographer.

Extraction and digestion of proteins

GCF samples were centrifuged at 10,000×g for 15 min at 4 °C and proteins were extracted. Protein content was estimated by Bradford Assay. Ammonium bicarbonate (50 mM) and dithiothreitol (100 mM) were put into 1.5 mL microcentrifuge tubes. A total of 100 ng of the GCF protein was then added to the tubes, and incubated at 95 °C for 5 min. The mixtures were allowed to cool at room temperature for 5 min. Iodoacetamide (100 mM) was added to the tubes and incubated in the dark at room temperature for 20 min. Trypsin (100 ng/μL) was then pipetted into the reaction tube and incubated at 37 °C for 3 h, followed by overnight incubation at 30 °C. Digested samples were dried using a SCANVAC Vacuum Concentrator (LaboGene, Denmark) at low speed. The samples were then dissolved in 0.1% formic acid in the reaction tubes and with vigorous vortexing.

Separation of digested peptides

Digested peptides were resolved using a Dionex UltiMate™ 3000 RSLCnano system (Thermo Fisher Scientific, USA) with an EASY-Spray Column Acclaim PepMap™ C18 100 A0 separation column, a 2 μm particle size, and a 50 μm internal diameter with 15 cm in length. The flow rate was 250 nL/min and mobile phase A and B contained 0.1% formic acid in water and 0.1% formic acid in acetonitrile, respectively. The gradient was the following: 5–40% mobile phase B for 91 min, 85% mobile phase B for 2 min, 85% mobile phase B for 3 min, 5% mobile phase B for 1 min and 5% mobile phase B for 4 min. The column temperature was maintained at 35 °C.

Mass spectrometry data acquisition

The mass spectrometric analysis was conducted using the Orbitrap Fusion Mass Spectrometer (Thermo Fisher Scientific, USA) in data-dependent mode. Full scan spectra were collected with orbitrap MS (OTMS) using the following parameters: scan range, 310–1800 m/z; resolving power of 120,000; automatic gain control (AGC) target of 4.0e5 (400,000); and maximum injection time of 50 ms. The method consisted of 3 s of top speed mode where the precursors were selected for a maximum of 3 s cycles. Only precursors with an assigned monoisotopic m/z and a charge state of 2–7 were further analysed for second stage mass spectrometry (MS2). All precursors were filtered using a 20 s dynamic exclusion window with intensity threshold of 5000. The MS2 spectra were analysed by ion trap mass spectrometry (ITMS) using the following parameters: rapid scan rate with a resolving power of 60,000, AGC target of 1.0e2 (100), 1.6 m/z isolation window, and a maximum injection time of 250 ms.

Protein identification

Proteins were identified using the Thermo Scientific™ Proteome Discoverer Software Version 2.1 where the mass spectra were searched against the Homo sapiens UniProt database [15]. All peptides were validated using the percolator® algorithm based on a q-value of less than a 5% false discovery rate (FDR). Following the identification of proteins, PANTHER software (Protein analysis through evolutionary relationships) was used to predict the biological processes affected during orthodontic treatment [16]. Search Tools or the Retrieval of Interesting Proteins/Genes (STRING) version 11.0 was then utilized to determine pathways involved during orthodontic treatment [17].

Statistical analysis

All analyses were performed using Microsoft Excel. The intra-examiner correlation coefficient (ICC) was calculated to determine the reliability of the measurements for the amount of root resorption at 6 months. Student’s t-test was used to determine the statistical signficance in PA during the first 6 months of orthodontic treatment. A difference was considered significant when p ≤ 0.05.


The differences in the root length ranged from 0.14–1.01 mm, while the crown length ranged from 0.15–0.91 mm (Table 1). The average difference in root resorption was 0.56 ± 0.30 mm, whereas the average difference in crown length was 0.37 ± 0.23 mm. Based on the radiographs, all subjects were shown to be having mild root resorption (< 2 mm). Since all subjects have similar root resorption rates, the samples were chosen randomly based on patient’s availability.

Table 1 Demographic, crown length and root length of orthodontic patients. C1 and R1 refer to the tooth length prior to orthodontic treatment, while C2 and R2 were taken after 6 months of orthodontic treatment. Subjects that were chosen labelled as*

Meanwhile, Fig. 1 shows the base-peak chromatogram for nano-LC monitored by mass spectrometry, representing the intensity of all peptide ions in the sample in a single scan. The GCF proteome from T0 to T6 showed consistent elution of proteins/peptides, ranging from 8 to 100 min. A total of 815 proteins were identified by LC-MS in human GCF. The total proteins at T0, T1, T3, and T6 were 439, 368, 354, and 502, respectively. Furthermore, 186 proteins were commonly found in all samples (Fig. 2). Following statistical analysis, 16 proteins were found to have changed significantly (Student’s t-test, p ≤ 0.05). They were S100A9, immunoglobulin J chain (IGJ), angiotensinogen (AGT), apolipoprotein D (APOD), beta-2-glycoprotein 1 (APOH), hemoglobin subunit beta (HBB), histidine-rich glycoprotein (HRG), immunoglobulin heavy variable 3–11 (IGHV3–11), immunoglobulin heavy constant gamma 1 (IGHG1), thymosin beta-4 (TMSB4X), zinc-alpha-2-glycoprotein (AZGP1), heat shock protein 1A (HSPA1A) (70 kDa), immunoglobulin heavy variable 4–34 (IGHV4–34), thymidine phosphorylase (TYMP), immunoglobulin kappa variable 3–20 (IGKV3–20), and vitronectin (VTN) (Table 2).

Fig. 1

Total ion current (TIC) chromatograms from three subjects (n = 3) with total of 12 samples and 3 injections/sample of human gingival crevicular fluid (GCF). a Prior to treatment (T0), b after 1 month of treatment (T1), c after 3 months of treatment (T3), and (d) after 6 months of treatment (T6). The graph shows the total intensity versus time (min). Two microlitres of the one microgram per microlitre sample was loaded for 110 min, with a flow rate of 250 nL/min in 0.1% mobile phase B (0.1% formic acid in acetonitrile)

Fig. 2

Venn diagram of a number of proteins found in the GCF of 3 individual subjects before orthodontic treatment (T0), after 1 month of treatment (T1), after 3 months of treatment (T3), and after 6 months of treatment (T6). The total protein counts at T0, T1, T3, and T6 were 436, 365, 351, and 498, respectively. One hundred eighty-six proteins (circled)-- were commonly found in all samples

Table 2 Biological process, protein class, and pathway of significant proteins. Cellular process (CP), metabolic process (MP), immune system process (ISP), localisation (L), response to stimulus (RS), and biological regulation (BR). The proteins found to be involved in the biological processes are represented with symbol X

These proteins were then analysed using the Protein Analysis Through Evolutionary Relationships (PANTHER, database. This led to the identification of six biological processes [18] that may be associated with these changes. Of the proteins, 24.4% changed significantly and they were involved in cellular processes (Fig. 3a). Other biological processes (more than 10%) were associated with metabolic processes, immune system processes, localisation, response to stimuli, and biological regulation, as shown in Fig. 3a. PANTHER found five protein classes that were involving during treatment (Fig. 3b). They were a calcium-binding proteins, defence or immunity proteins, enzyme modulators, signalling molecules, and transferases [18]. Two proteins, AGT and HRG, out of sixteen proteins were enzyme modulators. The SEQUEST score, pI/MW, peptide matched (% sequence coverage), fold change, and p-value of significant proteins are shown in Table 3.

Fig. 3

Pie chart representing the; a) biological process, b) protein class, and c) pathway of the proteins that changed significantly based on Protein Analysis Through Evolutionary Relationships (PANTHER) database. The biological processes involved were cellular processes, metabolic processes, immune system processes, localisation, responses to stimuli, and biological regulation. The protein classes involved were calcium-binding proteins, defence or immunity proteins, enzyme modulators, signalling molecules, and transferase. While, the identified proteins also involved in other pathways were angiotensin II-stimulated signalling through G proteins and beta-arrestin, the apoptosis signalling pathway, the gonadotropin-releasing hormone receptor pathway, and Parkinson pathway

Table 3 Differentially expressed proteins in GCF at T1, T3, and T6 during orthodontic treatment when compared to T0 (pre-treatment). There was a change in the expression of 16 proteins in the GCF. Increased and decreased protein abundance (PA) of significant proteins are represented in the form of ↑ and ↓ (Student’s T-test, p-value ≤0.05)

S100A9 is a calcium-binding protein and a signalling molecule. The identified proteins were also involved in other pathways, as suggested by PANTHER, and these were the angiotensinogen II-stimulated signalling through G proteins and β-arrestin [19], the apoptosis signalling pathway [18], the gonadotropin-releasing hormone receptor pathway [20], and Parkinson pathway [21] (Fig. 3c). Proteins that were significantly different in abundance were also analysed using STRING analysis. Thirteen proteins were predicted to be associated with six pathways (Table 4). The interactions between these proteins are shown in Fig. 4.

Table 4 Protein interaction network found in the GCF with an additional seven proteins in the Search Tools or the Retrieval of Interesting Proteins/Genes (STRING) analysis. A total of six pathways were found during orthodontic treatment
Fig. 4

Protein interaction networks of proteins that significantly changed in abundance, determined using STRING version 11.0. A total of 13 proteins were analysed. Seven additional proteins were added to determine protein interactions during orthodontic treatment. Each protein was represented with a node and its respective gene name. The colours of the node representations are as follows: red, positive regulation of biological process; dark green, response to stimulus; purple, metabolic process; dark blue, cellular process; light green, immune system process; and yellow, regulation of localisation. Meanwhile, interactions of proteins were identified using the edges connecting the nodes. The thickness of the edge is represented as follows: thinnest, low confidence (interaction score: 0.150); thin, medium confidence (interaction score: 0.400); thick, high confidence (interaction score: 0.700); and thickest, highest confidence (interaction score: 0.900)



Root resorption is an unavoidable consequence of orthodontic treatment. At present, root resorption can only be detected using radiographs. Periapical radiographs have previously been used to measure root resorption [1, 14, 22]. In this study, periapical radiographs were obtained on central incisors at T0 and T6 as evidence of root resorption. The roots and crowns of all subjects were measured three times. Although radiography is the gold standard for measurement of root resorption, it has limitations. Radiation exposure due to frequent radiography procedures can be detrimental to humans and may lead to radiation carcinogenesis [23]. The dosage of radiation should be as low as rationally obtainable, and the number of radiographs should be chosen wisely because the risk will increase as the frequency of X-ray exposure increases [23]. In this study, a periapical radiograph was obtained twice, before treatment and at 6 months of orthodontic treatment, and this was found to be acceptable [24].

Root resorption

All subjects were shown to have mild root resorption, most likely as they were in the early stage of orthodontic treatment (first 6 months). This has been previously observed in other reports (< 2 mm at sixth months of orthodontic treatment) [25]. GCF samples from the first sixth months of orthodontic treatment is therefore a useful source for the identification of proteins involved in root resorption. Prolonged orthodontic treatment may deteriorate the root, potentially causing severe root resorption [26]. It has been shown that severe root resorption (≥ 2 mm) occured in patients after 12 months of orthodontic treatment [25]. The long duration of orthodontic treatment caused the teeth to be exposed to a long-term jiggling force, causing an increase in root resorption [27]. In this study, we targeted proteins from the early phase of orthodontic treatment that can be useful root resorption early markers. Once identified and confirmed, in practice, detection of the identified biomarkers will be conducted using protein assay or ELISA techniques. When a certain threshold or panel of protein biomarkers have been detected, orthodontic force can be lowered or treatment can be stopped to prevent severe root resorption.

Gingival crevicular fluid (GCF) biomarker

To date, there are no reliable biomarkers from GCF that can be used in measuring root resorption during orthodontic treatment. Previous studies have been focused on identifying biomarkers related to periodontal disease, systemic diseases, and pubertal growth [28,29,30,31]. Furthermore, other previous studies used non-resorbing and resorbing deciduous molars to identify unique root resorption biomarkers [32].

Due to its ability to be site-specific, GCF can be an unparalleled source for identifying specific dental biomarkers of oral conditions [33]. GCF is known to be unique and different from saliva [34, 35]. Many researchers have classified GCF as an inflammatory exudate, and some have suggested that it is a transformed tissue transudate in a normal healthy state [7, 34, 36, 37], which consists of proteins, a distinct population of cells, exfoliated epithelial cells, and bacteria from neighbouring plaque [5]. In this study, we compared samples from T0, which is control (normal) to samples at T1, T3 and T6. T0 samples were acquired before subjects began wearing braces.


Proteomics has many advantages in the identification of potential disease biomarkers as it can detect low levels of a specific proteins [38]. In relation to our work, a biomarker can be used to detect the severity of root resorption. In this study, we used proteomics to determine the protein profile of GCF samples before treatment and at 6 months of orthodontic treatment. Previous studies have also used a proteomics approach to identify disease biomarkers [39, 40]. A separate study used proteomics to identify the protein expression of bacteria [41].

Potential early markers for root resorption

Table 3 shows that both S100A9 and IGJ had increased PA from T0 to T6. S100A9 is a calcium- and zinc-binding protein that plays a prominent role in the regulation of inflammatory processes and immune responses. It can induce neutrophil chemotaxis and adhesion, and it can increase the bactericidal activity of neutrophils by promoting phagocytosis [42, 43]. However, this protein was reported to be downregulated on day 14 compared to day 0 of orthodontic treatment [44]. This suggested that activation of inflammation occurred after 1 month as this protein has known roles in inflammation and bone resorption [44]. Another study stated that S100A9 was distinguished in osteoclasts and was involved in chondrocyte and osteoblast maturation [43].

IGJ (J Chain) has a prominent role in the binding immunoglobulin units - IgM or IgA [45]. It can induce the formation of larger polymers when binding to IgA. It also helps to bind immunoglobulins to secretory component [44]. A previous study have found significant differences in the mean serum levels of IgA, IgG, and IgM between well-controlled type 1 diabetic patients that wore orthodontic appliances and those who did not [46]. IGJ has been found to be involved in the inflammatory process and immune response during the early phase of orthodontic treatment. A separate study found increased expression of IGJ in type 2 diabetes mellitus patients with periodontitis compared to controls [47]. Periodontitis is a condition where the periodontium is inflamed, mainly induced by a bacterial infection [9]. This may provide an explanation for the high level of IGJ during the inflammatory process. Nonetheless, a past study reported that IGJ was downregulated after 14 days of orthodontic treatment [44]. This was likely due to the inflammation process not being active on day 14 of orthodontic treatment [48].

Conversely, AGT, APOD, APOH, HBB, HRG, IGHV3–11, IGHG1, TMSB4X, and AZGP1 showed decreased PA in all samples. A decrease in PA suggested that the proteins were indirectly involved during orthodontic treatment. An increased PA of HSPA1A and IGHV4–34 after 1 month of orthodontic treatment indicated a possible response to stress during the early phase of orthodontic treatment. HSPA1A was found to be associated with an apoptosis signalling pathway, a gonadotropin-releasing hormone receptor pathway, and Parkinson pathway (Fig. 3c). The acute inflammatory response has been shown to be involved during the early phase of orthodontic treatment, which may lead to apoptosis or cell death when induced by mechanical force [42]. Meanwhile, the increased PA of TYMP at T3 indicated growth promotion and angiogenic and chemotactic activity. TYMP has been found to have roles in maintaining the integrity of blood vessels, growth-promoting activity on endothelial cells, angiogenic activity in vivo, and chemotactic activity on endothelial cells in vitro [49]. A previous study reported that the production of angiogenic regulators in cultured human gingival and periodontal ligament fibroblasts was caused by mechanical stress [50].

IGKV3–20 showed increased PA after 6 months of orthodontic treatment, indicating that an immune response took place. VTN showed increased PA at T1 and T6, suggesting This that VTN was involved during inflammation and that it was affected by a mechanical force during the early phase of orthodontic treatment [51]. VTN is a cell adhesion and spreading factor found in serum and tissues. It interacts with glycosaminoglycans and proteoglycans [52]. It also acts as an inhibitor of the membrane-damaging effect of the terminal cytolytic complement pathway.

We found that there were significant differences in protein abundance pre-treatment and during orthodontic treatment. Of the proteins that showed differences in PA, only S100A9, IGJ, HSPA1A, IGHV4–34, TYMP, and VTN has the potential to be developed into potential early markers for root resorption, based on increases in PA during the early phase of orthodontic treatment. However, several limitations should be noted. Firstly, the sample size was small, and add to that the fact that the collected GCF samples were in a proportion of 9:1 female to male. Secondly, only three out of ten subjects were chosen randomly for further analyses using Orbitrap LC-MS/MS due to financial constraints. A total of 12 samples were collected at four time-points (T0, T1, T3 and T6) from three subjects. As all subjects have similar root resorption, the samples were selected randomly based on patient’s availability. With the current data, it was not possible to group based on low, medium, and high root resorption.


Early detection of root resorption is extremely essential and may help orthodontists in patient management. The use of biomarkers is a better alternative to radiography as it will reduce unnecessary exposure to radiation. Our research has identified significant differences in protein abundance between pre-treatment and during orthodontic treatment. Out of these differences, six proteins - S100A9, IGJ, HSPA1A, IGHV4–34, TYMP, and VTN - showed significant increases in PA during the early phase of orthodontic tooth movement, making them prime candidates as potential biomarkers that can be exploited for detection of early root resorption.

Availability of data and materials

The datasets generated and analysed during the current study are available in the Figshare repository,



Automatic gain control




Apolipoprotein D


Beta-2-glycoprotein 1




Biological regulation


Cellular process


Enzyme-linked immunosorbent assay


False discovery rate


Gingival crevicular fluid


Heat shock 70 kDa protein 1A


Hemoglobin subunit beta


Histidine-rich glycoprotein


Intra-examiner correlation coefficient


Immunoglobulin heavy constant gamma 1


Immunoglobulin heavy variable 3–11


Immunoglobulin kappa variable 3–20


Immunoglobulin heavy variable 4–34


Immunoglobulin J chain


Immune system process


Ion trap mass spectrometry






Liquid chromatography


Liquid chromatography-mass spectrometry


Liquid chromatography-tandem mass spectrometry










Metabolic process




Second stage mass spectrometry


Molecular weight


Mass-to-charge ratio




Orbitrap mass spectrometry


Protein abundance


Protein analysis through evolutionary relationships


Isoelectric point

RS :

Response to stimulus


Protein S100-A9


Search tools or the retrieval of interesting proteins/genes




One month of orthodontic treatment


Three months of orthodontic treatment


Six months of orthodontic treatment


Total ion current


Thymosin beta-4


Thymidine phosphorylase




Times gravity


Degree Celcius








  1. 1.

    Kaczor-Urbanowicz KE, Deutsch O, Zaks B, Krief G, Chaushu S, Palmon A. Identification of salivary protein biomarkers for orthodontically induced inflammatory root resorption. Proteomics Clin Appl. 2017;11(9–10):1–27.

    Google Scholar 

  2. 2.

    Hafiz A, Zahid D, Yushar A, Ibrahim Z, Mastura S, Octavianti F, et al. The psychosocial impact of dental aesthetics on orthodontics patients. J Int Dental Med Res. 2019;12(2):563–8.

    Google Scholar 

  3. 3.

    Machado R, Leoni GB, Yarid P, Cruz-Filho AM, Sousa-Neto MD. Idiopathic external apical root resorption: clinical report of 2 cases. Gen Dent. 2017;65(2):62–5.

    PubMed  Google Scholar 

  4. 4.

    Lopatiene K, Dumbravaite A. Risk factors of root resorption after orthodontic treatment. Stomatologija. 2008;10(3):89–95.

    PubMed  Google Scholar 

  5. 5.

    Khurshid Z, Mali M, Naseem M, Najeeb S, Zafar MS. Human gingival Crevicular fluids (GCF) proteomics: an overview. Dent J. 2017;5(12):1–8 Available from:

    Google Scholar 

  6. 6.

    Tsuchida S, Satoh M, Umemura H, Sogawa K, Kawashima Y, Kado S, et al. Proteomic analysis of gingival crevicular fluid for discovery of novel periodontal disease markers. Proteomics. 2012;12(13):2190–202.

    PubMed  Google Scholar 

  7. 7.

    Yuan C, Liu X, Zheng S. Matrix metalloproteinase-8 levels in oral samples as a biomarker for periodontitis in the Chinese population: an observational study. BMC Oral Health. 2018;18(1):1–6.

    Google Scholar 

  8. 8.

    Rody WJ, Wijegunasinghe M, Shannon Holliday L, McHugh KP, Wallet SM. Immunoassay analysis of proteins in gingival crevicular fluid samples from resorbing teeth. Angle Orthod. 2016;86(2):187–92.

    PubMed  Google Scholar 

  9. 9.

    Silva-Boghossian CM, Colombo APV, Tanaka M, Rayo C, Xiao Y, Siqueira WL. Quantitative proteomic analysis of gingival Crevicular fluid in different periodontal conditions. PLoS One. 2013;8(10):1–15.

    Google Scholar 

  10. 10.

    Kereshanan S, Stephenson P, Waddington R. Identification of dentine sialoprotein in gingival crevicular fluid during physiological root resorption and orthodontic tooth movement. Eur J Orthod. 2008;30(3):307–14.

    PubMed  Google Scholar 

  11. 11.

    Makedonas D, Lund H, Gröndahl K, Hansen K. Root resorption diagnosed with cone beam computed tomography after 6 months of orthodontic treatment with fixed appliance and the relation to risk factors. Angle Orthod. 2012;82(2):196–201.

    PubMed  Google Scholar 

  12. 12.

    Patel S, Durack C, Abella F, Shemesh H, Roig M, Lemberg K. Cone beam computed tomography in Endodontics - a review. Int Endod J. 2015;48(1):3–15.

    PubMed  Google Scholar 

  13. 13.

    Creanga AG, Geha H, Sankar V, Teixeira FB, McMahan CA, Noujeim M. Accuracy of digital periapical radiography and cone-beam computed tomography in detecting external root resorption. Imaging Sci Dent. 2015;45(3):153–8.

    PubMed  PubMed Central  Google Scholar 

  14. 14.

    Pelagio C, Maués R, Ramos Do Nascimento R, De O, Vilella V. Severe root resorption resulting from orthodontic treatment: prevalence and risk factors. Dent Press J Orthod. 2015;20(1):52–8.

    Google Scholar 

  15. 15.

    Bateman A. UniProt: a worldwide hub of protein knowledge. Nucleic Acids Res. 2019;47:D506–15.

    Google Scholar 

  16. 16.

    Thomas PD, Campbell MJ, Kejariwal A, Mi H, Karlak B, Daverman R, et al. PANTHER: a library of protein families and subfamilies indexed by function. Genome Res. 2003;13(9):2129–41.

    PubMed  PubMed Central  Google Scholar 

  17. 17.

    Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(Database Issue):D607–13.

    PubMed  PubMed Central  Google Scholar 

  18. 18.

    Mi H, Muruganujan A, Ebert D, Huang X, Thomas PD. PANTHER version 14: more genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools. Nucleic Acids Res. 2019;1:1–8.

    Article  Google Scholar 

  19. 19.

    Wei H, Ahn S, Barnes WG, Lefkowitz RJ. Stable interaction between β-Arrestin 2 and angiotensin type 1A receptor is required for β-Arrestin 2-mediated activation of extracellular signal-regulated kinases 1 and 2*. J Biol Chem. 2004;279(46):48255–61.

    PubMed  Google Scholar 

  20. 20.

    Fink MY, Pincas H, Choi G, Nudelman G, Sealfon SC. Research resource : gonadotropin-releasing hormone receptor- mediated signaling network in LβT2 cells : a pathway-based web- accessible knowledgebase. Mol Endocrinol. 2010;24(9):1863–71.

    PubMed  PubMed Central  Google Scholar 

  21. 21.

    Aridon P, Geraci F, Turturici G, D’amelio M, Savettieri G, Sconzo G. Protective role of heat shock proteins in Parkinson’s disease. Neurodegener Dis. 2011;8(4):155–68.

    PubMed  Google Scholar 

  22. 22.

    Megat Abdul Wahab R, Yamamoto Z, Sintian A, Abu Kasim N, Zainol Abidin IZ, Senafi S, et al. The effects of orthodontic forces during canine retraction using self-ligating brackets on gingival crevicular fluid enzyme activity, canine movement and root resorption. Sains Malaysiana. 2015;44(2):249–56.

    Google Scholar 

  23. 23.

    Okano T, Sur J. Radiation dose and protection in dentistry. Japanese Dent Sci Rev. 2010;46:112–21.

    Google Scholar 

  24. 24.

    Tsolakis AI, Kalavritinos M, Bitsanis E, Sanoudos M, Benetou V, Alexiou K, et al. Reliability of different radiographic methods for the localization of displaced maxillary canines. Am J Orthod Dentofac Orthop. 2018;153(2):308–14.

    Google Scholar 

  25. 25.

    Zainal Ariffin SH, Rus Din RD, Yamamoto Z, Mohamad Jaafar I, Senafi S, Megat Abdul Wahab R. Penyerapan akar gigi apeks luaran hasil rawatan ortodontik pada enam dan 12 bulan. Sains Malaysiana. 2017;46(8):1299–307.

    Google Scholar 

  26. 26.

    Wang Y, Min HY, Chen Z, Xie X, Qin C, Zheng Y, et al. Study the predisposing factors of root resorption during orthodontic uprighting treatment of impacted mandibular third molars. Int Orthod. 2019;17(2):249–55.

    PubMed  Google Scholar 

  27. 27.

    Nanekrungsan K, Patanaporn V, Janhom A, Korwanich N. External apical root resorption in maxillary incisors in orthodontic patients: associated factors and radiographic evaluation. Imaging Sci Dent. 2012;42(3):147–54.

    PubMed  PubMed Central  Google Scholar 

  28. 28.

    Luke R, Khan SN, Iqbal PS, Soman RR, Chakkarayan J, Krishnan V. Estimation of specific salivary enzymatic biomarkers in individuals with gingivitis and chronic periodontitis: a clinical and biochemical study. J Int Oral Health. 2015;7(9):54–7.

    PubMed  PubMed Central  Google Scholar 

  29. 29.

    Wen X, Franchi L, Chen F, Gu Y. Proteomic analysis of gingival crevicular fluid for novel biomarkers of pubertal growth peak. Eur J Orthod. 2018;40(4):414–22.

    PubMed  Google Scholar 

  30. 30.

    Tsuchida S, Satoh M, Takiwaki M, Nomura F. Current status of proteomic technologies for discovering and identifying gingival crevicular fluid biomarkers for periodontal disease. Int J Mol Sci. 2019;20(86):1–11.

    Google Scholar 

  31. 31.

    Carneiro LG, Venuleo C, Oppenheim FG, Salih E. Proteome data set of human gingival crevicular fluid from healthy periodontium sites by multidimensional protein separation and mass spectrometry. J Periodontal Res. 2012;47(2):248–62.

    PubMed  Google Scholar 

  32. 32.

    Rody WJ, Holliday LS, McHugh KP, Wallet SM, Spicer V, Krokhin O. Mass spectrometry analysis of gingival crevicular fluid in the presence of external root resorption. Am J Orthod Dentofac Orthop. 2014;145(6):787–98.

    Google Scholar 

  33. 33.

    Papagerakis P, Zheng L, Kim D, Said R, Ehlert AA, Chung KKM, et al. Chapter 41. Odontogenesis: methods and protocols, methods in molecular biology. Methods Mol Biol. 2019;1922:549–62.

    PubMed  Google Scholar 

  34. 34.

    Lamster IB, Ahlo JK. Analysis of gingival crevicular fluid as applied to the diagnosis of oral and systemic diseases. Ann N Y Acad Sci. 2007;1098:216–29.

    PubMed  Google Scholar 

  35. 35.

    Wahab RA, Kasim NA, Senafi S, Jemain AA, Abidin IZ, Shahidan MA, et al. Enzyme activity profiles and ELISA analysis of biomarkers from human saliva and gingival crevicular fluid during orthodontic tooth movement using self-ligating brackets. Oral Health Dent Manag. 2014;13(2):194–9.

    PubMed  Google Scholar 

  36. 36.

    Bostanci N, Heywood W, Mills K, Parkar M, Nibali L, Donos N. Application of label-free absolute quantitative proteomics in human gingival crevicular fluid by LC/MS E (gingival exudatome). J Proteome Res. 2010;9(5):2191–9 Available from:

    PubMed  Google Scholar 

  37. 37.

    Grant M, Wilson J, Rock P, Chapple I. Induction of cytokines, MMP9, TIMPs, RANKL and OPG during orthodontic tooth movement. Eur J Orthod. 2013;35(5):644–51.

    PubMed  Google Scholar 

  38. 38.

    Ghallab NA. Diagnostic potential and future directions of biomarkers in gingival crevicular fl uid and saliva of periodontal diseases: review of the current evidence. Arch Oral Biol. 2018;87:115–24.

    PubMed  Google Scholar 

  39. 39.

    Karsani SA, Saihen NA, Zain RB, Cheong S, Abdul RM. Comparative proteomics analysis of oral cancer cell lines : identification of cancer associated proteins. Proteome Sci. 2014;12(3):1–13.

    Google Scholar 

  40. 40.

    Lokamani I, Looi ML, Md Ali SA, Mohd Dali AZH, Ahmad Annuar MA, Jamal R. Gelsolin and ceruloplasmin as potential predictive biomarkers for cervical cancer by 2D-DIGE proteomics analysis. Pathol Oncol Res. 2014;20(1):119–29.

    PubMed  Google Scholar 

  41. 41.

    Yap TW, Rabu A, Abu Bakar FD, Abdul Rahim R, Mahadi NM, Illias RM, et al. Growth phase-dependent proteomes of the Malaysian isolated Lactococcus lactis dairy strain M4 using label-free qualitative shotgun proteomics analysis. Sci World J. 2014;2014:1–14.

    Google Scholar 

  42. 42.

    Simard J-C, Simon M-M, Tessier PA, Girard D. Damage-associated molecular pattern S100A9 increases bactericidal activity of human neutrophils by enhancing phagocytosis. J Immunol. 2011;186(6):3622–31.

    PubMed  Google Scholar 

  43. 43.

    Zreiqat H, Howlett CR, Gronthos S, Hume D, Geczy CL. S100A8/S100A9 and their association with cartilage and bone. J Mol Histol. 2007;38(5):381–91.

    PubMed  Google Scholar 

  44. 44.

    Ellias MF, Zainal Ariffin SH, Karsani SA, Abdul Rahman M, Senafi S, Megat Abdul Wahab R. Proteomic analysis of saliva identifies potential biomarkers for orthodontic tooth movement. Sci World J. 2012;2012:1–6.

    Google Scholar 

  45. 45.

    Max BYEE, Korsmeyer SJ, Human J. Chain gene. Structure and expression in B lymphoid cells. J Exp Med. 1985;161(4):832–49.

    PubMed  Google Scholar 

  46. 46.

    Taha AA. The correlation between root resorption and some immune parameters in well-controlled type Ι diabetic patients during orthodontic treatment. Tikrit J Dent Sci. 2015;3(1):95–104.

    Google Scholar 

  47. 47.

    Chan HH, Rahim ZH, Jessie K, Hashim OH, Taiyeb-Ali TB. Salivary proteins associated with periodontitis in patients with type 2 diabetes mellitus. Int J Mol Sci. 2012;13(4):4642–54.

    PubMed  PubMed Central  Google Scholar 

  48. 48.

    Ariffin SH, Faiz Ellias M, Wahab RM, Bakar Y, Senafi S. Profil aktiviti laktat dehidrogenase, asid fosfatase rintang tartarat dan alkalin fosfatase pada air liur semasa rawatan ortodontik. Sains Malaysiana. 2010;39(3):405–12.

    Google Scholar 

  49. 49.

    Usuki K, Saras J, Waltenberger J, Miyazono K, Pierce G, Thomason A, et al. Platelet-derived endothelial cell growth factor has thymidine phosphorylase activity. Biochem Biophys Res Commun. 1992;184(3):1311–6.

    PubMed  Google Scholar 

  50. 50.

    Yoshino H, Morita I, Murota SI, Ishikawa I. Mechanical stress induces production of angiogenic regulators in cultured human gingival and periodontal ligament fibroblasts. J Periodontal Res. 2003;38(4):405–10.

    PubMed  Google Scholar 

  51. 51.

    Atsawasuwan P, Shirazi S. Advances in orthodontic tooth movement: gene therapy and molecular biology aspect. Curr Approaches Orthod. 2018;1:51–72.

    Google Scholar 

  52. 52.

    Pettiette MT, Zhang S, Moretti AJ, Kim SJ, Naqvi AR, Nares S. MicroRNA expression profiles in external cervical Resorption. J Endod. 2019;45(9):1106–1113.e2.

    PubMed  PubMed Central  Google Scholar 

Download references


We would like to thank all the patients involved in this study. We also thank the Mass Spectrometry Technology Section, the Malaysia Genome Institute, and the National Institutes of Biotechnology Malaysia (NIBM) for the LC-MS/MS service. The authors declare that they have no conflicts of interest.


This study was funded by the Universiti Kebangsaan Malaysia grant GUP-2017-002.

Author information




FAMN contributed to conception, design, data acquisition, and interpretation, as well as the drafting and critical revision of the manuscript. SHZA contributed to the conception and design of the manuscript, as well as the critical revision of the manuscript. SAK contributed to conception, design, data interpretation, as well as the critical revision of the manuscript. RMAW contributed to conception and design, as well as critical revision of the manuscript. All authors had read and gave their final approval of the manuscript for submission.

Corresponding author

Correspondence to Rohaya Megat Abdul Wahab.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Research Committee, Universiti Kebangsaan Malaysia (UKM PPI/111/8/JEP-2018-268), and patients were briefed and written informed consent was obtained.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Mohd Nasri, F.A., Zainal Ariffin, S.H., Karsani, S.A. et al. Label-free quantitative proteomic analysis of gingival crevicular fluid to identify potential early markers for root resorption. BMC Oral Health 20, 256 (2020).

Download citation


  • Root resorption
  • Early biomarker
  • Gingival crevicular fluid
  • Proteomics
  • Orthodontic treatment