Abstract

Full-text study online at
https://spj.science.org/doi/10.34133/csbj.0003

Moderate dental fluorosis (MF) represents a pivotal yet understudied stage of enamel fluorosis, characterized by distinct opacities indicating substantial but incomplete enamel disruption. Despite its prevalence, the molecular basis underlying this transitional phenotype remains unclear. This study employed a system proteomic approach to delineate salivary proteomic alterations associated with MF and to elucidate their biological and mechanistic implications. Unstimulated whole saliva was collected from school-aged children with MF (n = 10) and age-matched controls without fluorosis (CF, n = 21) for comparative proteomic analysis using liquid chromatography–tandem mass spectrometry. Urinary and drinking water fluoride levels were assessed to confirm comparable exposure between groups. Differentially expressed proteins (DEPs) were identified and subjected to Gene Ontology enrichment analysis, pathway analysis, and protein–protein interaction network analysis. Among 101 shared salivary proteins, 12 DEPs were significantly different between MF and CF subjects. Up-regulated proteins (e.g., neutrophil defensin 3, protein LEG1 homolog, immunoglobulin kappa constant, pyruvate kinase, mucin-7, and alpha-enolase) converged on pathways related to immune activation and altered glycolytic metabolism. Conversely, down-regulated proteins (e.g., superoxide dismutase [Cu-Zn], neutrophil gelatinase-associated lipocalin, thymidine phosphorylase, metalloproteinase inhibitor 1, galectin-3-binding protein, and alpha-1B-glycoprotein) reflected compromised antioxidant defense and perturbed extracellular matrix remodeling. Network topology analysis revealed distinct MF-specific interactomes dominated by immune–epithelial clusters. Notably, the cystic fibrosis transmembrane conductance regulator emerged as a central hub connecting multiple DEPs and MF-exclusive proteins, suggesting that fluoride-mediated ion transport perturbation may underlie enamel hypomineralization and broader oral homeostatic imbalance. Collectively, these findings provide novel molecular insights into MF pathogenesis and establish salivary proteomic profiles as a promising noninvasive platform for biomarker discovery, disease monitoring, and understanding of fluorosis progression.

Introduction

… Advances in high-throughput proteomics now allow the systematic characterization of saliva, transforming it into a powerful molecular window for disease research [9,10]. Saliva is a particularly attractive biofluid: it is noninvasive to collect, contains proteins secreted directly from salivary glands and gingival crevicular fluid, and reflects early immune and metabolic responses [11]. Importantly, altered salivary protein profiles have already been linked to several oral diseases, including dental caries, oral cancers, periodontitis, amelogenesis imperfecta, and molar–incisor hypomineralization, underscoring their utility in capturing disease-specific molecular signatures [1214].

Investigating the salivary proteome in MF offers a unique opportunity to identify the early molecular signatures of fluoride toxicity. This approach promises to enhance our understanding of critical fluoride exposure thresholds and guide targeted strategies for prevention and early intervention, thereby mitigating the risk of irreversible enamel damage. Building on this rationale, our study focused on school-aged children from fluoride-endemic regions of Lamphun and Ratchaburi provinces, Thailand, where natural groundwater contamination poses a significant public health challenge. By employing comparative salivary proteomics, we aimed to uncover differentially expressed proteins and potential molecular biomarkers that reflect the specific biological responses to chronic fluoride exposure characteristic of MF. To robustly assess exposure, urinary fluoride concentrations were also quantified. This comprehensive approach was designed to provide novel molecular insights into the pathogenesis of dental fluorosis, bridging a critical knowledge gap between healthy enamel, moderate pathology, and severe disease stages.

Materials and Methods

Study design and participant recruitment

The cross-sectional study was conducted from December 2021 to January 2022 in fluoridated water regions of Lamphun and Ratchaburi provinces, Thailand. Participants, aged 6 to 16 years, were systematically recruited through collaboration with local schools and public health centers. Before enrollment, the study objectives, procedures, and potential benefits were described to children and their legal guardians during organized briefing sessions. Written informed consent from guardians and assent from the participants were subsequently obtained.

Participants were categorized into 2 groups based on their clinical dental fluorosis status: those with MF and those without fluorosis (control fluorosis [CF] group). Dental examinations were performed using the TF index, with scores of 3 and 4 designating the MF group and a score of 0 defining the CF group. Two calibrated evaluators independently performed all assessments, and a consensus diagnostic score was recorded on a standardized electronic form. The reliability of scoring was confirmed by calculating the intraclass correlation coefficient, which demonstrated excellent agreement: 0.989 for interobserver, 0.999 for intraobserver 1, and 1.000 for intraobserver 2. To ensure data validity, individuals were excluded if they presented with systemic disease, active caries, amelogenesis imperfecta, or any signs of oral mucosal lesions or periodontitis or if they had used antibiotics, antifungals, antivirals, corticosteroids, or mouthwash within the previous 3 months. Additional exclusion criteria include smoking, alcohol use, or narcotic use.

… Urine sample collection and fluoride measurement

Each participant provided a 24-h urine sample using 2.7-l polyethylene containers. The first morning void was discarded, and all subsequent urine in the next 24 h was collected into a wide-necked container and later transferred into screw-cap bottles. Samples were maintained at 4 °C during transport. Fluoride concentrations were measured using a fluoride-ion-selective electrode (Orion 4-Star Benchtop, Thermo Fisher, USA) in combination with Total Ionic Strength Adjustment Buffer III (TISAB III buffer) to stabilize pH and prevent ion interference. All analyses were conducted at the Intercountry Centre for Oral Health in Chiang Mai, Thailand.

Drinking water collection and analysis

Participants were instructed to collect samples of their most frequently consumed household drinking water in clean, labeled 60-ml containers (filling to approximately five-sixths capacity). Samples were transported under controlled conditions to the Intercountry Centre for Oral Health, Chiang Mai, Thailand, where fluoride concentrations were determined using a fluoride-ion-selective electrode in combination with TISAB III, following the same protocol used for urine fluoride quantification.

Saliva sample collection and processing

To minimize confounding variables, participants adhered to a strict regimen prior to specimen collection: an overnight fast of at least 12 h and abstinence from oral hygiene practices, food, and beverages on the morning of sampling. Unstimulated whole saliva was collected between 8:00 and 10:00 AM. Following a 30-s rinse with 10 ml of sterile water to remove oral debris and hydrate the mucosa, approximately 3 ml of saliva was passively drooled into sterile collection tubes. Samples were placed on ice immediately, rapidly transferred to dry ice for transport, and then stored long-term at ?80 °C until proteomic analysis.

…Statistical analysis

Demographic characteristics and fluoride concentrations in urine and water samples were summarized using descriptive statistics. Data normality was evaluated via the Shapiro–Wilk test. Comparisons of fluoride concentrations between the MF and CF groups were determined using the Mann–Whitney test. Differences in age and in up-regulated protein expression were analyzed with an unpaired t test with Welch’s correction. Sex distribution between groups was evaluated using the chi-square (Fisher’s exact) test. A 2-tailed P value <0.05 was considered indicative of statistical significance. All statistical analyses were conducted using GraphPad Prism v9.

Results

Participant characteristics and fluoride exposure

The comprehensive study workflow, encompassing participant recruitment, sample collection (drinking water, urine, and saliva), fluoride quantification, and salivary proteomic analysis, is summarized in Fig. 1A. Participants were stratified into 2 groups based on their TF fluorosis scores: the CF group (n = 21, TF score 0) and the MF group (n = 10, TF scores 3 and 4). Representative intraoral photographs of participants diagnosed in the MF group revealed enamel hypomineralization characterized by opaque white striations, mottled surface patterns, and diffuse opacities (Fig. 1B to D). Conversely, the CF group showed smooth, uniformly translucent enamel surfaces without visible defects (Fig. 1E to G)….

Fig. 1. Dental characteristics of study participants and fluoride concentration in urine and drinking water. (A) Schematic summary of the study design. (B to D) Moderate dental fluorosis (MF) representing white opacities, yellow-brown speckles, surface irregularities, and pitting. (E to G) Control participants (control fluorosis [CF]) showing smooth, glossy, creamy-white enamel. (H and I) Fluoride concentrations (ppm) in urine (yellow) and water (blue) for individual participants in the MF (H) and CF (I) groups. Yellow plots indicate urinary fluoride levels, and blue overlays depict the corresponding water fluoride levels of each participant. (Note: Due to the technical limitations of data acquisition, the water fluoride concentration for participant MF06 was undetectable, as were both the water and urinary fluoride concentrations for participants CF20 and CF21.)

Demographically, the participant ages were comparable between the CF (mean 10.74 ± 2.40 years) and MF (11.40 ± 0.97 years) groups (P = 0.302). However, a statistically significant difference was observed in sex distribution between the CF (male = 12; female = 9) and MF (male = 1; female = 9) groups (P = 0.02). Fluoride concentrations in drinking water and urine for each participant are illustrated in Fig. 1H and I. In the MF group, the mean urinary fluoride concentration was 2.56 ± 1.21 ppm (range: 0.43 to 4.01 ppm), while drinking water fluoride levels averaged 0.19 ± 0.06 ppm (range: <0.10 to 0.30 ppm). For the CF group, the mean urinary fluoride was 2.29 ± 1.60 ppm (range: 0.56 to 5.86 ppm), and water fluoride averaged 1.15 ± 2.07 ppm (range: <0.10 to 5.6 ppm). No statistically significant differences were detected between the groups for either urinary or drinking water fluoride concentrations. Notably, several CF participants exhibited transient elevations in urinary fluoride despite consistently low fluoride levels in their drinking water, suggesting the influence of additional or alternative exposure sources beyond the local water supply.

Differential salivary proteomic signatures and functional enrichment in MF

Quantitative proteomic profiling identified a total of 166 proteins in saliva. After excluding contaminants and poorly labeled proteins, a high-confidence dataset of 161 proteins remained. To minimize the plausible bias arising from sporadic detection, proteins present in fewer than 3 samples per group were removed, resulting in a refined dataset of 154 proteins. Of these, 101 proteins were common to both groups, whereas 37 were uniquely detected in the MF group and 16 in the CF group (Fig. 2A)…

Among the 101 shared proteins, 12 proteins exhibited statistically significant differences in abundance between groups: 6 were up-regulated and 6 were down-regulated in the MF group. The proteins found to be elevated in the MF group included neutrophil defensin 3 (DEFA3), protein LEG1 homolog (LEG1), immunoglobulin kappa constant (IGKC), pyruvate kinase (PKM), mucin-7 (MUC7), and alpha-enolase (ENO1). In contrast, the proteins with reduced abundance in the MF group were superoxide dismutase [Cu-Zn] (SOD1), neutrophil gelatinase-associated lipocalin (LCN2), thymidine phosphorylase (TYMP), metalloproteinase inhibitor 1 (TIMP1), galectin-3-binding protein (LGALS3BP), and alpha-1B-glycoprotein (A1BG) (Fig. 2B and Table S1). The key biological functions of these differentially expressed proteins, as annotated from the DAVID database, are listed in Table S1.

… Discussion

Dental fluorosis results from cumulative fluoride exposure through water, food, and oral care products, disrupting amelogenesis and leading to MF [7]. As a transitional stage marking key molecular and structural changes in enamel, MF provides a window into disease progression. This study identified salivary proteomic signatures that not only mirror enamel mineralization defects but also reveal systemic alterations in immunity, metabolism, and ECM regulation, offering potential biomarkers for early detection and prevention.

Interestingly, we observed comparable urinary and drinking water fluoride levels between our MF and control groups, echoing previous reports that conventional fluoride biomarkers do not always correlate directly with clinical severity [15]. This highlights that the biological impact of fluoride extends beyond a single-point concentration measurement, depending instead on cumulative intake from multiple sources and the host’s intrinsic ability to regulate fluoride metabolism [16]. Urinary fluoride was used as a biomarker of systemic exposure, reflecting cumulative fluoride exposure inferred rather than comprehensively quantified across all potential sources. Our findings thus reinforce the concept that salivary proteomics may provide a more dynamic and functionally relevant readout of biological responses to fluoride exposure than traditional methods.

Among the 12 significantly altered proteins, the up-regulation of DEFA3, MUC7, IGKC, LEG1, PKM, and ENO1 in the MF group collectively points to heightened immune activity and metabolic reprogramming. DEFA3, an antimicrobial peptide derived from neutrophils, suggests a compensatory defense mechanism aimed at protecting structurally compromised enamel from microbial adhesion and colonization [17]. The elevation of MUC7, being a major salivary mucin, reflects altered mucin composition, potentially contributing to chemical and antimicrobial defense [18]. Similarly, the increased abundance of IGKC indicates activation of adaptive immune pathways, consistent with its previously reported association with immune dysregulation in oral malignant disorders [19]. LEG1 was also elevated in the MF group. Although its role in oral biology has not been extensively characterized, LEG1 is a secreted protein expressed in epithelial and secretory tissues and has been implicated in mucosal defense and innate immune regulation [20,21]. Its detectability in saliva is therefore biologically plausible, and its increased abundance in MF may reflect broader mucosal or immune adaptations to fluoride-associated environmental stress rather than a fluorosis-specific effect. Metabolically, the rise in PKM and ENO1 underscores glycolytic reprogramming under fluoride stress. PKM, a pivotal regulator of glycolytic metabolism, influences immune cell infiltration, inflammation, and tissue remodeling, all of which are implicated in periodontal and neoplastic conditions [22,23]. The elevated ENO1 abundance may reflect a compensatory mechanism in response to fluoride’s competitive inhibition of ENO1, which can disrupt glycolytic flux and induce inflammatory stress [24,25]. Together, these findings align with experimental evidence that the up-regulation of glycolytic proteins can alter energy metabolism in ameloblasts and exacerbate fluorosis development in animal models [26]. Thus, the salivary proteome of the MF group reveals a pattern of immune priming and metabolic adjustment, reflecting both local enamel vulnerability and systemic biological responses to cumulative fluoride exposure.

In contrast to the up-regulated immune and metabolic proteins, a core feature of the MF group was the down-regulation of antioxidant and ECM-regulating proteins in saliva. The notable reduction of SOD1, an essential antioxidant enzyme expressed in the periodontal ligament and gingival fibroblasts, indicates a diminished capacity to neutralize superoxide anions [27]. Such depletion may heighten susceptibility to reactive oxygen species, disrupt mitochondrial redox homeostasis, and impair cell function [28], thereby exacerbating enamel pathology through compromised ameloblast activity [29]. Cytoscape-based PPI analysis further revealed SOD1 clustering with several MF-unique proteins, including GANAB, PFN1, and YWHAZ, which are associated with protein folding [30,31], cell migration [32,33], and cell proliferation and apoptosis [34,35]. Interestingly, SOD1 also showed marked interactions with several key mediators of tooth development, including mitogen-activated protein kinase (MAPK) signaling mediators (MAP3K7 and SMAD2), the anti-apoptotic regulator (BCL2 apoptosis regulator [BCL2]), and cell proliferation markers (minichromosome maintenance complex component 4 [MCM4] and minichromosome maintenance complex component 6 [MCM6]) [36,37]. This network also encompassed ubiquitin-like modifier-activating enzyme (UBA1; a ubiquitination enzyme) and isoform 3 of L-lactate dehydrogenase A chain (LDHA; a glycolytic enzyme), both previously reported to be up-regulated in severe fluorosis [10], suggesting that redox imbalance and metabolic adaptation in the MF group may converge with proteostasis and cell survival pathways that become more pronounced in advanced stages of the disease.

Similarly, the marked reduction in TIMP1, a key inhibitor of matrix metalloproteinases (MMPs), points to disturbed ECM turnover. Diminished TIMP1 may tip the balance toward unchecked MMP activity, which could weaken enamel-supporting structures and compromise tissue stability [38]. Network mapping linked TIMP1 to ECM remodeling proteins (MMP9, MMP10, MMP14, and collagen type V alpha 1 chain [COL5A1] [39]), as well as to growth factor receptors (fibroblast growth factor receptor 2 [FGFR2] [40] and Erb-B2 receptor tyrosine kinase 2 [ERBB2] [41]) and the cell–cell aggregation and cohesion mediator (kallikrein related peptidase 5 [KLK5] [42,43]). These associations implicate TIMP1 in broader processes of ECM remodeling, growth factor signaling, and cell–cell adhesion, all of which are highly vulnerable to fluoride-induced dysregulation in the oral environment.

The PPI network of several proteins exclusive to the MF cohort revealed reinforced interaction profiles, which may be attributable to chronic fluoride exposure. Beyond their associations with the CFTR channel, proteins such as YWHAZ and HSPA8 were notably linked to proteins previously documented as up-regulated in severe dental fluorosis [10]. This observation is further supported by the network’s inclusion of other proteins vital for tooth development and maintenance of oral health [44]. These included the stress-inducible chaperone heat shock protein 90 alpha family class A member 1 (HSP90AA1) [45], which is crucial for cell proteostasis, alongside key players in mineralization and ameloblast differentiation such as collagen type I alpha 1 chain (COL1A1) [46], estrogen receptor 1 (ESR1), estrogen receptor 2 (ESR2) [47], family with sequence similarity 20 C (FAM20C) [48], family with sequence similarity 83 member H (FAM83H) [49], signal transducers and activators of transcription 3 (STAT3) [50], and tryptophan-aspartate repeat domain 72 (WDR72) [51]. Notably, MF-exclusive salivary proteins converged within interaction networks centered on FAM20C, FAM83H, and COL1A1, which themselves interface with key enamel- and matrix-associated regulators, including family with sequence similarity 20, member A (FAM20A); enamelin (ENAM); and integrin beta-6 (ITGB6), implicated in enamel mineralization, ameloblast-related pathways, and ECM organization [52,53]. FAM20C, FAM20A, and FAM83H are involved in enamel matrix phosphorylation and maturation, whereas COL1A1 supports the structural integrity of mineralized tissues [54]. Given that the study cohort comprised children aged 6 to 16 years, encompassing developmental stages during which enamel formation of different tooth types may be ongoing or already complete, these salivary proteomic signatures are not interpreted as markers of active amelogenesis. Instead, the integration of enamel- and matrix-associated regulators within the MF-related PPI network suggests that chronic fluoride exposure may induce persistent molecular imprints linked to disrupted enamel biology, which remain detectable in saliva across developmental stages. Accordingly, salivary proteomics is employed as a noninvasive, hypothesis-generating approach, with careful consideration of its temporal and biological scope. These signals likely reflect a combination of enduring regulatory adaptations arising from earlier developmental exposure and secondary responses associated with altered enamel structure and its ongoing interaction with the oral environment, supporting the use of saliva as a noninvasive molecular window into fluorosis-associated biology [55,56]. Collectively, the dense functional interconnections within the MF-exclusive proteome underscore coordinated biological reprogramming processes that span oxidative stress, ECM integrity, and ameloblast function, distinguishing MF from mild forms and hinting at a trajectory toward more advanced pathology.

Network analysis provided further insight by revealing CFTR as a prominent bridging node connecting multiple functional modules within the inferred PPI network. Importantly, CFTR expression is highly cell type specific and temporally regulated, particularly during developmental processes such as amelogenesis [57,58]. Previous studies have shown that CFTR activity is enriched during discrete stages of enamel formation, most notably in maturation-stage ameloblasts, rather than being constitutively or continuously expressed in extracellular secretions such as saliva [59]. Accordingly, the absence of CFTR as a directly detected salivary protein in the present proteomic dataset does not preclude its functional or regulatory relevance within fluoride-associated molecular networks. Network analysis provided additional insight by identifying CFTR as a prominent bridging node connecting multiple functional modules within the inferred PPI network. Although CFTR was not experimentally detected as a differentially expressed salivary protein, its apparent centrality reflects inferred functional connectivity based on curated interaction databases. Within this network, CFTR was linked to proteins involved in glucose metabolism, immune modulation, protein folding, structural organization, protease regulation, and key cellular processes such as homeostasis, adhesion–migration, and proliferation–differentiation [19,24,25,3035]. Together, these network-based associations highlight a context-dependent, integrative role for CFTR in coordinating fluoride-associated biological responses, potentially through developmentally restricted ion transport and pH-regulatory mechanisms rather than via direct salivary expression.

This interpretation is supported by prior studies reporting functional correlations between CFTR and up-regulated salivary proteins in severe dental fluorosis (histone H4 [H4C3], plastin-2 [LCP1], LDHA, S100 calcium binding protein A9 [S100A9], and UBA1), underscoring its possible involvement in immune, metabolic, and survival pathways under fluoride stress [10]. Moreover, urinary proteomic studies in children with high fluoride exposure have identified systemic alterations in ECM, biomineralization, oxidative stress, and immune pathways, with several proteins linked to CFTR. The emergence of CFTR as a network-inferred integrative node suggests that CFTR-mediated ion transport and pH regulation may represent a shared molecular axis linking systemic fluoride exposure to localized salivary and oral responses across different exposure severities [60]. Mechanistically, CFTR is a gated anion channel localized not only at the apical plasma membrane of epithelial cells, including ameloblasts, but also within intracellular organelles such as mitochondria and lysosomes [61]. Beyond its canonical role in chloride transport, CFTR orchestrates epithelial and intracellular pH regulation and ion homeostasis by mediating the exchange of glutathione, halides, and polyatomic anions [61]. Disruption of CFTR function has been shown to impair secretory processes and intracellular signaling, particularly pathways related to chemokine production and glucose metabolism [62,63], and has been consistently linked to incomplete enamel mineralization across porcine, rat, mouse, and human models [57].

Since fluoride can be transported through CFTR, excessive exposure may competitively interfere with anion exchange, especially during the maturation stage of amelogenesis, thereby potentially compromising enamel deposition and contributing to enamel hypomineralization [57]. Taken together, the CFTR-centered interaction patterns observed in MF may reflect an early stage of molecular dysregulation that parallels, and may precede, the proteomic disturbances documented in severe disease. In this context, MF could represent a transitional state in which emerging CFTR-linked alterations foreshadow the broader immune, metabolic, and structural disruptions characteristic of advanced pathology.

Although a statistically significant imbalance in sex distribution was observed between the CF and MF groups, available epidemiological evidence indicates that biological sex is not a dominant determinant of dental fluorosis prevalence or severity in pediatric populations [6466]. While biological sex may influence certain immune and salivary parameters, enamel formation and amelogenesis during childhood are governed by conserved developmental processes that are not regarded as strongly sex dependent [67]. Accordingly, sex was not incorporated as a covariate in the proteomic analysis to avoid overparameterization in a high-dimensional dataset. This imbalance is therefore acknowledged as a methodological consideration, and future studies employing sex-balanced cohorts will be important to further delineate potential sex-specific contributions to fluoride-associated salivary proteomic alterations.

While strict inclusion criteria and standardized sampling procedures were applied to reduce biological and technical variability, this study was inherently exploratory and hypothesis generating, and the cohort size was necessarily limited. In addition, cumulative fluoride exposure was inferred rather than comprehensively quantified across all potential sources, with urinary fluoride serving as a biomarker of systemic exposure. Accordingly, the identified proteomic alterations should be interpreted as preliminary molecular signals associated with MF, reflecting potential inter-individual variability in biological responses, and warrant validation in larger, independent cohorts [60,68]. Importantly, this study establishes a foundational proteomic framework for MF while opening promising avenues for future research. Beyond protein abundance, incorporating posttranslational modifications and metabolomic data could uncover additional fluoride-related mechanisms. Longitudinal studies are crucial to clarify how these molecular patterns evolve during enamel maturation, and broader exposure assessments, integrating dietary and genetic factors, could refine susceptibility profiles. Together, these approaches would strengthen the translational potential of salivary proteomics as a noninvasive tool for early diagnosis and targeted prevention of fluorosis.

In conclusion, this study demonstrates that MF represents a biologically distinct stage marked by coordinated alterations in immunity, metabolism, antioxidant defense, and ECM regulation (Fig. 5). The up-regulation of DEFA3, MUC7, IGKC, PKM, and ENO1 suggests intensified defense and glycolytic reprogramming, while the down-regulation of SOD1 and TIMP1 reflects increased susceptibility to oxidative stress and defective ECM remodeling. CFTR emerged as a network-inferred central hub linking immune, metabolic, and structural pathways, consistent with its role in enamel mineralization and ion homeostasis. These proteomic shifts, observed despite similar urinary and water fluoride levels between groups, highlight the role of cumulative fluoride exposure from multiple sources. Importantly, the molecular alterations identified in MF may foreshadow the transition toward severe disease, underscoring the value of salivary proteomics as a noninvasive tool for early detection, monitoring, and prevention of this prevalent condition….

Fig. 5. Summary schematic demonstrating the key findings of the study. (A) Clinical relevance: The findings suggest biological transition between clinical stages. Salivary proteomics offers a noninvasive platform for early detection and monitoring of fluorosis progression in response to cumulative fluoride exposure. ROS, reactive oxygen species.(B) Molecular signature: Salivary proteomics revealed a distinct profile in moderate fluorosis, by the up-regulation of immune (DEFA3, MUC7, and IGKC) and metabolic (PKM and ENO1) proteins and down-regulation of antioxidant (SOD1) and extracellular-matrix (ECM)-related (TIMP1) proteins. The cystic fibrosis transmembrane conductance regulator (CFTR) acted as a central hub coordinating these changes.

Ethical Approval

Written informed consent was obtained from all participants involved in the study.

Acknowledgments

Funding: This research was supported by the Second Century Fund (C2F), Chulalongkorn University, Thailand. K.S. was supported by the Faculty of Dentistry (DRF 68_006); the Ratchadaphiseksomphot Endowment Fund, Chulalongkorn University (The Exchange Faculty Travel Grant; Grant No. CTG168039); and the Thailand Science Research and Innovation Fund Chulalongkorn University (HEA_FF_69_051_3200_004). T.P. was supported by the Health Systems Research Institute (69-143 and 69-145), Thailand Science Research and Innovation Fund Chulalongkorn University (HEA_FF_69_036_3200_003), Chulalongkorn University. We acknowledge the Second Century Fund (C2F), Chulalongkorn University, for support of the Future Dent Digital Center.

Author contributions: A.P. contributed to data analysis, interpretation, and manuscript preparation. P.G. and T.K. conducted the experimental work and contributed to methodology development. H.-S.J. provided scientific advice and technical expertise. K.S. and T.P. conceptualized and supervised the study, oversaw project administration, acquired funding, and critically revised the manuscript. All authors have read and approved the final version of the manuscript.

Competing interests: The authors declare that they have no competing interests.
Data Availability

The datasets generated and/or analyzed during this study are available from the corresponding author upon reasonable request.

Supplementary Materials

Tables S1 to S11
References

1. Zohoori FV, Buzalaf M. Fluoride. Adv Nutr. 2022;13(6):2679–2680.
2. Zhang K, Lu Z, Guo X. Advances in epidemiological status and pathogenesis of dental fluorosis. Front Cell Dev Biol. 2023;11:1168215.
3. Kawasaki MPT, Kawasaki K, Oommen S, Otsuka-Tanaka Y, Hishinuma M, Nomoto T, Maeda T, Takubo K, Suda T, Sharpe PT, et al. R-spondins/Lgrs expression in tooth development. Dev Dyn. 2014;243(6):844–851.
4. Srivastava R, Tangade P, Priyadarshi S. Managing dental fluorosis: A guide to improving appearance of the teeth. Arch Dent Res. 2023;13(1):25–29.
5. Rosas López Portillo F, Rosas Alcaraz VM, Gámez Valenzuela J, Molina Ángulo IF, Gonzáles Jiménez MV, Olivas Velázquez AK, Reyna Solis CL, Rocha López A, Cortina Aguilar JA, Villalobos Rodelo JJ, et al. Dental fluorosis: A review. Int J Res Med Sci. 2024;12(4):1306–1310.
6. Samaranayake L, Porntaveetus T, Tsoi J, Tuygunov N. Facts and fallacies of the fluoride controversy: A contemporary perspective. Int Dent J. 2025;75(4): Article 100833.
7. EFSA Scientific Committee, Bennekou SH, Allende A, Bearth A, Casacuberta J, Castle L, Coja T, Crepet A, Hoogenboom R, Knutsen H, et al. Updated consumer risk assessment of fluoride in food and drinking water including the contribution from other sources of oral exposure. EFSA J. 2025;23(7): Article e9478.
8. Blackburn JKK, Porntaveetus T, Kawasaki M, Otsuka-Tanaka Y, Miake Y, Ota MS, Watanabe M, Hishinuma M, Nomoto T, Oommen S, et al. Excess NF-?B induces ectopic odontogenesis in embryonic incisor epithelium. J Dent Res J Dent Res. 2015;91(1):121–128.
9. Ajrithirong P, Krasaesin A, Sriarj W, Gavila P, Chetruengchai W, Sriwattanapong K, Manaspon C, Samaranayake L, Porntaveetus T. The metagenome and metabolome signatures of dental biofilms associated with severe dental fluorosis. J Oral Microbiol. 2025;17(1):2560591.
10. Gavila P, Ajrithirong P, Chumnanprai S, Kalpongnukul N, Pisitkun T, Chantarangsu S, Sriwattanapong K, Tagami J, Porntaveetus T. Salivary proteomic signatures in severe dental fluorosis. Sci Rep. 2024;14(1):18372.
11. Umapathy VR, Natarajan PM, Swamikannu B. Review insights on salivary proteomics biomarkers in oral cancer detection and diagnosis. Molecules. 2023;28(13):5283.
12. Nowwarote N, Theerapanon T, Osathanon T, Pavasant P, Porntaveetus T, Shotelersuk V. Amelogenesis imperfecta: A novel FAM83H mutation and characteristics of periodontal ligament cells. Oral Dis. 2018;24(8):1522–1531.
13. Porntaveetus T, Srichomthong C, Ohazama A, Suphapeetiporn K, Shotelersuk V. A novel GJA1 mutation in oculodentodigital dysplasia with extensive loss of enamel. Oral Dis. 2017;23(6):795–800.
14. Sriwattanapong K, Theerapanon T, Khamwachirapitak C, Sae-ear P, Srijunbarl A, Porntaveetus T, Shotelersuk V. Deep dental phenotyping and a novel FAM20A variant in patients with amelogenesis imperfecta type IG. Oral Dis. 2024;30(2):537–550.
15. Singhal R, Namdev R, Kumar A, Bhagol A, Supriya S. Correlation of fluoride intake with haemoglobin level and intelligence quotient in 8–12 year aged children: An observational study from India. BMC Public Health. 2025;25(1):788.
16. García Pérez A, Villanueva Gutiérrez T, González-Aragón Pineda AE, Murillo Santos KL, Pérez Pérez NG. Molar–incisor hypomineralization is associated with the prevalence of thinness among schoolchildren in communities with different fluoride levels in the drinking water. Int J Dent. 2024;2024:6212877.
17. Teutle-Coyotecatl B, Contreras-Bulnes R, Rodríguez-Vilchis LE, Scougall-Vilchis RJ, Velazquez-Enriquez U, Almaguer-Flores A, Arenas-Alatorre JA. Effect of surface roughness of deciduous and permanent tooth enamel on bacterial adhesion. Microorganisms. 2022;10(9):1701.
18. Frenkel ES, Ribbeck K. Salivary mucins in host defense and disease prevention. J Oral Microbiol. 2015;7:29759.
19. Kurkalang S, Roy S, Acharya A, Mazumder P, Mazumder S, Patra S, Ghosh S, Sarkar S, Kundu S, Biswas NK, et al. Single-cell transcriptomic analysis of gingivo-buccal oral cancer reveals two dominant cellular programs. Cancer Sci. 2023;114(12):4732–4746.
20. Foratori-Junior GAVT, Grizzo LT, Jesuino BG, Castilho AVSS, Buzalaf MAR, Sales-Peres SHC. Is there a difference in the proteomic profile of stimulated and unstimulated saliva samples from pregnant women with/without obesity and periodontitis? Cells. 2023;12(10):1389.
21. Dang YWJ, Liu C, Zhang K, Jinrong P, He J. Evolutionary and molecular characterization of liver-enriched gene 1. Sci Rep. 2020;10(1):4262.
22. Gao H. PKM2’s dual role in periodontitis: Regulating inflammation and bone metabolism imbalance. Oral Dis. 2025.
23. Xue S, Luo Z, Mao Y, Liu S. A comprehensive analysis of the pyruvate kinase M1/2 (PKM) in human cancer. Gene. 2025;937: Article 149155.
24. Ji H, Wang J, Guo J, Li Y, Lian S, Guo W, Yang H, Kong F, Zhen L, Guo L, et al. Progress in the biological function of alpha-enolase. Anim Nutr. 2016;2(1):12–17.
25. Strunecka A, Strunecky O. Mechanisms of fluoride toxicity: From enzymes to underlying integrative networks. Appl Sci. 2020;10(20):7100.
26. Ba Y, Yang S, Yu S, Hou X, Du Y, Gao M, Zuo J, Fu X, Li Z, Huang H, et al. Role of glycolysis/gluconeogenesis and HIF-1 signaling pathways in rats with dental fluorosis integrated proteomics and metabolomics analysis. Int J Mol Sci. 2022;23(15):8266.
27. Ramesh KSV, Swetha P, Mohan Kumar P, Sruthima NVS, Naresh C. Estimation of superoxide dismutase levels in saliva and gingival crevicular fluid among smokers and non-smokers in periodontitis patients—An observational study. Niger Med J. 2019;60(3):133–137.
28. Pal P, Jha NK, Pal D, Jha SK, Anand U, Gopalakrishnan AV, Dey A, Mukhopadhyay PK. Molecular basis of fluoride toxicities: Beyond benefits and implications in human disorders. Genes Dis. 2023;10(4):1470–1493.
29. Suzuki M, Bandoski C, Bartlett JD. Fluoride induces oxidative damage and SIRT1/autophagy through ROS-mediated JNK signaling. Free Radic Biol Med. 2015;89:369–378.
30. Chiu CC, Lin CY, Lee LY, Chen YJ, Lu YC, Wang HM, Liao CT, Chang JT, Cheng AJ. Molecular chaperones as a common set of proteins that regulate the invasion phenotype of head and neck cancer. Clin Cancer Res. 2011;17(14):4629–4641.
32. Kolokythas A, Haupt DB, Haupt AM, Adami G. Decreased expression levels of pfn1, tmsb4x and krt13 in advanced stage oral squamous cell carcinoma. Is there an association with disease progression? J Oral Maxillofac Surg. 2015;73(9): Article e1.
33. Adami GR, O’Callaghan TN, Kolokythas A, Cabay RJ, Zhou Y, Schwartz JL. A loss of profilin-1 in late-stage oral squamous cell carcinoma. J Oral Pathol Med. 2017;46(7):489–495.
34. Hu M, Song HY, Chen L. Quercetin acts via the G3BP1/YWHAZ axis to inhibit glycolysis and proliferation in oral squamous cell carcinoma. Toxicol Mech Methods. 2023;33(2):141–150.
35. Hu S, Chen C, Chen H, Yu X, Li X, Bai Y, Shao C. YWHAZ gene contributes to the radioresistance of oral squamous cell carcinoma cells. Radiat Med Prot. 2024;5(1):30–36.
36. Apellániz D, Pereira-Prado V, Sicco E, Vigil-Bastitta G, González-González R, Mosqueda-Taylor A, Molina-Frechero N, Hernandez M, Sánchez-Romero C, Bologna-Molina R. Comparative study of the minichromosome maintenance proteins complex (MCM 4/5/6) in ameloblastoma and unicystic ameloblastoma. Int J Surg Pathol. 2018;26(8):714–720.
37. Yamashita S, Okamoto M, Mendonca M, Fujiwara N, Kitamura E, Chang C-SS, Brueckner S, Shindo S, Kuriki N, Cooley MA, et al. Fluoride alters gene expression via histone H3K27 acetylation in ameloblast-like LS8 cells. Int J Mol Sci. 2024;25(17):9600.
38. Grünwald B, Schoeps B, Krüger A. Recognizing the molecular multifunctionality and interactome of TIMP-1. Trends Cell Biol. 2019;29(1):6–19.
39. Shaikh S, Dhar H, Moorthy M, Bhat V, Basu S, Banerjee D, Mishra DK, Datta S, Mukherjee G. The spatial distribution of intermediate fibroblasts and myeloid-derived cells dictate lymph node metastasis dynamics in oral cancer. J Transl Med. 2024;22(1):759.
40. Du W, Du W, Yu H. The role of fibroblast growth factors in tooth development and incisor renewal. Stem Cells Int. 2018;2018:7549160.
41. Meng Y, Yang P, Ma L. Prognostic and clinical implications of c-erbB-2 expression in patients with oral cancer: A meta-analysis. Medicine. 2020;99(23): Article e20575.
42. Jiang R, Shi Z, Johnson JJ, Liu Y, Stack MS. Kallikrein-5 promotes cleavage of desmoglein-1 and loss of cell-cell cohesion in oral squamous cell carcinoma. J Biol Chem. 2011;286(11):9127–9135.
43. Alves MG, Kodama MH, da Silva EZM, Gomes BBM, da Silva RAA, Vieira GV, Alves VM, da Fonseca CK, Santana AC, Cecílio NT, et al. Relative expression of KLK5 to LEKTI is associated with aggressiveness of oral squamous cell carcinoma. Transl Oncol. 2021;14(1): Article 100970.
44. Thaweesapphithak SSJ, Kamolvisit W, Theerapanon T, Porntaveetus T, Shotelersuk V. Cleidocranial dysplasia and novel RUNX2 variants: Dental, craniofacial, and osseous manifestations. J Appl Oral Sci. 2022;30: Article e20220028.
45. Zhang H, Huang J, Fan X, Miao R, Wang Y. HSP90AA1 promotes the inflammation in human gingival fibroblasts induced by Porphyromonas gingivalis lipopolysaccharide via regulating of autophagy. BMC Oral Health. 2022;22(1):366.
46. Budsamongkol TIN, Theerapanon T, Yodsanga S, Porntaveetus T, Shotelersuk V. A novel mutation in COL1A2 leads to osteogenesis imperfecta/Ehlers-Danlos overlap syndrome with brachydactyly. Genes Dis. 2019;6(2):138–146.
47. Madalena IR, Reis CLB, Matsumoto MAN, Stuani MBS, Mattos NHR, Oliveira DSB, de Menezes Oliveira MAH, Roskamp L, Kuchler EC, Baratto-Filho F. Investigating the association between dental age and polymorphisms in genes encoding estrogen receptors. J Appl Oral Sci. 2023;31: Article e20230184.
48. Zhang R, Ren Y, Ju Y, Zhang Y, Zhang Y, Wang Y. FAM20C: A key protein kinase in multiple diseases. Genes Dis. 2025;12(2): Article 101179.
49. Xin W, Wenjun W, Man Q, Yuming Z. Novel FAM83H mutations in patients with amelogenesis imperfecta. Sci Rep. 2017;7(1):6075.
50. Chan L, Lu J, Feng X, Lin L, Yao Y, Zhang X. Loss of Stat3 in Osterix+ cells impairs dental hard tissues development. Cell Biosci. 2023;13(1):75.
51. Katsura K, Nakano Y, Zhang Y, Shemirani R, Li W, Den Besten P. WDR72 regulates vesicle trafficking in ameloblasts. Sci Rep. 2022;12(1):2820.
52. Sriwattanapong KTS, Khamwachirapitak C, Sae-Ear P, Prommanee S, Sa-Ard-Iam N, Phothichailert S, Jung HS, Shotelersuk V, Porntaveetus T. FAM20A deficiency drives transcriptomic dysregulation and functional impairment in gingival fibroblasts. Cell Prolif. 2025;22: Article e70096.
53. Sriwattanapong KTT, Boonprakong L, Srijunbarl A, Porntaveetus T, Shotelersuk V. Novel ITGB6 variants cause hypoplastic-hypomineralized amelogenesis imperfecta and taurodontism: Characterization of tooth phenotype and review of literature. BDJ Open. 2023;9(1):15.
54. Sriwattanapong KTT, Khamwachirapitak C, Sae-Ear P, Sa-Ard-Iam N, Shotelersuk V, Porntaveetus T. In-depth investigation of FAM20A insufficiency effects on deciduous dental pulp cells: Altered behaviours, osteogenic differentiation, and inflammatory gene expression. Int Endod J. 2024;57(6):745–758.
55. Constantin VLI, Goriuc A, Budala DG, Bida FC, Cojocaru C, Butnaru O-M, Virvescu DI. Salivary biomarkers identification: Advances in standard and emerging technologies. Oral Dis. 2025;5(2):26.
56. Ye JZF, Luo Z, Ou X. Comparative salivary proteomics analysis of children with and without early childhood caries using the DIA approach: A pilot study. Proteomics Clin Appl. 2024;18(6): Article e2400006.
57. Bronckers AL, Lyaruu DM, Guo J, Bijvelds MJ, Bervoets TJ, Zandieh-Doulabi B, Medina JF, Li Z, Zhang Y, DenBesten P. Composition of mineralizing incisor enamel in cystic fibrosis transmembrane conductance regulator-deficient mice. Eur J Oral Sci. 2015;123(1):9–16.
58. Nitayavardhana ITT, Srichomthong C, Piwluang S, Wichadakul D, Porntaveetus T, Shotelersuk V. Four novel mutations of FAM20A in amelogenesis imperfecta type IG and review of literature for its genotype and phenotype spectra. Mol Gen Genomics. 2020;295(4):923–931.
59. Ghigo ADSC, Hart M, Mitash N, Swiatecka-Urban A. Cell signaling and regulation of CFTR expression in cystic fibrosis cells in the era of high efficiency modulator therapy. J Cyst Fibros. 2023;1:S12–S16.
60. Pongsapipatana NGG, Cho SD, El-Tanani M, Rabbani SA, Islam S, Essigmann JM, Sriwattanapong K, Porntaveetus T. Urinary proteomic biomarkers of systemic fluoride exposure in children. Ecotoxicol Environ Saf. 2026;309: Article 119643.
61. Lukasiak A, Zajac M. The distribution and role of the CFTR protein in the intracellular compartments. Membranes. 2021;11(11):804.
62. Koivula FNM, McClenaghan NH, Harper AGS, Kelly C. Islet-intrinsic effects of CFTR mutation. Diabetologia. 2016;59(7):1350–1355.
63. Cabrini G, Rimessi A, Borgatti M, Lampronti I, Finotti A, Pinton P, Gambari R. Role of cystic fibrosis bronchial epithelium in neutrophil chemotaxis. Front Immunol. 2020;11:1438.
64. Liu JZBR, Chen C, Wang R. The occurrence, severity degree, and associated risk factors of dental fluorosis among the 12-year-old schoolchildren in Jilin, China. Medicine. 2021;100(22): Article e23820.
65. Arheiam AAA, Gaber A, Fakron S. Dental fluorosis and its associated factors amongst Libyan schoolchildren. Int Dent J. 2022;72(6):853.
66. Irigoyen-Camacho ME, Perez-Perez N, Zepeda-Zepeda MA, Velazquez-Alva MC, Castaño-Seiquer A, Barbero-Navarro I, Sanchez-Perez L. Relationships between dental fluorosis and fluoride concentrations in bottled water and groundwater in low-income children in Mexico. Front Oral Health. 2023;4:1187463.
67. Park SYJS, Han JH, Shin JE, Lee JH, Kang CM. Natal factors affecting developmental defects of enamel in preterm infants: A prospective cohort study. Sci Rep. 2024;14(1):2089.

68. Dakna MHK, Kalousis A, Carpentier S, Kolch W, Schanstra JP, Haubitz M, Vlahou A, Mischak H, Girolam M. Addressing the challenge of defining valid proteomic biomarkers and classifiers. BMC Bioinformatics. 2010;11:594.