Skip to main content

Salivary microbiome in patients undergoing hemodialysis and its associations with the duration of the dialysis



Chronic kidney disease (CKD) patients, especially those with end stage renal disease (ESRD) undergoing hemodialysis (HD), exhibit high prevalence of periodontitis. This cross-sectional study aimed to investigate the periodontal status of HD patients and its relationship with salivary microbiome.


One hundred eight HD patients and one hundred healthy control individuals were recruited. They were subjected to periodontal examination followed by saliva samples collection for 16S rRNA gene sequencing.


The HD patients were with worse periodontal health status, and exhibited higher salivary microbial diversity and lower richness. The periodontal pathogens were significantly enriched in the HD patients. The inferred functional analyze showed microbes enriched in the HD patients were mainly related to metabolism. Despite the periodontal status and overall structure of the microbiome were not significantly altered as the HD duration prolonged, the abundance of Lachnospiraceae [G-2] sp. |HMT_096| is positively correlated with the duration of HD and the community periodontal index (CPI). Five OTUs (operational taxonomic units) belonging to the phyla Firmicutes were enriched as the duration prolonged, and four OTUs originated from the phyla Proteobacteria were negatively related with the CPI index. ESRD patients undergoing HD exhibited microbiota structural, compositional and functional differences compared with the healthy controls. And the species changed as the duration of hemodialysis prolonged.


End stage renal disease changes salivary microbiome and is a risk factor for oral dysbiosis.

Peer Review reports


Chronic kidney disease (CKD) is characterized by a presence of injury and/or a progressive loss of renal function over a period of months or years culminating with end stage kidney disease (ESKD) [1]. Due to serious systemic sequelae, CKD has emerged as a major public health, the global burden of which has increased substantially in recent years [2]. When the functional capacity of the kidney decreases below 5–10% of the normal efficiency, hemodialysis (HD), the most common form of the renal replacement therapy, can improve the long-term survival of patients with ESKD [3].

CKD is a risk factor for chronic periodontitis [4,5,6]. For patients with ESRD, especially, worse oral hygiene and higher prevalence of periodontal diseases were reported [7,8,9]. Moreover, the risks of oral problem increase as the renal failure progresses [10].

A complex community of bacteria, archaea, protozoa, fungi and viruses colonize our human body, which encompasses 10 times more cells and 100 times more genes than the host [11]. The oral microbiome can damage the dynamic balance between health and disease, locally and systemically [12]. The microbiome in saliva sheds from adhering microbial niches on various intraoral surfaces and appears to be the representative of the oral microbiome integrity. Salivary microbiome has a critical role on the genesis of oral diseases, mainly including dental caries, periodontitis and peri-implantitis [13,14,15]; non-oral diseases such as cardiovascular disease [16], cancer [17, 18] and low birth weight neonate [19] are also somewhat related to the dysbiosis of salivary microbiome.

Gut dysbiosis has already been well described in CKD patients [20, 21] and ESRD patients undergoing HD [22]. The oral cavity is the primary entry site of the gastrointestinal, additionally, and the oral pathogenic microbes has been described to induce the dysbiosis of gut microbiota in an animal model [23]. Previous studies, using microbiological methods with limited scope, have found that the contents of periodontal pathogens increased in ESRD patients as compared to non-CKD controls [24, 25]. However, until now, very few studies analyzed the salivary microbiota compositions in HD patients.

In this study, we aimed to compare the periodontal status and saliva microbiota between HD patients and healthy persons, and further investigated the association between the salivary microbiome and the duration of HD treatment. These results enrich our understanding of the relationship between general disease state, oral microbiome changes and oral health status.


Study population

This study was approved by the Institutional Review Board of the Sichuan University Hospital. The case group comprised 108 patients undergoing HD three times per week at the West China Hospital, Sichuan University, Chengdu, China. They were in stable clinical conditions for at least 3 months. Patients with diabetic nephropathy or kidney transplantation, or patients who had been on absorbable antibiotics or undergoing periodontal disease treatment within the preceding 3 months or have smoking habits were excluded. 100 healthy matched subjects were recruited as controls. All participants gave a written informed consent before participation, and the study was carried out in accordance with those ethical standards in the Declaration of Helsinki.

Oral examination

One trained dentist performed all the clinical measurements. The Community Periodontal Index (CPI) was used to determine the oral periodontal health status, with the score of 6 sites for 10 World Health Organization index teeth [26]. The number of decayed, missing, and filled teeth (DMFT) was assessed according WHO guidelines [27]. To achieve a high quality of the data, an examiner repeatability exercise was undertaken prior to the initiation of the study. The kappa score for intra-examiner reliability was 0.94.

Saliva collection and DNA isolation

After the oral examination, saliva samples were collected by the same dentist. All participates were instructed not to eat or brush their teeth at least 2 h before sampling. A Salivette (Sarstedt, GER) collection swab with hard pack was used for the saliva collection and storage. Whole unstimulated saliva was collected using the swab in floor of the mouth. The swab was placed back into the collection tube, and stored at − 80 °C before subsequent analysis. Total DNA was isolated from saliva using the QIAamp DNA micro Kit (QIAGEN Sciences, MD, USA) per the manufacturer’s recommendations, with minor modification by adding an extra lysozyme (3 mg ml− 1, 1.5 h) treatment step to lyse the bacterial cell.

16S rRNA gene amplification and sequencing

Total DNA was extracted from saliva using the QIAamp DNA micro Kit (QIAGEN Sciences, MD, USA) per the manufacturer’s recommendations, with minor modification by adding an extra lysozyme (3 mg ml− 1, 1.5 h) treatment step for lysing the bacterial cell. The V3-V4 regions of the bacterial 16S rRNA gene were PCR amplified using barcoded 338F and 806R. Amplicons were extracted from 2% agarose gels and purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) and quantified using QuantiFluor™ -ST (Promega, USA). Purified amplicons were pooled in equimolar and paired-end sequenced (2 × 250) on an Illumina MiSeq platform (Illumina, San Diego, CA) according to the standard protocols.

Sequence information and data analysis

The raw fastq files were first quality-filtered by Trimmomatic and then merged by FLASH. OTUs were clustered by 97% cutoff using UPARSE (version 7.1). To determine the species name and HOT (Oral Taxon ID), sequences were characterized based on the HOMD database. The relative abundance of each bacterial taxon was calculated and typically presented as “Mean ± SD”. The core microbiome of salivary community was analyzed based on previous study [28]. Mann-Whitney U-test was used to compare the differences between the two groups (P < 0.05), and Kruskal-Wallis test was used for comparation among three groups (P < 0.05). X2 test or Fisher’s exact test was used to compare the categorical variables in the contingency tables (P < 0.05). Principal co-ordinates analysis (PCoA) and supervised Partial least squares-discriminant analysis (PLS­DA) analysis were assessed using the Simca­P software (version 13.0, Umetrics, Umea, Sweden). Linear discriminant analysis (LDA) effect size (LEfSe) method was implemented with LEfSe (version 1.0). The functionality of the microbiota was performed using PiCRUST and compared between the two groups. Cytoscape (version 3.6.0) was used to visualize the correlations between the microbiota and duration of HD and other clinical index, with a Spearman’s correlation coefficient (|correlation ρ| > 0.3 or 0.4, and P < 0.05).


Sociodemographic and Oral health characteristics of the study population

One hundred eight HD patients and one hundred healthy subjects were recruited. The two groups matched in age, gender and body mass indices (BMIs) (Table 1). For the HD group, individuals with CPI 1 (23.0% v. 0.9%) and CPI 2 (39.0% v. 10.2%) were less common, while the percentages of CPI 3 (27.0% v. 41.7%) and CPI 4 (3.0% v. 47.2%) individuals were all significantly higher than those in the healthy controls. Both groups have no individuals with CPI X. There were 45 (41.7%) individuals with moderate periodontitis and 51 (47.2%) individuals with severe periodontitis in the HD group, compared with 27 (27.0%) and 3 (3.0%) in the control group. Namely, the percentage of individuals with CPI 3 and CPI 4 in the HD group was nearly 3 times that in the control group. Meanwhile, the HD group demonstrated greater number of missing teeth compared with healthy group (P = 0.013).

Table 1 Demographics and Clinical Parameters of all Patients

Overview of salivary microbial differences among the patients

In total, 29 known phyla and 631 genera were identified, and 1430 OTUs were detected at 3% dissimilarity. To study the richness and diversity of the microbiome from the samples, we compared a series of alpha diversity metrics between the HD patients and the healthy controls, including the Chao index, Shannon index and observed OTUs (Fig. 1a). The HD patients had a lower richness estimator (Chao, P = 0.012), and a higher diversity index (Shannon, P = 0.016). Uncultivated phylotypes did not significantly differ between the two groups (P = 0.610). We generated PCoA plots to investigate the relationship of the composition. The unweighted UniFrac plot showed a distinct separation of the HD samples and the healthy controls based on PC1 (13.28%) and PC2 (8.91%) (Fig. 1b). Consistently, the dissimilarity tests, Adonis and ANOSIM, also revealed significant difference between the two microbial communities (supplementary Table S1). We constructed a Venn diagram and found that 696 OTUs were common in the two communities, while 298 and 436 OTUs were unique to the salivary microbiome of healthy controls and that of HD patients, respectively (Fig. 1c).

Fig. 1

Overview of the salivary microbiota. a: Detailed characteristics of alpha diversity (Mean ± SD); b: PCoA discriminant analysis plot; c: Venn diagram based on OTUs. PCoA: Principal co-ordinates analysis; OTUs: Operational taxonomic units

Next, we evaluated the oral microbial community characteristics of the two niches. The taxonomic compositions at the phylum level varied widely between groups (Fig. 2a). Firmicutes, Proteobacteria, Bacteroidetes, Actinobacteria and Fusobacteria phyla were the dominant phylum (97.66% v. 95.91%) in both groups. The relative abundances of Firmicutes (P = 0.008), Bacteroidetes (P = 0.002), Spirochaetae (P = 0.001), Synergistetes (P < 0.001), Tenericutes (P < 0.001), Gracilibacteria (P < 0.001), and the total relative abundance of the other 21 phyla (P < 0.001) were significantly higher in the HD patients, whereas the abundances of Proteobacteria (P = 0.074) and Actinobacteria (P < 0.001) decreased (Fig. 2b). The microbial shift was illustrated in more detail at the genus level. We analyzed the top 30 genera (supplementary Table S2), and the taxa with significant difference were shown in supplementary Table S3, where more than half of these 30 genera were listed. Compared with the healthy cohort, the relative abundances of Ruminococcaceae (P < 0.001), Capnocytophaga (P < 0.001), Porphyromonas (P = 0.001), Veillonellaceae (P < 0.001), Granulicatella (P < 0.001) and other 9 genera were significantly enriched in the HD samples. In contrast, Lautropia (P = 0.031), Prevotella [G-7] (P < 0.001), Actinomyces (P < 0.001), Veillonella (P < 0.001), Rothia (P = 0.031) and Leptotrichia (P = 0.027) showed significant overabundance in the healthy controls.

Fig. 2

Taxonomic classification of the salivary microbiota at the phylum level. a: Relative taxa abundance comparison of the salivary microbiota at the phylum level; b: Relative abundance of the top 10 bacterial phyla. The significance of differences between groups was assessed using Mann–Whitney U-test (*P < 0.05, **P < 0.01, ***P < 0.001)

Changes of the Core microbiota

The core microbiome with average relative abundance > 1% or species with statistically difference were charted in Fig. 3. Although many bacteria were prevalent in both communities, Rothia aeria (P < 0.001), Actinomyces odontolyticus (P < 0.001), Fusobacterium periodonticum (P = 0.001), Oribacterium asaccharolyticum (P = 0.033) were unique with relative abundances > 1% to the healthy group. On the contrary, Porphyromonas gingivalis (P < 0.001), Neisseria elongate (P < 0.001), Catonella morbi (P = 0.008), Porphyromonas endodontalis (P < 0.001) and Capnocytophaga leadbetteri (P < 0.001) were only detected in the HD group as core salivary microbiome with a mean relative abundance > 1%.

Fig. 3

Core salivary microbiome with relative abundance > 1% or relative abundance of < 1% while presenting statistically difference in the HD patients or the healthy controls. The core microbiome was defined as species present in 75% or greater of samples. The inner box labeled with 1 indicates the species is more abundant in the healthy controls than in the HD patients. The inner box labeled with 2 indicates the species is more abundant in the HD patients than in the controls. The significance of differences between groups was assessed using Mann-Whitney U-test (*P < 0.05, **P < 0.01, ***P < 0.001)

Species associated with periodontal diseases

Considering CKD as a risk factor for periodontal health, we evaluated the microbes that are the important periodontal pathogens belonging to the “red complex” and “orange complex” [29]. HD patients appeared to have higher levels of periodontal pathogens except for Prevotella nigrescens (P = 0.147) (supplementary Table S4).

Predicted microbial function

We generated PICRUSt to predict the gene functions in the salivary microbiome (Fig. 4). Microbes with greater relative abundances in the healthy community had functions that were most related to cellular processes and signaling (ABC transporters, secretion system and pores ion channels) and energy metabolism (oxidative phosphorylation). In contrast, the microbes in the samples of the HD patients were more likely to have functionality related to metabolism including carbohydrate metabolism (fructose and mannose metabolism), amino acid metabolism (valine leucine and isoleucine degradation, beta alanine metabolism), xenobiotics biodegradation and metabolism, as well as lipid metabolism. Moreover, microbes with the function of sporulation also increased in the HD patients.

Fig. 4

Linear discriminant analysis scores of the enriched microbial functions in the healthy controls (Green) and the HD patients (Red). LDA: Linear discriminant analysis

Associations of the salivary microbiome with duration of HD therapy and clinical indices

Next, we performed a correlation analysis of the differential microbiota at the OTU level with the duration of HD therapy and oral clinical indices. Figure 5 showed the constructed related network (|correlation ρ| > 0.3 or 0.4, and P < 0.05), and the supplementary Table S5 provided the detailed description. The abundance of Lachnospiraceae [G-2] sp. | HMT_096 | was positively correlated with both the duration of HD and the CPI index of the patients. Agrobacterium tumefaciens, Pseudomonas marincola, Vibrio gigantis, Streptococcus anginosus, Pseudomonas orientalis and Arthrospira platensi showed moderate correlations (|correlation ρ| > 0.4) with CPI index, were negatively correlated with that periodontal indicator. In addition, there were another 6 species exhibiting weak positive correlations with the duration of HD therapy (0.3 < correlation ρ < 0.4, P < 0.05).

Fig. 5

Associations of the salivary microbiotas with duration of HD and CPI index. Correlations were detected, and the indices with Spearman’s correlation (|correlation ρ| coefficients > 0.3 or 0.4, P < 0.05) were visualized in Cytoscape. If the correlations are negative, the connecting line is blue; if the correlation is positive, the connecting line is red. The line thickness represents the size of the correlation coefficient

To determine whether the duration of HD therapy could be a potential effect on the ESRD patients’ salivary microbiome, we divided all the patients into three groups: (i) HD1 group (n = 23): 3 m < duration of HD ≤1y; (ii) HD2 group (n = 61): 1y < duration of HD ≤5y and (iii) HD3 group (n = 24): duration of HD >5y. All the patients in the HD3 group were with periodontitis and the percentage of individuals with CPI 3 and CPI 4 was slightly higher than the other two groups (supplementary Table S6). The overall structure of the microbial community (supplementary Figure S1, supplementary Table S7) were not altered significantly. No significant difference was observed at the phylum level (supplementary Table S8) and among periodontal pathogens (supplementary Table S9).


Our previous study have found that the ESRD patients undergoing HD suffered a high prevalence and severity of periodontitis [8], further related with the survival of this population [30, 31]. Periodontitis has been identified as a novel risk factor for morbidity and mortality in patients undergoing long-term hemodialysis treatment [32,33,34,35]. Moreover, nonsurgical periodontal treatment, as a relatively simple intervention, has beneficial impact on periodontal status, circulating inflammatory markers and nutritional status in ESRD patients [36,37,38] that underlines the importance of oral health and efficient periodontal treatment to maintain prolong survival and quality of life in these patients. However, very few patients have regular dental care [39, 40] suggesting that more effort should be paid to provide reasonable clinical guideline for dental treatment [41] and to improve the oral health status of this population.

ESRD with HD alters the gut microbiome and may have an impact on the human oral microbiome. Alterations in microbial ecology are related to both oral and systemic diseases, including periodontitis. A recent study described the profile of the subgingival microbiome of peritoneal dialysis patients with periodontitis [42]. In this study, we aimed to investigate the microbial salivary profiles of HD patients and further explore the association between salivary microbiome, periodontal health status and duration of HD. Our finding of variations in oral health status and in salivary microbiome between the HD patients and the healthy individuals suggests that ESRD with HD obviously disturbed the salivary microbiome and periodontal health; duration of HD may further alter the abundance of some specific species, meanwhile, alterations of species may also influence the periodontal health states of the HD patients.

Saliva is a microbial reservoir that contains microbes shedding from both supra and niches. Hence, intraoral transmissions of pathogenic bacteria are likely to be mediated through saliva [43]. Salivary microbiome can enter the tract by swallowing, which may cause imbalance of the intestinal flora [23]. To elucidate the changes of salivary microbiome in ESRD patients can expand the understanding of the interactions between oral microbiome and systemic health. A number of previous studies have used saliva sample as an easy, noninvasive and inexpensive diagnostic tool to assess healthy or disease conditions over the past decades [44,45,46], including periodontitis [47]. Salivary changes have already been reported by previous studies [48], and understanding of the changes in salivary microbiome may present a new point to explore the relationship between ESRD undergoing HD and the other oral manifestations caused by this status, such as taste change, calculus formation and fungal infection. Moreover, previous study reported the salivary microbiome of immunoglobulin A nephropathy (IgAN) differed from that of healthy controls [49]. Inspired by these findings, in this study, we collected saliva samples from the 208 individuals.

Taxonomic differences were detected between the HD patients and the healthy controls at phylum, genus and OTU level. Besides “the others” sort, 7 of the top 10 phyla show significant difference compare with the healthy controls, indicating obvious overall changes in the compositions of the two microbial communities. Changes at phylum level in this study are consistent with previous finding, suggesting that the development of periodontal diseases is intimately linked to shift from a symbiotic dental biofilm, composed mostly of facultative anaerobic bacteria (e.g. Streptococci and Actinomyces), to a dysbiotic microbial community of anaerobic microbes belonging to the phyla essentially comprised of Firmicutes, Spirochaetes, Synergistetes, and Bacteroidete [50]. An increase in Firmicutes has also been reported in the gut microbiome of the ESRD patients during HD [22].

The core salivary microbiome of both HD patients and healthy controls comprised of 104 species (58 species shared by the two groups, while 46 differed). The data indicate that there exists a core salivary microbiome composed of microbes that are most suited to healthy individuals, whereas ESRD with HD therapy might modify this environment. Actinomyces odontolyticus, core microbiome of the healthy community with relative abundance > 1%, is a gram-positive facultative aerobic bacterium, which decreased significantly in the HD community. Previous studies have reported that A. odontolyticus was the host of the candidate phylum Saccharibacteria, being associated with human inflammatory mucosal diseases, thereby killed by a parasitic phase [51, 52]. By contrast, Porphyromonas gingivalis, core microbiome of the HD community with relative abundance > 1%, is a convincedly proposed keystone bacterium in the genesis and development of periodontitis, which was at significantly higher level in the HD group. It is worth mentioning that all the three important periodontal pathogens belonging to the “red complex” were significantly elevated in HD patients, suggesting that ESRD with HD creates a pathogen-rich oral community that could further inducing a high prevalence and a serious degree of periodontitis in this cohort. On the contrary, early study has failed to find the periodontitis-associated taxa differed between CKD and control group [53]. The reasons might include that participants were in different stages of CKD and different primers and sequencing platform were used. Notably, there is a bidirectional relationship between oral microbiome and kidney disease. Evidence suggests that periodontal pathogens and those metabolic products have adverse effects on the survival of ESKD patients who are undergoing dialysis [30]. Moreover, moderate-to-severe periodontal disease also may increase the risk of cardiovascular mortality in the dialysis patients [54]. Though nephropathies also favor the appearance of oral diseases, periodontal examination is still not part of their daily standard medical assessments. Hence, periodontitis in this population must be carefully evaluated by dental professionals, and clinical dental treatment should be taken to reverse this amendable inflammation status.

According to the LDA score, it inferred that the metagenomic pathways that were enriched in the patients undergoing HD were mainly those involved in metabolism, ranging widely from carbohydrate, amino acid and lipid. Porphyromonas, showing significantly higher abundance in the HD group, were amino acid degrading species, which can break down peptides and proteins into amino acids and then degrade them to produce short-chain fatty acids, sulfur compounds, ammonia, and indole/skatole, which all act as virulent influences in periodontitis [55]. The relationship between cellular processes and signaling and the healthy population is unclear.

The severity of periodontal disease slightly worsened as the HD duration prolonged. Further, we found that OTUs, positively related with the duration of HD, were mainly originated from the phyla Firmicutes, while OTUs being negatively related with CPI index of the HD patients were principally from the phyla Proteobacteria, which is similar with the higher ratio of Firmicutes and Proteobacteria in the HD community. Interestingly, Lachnospiraceae [G-2] sp. | HMT_096 |, from the genera Firmicutes, was enriched as the extension of the duration of HD and was more abundant in the patients with severe periodontitis. The duration may not change the overall composition and structure of the salivary microbial community significantly, but contribute to changes of some specific species that may further stimulate a bad effect on the patients’ periodontal status. Considering that the available data about oral microbiota of HD patients remain limited, all these correlated nods may be helpful in finding potential periodontal pathogens or probiotic bacteria and new reasonable periodontal treatment with the HD population in the future.


ESRD patients undergoing HD exhibited more severe periodontal health status, and had microbiota structural, compositional and functional differences compared with the healthy controls. Furthermore, we identified the species changed as the duration of hemodialysis prolonged. End stage renal disease changes salivary microbiome and is a risk factor for oral dysbiosis.

Availability of data and materials

The raw sequence data and the patient data analyzed during the current study are available from the corresponding author on reasonable request.



Chronic kidney disease


Community periodontal index


The number of decayed, missing, and filled teeth


Deoxyribonucleic acid


Especially end stage renal disease


Linear discriminant analysis


Linear discriminant analysis effect size


Immunoglobulin A nephropathy


Undergoing hemodialysis


Operational taxonomic units


Principal co-ordinates analysis


Partial least squares-discriminant analysis


Polymerase chain reaction


Ribosomal ribonucleic acid


  1. 1.

    Daugirdas JT. Handbook of chronic kidney disease management: Wolters Kluwer/Lippincott Williams & Wilkins Health; 2011.

    Google Scholar 

  2. 2.

    Mortality GBD. Causes of death C: global, regional, and national age-sex specific all-cause and cause-specific mortality for 240 causes of death, 1990-2013: a systematic analysis for the global burden of disease study 2013. Lancet. 2015;385(9963):117–71.

    Article  Google Scholar 

  3. 3.

    Himmelfarb J, Ikizler TA. Hemodialysis. N Engl J Med. 2010;363(19):1833–45.

    PubMed  Article  CAS  Google Scholar 

  4. 4.

    Nylund KM, Meurman JH, Heikkinen AM, Furuholm JO, Ortiz F, Ruokonen HM. Oral health in patients with renal disease: a longitudinal study from predialysis to kidney transplantation. Clin Oral Investig. 2018;22(1):339–47.

    PubMed  Article  Google Scholar 

  5. 5.

    Ruokonen H, Nylund K, Furuholm J, Meurman JH, Sorsa T, Kotaniemi K, et al. Oral health and Mortality in patients with chronic kidney disease. J Periodontol. 2017;88(1):26–33.

    PubMed  Article  Google Scholar 

  6. 6.

    Ruospo M, Palmer SC, Craig JC, Gentile G, Johnson DW, Ford PJ, et al. Prevalence and severity of oral disease in adults with chronic kidney disease: a systematic review of observational studies. Nephrol Dial Transplant. 2014;29(2):364–75.

    PubMed  Article  Google Scholar 

  7. 7.

    Sharma P, Dietrich T, Ferro CJ, Cockwell P, Chapple IL. Association between periodontitis and mortality in stages 3-5 chronic kidney disease: NHANES III and linked mortality study. J Clin Periodontol. 2016;43(2):104–13.

    PubMed  PubMed Central  Article  Google Scholar 

  8. 8.

    Zhao D, Zhang S, Chen X, Liu W, Sun N, Guo Y, et al. Evaluation of periodontitis and bone loss in patients undergoing hemodialysis. J Periodontol. 2014;85(11):1515–20.

    PubMed  Article  Google Scholar 

  9. 9.

    de Souza CM, Braosi AP, Luczyszyn SM, Olandoski M, Kotanko P, Craig RG, et al. Association among oral health parameters, periodontitis, and its treatment and mortality in patients undergoing hemodialysis. J Periodontol. 2014;85(6):e169–78.

    PubMed  Article  Google Scholar 

  10. 10.

    Cengiz MI, Sumer P, Cengiz S, Yavuz U. The effect of the duration of the dialysis in hemodialysis patients on dental and periodontal findings. Oral Dis. 2009;15(5):336–41.

    PubMed  Article  CAS  Google Scholar 

  11. 11.

    Human Microbiome Project C. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486(7402):207–14.

    Article  CAS  Google Scholar 

  12. 12.

    Pflughoeft KJ, Versalovic J. Human microbiome in health and disease. Annu Rev Pathol. 2012;7:99–122.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  13. 13.

    Takeshita T, Matsuo K, Furuta M, Shibata Y, Fukami K, Shimazaki Y, et al. Distinct composition of the oral indigenous microbiota in south Korean and Japanese adults. Sci Rep. 2014;4:6990.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  14. 14.

    Crielaard W, Zaura E, Schuller AA, Huse SM, Montijn RC, Keijser BJ. Exploring the oral microbiota of children at various developmental stages of their dentition in the relation to their oral health. BMC Med Genet. 2011;4:22.

    CAS  Google Scholar 

  15. 15.

    Duan XB, Wu TX, Guo YC, Zhou XD, Lei YL, Xu X, et al. Marginal bone loss around non-submerged implants is associated with salivary microbiome during bone healing. Int J Oral Sci. 2017;9(2):95–103.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  16. 16.

    Fak F, Tremaroli V, Bergstrom G, Backhed F. Oral microbiota in patients with atherosclerosis. Atherosclerosis. 2015;243(2):573–8.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  17. 17.

    Fan X, Alekseyenko AV, Wu J, Peters BA, Jacobs EJ, Gapstur SM, et al. Human oral microbiome and prospective risk for pancreatic cancer: a population-based nested case-control study. Gut. 2018;67(1):120–7.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  18. 18.

    Furquim CP, Soares GM, Ribeiro LL, Azcarate-Peril MA, Butz N, Roach J, et al. The salivary microbiome and Oral cancer risk: a pilot study in Fanconi anemia. J Dent Res. 2017;96(3):292–9.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  19. 19.

    Costello EK, Carlisle EM, Bik EM, Morowitz MJ, Relman DA. Microbiome assembly across multiple body sites in low-birthweight infants. MBio. 2013;4(6):e00782–13.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  20. 20.

    De Angelis M, Francavilla R, Piccolo M, De Giacomo A, Gobbetti M. Autism spectrum disorders and intestinal microbiota. Gut Microbes. 2015;6(3):207–13.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  21. 21.

    Simoes-Silva L, Araujo R, Pestana M, Soares-Silva I, Sampaio-Maia B. The microbiome in chronic kidney disease patients undergoing hemodialysis and peritoneal dialysis. Pharmacol Res. 2018;130:143–51.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  22. 22.

    Vaziri ND, Wong J, Pahl M, Piceno YM, Yuan J, DeSantis TZ, et al. Chronic kidney disease alters intestinal microbial flora. Kidney Int. 2013;83(2):308–15.

    PubMed  Article  PubMed Central  Google Scholar 

  23. 23.

    Nakajima M, Arimatsu K, Kato T, Matsuda Y, Minagawa T, Takahashi N, et al. Oral Administration of P. gingivalis induces Dysbiosis of gut microbiota and impaired barrier function leading to dissemination of Enterobacteria to the liver. PLoS One. 2015;10(7):e0134234.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  24. 24.

    Araujo MV, Hong BY, Fava PL, Khan S, Burleson JA, Fares G, et al. End stage renal disease as a modifier of the periodontal microbiome. BMC Nephrol. 2015;16:80.

    PubMed  PubMed Central  Article  Google Scholar 

  25. 25.

    Schmalz G, Kauffels A, Kollmar O, Slotta JE, Vasko R, Muller GA, et al. Oral behavior, dental, periodontal and microbiological findings in patients undergoing hemodialysis and after kidney transplantation. BMC Oral Health. 2016;16(1):72.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  26. 26.

    Cutress TW, Ainamo J, Sardo-Infirri J. The community periodontal index of treatment needs (CPITN) procedure for population groups and individuals. Int Dent J. 1987;37(4):222–33.

    PubMed  CAS  PubMed Central  Google Scholar 

  27. 27.

    Organization WH. Oral health surveys: basic methods: World Health Organization; 2013.

    Google Scholar 

  28. 28.

    Duan X, Wu T, Xu X, Chen D, Mo A, Lei Y, et al. Smoking may Lead to marginal bone loss around non-submerged implants during bone healing by altering salivary microbiome: a prospective study. J Periodontol. 2017;88(12):1297–308.

    PubMed  Article  PubMed Central  Google Scholar 

  29. 29.

    Haffajee AD, Socransky SS, Patel MR, Song X. Microbial complexes in supragingival plaque. Oral Microbiol Immunol. 2008;23(3):196–205.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  30. 30.

    Ricardo AC, Athavale A, Chen J, Hampole H, Garside D, Marucha P, et al. Periodontal disease, chronic kidney disease and mortality: results from the third National Health and nutrition examination survey. BMC Nephrol. 2015;16:97.

    PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Palmer SC, Ruospo M, Wong G, Craig JC, Petruzzi M, De Benedittis M, et al. Dental health and Mortality in people with end-stage kidney disease treated with hemodialysis: a multinational cohort study. Am J Kidney Dis. 2015;66(4):666–76.

    PubMed  Article  Google Scholar 

  32. 32.

    Fisher MA, Taylor GW, Shelton BJ, Jamerson KA, Rahman M, Ojo AO, et al. Periodontal disease and other nontraditional risk factors for CKD. Am J Kidney Dis. 2008;51(1):45–52.

    PubMed  Article  Google Scholar 

  33. 33.

    Deschamps-Lenhardt S, Martin-Cabezas R, Hannedouche T, Huck O. Association between periodontitis and chronic kidney disease: systematic review and meta-analysis. Oral Dis. 2019;25(2):385–402.

    PubMed  Article  CAS  Google Scholar 

  34. 34.

    Hickey NA, Shalamanova L, Whitehead KA, Dempsey-Hibbert N, van der Gast C, Taylor RL. Exploring the putative interactions between chronic kidney disease and chronic periodontitis. Crit Rev Microbiol. 2020;1:1–17.

    Google Scholar 

  35. 35.

    Kapellas K, Singh A, Bertotti M, Nascimento GG, Jamieson LM. Perio CKDc: periodontal and chronic kidney disease association: a systematic review and meta-analysis. Nephrology (Carlton). 2019;24(2):202–12.

    Article  Google Scholar 

  36. 36.

    Fang F, Wu B, Qu Q, Gao J, Yan W, Huang X, et al. The clinical response and systemic effects of non-surgical periodontal therapy in end-stage renal disease patients: a 6-month randomized controlled clinical trial. J Clin Periodontol. 2015;42(6):537–46.

    PubMed  Article  CAS  Google Scholar 

  37. 37.

    Tonetti MS, D'Aiuto F, Nibali L, Donald A, Storry C, Parkar M, et al. Treatment of periodontitis and endothelial function. N Engl J Med. 2007;356(9):911–20.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  38. 38.

    Miyata Y, Obata Y, Mochizuki Y, Kitamura M, Mitsunari K, Matsuo T, et al. Periodontal disease in patients receiving dialysis. Int J Mol Sci. 2019;20:15.

    Article  CAS  Google Scholar 

  39. 39.

    Grubbs V, Plantinga LC, Tuot DS, Powe NR. Chronic kidney disease and use of dental services in a United States public healthcare system: a retrospective cohort study. BMC Nephrol. 2012;13:16.

    PubMed  PubMed Central  Article  Google Scholar 

  40. 40.

    Grubbs V, Plantinga LC, Crews DC, Bibbins-Domingo K, Saran R, Heung M, et al. Vulnerable populations and the association between periodontal and chronic kidney disease. Clin J Am Soc Nephrol. 2011;6(4):711–7.

    PubMed  PubMed Central  Article  Google Scholar 

  41. 41.

    Yuan Q, Xiong QC, Gupta M, Lopez-Pintor RM, Chen XL, Seriwatanachai D, et al. Dental implant treatment for renal failure patients on dialysis: a clinical guideline. Int J Oral Sci. 2017;9(3):125–32.

    PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Zhang X, Chen H, Lu W, Zhu L, Zhou W, Song Z. Characterization of the subgingival microbiota in the peritoneal dialysis patients with periodontitis. Arch Oral Biol. 2020;115:104742.

    PubMed  Article  CAS  Google Scholar 

  43. 43.

    Greenstein G, Lamster I. Bacterial transmission in periodontal diseases: a critical review. J Periodontol. 1997;68(5):421–31.

    PubMed  Article  CAS  Google Scholar 

  44. 44.

    Acharya A, Chan Y, Kheur S, Jin LJ, Watt RM, Mattheos N. Salivary microbiome in non-oral disease: a summary of evidence and commentary. Arch Oral Biol. 2017;83:169–73.

    PubMed  Article  Google Scholar 

  45. 45.

    Prasad S, Tyagi AK, Aggarwal BB. Detection of inflammatory biomarkers in saliva and urine: potential in diagnosis, prevention, and treatment for chronic diseases. Exp Biol Med (Maywood). 2016;241(8):783–99.

    Article  CAS  Google Scholar 

  46. 46.

    Liebsch C, Pitchika V, Pink C, Samietz S, Kastenmuller G, Artati A, et al. The saliva metabolome in association to Oral health status. J Dent Res. 2019;98(6):642–51.

    PubMed  Article  CAS  Google Scholar 

  47. 47.

    Mikkonen JJ, Singh SP, Herrala M, Lappalainen R, Myllymaa S, Kullaa AM. Salivary metabolomics in the diagnosis of oral cancer and periodontal diseases. J Periodontal Res. 2016;51(4):431–7.

    PubMed  Article  CAS  Google Scholar 

  48. 48.

    Kaushik A, Reddy SS, Umesh L, Devi BK, Santana N, Rakesh N. Oral and salivary changes among renal patients undergoing hemodialysis: a cross-sectional study. Indian J Nephrol. 2013;23(2):125–9.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  49. 49.

    Piccolo M, De Angelis M, Lauriero G, Montemurno E, Di Cagno R, Gesualdo L, et al. Salivary microbiota associated with immunoglobulin a nephropathy. Microb Ecol. 2015;70(2):557–65.

    PubMed  Article  CAS  Google Scholar 

  50. 50.

    Matsubara VH, Bandara HM, Ishikawa KH, Mayer MP, Samaranayake LP. The role of probiotic bacteria in managing periodontal disease: a systematic review. Expert Rev Anti-Infect Ther. 2016;14(7):643–55.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  51. 51.

    He X, McLean JS, Edlund A, Yooseph S, Hall AP, Liu SY, et al. Cultivation of a human-associated TM7 phylotype reveals a reduced genome and epibiotic parasitic lifestyle. Proc Natl Acad Sci U S A. 2015;112(1):244–9.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  52. 52.

    Shen M, Yang Y, Shen W, Cen L, McLean JS, Shi W, et al. A linear plasmid-like prophage of Actinomyces odontolyticus promotes biofilm assembly. Appl Environ Microbiol. 2018;84:17.

    Google Scholar 

  53. 53.

    Hu J, Iragavarapu S, Nadkarni GN, Huang R, Erazo M, Bao X, et al. Location-specific Oral microbiome possesses features associated with CKD. Kidney Int Rep. 2018;3(1):193–204.

    PubMed  Article  PubMed Central  Google Scholar 

  54. 54.

    Kshirsagar AV, Craig RG, Moss KL, Beck JD, Offenbacher S, Kotanko P, et al. Periodontal disease adversely affects the survival of patients with end-stage renal disease. Kidney Int. 2009;75(7):746–51.

    PubMed  Article  PubMed Central  Google Scholar 

  55. 55.

    Takahashi N. Oral microbiome metabolism: from “who are they?” to “what are they doing?”. J Dent Res. 2015;94(12):1628–37.

    PubMed  Article  CAS  Google Scholar 

Download references


We wish to thank all participants in this study.


The corresponding author QY is grateful to the National Natural Science Foundation of China and Sichuan Province Science and Technology Innovation Team Program for financially supporting this research (NSFC81722014 and 2017TD0016), respectively. The co-corresponding author YQL is grateful to the National Natural Science Foundation of China for financially supporting this research (NSFC31870065). Author XDZ thanks the National Key R & D Program of China during the 13th Five-Year Plan for providing financial support for this study (2016YFC1102700). The funding body had no role in study design, data collection and analyses, decision to publish or preparation of the manuscript.

Author information




XBD contributed to data acquisition, analysis and interpretation, drafted and critically revised the manuscript; QY and YQL contributed to design and critically revised the manuscript; XLC, MG and DS contributed to conception and critically revised the manuscript; HXX, QCX, TX and DL contributed to data acquisition; ACM, XDZ and XT contributed to data analysis and interpretation. All authors gave final approval and agree to be accountable for all aspects of the work.

Corresponding authors

Correspondence to Yuqing Li or Quan Yuan.

Ethics declarations

Ethics approval and consent to participate

The design and procedures of this study were approved by the Institutional the West China Hospital of Stomatology, Sichuan University. Participants provided informed consent via a signed statement before participation.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1: Table S1.

Dissimilarity Tests between the Healthy Controls and the HD Patients. Table S2. Relative Abundance of the Top 30 Bacterial Genera. Table S3. Relative abundance of the top 30 bacterial genera with significant difference between the groups. Table S4. Relative Abundance of the Species Associated with Diseases in the Healthy Controls and the HD Patients. Table S5. The detailed information of the Cytoscape analyzation. Table S6. Demographics and Clinical Parameters of Patients in the Three Groups. Table S7. Comparison of Bacterial Diversity, Richness, Observed Operational Taxonomic Units (OTUs) and Uncultured Species among the Three Groups. Table S8. Relative Abundance of the Top 10 Bacterial Phyla. Table S9. Relative Abundance of the Species Associated with Diseases in the Three Groups. Figure S1. PLS-DA discriminant analysis plot. PLS-DA: Partial least squares-discriminant analysis; HD1 group (N = 23): 3 m < duration of HD ≤1y; HD2 group (N = 61): 1y < duration of HD ≤5y; HD3 group (N = 24):5y < duration of HD.

Rights and permissions

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

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Duan, X., Chen, X., Gupta, M. et al. Salivary microbiome in patients undergoing hemodialysis and its associations with the duration of the dialysis. BMC Nephrol 21, 414 (2020).

Download citation


  • Saliva
  • Oral microbiome
  • Hemodialysis
  • Periodontal disease