Abstract
Background:
Chlorhexidine is a standard antiseptic in periodontal therapy, yet its long-term impact on the commensal oral microbiota and antimicrobial resistance remains unclear.
Objective:
This study evaluated the ecological effects of 0.06% chlorhexidine toothpaste on the salivary microbiota of healthy individuals, focusing on taxonomic composition, functional shifts, and prevalence of predicted antimicrobial resistance pathways.
Methods:
This study was designed as a prospective, single-arm interventional trial to evaluate the temporal dynamics of the oral microbiota. Saliva samples were collected from eleven healthy participants at three specific time points: baseline (following a 4-week run-in period using a reference fluoride toothpaste), intervention (after 12 weeks of using 0.06% chlorhexidine digluconate toothpaste), and washout (after a subsequent 12-week period using the reference fluoride toothpaste). Sequencing targeted the 16S rRNA V3–V4 region.
Results:
Alpha diversity remained stable throughout the study. However, beta diversity analysis revealed significant shifts in community structure, which were most pronounced between the intervention and washout phases. We identified forty-one differentially abundant metabolic pathways, including a significant downregulation of predicted beta-lactam resistance pathways during the washout period. While chlorhexidine did not increase predicted resistance-related pathways, delayed compositional shifts toward gram-negative taxa were observed at the end of the 12-week washout period.
Conclusion:
These findings suggest that 0.06% chlorhexidine is not ecologically neutral, and dental professionals should carefully consider the appropriateness of chlorhexidine for long-term use.
Clinical Trial Registration:ClinicalTrials.gov, NCT07570927.
Introduction
Oral health depends critically on the ecological balance of the oral microbiota. While a stable oral microbiota supports mucosal defense and homeostasis, its disruption or dysbiosis precipitates common oral diseases, including periodontitis and dental caries (, ). This ecosystem is significantly influenced by oral hygiene habits and chemical agents applied. Antiseptics are widely used in dental products to suppress oral pathogens; however, it is crucial to determine whether these agents act selectively or if they also indiscriminately reduce beneficial commensal bacteria (, ).
Chlorhexidine (CHX) is a broad-spectrum antiseptic widely used in clinical dentistry for its antimicrobial properties (). Dental products containing CHX that are currently available to the public feature varying concentrations of this antiseptic and are marketed in diverse formulations, including mouthwashes, toothpastes, and gels (). For routine or full-mouth applications, the highest concentrations in commercially available over-the-counter products reach approximately 0.2% CHX; formulations containing 0.12% CHX or above are recommended for short-term use, typically limited to a maximum of 28 days (with higher concentrations restricted even further). Conversely, lower concentrations, such as 0.06% CHX, are clinically indicated for extended periods of up to 6 months, or are marketed for long-term use without explicit time restrictions (, ).
CHX's effectiveness in reducing oral bacterial load has been well-documented (, ), but its long-term impact on the salivary microbiota and its potential implications for antibiotic resistance remain areas of active research (, , ).
CHX exerts its antimicrobial activity by disrupting bacterial cell membranes. It generally displays higher effectiveness against gram-positive bacteria due to their higher negative surface charge, which facilitates strong binding with the cationic antiseptic. In contrast, gram-negative bacteria are less sensitive due to the complexity of their cell envelope, which generally impedes the rapid penetration of the antiseptic (). This differential sensitivity suggests a potential for inducing ecological dysbiosis by altering the proportional representation of species. Long-term use of CHX could potentially lead to undesirable ecological changes that favor resistant strains or disrupt the diversity of commensal species (, ). Understanding these changes is essential for assessing the broader impacts of CHX use on oral health.
Antimicrobial resistance (AMR) represents a significant global health challenge, with implications that extend beyond individual treatments (). Consequently, clarifying the relationship between CHX use and AMR is of paramount importance for the dental community. Recent evidence suggests a possible link between exposure to CHX and the development of AMR. Specifically, in vitro studies have indicated reduced CHX sensitivity in bacterial strains that already exhibit resistance to certain antibiotics (, ). Furthermore, emerging evidence suggests that CHX exposure may drive the development of cross-resistance to clinically relevant antibiotics (). However, it is important to note that these findings are predominantly based on in vitro models utilizing sub-lethal concentrations of CHX, which differ from those applied in clinical practice. To date, there is a lack of in vivo studies investigating this relationship in the context of commercial dental products and standard usage protocols. Therefore, investigating the effects of CHX on the emergence and persistence of antibiotic-resistant strains in the oral microbiota is crucial for guiding safe and effective clinical usage.
Therefore, the aim of this study was to investigate the longitudinal effects of chlorhexidine on the composition and functional potential of the salivary microbiota in healthy individuals. The study was designed to address the research question of whether daily use of 0.06% chlorhexidine digluconate induces ecological shifts in the salivary microbiota and alters the prevalence of predicted AMR markers.
Material and methods
Study and participants
This single-arm, prospective, longitudinal interventional study was conducted in strict adherence to the pre-established study protocol, which is publicly archived and accessible at https://dx.doi.org/10.17504/protocols.io.ewov1244ogr2/v1. To ensure maximum transparency, reproducibility, and methodological rigor, the study design and experimental workflow were structured in accordance with the STORMS (Strengthening the Reporting of Microbiome Studies (), development framework, and followed the specialized clinical guidelines for oral microbiome research proposed by Zaura et al. ().
The study was approved by the Ethics Committee of the Faculty of Medicine, Masaryk University (Approval No. 68/2020). All participants were fully informed about the nature, scope, and potential risks of the research, and provided written informed consent prior to their formal inclusion in the study cohort.
Participants were recruited from the local community in Brno, Czech Republic, between March and April 2022, through direct personal communication by the principal investigator. To ensure a final target sample size of eleven despite anticipated attrition, an initial cohort of 14 young, non-smoking, systemically healthy adults from a single geographical region was enrolled. Selection followed strict criteria detailed in the Study Protocol; briefly, eligible participants had to be unmedicated, with no antibiotic therapy in the preceding 6 months and no special dietary regulations.
Dental eligibility was confirmed via clinical examination. Inclusion required an intact natural dentition free of untreated caries, inflammatory periodontal diseases, dental implants, or orthodontic appliances. To minimize confounding effects of biofilm variability and ensure uniform antiseptic exposure, only participants with excellent baseline oral hygiene were selected [Plaque Index < 0.2; Gingival Index < 0.143; Silness and Löe ()]. Regular use of interdental cleaning aids was required, and hygiene habits were monitored at every sampling point to ensure consistency. During the study, three participants were excluded: one due to relocation, one due to antibiotic treatment, and one who discontinued the intervention due to extrinsic tooth staining. Baseline characteristics are summarized in Table 1.
Table 1
| Gender | Female:Male | 6:5 |
|---|---|---|
| Age | 30.4 ± 5.8 | Range: 25–43 |
| PI | 0.023 ± 0.06 | Range: 0–0.2 |
| GI | 0.052 ± 0.04 | Range: 0–0.14 |
| Toothbrush type | Manual only | 8 (72.7%) |
| Manual + Sonic combined | 3 (27.3%) | |
| Interdental cleaning | Interdental brushes (IDB) only | 8 (72.7%) |
| IDB + Floss combined | 3 (27.3%) |
Baseline and clinical characteristics of the study population.
Data are presented as mean ± standard deviation (range) or as absolute frequency (percentage). PI, plaque index; GI, gingival index; IDB, interdental brush. Clinical indices were assessed according to the methodology described in the study protocol.
Study design and sampling timeline
The trial comprised three consecutive phases, during which participants were restricted to the dental preparations supplied by the investigators; the use of other treatments like oral probiotics, natural rinses, or additional antiseptic products was prohibited.
During the initial 4-week run-in period (Phase 1; baseline), participants maintained their habitual mechanical oral hygiene practices, including manual or electric toothbrushes and interdental cleaning aids, but were required to use only a provided standard non-medicated toothpaste (Elmex® JUNIOR, 1,400 ppm amine fluoride) to establish a stable microbiological baseline. This was followed by a 12-week intervention period (Phase 2; CHX intervention), during which participants replaced the reference fluoride toothpaste with a test toothpaste containing 0.06% chlorhexidine digluconate, used twice daily according to their established hygiene routine. Finally, during the 12-week washout period (Phase 3; washout), participants reverted to the original reference toothpaste to evaluate the reversibility of the expected microbial shifts. The exact formulation and active ingredients of both dentifrices are detailed in Table 2. Throughout all phases, participants were required to report any antibiotic treatment, which resulted in immediate exclusion from the study. Furthermore, clinical health status and dental indices were monitored at every sampling point to ensure that oral health and hygiene performance remained consistent throughout the entire 28-week duration.
Table 2
| Toothpaste type | Complete ingredient list |
|---|---|
| Reference (baseline & washout) | Aqua, Hydrated Silica, Sorbitol, Hydroxyethylcellulose, Olaflur 1,400 ppm, Aroma, Saccharin, Limonene. * Commercially available as Elmex JUNIOR, CP GABA, Świdnica, Poland |
| Intervention (0.06% CHX) | Aqua, Sorbitol, Hydrated Silica, Isomalt, Panthenol, Hydroxyethylcellulose, PEG-40 Hydrogenated Castor Oil, Cocamidopropyl Betaine, Aroma, Titanium Dioxide, Sodium Fluoride 1,450 ppm, Sodium Saccharin, Sodium Chloride, Aloe Barbadensis Leaf Juice, Sodium Methylparaben, Glycerin, Chlorhexidine Digluconate 0,06%, Cetylpyridinium Chloride, Tocopherol, Sodium Benzoate, Limonene. * Commercially available as GUM Paroex, Sunstar, Barcelona, Spain |
Summary of the toothpaste types used in the study.
Sample collection, DNA extraction, and sequencing
Unstimulated whole saliva samples were collected at the end of each study phase into sterile tubes. To minimize diurnal variability, sampling was consistently performed in the morning before brushing teeth and consuming any food or drink. Participants provided the samples via active expectoration until a standardized volume of 5 mL was reached, with the entire collection process taking no longer than 9 min. Samples were immediately frozen at −20 °C until further processing.
Total bacterial DNA was extracted using the DNeasy® PowerSoil® Pro Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. The V3–V4 hypervariable regions of the 16S ribosomal RNA (16S rRNA) gene were amplified using the standard Illumina 16S metagenomic sequencing library preparation protocol. The specific primers used were: Forward 5′-CCTACGGGNGGCWGCAG-3′ and Reverse 5′-GACTACHVGGGTATCTAATCC-3′ (). Paired-end sequencing (2 × 300 bp) was performed on the Illumina MiSeq platform (Illumina Inc., San Diego, CA, USA) at BioVendor—Laboratorní medicína a.s. (Brno, Czech Republic).
Bioinformatic processing and statistical analysis
Demultiplexed raw reads underwent quality control, adapter trimming, denoising, and consensus chimera removal using the DADA2 pipeline (implemented within QIIME 2) to generate Amplicon Sequence Variants (ASVs) (). Taxonomic assignment was performed using a Naive Bayes classifier trained against the SILVA reference database (version 138.1) with a confidence threshold of 0.7. To remove rare, low-frequency, and potentially spurious features, ASVs were filtered during the bioinformatic pipeline to retain only robust, annotated features for downstream analysis.
Downstream statistical evaluation and data visualization were performed using MicrobiomeAnalyst 3.0 with further data filtering disabled to preserve the pre-filtered ASV structure (). Alpha diversity (Shannon and Simpson indices) was compared using Mann–Whitney U or Kruskal–Wallis tests, as appropriate. Beta diversity was visualized via Principal Coordinates Analysis (PCoA) and calculated using two complementary metrics: Aitchison distance [Euclidean distance on Centered Log-Ratio (CLR) normalized data] and Bray–Curtis dissimilarity on Total Sum Scaling (TSS) normalized data. Differences in overall microbial community composition between timepoints were formally tested via Permutational Multivariate Analysis of Variance (PERMANOVA) with 999 permutations, followed by post-hoc pairwise comparisons. Functional prediction of the salivary microbiota, including metabolic potential and AMR pathways, was conducted using PICRUSt2 (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2) ().
Differential abundance analysis was performed at the phylum, family, genus, and species levels to identify shifts between the baseline, intervention, and washout phases. DESeq2 was employed as the primary statistical tool. To ensure robustness, a sensitivity analysis was conducted using edgeR (Empirical analysis of digital gene expression data in R) with identical threshold settings, and concordance between the two methods was evaluated. Taxa showing nominal significance (p < 0.05) were also recorded to identify potentially biologically relevant features warranting further investigation.
Multiple testing corrections were performed using the Benjamini–Hochberg method to control the false discovery rate (FDR). For functional pathways predicted by PICRUSt2, significance was assessed using adjusted p-values (p-adj). In both cases, a threshold of <0.05 was considered statistically significant.
Generative artificial intelligence tools were utilized in the preparation of this manuscript for visual and linguistic refinement. Specifically, Figure 6 (Longitudinal dynamics of predicted antimicrobial resistance pathways) was generated and visually formatted using ChatGPT (GPT-4 architecture, OpenAI, San Francisco, CA, USA) based on analyzed data from PICRUSt2. Additionally, language refinement of the manuscript text was assisted by the Gemini large language model (Google LLC, Mountain View, CA, USA).
Results
A total of 33 sequenced samples from 11 participants were analyzed. To ensure consistent oral health status, dental indices were determined at each of the three sampling time points (baseline, post-intervention, and post-washout). Throughout all three phases of the project, no significant changes in dental index values were recorded, adhering to the inclusion criteria defined in the Study Protocol. All salivary samples yielded high-quality raw reads ranging from 61,000 to 142,000 reads per sample. Following bioinformatic processing and quality filtering through the DADA2 pipeline, a total of 233 ASVs were identified.
Taxonomic classification using the SILVA reference database (v138.1) identified 8 bacterial phyla, 14 classes, 33 orders, 47 families, 87 genera, and 117 species.
Microbial diversity
To assess the impact of chlorhexidine intervention on salivary microbiota ecology, alpha and beta diversity across the study phases were evaluated. No statistically significant differences in species richness or evenness were detected between phases via alpha diversity analysis using Shannon and Simpson indices (p > 0.05; Figures 1A,B). In contrast, beta diversity analysis revealed a statistically significant shift in salivary microbial community structure across the study phases, consistent across both evaluated metrics (PERMANOVA, p < 0.05). Specifically, overall community differences were significant when assessed by both Bray–Curtis dissimilarity (p = 0.03, Figure 1C) and Aitchison distance (p = 0.022). Pairwise comparisons confirmed that this temporal shift was driven specifically by a significant difference between Phase 2 and Phase 3 (Bray–Curtis, FDR = 0.018; Aitchison, FDR = 0.006).
Figure 1
Taxonomic composition and differential abundance
To identify the drivers of the observed community shifts, we further analyzed changes in bacterial composition at the taxonomic levels.
Taxonomic profiling at the phylum level identified Firmicutes as the predominant phylum across all study phases (Figure 2A). However, significant compositional shifts were detected, particularly at the end of the washout period. Compared to baseline, the relative abundances of Proteobacteria and Campylobacterota significantly increased (FDR < 0.05), whereas Actinobacteriota significantly decreased (FDR < 0.05). Furthermore, comparisons between the intervention and washout phases confirmed the significant expansion of Proteobacteria and Campylobacterota (FDR < 0.05). A reduction in Actinobacteriota and Bacteroidota was also observed in the washout period relative to the intervention period; however, these changes were not significant after correction for multiple testing (FDR > 0.05).
Figure 2
At the family level, the relative abundances of Pasteurellaceae and Campylobacteraceae were significantly higher in the washout period compared to baseline (FDR < 0.05). Significant increases in these families were also observed in washout period compared to the intervention period (FDR < 0.05).
Taxonomic profiles displaying bacterial genera and species across all study time points are provided in Supplementary Figures 1 and 2. To provide a comprehensive overview, the 25 most abundant genera were visualized in a heatmap (Figure 3). Consistent with these analyses, the salivary microbiota was dominated by Streptococcus, followed by Haemophilus, Veillonella, Prevotella (including the Prevotella 7 group), and Actinomyces. Detailed analysis revealed that the predominant fraction of streptococci (mean 27.1% of total abundance) consisted of sequences unassigned at the species level. Phylogenetically, this cluster corresponds largely to the Mitis group (Streptococcus mitis, Streptococcus oralis) and exhibited high stability throughout the study, with no statistically significant fluctuations between phases. Other health-associated species (Streptococcus parasanguinis, Streptococcus salivarius) were detected at substantially lower abundances (<1%), while the cariogenic species Streptococcus mutans was present only in negligible quantities (<0.04%).
Figure 3
To ensure robustness, differentially represented taxa were identified using the intersection of two statistical models (EdgeR and DESeq2), as illustrated in Figure 4A. This consensus approach confirmed that the shift in washout phase was primarily driven by the significant enrichment of Haemophilus, Aggregatibacter, Actinobacillus, and Campylobacter (FDR < 0.05). In contrast, Solobacterium and Atopobium were identified as taxa with significantly reduced abundance compared to baseline (Figure 4B).
Figure 4
At the species level, unassigned Haemophilus spp. were significantly more abundant in the washout phase compared to both baseline and the intervention phase (FDR < 0.05). Additionally, a significant increase in unassigned Actinobacillus spp. was observed in washout period when compared to intervention period (FDR < 0.05).
Changes in salivary microbiota function and prediction of AMR potential
Functional prediction analysis using PICRUSt2 and LinDA (Linear models for differential abundance analysis) identified distinct clustering of microbial communities across study phases (Figure 5A), indicating a significant overall shift in metabolic potential. Specifically, a total of 41 metabolic pathways with differing abundances between the intervention and washout phases were identified (p-adj < 0.05). These differentially abundant pathways are visualized in the heatmap (Figure 5B), while a subset of those with the highest statistical significance (p-adj < 0.01) is further detailed in Supplementary Figure 3.
Figure 5
The enriched pathways primarily represented functional categories related to the biosynthesis of cellular membrane components and essential cofactors. Specifically, pathways related to the superpathway of fatty acid biosynthesis initiation (FASYN-INITIAL-PWY), fatty acid elongation—saturated (FASYN-ELONG-PWY), biotin biosynthesis I (BIOTIN-BIOSYNTHESIS-PWY), and 8-amino-7-oxononanoate biosynthesis I (PWY-6519) were significantly enriched during the washout period compared to the intervention period (p < 0.01, p-adj < 0.04; Figure 5B). Similarly, pathways involved in lipopolysaccharide biosynthesis, such as the superpathway of (Kdo)_2-lipid A biosynthesis (KDO-NAGLIPASYN-PWY) and ADP-L-glycero-beta-D-manno-heptose biosynthesis (PWY0-1241) showed significant enrichment. In contrast, metabolic pathways involved in the nitrogen cycle (e.g., DENITRIFICATION-PWY, PWY-490-3) remained stable throughout the study period (p > 0.05).
Out of 310 predicted metabolic pathways, a targeted analysis was performed on pathways associated with antimicrobial resistance mechanisms. Statistical analysis revealed significant alterations in the abundance of peptidoglycan biosynthesis pathways linked to resistance. Specifically, pathways PWY-6470 and PWY-6471 associated with beta-lactam resistance were significantly downregulated in the washout phase compared to the chlorhexidine intervention phase (p-adj < 0.002, Supplementary Table 1). Longitudinal analysis of standardized abundance (Z-score) highlights the distinct dynamics of these markers against the global metabolic background. While the overall community profile showed variable trends, both identified resistance pathways followed a specific pattern of upregulation during the intervention, followed by a sharp decline after discontinuation (Figure 6).
Figure 6
Discussion
Given the widespread long-term use of lower-concentration antiseptics, this study investigated the modulation of salivary microbiota by 0.06% chlorhexidine digluconate over a 12-week period and assessed the reversibility of these changes during a subsequent 12-week washout.
Since its introduction in the 1970s, the efficacy of chlorhexidine in reducing bacterial load and inhibiting dental plaque formation has been well established (–). However, this antiseptic also has associated drawbacks. While the most discussed side effect is tooth discoloration, which is well-documented particularly with long-term use of CHX mouthwashes (, , ), only one out of 14 participants in our cohort exhibited visible pigmentation after CHX toothpaste. Therefore, a biologically more significant concern is the alteration of the oral microbiota (, ).
Taxonomic shifts and recolonization dynamics
In the present study, our analysis of the salivary microbiota revealed a significant shift from gram-positive to gram-negative taxa, which was primarily caused by a marked increase in Proteobacteria during the washout phase. Firmicutes remained the dominant phylum across all three phases of the study. However, the stable relative abundance of Firmicutes during the CHX intervention (Figure 2A) requires cautious interpretation. Given that CHX reduces the total bacterial load (, ), the observation that phylum proportions remained comparable to baseline implies a uniform suppression of the entire community. Consequently, while Figure 2B depicts a stable relative profile for Firmicutes in the intervention phase, comparable to baseline, absolute quantification would undoubtedly demonstrate a significant decline in Firmicutes abundance (and all other phyla). The ecological dynamics shifted only during the subsequent washout period, where we observed distinct compositional changes that moved away from the initial homeostatic state.
The impact of CHX on the oral microbiota as reported in current literature varies significantly, suggesting that the microbial response is influenced by both CHX concentration and study duration. Interestingly, our observation of a non-significant trend toward Proteobacteria reduction during the intervention phase aligns with similar findings by Wattanawongwan et al. (), who also reported a decrease that did not reach statistical significance using a 0.12% CHX formulation. This shared observation of a downward trend at lower clinical concentrations stands in contrast to studies utilizing 0.2% CHX, where a significant and immediate enrichment of Proteobacteria was documented (, ). These differences suggest that while high doses of CHX can rapidly induce dysbiotic changes during the active phase of treatment, lower concentrations may have a more subtle suppressive effect on the salivary microbiota.
Interestingly, there were no statistically significant differences at the genus level between Phase 1 (baseline) and Phase 2 (CHX intervention) of our study. This finding contrasts with previous reports utilizing higher CHX concentrations (e.g., 0.12% or 0.2%), which have documented significant bacterial shifts after CHX treatment. However, consensus among these studies is lacking, with varying reports on which specific genera increase or decrease (, , ). In our study, changes in genus abundance were only apparent during the washout phase. The statistically significant genera enriched in Phase 3—Haemophilus, Actinobacillus, and Aggregatibacter—all belong to the phylum Proteobacteria, specifically the family Pasteurellaceae. Although the precise drivers of this synchronized surge remain to be fully elucidated, it is essential to interpret this increased abundance of gram-negative taxa with caution. Given that alpha diversity and clinical indices remained stable throughout our study, these microbial changes likely represent a compositional adaptation rather than a clinically relevant dysbiosis, particularly in individuals with a healthy oral microbiota.
To validate the composition of the baseline salivary microbiota, we compared our data with the core microbiome profile defined by Relvas et al. (). We observed a high degree of concordance, with 14 of the 15 most abundant genera in our cohort matching the core taxa reported in their study. The only distinction was the presence of Aggregatibacter among the top 15 genera in our dataset (Supplementary Figure 1). Importantly, this minor discrepancy occurred despite both studies targeting the same 16S rRNA variable region (V3–V4) and utilizing the SILVA reference database, and it likely stems from database versioning. In earlier versions, sequences currently assigned to Aggregatibacter may have been categorized within the genera Actinobacillus or Haemophilus.
Functional metabolic shifts, predicted resistance-related pathways, and systemic implications
Previous studies indicate that CHX impairs the ability of oral bacteria to reduce nitrates to nitrites (), based on short-term (7 days) exposures to higher-concentration (0.2% CHX) mouthwashes (, ). In contrast, our results suggest that toothpaste with a lower concentration (0.06% CHX) does not affect this ability to reduce nitrates.
Significant shifts in metabolic pathways were most prominent during the washout phase, characterized by an enrichment of pathways associated with lipopolysaccharide (LPS) and cell envelope biosynthesis. Unlike peptidoglycan biosynthesis, which is common to most bacteria, LPS production is a hallmark of gram-negative taxa (). Therefore, these functional increases directly align with the expanded abundance of gram-negative Proteobacteria observed during this period, providing robust functional evidence for the selective proliferation of this group.
Currently, there is a lack of studies investigating the impact of chlorhexidine on predicted AMR pathways in vivo. The predicted abundance of specific peptidoglycan biosynthesis pathways (PWY-6470 and PWY-6471) exhibited a minor, non-significant increase after CHX treatment, followed by a substantial and significant decrease during the washout phase. These dynamics correspond with the significant fluctuations observed in the abundance of the family Atopobiaceae throughout the study phases (p < 0.01). To the best of our knowledge, the only study examining this topic is the work by Bartsch et al. (), who investigated the oral microbiota in patients exposed to 0.2% CHX. They examined both dental plaque and saliva and, with regard to AMR, focused specifically on the presence of resistance genes. Although they observed a non-significant increase in a tetracycline resistance gene following the CHX intervention, they did not find any statistically significant changes in the prevalence of AMR genes after FDR correction in either dental plaque or saliva, either after CHX treatment (4 weeks) or during the subsequent washout period (4 weeks).
The delayed shift toward gram-negative bacteria observed during the washout phase is noteworthy not only because of its importance for local oral health but also in the broader context of the oral-gut axis. Since the oral cavity serves as a continuous reservoir for microbes capable of translocating to the gastrointestinal tract (), even these temporary compositional changes might have broader implications for systemic microbial homeostasis.
Methodological strengths, limitations, and clinical relevance
The advantage of our study lies in the long-term monitoring of the salivary microbiota exposed to CHX and especially wash out, which provides a unique insight into the state of the salivary microbiota after CHX exposure. The selection of participants was also adapted to this monitoring. The selection of a group with excellent oral hygiene was a strategic decision to minimize microbial shifts associated with mature dental biofilm and gingivitis, which typically involve an increase in gram-negative and anaerobic bacteria (, ). Maintaining a consistent level of oral hygiene and dental plaque control throughout the entire 28 weeks was therefore crucial for this study, and this requirement was successfully met.
Our strict inclusion criteria limited participation to systemically healthy, nonsmoking young adults with regular interdental hygiene habits. While a healthy baseline oral microbiota minimizes confounding factors and significantly improves internal validity, it inherently limits the direct generalizability of our findings to other population groups.
Assuming that CHX non-selectively reduces oral bacteria, the recolonization dynamics in patients with gingivitis or periodontitis could be different, since the initial state of the oral microbiota in these patients is distinct from that of healthy individuals; therefore, the recovery of the oral microbiota following exposure to CHX may also vary (). Similarly, older populations or smokers, who have a different oral microenvironment and altered host immune responses, may exhibit different microbial changes following exposure to CHX (). Therefore, future longitudinal studies are warranted to assess whether the changes in oral microbiota composition observed in our healthy cohort also occur in at-risk groups.
It is essential to consider the exact formulations of the dentifrices used in this study, as they were explicitly selected to reflect typical home oral care conditions. The reference toothpaste, used during the baseline and washout phases, contained amine fluoride, which is not ecologically neutral. While direct evidence regarding its specific effect on the salivary microbiota is limited, studies on dental plaque demonstrate that fluoride formulations actively influence bacterial composition and metabolic pathways (–). Given these distinct ecological dynamics, the baseline microbiota in our study represents a state already conditioned by habitual fluoride exposure. Consequently, the observed microbial shifts should be interpreted as the result of sequential exposure (fluoride at baseline, the test formulation during the CHX intervention, and a return to fluoride in the washout phase) rather than the isolated effect of CHX alone. Furthermore, during the intervention phase, participants were exposed to a home-care toothpaste rather than a pure chlorhexidine solution. As detailed in Table 2, this formulation also contained fluoride and cetylpyridinium chloride, a substance known to synergistically enhance the antiseptic effect ().
A notable limitation of this study is the reliance on 16S rRNA gene sequencing targeting the V3–V4 region. While this approach is a widely established standard for profiling the oral microbiota, it presents inherent methodological constraints compared to whole-genome shotgun sequencing. Specifically, it lacks sufficient taxonomic resolution to accurately identify bacteria at the species level, restricting our findings primarily to broader taxonomic ranks (Supplementary Figure 2). Furthermore, our analysis is based exclusively on relative abundance. Because CHX is an antiseptic known to significantly reduce the overall bacterial load, changes in relative abundance may not accurately reflect the absolute biological expansion or reduction of specific taxa. The absence of total bacterial quantification is a limitation that must be considered.
From a clinical perspective, the use of 0.06% chlorhexidine digluconate toothpaste induced statistically significant changes in the oral microbiota, which became most evident at the end of the 12-week washout phase. While long-term exposure to this low-concentration antiseptic did not increase the prevalence of predicted antimicrobial resistance pathways, taxonomic changes in the oral bacterial community were observed. Clinically, dental professionals should carefully weigh the therapeutic benefits of CHX products against the risk of inducing these delayed compositional shifts. Our results indicate that the appropriateness of CHX for long-term use should be reconsidered, as its ecological impact on the oral microbiota is not neutral even at lower concentrations.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ebi.ac.uk/ena, PRJEB108494 https://doi.org/10.5281/zenodo.19655210, Zenodo.
Ethics statement
The studies involving humans were approved by Ethics Committee of the Faculty of Medicine, Masaryk University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
VC: Conceptualization, Methodology, Writing – original draft, Writing – review & editing, Formal analysis, Visualization. KB: Methodology, Writing – original draft. MZ: Formal analysis, Methodology, Writing – original draft, Writing – review & editing. ML: Writing – review & editing. MK: Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the Internal Grant Agency of Masaryk University, project no. MUNI/A/1807/2025.
Acknowledgments
The authors would like to thank all the volunteers who participated in this study for their time and cooperation. We would like to thank the BioVendor—Laboratorní medicína a.s. (Brno, Czech Republic) for performing the 16S rRNA gene sequencing. The authors acknowledge the use of generative artificial intelligence (AI) tools, specifically ChatGPT (OpenAI) and Gemini (Google), for visual data formatting and language refinement during the preparation of this manuscript.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was used in the creation of this manuscript. Generative artificial intelligence tools were utilized in the preparation of this manuscript for visual and linguistic refinement. Specifically, Figure 6 (Longitudinal dynamics of predicted antimicrobial resistance pathways) was generated and visually formatted using ChatGPT (GPT-4 architecture, OpenAI, San Francisco, CA, USA) based on analyzed data from PICRUSt2. Additionally, language refinement of the manuscript text was assisted by the Gemini large language model (Google LLC, Mountain View, CA, USA).
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/froh.2026.1887909/full#supplementary-material
Supplementary Figure 1Relative abundance of the top 15 bacterial genera across the three study phases. The stacked bar plots represent the salivary microbiota composition at the genus level for each participant (n = 11). The study phases are defined as follows: (1) Baseline: samples collected following a 4-week run-in period; (2) Intervention: samples collected after 12 weeks of using 0.06% chlorhexidine digluconate toothpaste; and (3) Washout: samples collected after a subsequent 12-week period using the reference fluoride toothpaste. Only the top 15 genera by mean relative abundance are displayed.
Supplementary Figure 2Relative abundance of the top 15 bacterial species across the three study phases. The stacked bar plots display the salivary microbiota composition at the species level for each participant (n = 11). The study phases are defined as follows: (1) Baseline: samples collected following a 4-week run-in period; (2) Intervention: samples collected after 12 weeks of using 0.06% chlorhexidine digluconate toothpaste; and (3) Washout: samples collected after a subsequent 12-week period using the reference fluoride toothpaste. Only the top 15 species by mean relative abundance are displayed.
Supplementary Figure 3Statistical validation of metabolic shifts. Error plot of highly significant predicted metabolic pathways (p-adj < 0.01) comparing the intervention phase (Phase 2) and the washout phase (Phase 3). The plot displays mean differences with 95% confidence intervals.
Abbreviations
AMR, antimicrobial resistance; ASV, amplicon sequence variant; CHX, chlorhexidine; CLR, centered log-ratio; edgeR, empirical analysis of digital gene expression data in R; EMBL-EBI, European Molecular Biology Laboratory—European Bioinformatics Institute; ENA, European Nucleotide Archive; FDR, false discovery rate; LinDA, linear models for differential abundance analysis; p-adj, adjusted p-value; PCoA, principal coordinates analysis; PERMANOVA, permutational multivariate analysis of variance; PICRUSt2, phylogenetic investigation of communities by reconstruction of unobserved states 2; TSS, total sum scaling; 16S rRNA, 16S ribosomal ribonucleic acid.
References
1.
KilianMChappleILCHannigMMarshPDMeuricVPedersenAMLet al. The oral microbiome – an update for oral healthcare professionals. Br Dent J. (2016) 221:657–66. 10.1038/sj.bdj.2016.865
2.
SharmaNBhatiaSSingh SodhiABatraN. Oral microbiome and health. AIMS Microbiol. (2018) 4:42–66. 10.3934/microbiol.2018.1.42
3.
YamashitaYTakeshitaT. The oral microbiome and human health. J Oral Sci. (2017) 59:201–6. 10.2334/JOSNUSD.16-0856
4.
JayanandanMVeeraraghavanVPGovindarajanSMariyappa SubramaniS. Development of oral dysbiosis following use of antimicrobial mouthwashes: a systematic review. Odontology. (2025) 114:1069–80. 10.1007/s10266-025-01269-0
5.
BrookesZLSBelfieldLAAshworthACasas-AgustenchPRajaMPollardAJet al. Effects of chlorhexidine mouthwash on the oral microbiome. J Dent. (2021) 113:103768. 10.1016/j.jdent.2021.103768
6.
BrookesZLSBescosRBelfieldLAAliKRobertsA. Current uses of chlorhexidine for management of oral disease: a narrative review. J Dent. (2020) 103:103497. 10.1016/j.jdent.2020.103497
7.
JamesPWorthingtonHVParnellCHardingMLamontTCheungAet al. Chlorhexidine mouthrinse as an adjunctive treatment for gingival health. Cochrane Database Syst Rev. (2017) 2017:CD008676. 10.1002/14651858.CD008676.pub2
8.
VekslerAEKayrouzGANewmanMG. Reduction of salivary bacteria by pre-procedural rinses with chlorhexidine 0.12%. J Periodontol. (1991) 62:649–51. 10.1902/jop.1991.62.11.649
9.
LöeHRindom SchiøttC. The effect of mouthrinses and topical application of chlorhexidine on the development of dental plaque and gingivitis in man. J Periodontal Res. (1970) 5:79–83. 10.1111/J.1600-0765.1970.TB00696.X
10.
Poppolo DeusFOuanounouA. Chlorhexidine in dentistry: pharmacology, uses, and adverse effects. Int Dent J. (2022) 72:269–77. 10.1016/J.IDENTJ.2022.01.005
11.
National Center for Biotechnology Information (NCBI). Chlorhexidine Gluconate | C34H54Cl2N10O14 | CID 9552081 - PubChem. Available online at:https://pubchem.ncbi.nlm.nih.gov/compound/chlorhexidine%20gluconate(Accessed January 5, 2026).
12.
BescosRAshworthAClarkeCBrookesZLBelfieldLRodilesAet al. Effects of chlorhexidine mouthwash on the oral microbiome. Sci Rep. (2020) 10:5254. 10.1038/S41598-020-61912-4
13.
SedghiLDiMassaVHarringtonALynchSVKapilaYL. The oral microbiome: role of key organisms and complex networks in oral health and disease. Periodontol 2000. (2021) 87:107–31. 10.1111/PRD.12393
14.
World Health Organization. Global antibiotic resistance surveillance report 2025: WHO Global Antimicrobial Resistance and Use Surveillance System (GLASS) (2025). p. 1–114.
15.
RoyerGOrtiz de la RosaJ-MVuilleminXLacombeBChauFClermontOet al. Reduced chlorhexidine susceptibility is associated with tetracycline resistance tet genes in clinical isolates of Escherichia coli. Antimicrob Agents Chemother. (2022) 66(3). 10.1128/AAC.01972-21
16.
Van den PoelBSaegemanVSchuermansA. Increasing usage of chlorhexidine in health care settings: blessing or curse? A narrative review of the risk of chlorhexidine resistance and the implications for infection prevention and control. Eur J Clin Microbiol Infect Dis. (2022) 41:349–62. 10.1007/S10096-022-04403-W
17.
Tag ElDeinMAYassinASEl-TayebOKashefMT. Chlorhexidine leads to the evolution of antibiotic-resistant Pseudomonas aeruginosa. Eur J Clin Microbiol Infect Dis. (2021) 40:2349–61. 10.1007/S10096-021-04292-5
18.
MirzayiCRensonAFurlanelloCSansoneS-AZohraFElsafourySet al. Reporting guidelines for human microbiome research: the STORMS checklist. Nat Med. (2021) 27(11):1885–92. 10.1038/s41591-021-01552-x
19.
ZauraEPappalardoVYBuijsMJVolgenantCMCBrandtBW. Optimizing the quality of clinical studies on oral microbiome: a practical guide for planning, performing, and reporting. Periodontol 2000. (2021) 85:210–36. 10.1111/PRD.12359
20.
LöeH. The gingival index, the plaque index and the retention index systems. J Periodontol. (1967) 38:610–6. 10.1902/JOP.1967.38.6.610
21.
KlindworthAPruesseESchweerTPepliesJQuastCHornMet al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. (2013) 41:e1. 10.1093/nar/gks808
22.
CallahanBJMcMurdiePJRosenMJHanAWJohnsonAJAHolmesSP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. (2016) 13:581–3. 10.1038/NMETH.3869
23.
LuYZhouGEwaldJPangZShiriTXiaJ. Microbiomeanalyst 2.0: comprehensive statistical, functional and integrative analysis of microbiome data. Nucleic Acids Res. (2023) 51:W310–8. 10.1093/nar/gkad407
24.
DouglasGMMaffeiVJZaneveldJRYurgelSNBrownJRTaylorCMet al. PICRUSt2 for prediction of metagenome functions. Nat Biotechnol. (2020) 38:685–8. 10.1038/s41587-020-0548-6
25.
WattsAAddyM. Tooth discolouration and staining: a review of the literature. Br Dent J. (2001) 190:309–16. 10.1038/sj.bdj.4800959
26.
WattanawongwanWKrasaesinAKhieotaTThongchotchatVPorntaveetusTWiriyakijjaP. Effects of chlorhexidine and a polyherbal mouthwash on the oral microbiome and user satisfaction: a randomized controlled trial. Clin Oral Investig. (2025) 29(12):555. 10.1007/s00784-025-06627-4
27.
AshworthAClarkeCFarnhamGLiddleLBurleighMRodilesAet al. Dietary intake of inorganic nitrate in vegetarians and omnivores and its impact on blood pressure, resting metabolic rate and the oral microbiome. Free Radic Biol Med. (2019) 138:63–72. 10.1016/j.freeradbiomed.2019.05.010
28.
AnnavajhalaMKMayMCompresGFreedbergDEGrahamRStumpSet al. Relationship of the esophageal microbiome and tissue gene expression and links to the oral microbiome: a randomized clinical trial. Clin Transl Gastroenterol. (2020) 11:e00235. 10.14309/ctg.0000000000000235
29.
RelvasMRegueira-IglesiasABalsa-CastroCSalazarFPachecoJJCabralCet al. Relationship between dental and periodontal health status and the salivary microbiome: bacterial diversity, co-occurrence networks and predictive models. Sci Rep. (2021) 11(1):929. 10.1038/s41598-020-79875-x
30.
LimaLGasparSRochaBSAlvesRAlmeidaMG. Current clinical framework on nitric oxide role in periodontal disease and blood pressure. Clin Oral Investig. (2024) 28(10):521. 10.1007/s00784-024-05913-x
31.
KapilVHaydarSMAPearlVLundbergJOWeitzbergEAhluwaliaA. Physiological role for nitrate-reducing oral bacteria in blood pressure control. Free Radic Biol Med. (2013) 55:93–100. 10.1016/j.freeradbiomed.2012.11.013
32.
BertaniBRuizN. Function and biogenesis of lipopolysaccharides. EcoSal Plus. (2018) 8. 10.1128/ecosalplus.esp-0001-2018
33.
BartschSKohnertEKreutzCWoelberJPAndersonABurkhardtA-Set al. Chlorhexidine digluconate mouthwash alters the oral microbial composition and affects the prevalence of antimicrobial resistance genes. Front Microbiol. (2024) 15:1429692. 10.3389/fmicb.2024.1429692
34.
KunathBJDe RudderCLacznyCCLetellierEWilmesP. The oral–gut microbiome axis in health and disease. Nat Rev Microbiol. (2024) 22(12):791–805. 10.1038/s41579-024-01075-5
35.
Van Der PloegGRBrandtBWKeijserBJFVan Der VeenMHVolgenantCMCZauraEet al. Multi-way modelling of oral microbial dynamics and host-microbiome interactions during induced gingivitis. NPJ Biofilms Microbiomes. (2024) 10:89. 10.1038/S41522-024-00565-X
36.
AbuslemeLHoareAHongBYDiazPI. Microbial signatures of health, gingivitis, and periodontitis. Periodontol 2000. (2021) 86:57–78. 10.1111/PRD.12362
37.
IniestaMChamorroCAmbrosioNMarínMJSanzMHerreraD. Subgingival microbiome in periodontal health, gingivitis and different stages of periodontitis. J Clin Periodontol. (2023) 50:905–20. 10.1111/JCPE.13793
38.
MasonMRPreshawPMNagarajaHNDabdoubSMRahmanAKumarPS. The subgingival microbiome of clinically healthy current and never smokers. ISME J. (2015) 9:268–72. 10.1038/ISMEJ.2014.114
39.
Carda-DiéguezMMoazzezRMiraA. Functional changes in the oral microbiome after use of fluoride and arginine containing dentifrices: a metagenomic and metatranscriptomic study. Microbiome. (2022) 10(1):159. 10.1186/S40168-022-01338-4
40.
ZhangWXieQYHuangSChuCHLoECM. Changes in oral microbiome in preschool children after using toothpaste with different fluoride concentrations. J Dent. (2025) 161:105941. 10.1016/J.JDENT.2025.105941
41.
ZhangQGuanLGuoJChuanATongJBanJet al. Application of fluoride disturbs plaque microecology and promotes remineralization of enamel initial caries. J Oral Microbiol. (2022) 18(1). 10.1080/20002297.2022.2105022
42.
López-LópezAMiraA. Shifts in composition and activity of oral biofilms after fluoride exposure. Microb Ecol. (2020) 80:729–38. 10.1007/S00248-020-01531-8
43.
ZhuXLiYZhangXZhangPTianQMaCet al. Combination of cetylpyridinium chloride and chlorhexidine acetate: a promising candidate for rapid killing of gram-positive/gram-negative bacteria and fungi. Curr Microbiol. (2023) 80:97. 10.1007/S00284-023-03198-Z
Summary
Keywords
antimicrobial resistance, chlorhexidine digluconate, oral bacteria, oral microbiota, prediction of resistance, salivary microbiota, antiseptics
Citation
Chuchmova V, Brodíková K, Zapletalova M, Lindholm M and Krsek M (2026) Delayed shifts in salivary microbiota composition following long-term use of 0.06% chlorhexidine toothpaste. Front. Oral Health 7:1887909. doi: 10.3389/froh.2026.1887909
Received
21 May 2026
Revised
17 July 2026
Accepted
28 July 2026
Published
11 August 2026
Volume
7 - 2026
Edited by
Gabriel Nunes, State University of Campinas, Brazil
Reviewed by
Denis Bourgeois, Université Claude Bernard Lyon 1, France
Ali Al-Ahmad, University of Freiburg Medical Center, Germany
Updates
Copyright
© 2026 Chuchmova, Brodíková, Zapletalova, Lindholm and Krsek.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Veronika Chuchmova veronika.chuchmova@med.muni.cz
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.