Xuyu Xiang
Bo Peng
Kai Liu1,2
Hao Li
Yingzi Ming- 1The Transplantation Center of the Third Xiangya Hospital, Central South University, Changsha, China
- 2Engineering and Technology Research Center for Transplantation Medicine of National Health Commission, Changsha, China
Introduction: Renal transplantation is an effective treatment for the end stage renal disease (ESRD). However, how salivary microbiota changes during perioperative period of renal transplant recipients (RTRs) has not been elucidated.
Methods: Five healthy controls and 11 RTRs who had good recovery were enrolled. Saliva samples were collected before surgery and at 1, 3, 7, and 14 days after surgery. 16S rRNA gene sequencing was performed.
Results: There was no significant difference in the composition of salivary microbiota between ESRD patients and healthy controls. The salivary microbiota of RTRs showed higher operational taxonomic units (OTUs) amount and greater alpha and beta diversity than those of ESRD patients and healthy controls, but gradually stabilized over time. At the phylum level, the relative abundance of Actinobacteria, Tenericutes and Spirochaetes was about ten times different from ESRD patients or healthy controls for RTRs overall in time. The relative abundance of Bacteroidetes, Fusobacteria, Patescibacteria, Leptotrichiaceae and Streptococcaceae was correlated with serum creatinine (Scr) after renal transplantation.
Discussion: In short, salivary microbiota community altered in the perioperative period of renal transplantation and certain species of salivary microbiota had the potential to be a biomarker of postoperative recovery.
1. Introduction
End-stage renal disease (ESRD) represents a serious public health problem fueled by aging populations and a pandemic of chronic non-communicable diseases, which is characterized by high mortality and economic burden. Renal transplantation is one of the effective treatments, with the hope of recovery for patients to normal life. However, there is still a lack of highly sensitive and specific biomarkers with minimal invasion and cost to assess recovery or rejection during the perioperative period.
The oral cavity consists of teeth, gingival groove, tongue, soft and hard palates, buccal mucosa, and tonsils. All the above areas are inhabited by microbiota and soaked in saliva all the time. Each salivary gland is highly permeable and surrounded by capillaries, a feature that allows for a freer exchange of substances between the salivary glands and blood (Wilson et al., 2014). Therefore, the salivary microbiota has the potential to be a bridge between oral (Belstrøm et al., 2018a) and systemic conditions.
Indeed, the relationship between chronic kidney disease (CKD) and gut microbiota has been widely investigated, both regarding changes in the floras of patients with CKD (Crespo-Salgado et al., 2016; Meijers et al., 2019; Ren et al., 2020) and regarding the mechanisms of gut microbiota in the development of CKD (Wang X. et al., 2020; Zhu et al., 2021; Wang et al., 2023). Saliva, one of the largest sources of gut microbiota, may play an important role in kidney disease that salivary microbiota ectopically colonizing the gut may be closely associated with the development of kidney disease and renal function. At the same time, several studies have discussed changes in salivary flora in patients with CKD (Hu et al., 2018; Duan et al., 2020; Guo et al., 2022; Liu et al., 2022). The overall composition of the salivary microbiota in CKD patients is significantly different from that of the healthy population, although the variation in individual flora or individual indicators is not entirely consistent across studies. Hence, the possibility of salivary microbiota functioning at CKD in situ cannot be ruled out. In summary, salivary microbiota has the potential as a diagnostic and therapeutic target for ESRDs or renal transplant recipients (RTRs).
Based on previous studies, we speculate that salivary microbiota in patients after renal transplantation will be significantly different from the preoperative flora and this change may be associated with renal function. Although the alteration of salivary floras in patients with ESRD has been studied, how salivary microbiota dynamic changes during the perioperative period of RTRs and the association between salivary microbiota and postoperative recovery have not been elucidated. Therefore, our study is the first to examine the variations of salivary microbiota during the perioperative period of renal transplantation and the relationship between salivary microbiota and renal function. We aimed to find some special floras associated with the return of renal function as clinical biomarkers.
2. Materials and methods
2.1. Subjects and sample collection
From 1 October 2022 to 18 October 2022, a total of 11 consecutive ESRD patients received renal transplantation in our center and were enrolled. Saliva samples were collected before surgery and at 1, 3, 7, and 14 days after surgery. Saliva samples from five healthy people were also collected as healthy controls. None of the above subjects had oral antibiotics, cortisol, smoking, or drinking history within 6 months.
Before collection, patients fasted for half an hour and rinsed their mouths. Patients spit the saliva into a sterile tube until it reaches 2 ml. Saliva was stored at −80°C immediately after collection.
The study protocol was approved (22207) by the Ethics Committee of the Third Xiangya Hospital of Central South University (Changsha, China). Written informed consent was obtained from all study participants. Experiments were carried out in accordance with the ethical guidelines set by the Declaration of Helsinki 1964 and its later amendments.
2.2. Sequencing
2.2.1. Sampling and DNA extraction
Total genome DNA from samples was extracted using the CTAB/SDS method. DNA concentration and purity were monitored on 1% agarose gel. According to the concentration, DNA was diluted to 1 ng/μl using sterile water.
2.2.2. Amplicon generation
16S rRNA genes were amplified using the specific primer 341F (CCTAYGG-GRBGCASCAG) and 806R (GGACTACNNGGGTATCTAAT) with the barcode. All PCR reactions were carried out in 30 μl of reactions with 15 μl of Phusion® High-Fidelity PCR Master Mix (New England Biolabs); 0.2 μM of forward and reverse primers, and about 10 ng of template DNA. Thermal cycling consisted of initial denaturation at 98°C for 1 min, followed by 30 cycles of denaturation at 98°C for 10 s, annealing at 50°C for 30 s and elongation at 72°C for 60 s. And final extension at 72°C for 5 min.
2.2.3. PCR products quantification and qualification
The same volume of 1X loading buffer (containing SYB green) with the PCR products and operate electrophoresis was mixed on a 2% agarose gel for detection. Samples with a bright main strip between 400 and 450 bp were chosen for further experiments.
2.2.4. PCR products mixing and purification
PCR products were mixed in equidensity ratios. Then, the mixture of PCR products was purified with AxyPrep DNA Gel Extraction Kit (AXYGEN).
2.2.5. Library preparation and sequencing
Sequencing libraries were generated using NEBNext® Ultra™ DNA Library Prep Kit for Illumina (NEB, USA) following the manufacturer's recommendations, and index codes were added. The library quality was assessed on the Qubit@ 2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 system. Finally, the library was sequenced on an Illumina NovaSeq 6000 platform, and 250 bp paired-end reads were generated. Sequences are deposited under SRA PRJNA904953.
2.3. Data analysis
2.3.1. OTU cluster and species annotation
Paired-end reads from the original DNA fragments were merged using FLASH. Sequences analysis was performed by the UPARSE software package using the UPARSE-OTU and UPARSE-OTUref algorithms. In-house Perl scripts were used to analyze alpha (within samples) and beta (among samples) diversity. Sequences with ≥97% similarity were assigned to the same OTU. We picked a representative sequence for each OTUs and used the RDP classifier to annotate taxonomic information for each representative sequence based on Silva 132 database.
2.3.2. Phylogenetic distance and community distribution
Graphical representation of the relative abundance of bacterial diversity from phylum to species can be visualized using the Krona chart. Cluster analysis was preceded by principal component analysis (PCA), which was applied to reduce the dimension of the original variables using the QIIME software package. We used weighted UniFrac distance for principal coordinate analysis (PCoA) and Unweighted Pair Group Method with Arithmetic mean for the abbreviation (UPGMA) Clustering.
2.4. Statistical analysis
Linear discriminant analysis Effect Size (LEfSe) was used for the quantitative analysis of biomarkers within different groups. To identify differences in microbial communities between the two groups, ANOSIM and ADONIS were performed based on the Bray–Curtis dissimilarity distance matrices. A Wilcoxon rank-sum test and unpaired t-test were performed to evaluate differences between the two groups in alpha diversity, principal coordinates, and community difference analysis. Pearson correlation analysis was used to assess the correlation between microbiota and creatinine. A p-value of < 0.05 was required for the results to be considered statistically significant.
3. Results
3.1. Study population
The clinical information of RTRs (age 44.8 ± 12.2 years; 63.6% males) is shown in Table 1. The mean body weight index (BMI) was 20.5 ± 3.8 kg/m2 for RTRs. All RTRs received antihuman thymocyte globulin (ATG) for induction, the same triple immunosuppressive therapy, FK506, mycophenolate mofetil (MMF) plus steroids, and meropenem as a primary antibiotic. The saliva samples were collected before surgery (ESRD, n = 11) and at 1 day (RTR1, n = 9), 3 days (RTR3, n = 11), 7 days (RTR7, n = 11), and 14 days (RTR14, n = 8) after surgery. Generally, the specimens were divided into two groups, namely the ESRD group and the RTR group.
Table 1. Characteristics of renal transplant patients and control subjects.
Healthy controls (HCs, n = 5) ranged in age from 30 to 56 years and consisted of three men and two women. The mean BMI was 24.0 ± 3.1 kg/m2, and the mean serum creatinine (Scr) was 71.3 ± 16.7 umol/L for HCs. We collected saliva samples (n = 9) from each of them two times, 7 days apart, to form the HC group. They all reported no history of chronic diseases or medication.
3.2. Impact of renal transplantation on salivary microbiota in individuals
First, according to the rarefaction curve, Shannon curve, and rank-abundance curve (Supplementary Figure 1), we found that the number of reads of most samples is reasonable. The curves tended to be flat, which indicated that the number of reads was relatively large enough to reflect species richness.
The Venn graph demonstrated the shared and unique OTUs between the three groups (Figure 1A). Overall, the RTR group has more OTUs. However, there was no significant difference in the number of OTUs in a single sample between the three groups. Figures 1B, C show the species composition of each group and individual samples at the phylum level. The ESRD and HC groups had relatively similar species composition, whereas the RTR group was quite different, especially in the relative abundance of Actinobacteria, Tenericutes, and Spirochaetes.
Figure 1. Composition characteristics of salivary microbiota in RTR, ESRD, and HC groups: (A) Venn graph for the OTUs of RTR, ESRD, and HC groups; (B) Salivary species composition of RTR, ESRD, and HC groups in the phylum level; (C) Salivary species composition of each individual in the phylum level; (D) Shannon index of RTR, ESRD, and HC groups; (E) PCoA graph of RTR and ESRD groups; (F) PCoA graph of ESRD and HC groups.
Alpha diversity, including Shannon, Simpson, and so on, provided a measurement of species diversity within a sample. The larger the Shannon index, the greater the diversity. The ESRD group was close to the HC group in alpha diversity. Group RTR always had larger intra-group differences (Figure 1D, RTR vs. ESRD: p < 0.05).
Beta diversity was used to study the intrinsic composition of the microbial structure. The closer the samples were to each other, the more similar the species' composition was. The PCoA was analyzed based on weighted UniFrac distance. According to the PCoA, the microbial composition of the RTR group was significantly different from those of the ESRD (Figure 1E) and HC (Supplementary Figure 2) groups, which was proved by the ADONIS analysis (RTR vs. ESRD: p = 0.001, RTR vs. HC: p = 0.001). On the contrary, the beta diversity of the ESRD group was not significantly different from that of the HC group (Figure 1F).
LEfSe analysis was used to describe differential species between different groups. When the LDA value was >2, the species was a statistically significant biomarker between groups. The results showed that when compared with the ESRD (Figure 2A) and HC (Figure 2B) groups, the relative abundances of Burkholderiaceae, Lautropia, and Actinobacteria of the RTR group were significantly increased, and Neisseriaceae and Neisseria were significantly decreased. The composition of OTU sequences was further transformed into KEGG orthodontics to analyze the differences in predicted function. The pathways related to membrane transport, carbohydrate metabolism, and signal transduction were significantly enriched in the RTR group (Figures 2C, D).
Figure 2. Differential species and KEGG analysis for RTR, ESRD, and HC groups: (A) LEfSe analysis of the salivary microbiota composition between RTR and ESRD groups; (B) LEfSe analysis of the salivary microbiota composition between RTR and HC groups; (C) LEfSe analysis of the predicted function between RTR and ESRD groups; (D) LEfSe analysis of the predicted function between RTR and HC groups.
3.3. The dynamic change in salivary microbiota during the early stage post-renal transplantation
The unique community occupied most OTUs of the RTR1, RTR3, and RTR7 and differed between every two adjacent time points (Figure 3A, RTR1 vs. RTR14: p < 0.05). Over time, the OTUs of RTRs gradually decreased and the shared OTUs with ESRD or HC groups increased, and RTR1 or RTR7 was significantly different from ESRD or HC (Figure 3B, p < 0.05). Figures 3C, D show that the Ace index and intra-group differences of the RTR group descended and approached ESRD and HC groups over time (RTR1 vs. ESRD: p < 0.05, RTR1 vs. RTR14: p < 0.05). Other alpha diversity indexes also showed the same changes in the salivary microbiota of RTRs at different time points (Supplementary Figure 3). At the phylum level, the microbial composition of the RTR group was constantly changing (Figure 3E) but always differed from ESRD or HC groups (Figure 3F). From PCoA, the intrinsic microbial composition of RTR1 and RTR14 groups significantly differed (p = 0.001) from ESRD or HC groups (Figure 3H and Supplementary Figure 4) and the intrinsic microbial composition significantly changed between 3 and 7 days after surgery (Figure 3G and Supplementary Figure 4).
Figure 3. Composition characteristics of salivary microbiota at different time points and states: (A) the Venn graph for OTUs of the RTR group at different time points; (B) the Venn graph for OTUs of the RTR1, RTR14, ESRD, and HC; (C) the Ace index of the RTR group at different time points; (D) the Ace index of the RTR1, RTR14, ESRD, and HC; (E) salivary species composition of RTR group at different time points at the phylum level; (F) salivary species composition of RTR1, RTR14, ESRD, and HC at the phylum level; (G) PCoA graph of RTR3 and RTR7; (H) PCoA graph of RTR14 and ESRD.
3.4. Certain species of salivary microbiota were associated with the recovery of renal function
The correlation between the predominant species of salivary microbiota in RTRs and ESRD and the corresponding Scr on the day of saliva collection was analyzed (Table 2 and Figure 4). In the RTRs group, Bacteroidetes, Fusobacteria, Patescibacteria, and Leptotrichiaceae were positively correlated with Scr, whereas Streptococcaceae was negatively correlated with Scr after renal transplantation. However, these floras were not significantly associated with Scr in the ESRD group.
Table 2. Pearson correlation between the salivary microbiota and Scr after renal transplantation.
Figure 4. Correlation between Bacteroidetes, Fusobacteria, Patescibacteria, and Scr: (A) the scatter diagram of relative abundance of Bacteroidetes and Scr concentration; (B) the scatter diagram of relative abundance of Fusobacteria and Scr concentration; (C) the scatter diagram of relative abundance of Patescibacteria and Scr concentration.
4. Discussion
Salivary microbiota is more stable than gut microbiota, and factors that alter gut microbiota may not significantly alter salivary microbiota (David et al., 2014; Tuganbaev et al., 2022). ESRD patients have changes in the composition of their gut microbiota compared with healthy people (Rysz et al., 2021; Shivani et al., 2022). However, as shown in this research, the salivary microbiota of ESRD patients was similar to that of HCs in terms of the number of individual species, the relative abundance of dominant flora, alpha diversity, and beta diversity, and did not alter significantly due to chronic renal impairment, different long-term treatments, or accompanying changes in life habits. In contrast, the salivary microbiota of RTRs showed huge differences compared with ESRD and HC groups. RTRs contained nearly 10 times as many species of unique salivary microbiota. From both alpha diversity and beta diversity, the RTR group showed higher richness and intra-group differences than ESRD or HC groups. At the phylum level, the relative abundance of Actinobacteria, Tenericutes, and Spirochaetes was about 10 times higher than that of ESRD or HC groups. Actinobacteria is a ubiquitous gram-positive phylum, which has attracted much attention as a rich source of bioactive substances and a complex evolution and diversification process (Miao and Davies, 2010; Barka et al., 2016). As an oral bacterium, Actinobacteria may play a role in the etiology of diabetes (Long et al., 2017; Matsha et al., 2020). The Tenericutes were composed of bacteria that lack a peptidoglycan cell wall. The most well-studied branch of this phylum was Mollicutes, including Mycoplasma. To date, most studies had focused on pathogenic strains of the Mycoplasma order (Wang Y. et al., 2020). As reported, Mycoplasma was associated with oral leukoplakia (Mizuki et al., 2015, 2017), mucositis (Morand and Hatami, 2018), and Fanconi anemia-associated oral carcinoma (Henrich et al., 2014). Spirochaetes were important pathogenic bacteria in the clinic, but they were not well-understood. Some of these caused Lyme disease, leptospirosis, syphilis, and other human diseases. Moreover, Spirochaetes were closely related to periodontal disease and gingivitis (Reed et al., 2018; Yousefi et al., 2020; Zeng et al., 2021), and, in turn, periodontal disease impacted the risk of systemic diseases such as diabetes (Deng et al., 2018). Taken together, these changes occurring in the salivary microbiota of RTRs appeared to be associated with the new onset diabetes, periodontal disease, and gingivitis after renal transplant.
During the early stage (<14 days) after surgery, the salivary microbiota of RTRs was not static. From the Venn graph, we could see that the total number of species decreased over time and most of the salivary microbiota were species unique at each time point. At the phylum level, the relative abundance of Actinobacteria, Cyanobacteria, Epsilonbacteraeota, Tenericutes, and Spirochaetes changed incrementally with time. Among them, Actinobacteria, Tenericutes, and Spirochaetes in RTRs had changed the most. From the point of alpha diversity, including kinds of evaluation indexes such as Ace, Chao1, Shannon, and Simpson, the composition of the salivary microbiota generally moved toward less richness and less variation within groups for each index. As shown in the PCoA figure, RTR1, RTR3, RTR7, and RTR14 groups all had different flora structures. However, the microbiota structure had significant differences only between RTR3 and RTR7 but not RTR1 and RTR3 or RTR7 and RTR14 which may be limited by the insufficient sample size. Moreover, the Chao1 and Ace, two alpha diversity evaluation indexes, described the distribution of bacteria with low abundances and decreased as the number of OTUs decreased during the perioperative period. Hence, we speculated that these low-abundance floras, which occurred in huge changes at different time points after renal transplant, occupy the majority of OTUs. As discussed, RTR14 was closer to HC in terms of OTU number and alpha diversity of salivary microbiota than RTR1. Hence, we considered that changes in the composition of the salivary microbiota of RTRs were a process of stabilization during the early stage after renal transplant.
All renal transplant patients received ATG preoperatively, a drug that inhibited thymocyte activity and helped patients fight off rejection. ATG was also often used to treat severe aplastic anemia and altered patients' salivary microbiota but did not lead to a clear change in diversity over time (Ames et al., 2019). Similarly, RTRs were also all treated with FK506 and MMF which were related to oral cancer (Li et al., 2021) or oral ulcer (Asare and Gatzke, 2020), and these oral diseases were also related to salivary floras (Lin et al., 2021; Bai et al., 2022). Not only salivary microbiota but these immunosuppressants had been linked to the altered intestinal microbiota. Intaking a moderate dose of FK506 maintained immunosuppression, induced normal graft function of the liver, maintained gut barrier integrity, and low plasma endotoxin levels. In addition, it also led to increased species richness and rare species abundance which was consistent with our findings (Jiang et al., 2018). FK506 treatment significantly improved the relative abundance of Bacteroides (Zhang et al., 2018). In our study, the relative abundance of Bacteroides decreased than increased. As an ongoing drug, the effect of FK506 on the elevating relative abundance of Bacteroides may not manifest until the latter part of the early stage after kidney transplantation. MMF enhanced colonic integrity and decreased sympathetic drive in the gut which was associated with the improvement of gut dysbiosis, including the increased abundance of Proteobacteria and Bacteroidetes and decreased abundance of Firmicutes (Robles-Vera et al., 2021) such as the results of our study that the abundance of Firmicutes continuously reduced. MMF increased the alpha diversity of gut microbiota embodied in the first postoperative day of our research (Llorenç et al., 2022). Based on the induction therapy and immunosuppressive therapy, RTRs had immune dysfunction such as acquired immunodeficiency disease (AIDS) patients, according to which we speculated that salivary microbiota changes were similar to those in AIDS patients. In the research of Perez Rosero et al., the significant reduction in the frequency of oral neutrophils in the oral cavity of AIDS individuals was positively related to their CD4+ T cell count and observed OTUs indexes raised in AIDS individuals as alpha diversity of salivary microbiota (Perez Rosero et al., 2021). Interestingly, Alpha diversity altered as the disease progresses (Guo et al., 2021). Compared with healthy people, AIDS patients exhibited a lower abundance of salivary Fusobacteria resembles our study (Yang et al., 2020). For AIDS patients, antiretroviral therapy was an effective treatment. After the treatment, the patient's immune function would be restored to some extent, which was due to the changes in immune function during the perioperative period in RTRs as the gradual recovery of T cell abundance occurred (Bouteloup et al., 2017). These two statuses were all accompanied by decreased salivary alpha diversity (Imahashi et al., 2021). Although RTRs were essentially on constant antibiotics, previous studies had demonstrated that antibiotic use appears to have little effect on salivary flora composition (Tuganbaev et al., 2022). Coincidentally, some external factors which may affect salivary microbiota for patients, such as diet habits (Marsh et al., 2016), drinking water (Sinha et al., 2021), oral hygiene (Belstrøm et al., 2018b), and living environment (David et al., 2014), changed between these time points of saliva collection. In conclusion, we speculated that the factors mentioned earlier functioned together and led to the alternation of salivary microbiota in RTRs like increasing and then gradually decreasing the number of OTUs and alpha diversity index and changing the composition of species and relative abundance of dominant flora with various trends.
Finally, we analyzed the relationship between the dominant flora in saliva and Scr. We found that Bacteroidetes, Fusobacteria, Patescibacteria, and Leptotrichiaceae were positively correlated with Scr, and Streptococcaceae was negatively correlated with Scr after renal transplant. Therefore, these strains could be biomarkers of postoperative recovery of RTRs.
Although the presence of the floras Bacteroidetes, Fusobacteria, Patescibacteria, Leptotrichiaceae, and Streptococcaceae in saliva and their potential correlation with renal function have barely been researched, several studies have elucidated the relation of some of them in the gut and renal dysfunction. Studies have shown an increase in the relative abundance of gut Bacteroidetes in patients with stage 4–5 chronic kidney disease or patients with ESRD receiving hemodialysis (Crespo-Salgado et al., 2016; Wu et al., 2021). Although urinary stones are unlikely to cause kidney damage, urolithiasis patients had significantly lower microbial abundance and higher proportions of Bacteroidetes (Zhou et al., 2020). In a study by Li et al., uremic clearance granules enhanced renal function and decreased levels of Scr, blood urea nitrogen, inflammatory responses, and NF-κB and MAPK expressions in renal tissues of ESRD rats. At the same time, the relative abundances of gut Bacteroidetes descended in response to uremic clearance granules (Li et al., 2022). As a prescription of traditional Chinese medicine for treating chronic kidney disease, the Shenyan Kangfu tablet alleviated renal dysfunction, glomerular and tubular damage, and renal inflammation and reduced the relative abundances of gut Bacteroidetes in the mouse with diabetic kidney disease (Chen et al., 2021). Accompanied by the fecal microbiota transplant, a significant increase of gut Bacteroidetes had the closest correlation with worse response to high salt of salt-sensitive rats, evidenced by increased albuminuria, systolic arterial pressure, and renal T-cell infiltration (Abais-Battad et al., 2021). By contrast, SGL5213 and Daphnetin, two proven renoprotectants, saved kidney function in mice or rats with renal injury and elevated the relative abundances of gut Bacteroidetes (Ho et al., 2021; Zhou et al., 2022).
The relative abundance of Fusobacteria in patients with immunoglobulin A nephropathy or membranous nephropathy exhibited significant elevation when compared with healthy controls (Hu et al., 2020; Zhang et al., 2020; Sugurmar et al., 2021). The microbiota structure showed the same change in type 2 diabetes mellitus, chronic kidney disease, and renal uric acid stone patients (Salguero et al., 2019; Cao et al., 2022). Deltamethrin, as a widely used pyrethroid insecticide, had brought serious problems to the healthy breeding of aquatic animals. A high concentration of deltamethrin damaged the intestine and trunk kidney of goldfish or channel catfish in the early stage with a significant increase or decrease in the abundance of Fusobacteria (Zhou et al., 2021; Yang et al., 2022).
In summary, gastrointestinal Bacteroidetes and Fusobacteria in humans and mice were positively correlated with renal dysfunction which was consistent with our results. Hence, we speculated that these two floras and even more flora may have some connection with renal dysfunction. However, whether the changes in the digestive tract environment brought by renal dysfunction favored their colonization of the digestive tract or their colonization of the digestive tract promoted renal dysfunction remained to be proven.
5. Conclusion
This study has illustrated differences in salivary microbiota communities among RTRs, ESRD patients, and HCs, examining changes in the salivary microbiota community during the short period after renal transplantation. We speculated that changes in the salivary microbiota were a process of stabilization during the early stage after renal transplant, and certain species of salivary microbiota had the potential to be a biomarker of postoperative recovery. Our study first discussed the salivary microbiota variations associated with renal transplantation and the relationship between salivary microbiota and renal function.
Data availability statement
The data presented in the study are deposited in the NCBI repository, accession number PRJNA904953, https://www.ncbi.nlm.nih.gov/bioproject/PRJNA904953.
Ethics statement
The study protocol was approved (22207) by the Ethics Committee of the Third Xiangya Hospital of Central South University (Changsha, China). Written informed consent was obtained from all study participants. Experiments were carried out in accordance with the ethical guidelines set by the Declaration of Helsinki 1964 and its later amendments. The patients/participants provided their written informed consent to participate in this study.
Author contributions
XX, BP, YZ, and YM conceived, designed, and directed the manuscript. XX, BP, and KL wrote and revised the manuscript. XX, PD, and HL participated in the performance of the research. XX, TW, and BP analyzed data. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the National Natural Science Foundation of China (Grant Number: 81771722).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
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/fmicb.2023.1122101/full#supplementary-material
References
Abais-Battad, J. M., Saravia, F. L., Lund, H., Dasinger, J. H., Fehrenbach, D. J., Alsheikh, A. J., et al. (2021). Dietary influences on the Dahl SS rat gut microbiota and its effects on salt-sensitive hypertension and renal damage. Acta Physiol. 232, e13662. doi: 10.1111/apha.13662
Ames, N. J., Barb, J. J., Ranucci, A., Kim, H., Mudra, S. E., Cashion, A. K., et al. (2019). The oral microbiome of patients undergoing treatment for severe aplastic anemia: A pilot study. Ann. Hematol. 98, 1351–1365. doi: 10.1007/s00277-019-03599-w
Asare, K., and Gatzke, C. B. (2020). Mycophenolate-induced oral ulcers: Case report and literature review. Am. J. Health Syst. Pharm. 77, 523–528. doi: 10.1093/ajhp/zxz358
Bai, H., Yang, J., Meng, S., and Liu, C. (2022). Oral microbiota-driven cell migration in carcinogenesis and metastasis. Front. Cell. Infect. Microbiol. 12, 86f4479. doi: 10.3389/fcimb.2022.864479
Barka, E. A., Vatsa, P., Sanchez, L., Gaveau-Vaillant, N., Jacquard, C., Meier-Kolthoff, J. P., et al. (2016). Taxonomy, physiology, and natural products of actinobacteria. Microbiol. Mol. Biol. Rev. 80, 1–43. doi: 10.1128/MMBR.00019-15
Belstrøm, D., Grande, M. A., Sembler-Møller, M. L., Kirkby, N., Cotton, S. L., Paster, B. J., et al. (2018a). Influence of periodontal treatment on subgingival and salivary microbiotas. J. Periodontol. 89, 531–539. doi: 10.1002/JPER.17-0377
Belstrøm, D., Sembler-Møller, M. L., Grande, M. A., Kirkby, N., Cotton, S. L., Paster, B. J., et al. (2018b). Impact of oral hygiene discontinuation on supragingival and salivary microbiomes. JDR Clin. Transl. Res. 3, 57–64. doi: 10.1177/2380084417723625
Bouteloup, V., Sabin, C., Mocroft, A., Gras, L., Pantazis, N., Le Moing, V., et al. (2017). Reference curves for CD4 T-cell count response to combination antiretroviral therapy in HIV-1-infected treatment-naïve patients. HIV Med. 18, 33–44. doi: 10.1111/hiv.12389
Cao, C., Fan, B., Zhu, J., Zhu, N., Cao, J. Y., Yang, D. R., et al. (2022). Association of gut microbiota and biochemical features in a Chinese population with renal uric acid stone. Front. Pharmacol. 13, 888883. doi: 10.3389/fphar.2022.888883
Chen, Q., Ren, D., Wu, J., Yu, H., Chen, X., Wang, J., et al. (2021). Shenyan Kangfu tablet alleviates diabetic kidney disease through attenuating inflammation and modulating the gut microbiota. J. Nat. Med. 75, 84–98. doi: 10.1007/s11418-020-01452-3
Crespo-Salgado, J., Vehaskari, V. M., Stewart, T., Ferris, M., Zhang, Q., Wang, G., et al. (2016). Intestinal microbiota in pediatric patients with end stage renal disease: A Midwest Pediatric Nephrology Consortium study. Microbiome 4, 50. doi: 10.1186/s40168-016-0195-9
David, L. A., Materna, A. C., Friedman, J., Campos-Baptista, M. I., Blackburn, M. C., Perrotta, A., et al. (2014). Host lifestyle affects human microbiota on daily timescales. Genome Biol. 15, R89. doi: 10.1186/gb-2014-15-7-r89
Deng, Z. L., Sztajer, H., Jarek, M., Bhuju, S., and Wagner-Döbler, I. (2018). Worlds apart—Transcriptome profiles of key oral microbes in the periodontal pocket compared to single laboratory culture reflect synergistic interactions. Front. Microbiol. 9, 124. doi: 10.3389/fmicb.2018.00124
Duan, X., Chen, X., Gupta, M., Seriwatanachai, D., Xue, H., Xiong, Q., et al. (2020). Salivary microbiome in patients undergoing hemodialysis and its associations with the duration of the dialysis. BMC Nephrol. 21, 414. doi: 10.1186/s12882-020-02009-y
Guo, S., Wu, G., Liu, W., Fan, Y., Song, W., Wu, J., et al. (2022). Characteristics of human oral microbiome and its non-invasive diagnostic value in chronic kidney disease. Biosci. Rep. 42, 10694. doi: 10.1042/BSR20210694
Guo, Y., Xia, W., Wei, F., Feng, W., Duan, J., Sun, X., et al. (2021). Salivary microbial diversity at different stages of human immunodeficiency virus infection. Microbial Pathog. 155, 104913. doi: 10.1016/j.micpath.2021.104913
Henrich, B., Rumming, M., Sczyrba, A., Velleuer, E., Dietrich, R., Gerlach, W., et al. (2014). Mycoplasma salivarium as a dominant coloniser of Fanconi anaemia associated oral carcinoma. PLoS ONE 9, e92297. doi: 10.1371/journal.pone.0092297
Ho, H. J., Kikuchi, K., Oikawa, D., Watanabe, S., Kanemitsu, Y., Saigusa, D., et al. (2021). SGLT-1-specific inhibition ameliorates renal failure and alters the gut microbial community in mice with adenine-induced renal failure. Physiol. Rep. 9, e15092. doi: 10.14814/phy2.15092
Hu, J., Iragavarapu, S., Nadkarni, G. N., Huang, R., Erazo, M., Bao, X., et al. (2018). Location-specific oral microbiome possesses features associated with CKD. Kidney Int. Rep. 3, 193–204. doi: 10.1016/j.ekir.2017.08.018
Hu, X., Du, J., Xie, Y., Huang, Q., Xiao, Y., Chen, J., et al. (2020). Fecal microbiota characteristics of Chinese patients with primary IgA nephropathy: A cross-sectional study. BMC Nephrol. 21, 97. doi: 10.1186/s12882-020-01741-9
Imahashi, M., Ode, H., Kobayashi, A., Nemoto, M., Matsuda, M., Hashiba, C., et al. (2021). Impact of long-term antiretroviral therapy on gut and oral microbiotas in HIV-1-infected patients. Sci. Rep. 11, 960. doi: 10.1038/s41598-020-80247-8
Jiang, J. W., Ren, Z. G., Lu, H. F., Zhang, H., Li, A., Cui, G. Y., et al. (2018). Optimal immunosuppressor induces stable gut microbiota after liver transplantation. World J. Gastroenterol. 24, 3871–3883. doi: 10.3748/wjg.v24.i34.3871
Li, X., Zheng, J., Wang, J., Tang, X., Zhang, F., Liu, S., et al. (2022). Effects of uremic clearance granules on p38 MAPK/NF-κB signaling pathway, microbial and metabolic profiles in end-stage renal disease rats receiving peritoneal dialysis. Drug Design Dev. Ther. 16, 2529–2544. doi: 10.2147/DDDT.S364069
Li, Y., Wang, Y., Li, J., Ling, Z., Chen, W., Zhang, L., et al. (2021). Tacrolimus inhibits oral carcinogenesis through cell cycle control. Biomed. Pharmacother. 139, 111545. doi: 10.1016/j.biopha.2021.111545
Lin, D., Yang, L., Wen, L., Lu, H., Chen, Q., Wang, Z., et al. (2021). Crosstalk between the oral microbiota, mucosal immunity, and the epithelial barrier regulates oral mucosal disease pathogenesis. Mucosal Immunol. 14, 1247–1258. doi: 10.1038/s41385-021-00413-7
Liu, F., Sheng, J., Hu, L., Zhang, B., Guo, W., Wang, Y., et al. (2022). Salivary microbiome in chronic kidney disease: What is its connection to diabetes, hypertension, and immunity? J. Transl. Med. 20, 387. doi: 10.1186/s12967-022-03602-5
Llorenç, V., Nakamura, Y., Metea, C., Karstens, L., Molins, B., Lin, P., et al. (2022). Antimetabolite drugs exhibit distinctive immunomodulatory mechanisms and effects on the intestinal microbiota in experimental autoimmune uveitis. Investig. Ophthalmol. Vis. Sci. 63, 30. doi: 10.1167/iovs.63.3.30
Long, J., Cai, Q., Steinwandel, M., Hargreaves, M. K., Bordenstein, S. R., Blot, W. J., et al. (2017). Association of oral microbiome with type 2 diabetes risk. J. Periodontal Res. 52, 636–643. doi: 10.1111/jre.12432
Marsh, P. D., Do, T., Beighton, D., and Devine, D. A. (2016). Influence of saliva on the oral microbiota. Periodontology 70, 80–92. doi: 10.1111/prd.12098
Matsha, T. E., Prince, Y., Davids, S., Chikte, U., Erasmus, R. T., Kengne, A. P., et al. (2020). Oral microbiome signatures in diabetes mellitus and periodontal disease. J. Dental Res. 99, 658–665. doi: 10.1177/0022034520913818
Meijers, B., Evenepoel, P., and Anders, H. J. (2019). Intestinal microbiome and fitness in kidney disease. Nat. Rev. Nephrol. 15, 531–545. doi: 10.1038/s41581-019-0172-1
Miao, V., and Davies, J. (2010). Actinobacteria: the good, the bad, and the ugly. Antonie van Leeuwenhoek 98, 143–150. doi: 10.1007/s10482-010-9440-6
Mizuki, H., Abe, R., and Mikami, T. (2017). Ultrastructural changes during the life cycle of Mycoplasma salivarium in oral biopsies from patients with oral leukoplakia. Front. Cell. Infect. Microbiol. 7, 403. doi: 10.3389/fcimb.2017.00403
Mizuki, H., Kawamura, T., and Nagasawa, D. (2015). In situ immunohistochemical detection of intracellular Mycoplasma salivarium in the epithelial cells of oral leukoplakia. J. Oral Pathol. Med. 44, 134–144. doi: 10.1111/jop.12215
Morand, M., and Hatami, A. (2018). Multiple superficial oral mucoceles after Mycoplasma-induced mucositis. Pediatr. Dermatol. 35, e210–e1. doi: 10.1111/pde.13515
Perez Rosero, E., Heron, S., Jovel, J., O'Neil, C. R., Turvey, S. L., Parashar, P., et al. (2021). Differential signature of the microbiome and neutrophils in the oral cavity of HIV-infected individuals. Front. Immunol. 12, 780910. doi: 10.3389/fimmu.2021.780910
Reed, L. A., O'Bier, N. S., Oliver, L. D. Jr., Hoffman, P. S., and Marconi, R. T. (2018). Antimicrobial activity of amixicile against Treponema denticola and other oral spirochetes associated with periodontal disease. J. Periodontol. 89, 1467–1474. doi: 10.1002/JPER.17-0185
Ren, Z., Fan, Y., Li, A., Shen, Q., Wu, J., Ren, L., et al. (2020). Alterations of the human gut microbiome in chronic kidney disease. Adv. Sci. 7, 2001936. doi: 10.1002/advs.202001936
Robles-Vera, I., de la Visitación, N., Toral, M., Sánchez, M., Gómez-Guzmán, M., Jiménez, R., et al. (2021). Mycophenolate mediated remodeling of gut microbiota and improvement of gut-brain axis in spontaneously hypertensive rats. Biomed. Pharmacother. 135, 111189. doi: 10.1016/j.biopha.2020.111189
Rysz, J., Franczyk, B., Ławiński, J., Olszewski, R., Ciałkowska-Rysz, A., Gluba-Brzózka, A., et al. (2021). The impact of CKD on uremic toxins and gut microbiota. Toxins 13, 40252. doi: 10.3390/toxins13040252
Salguero, M. V., Al-Obaide, M. A. I., Singh, R., Siepmann, T., and Vasylyeva, T. L. (2019). Dysbiosis of Gram-negative gut microbiota and the associated serum lipopolysaccharide exacerbates inflammation in type 2 diabetic patients with chronic kidney disease. Exp. Therapeut. Med. 18, 3461–3469. doi: 10.3892/etm.2019.7943
Shivani, S., Kao, C. Y., Chattopadhyay, A., Chen, J. W., Lai, L. C., Lin, W. H., et al. (2022). Uremic toxin-producing bacteroides species prevail in the gut microbiota of Taiwanese CKD patients: An analysis using the new Taiwan microbiome baseline. Front. Cell. Infect. Microbiol. 12, 726256. doi: 10.3389/fcimb.2022.726256
Sinha, R., Zhao, N., Goedert, J. J., Byrd, D. A., Wan, Y., Hua, X., et al. (2021). Effects of processed meat and drinking water nitrate on oral and fecal microbial populations in a controlled feeding study. Environ. Res. 197, 111084. doi: 10.1016/j.envres.2021.111084
Sugurmar, A. N. K., Mohd, R., Shah, S. A., Neoh, H. M., and Cader, R. A. (2021). Gut microbiota in immunoglobulin A nephropathy: A Malaysian perspective. BMC Nephrol. 22, 145. doi: 10.1186/s12882-021-02315-z
Tuganbaev, T., Yoshida, K., and Honda, K. (2022). The effects of oral microbiota on health. Science 376, 934–936. doi: 10.1126/science.abn1890
Wang, H., Ainiwaer, A., Song, Y., Qin, L., Peng, A., Bao, H., et al. (2023). Perturbed gut microbiome and fecal and serum metabolomes are associated with chronic kidney disease severity. Microbiome 11, 3. doi: 10.1186/s40168-022-01443-4
Wang, X., Yang, S., Li, S., Zhao, L., Hao, Y., Qin, J., et al. (2020). Aberrant gut microbiota alters host metabolome and impacts renal failure in humans and rodents. Gut 69, 2131–2142. doi: 10.1136/gutjnl-2019-319766
Wang, Y., Huang, J. M., Zhou, Y. L., Almeida, A., Finn, R. D., Danchin, A., et al. (2020). Phylogenomics of expanding uncultured environmental Tenericutes provides insights into their pathogenicity and evolutionary relationship with Bacilli. BMC Genomics 21, 408. doi: 10.1186/s12864-020-06807-4
Wilson, K. F., Meier, J. D., and Ward, P. D. (2014). Salivary gland disorders. Am. Fam. Physician 89, 882–888. Available online at: https://www.aafp.org/pubs/afp/issues/2014/0601/p882.html
Wu, R., Ruan, X. L., Ruan, D. D., Zhang, J. H., Wang, H. L., Zeng, Q. Z., et al. (2021). Differences in gut microbiota structure in patients with stages 4-5 chronic kidney disease. Am. J. Transl. Res. 13, 10056–10074.
Yang, L., Dunlap, D. G., Qin, S., Fitch, A., Li, K., Koch, C. D., et al. (2020). Alterations in oral microbiota in HIV are related to decreased pulmonary function. Am. J. Respirat. Crit. Care Med. 201, 445–457. doi: 10.1164/rccm.201905-1016OC
Yang, Y., Zhu, X., Huang, Y., Zhang, H., Liu, Y., Xu, N., et al. (2022). RNA-seq and 16S rRNA analysis revealed the effect of deltamethrin on channel catfish in the early stage of acute exposure. Front. Immunol. 13, 916100. doi: 10.3389/fimmu.2022.916100
Yousefi, L., Leylabadlo, H. E., Pourlak, T., Eslami, H., Taghizadeh, S., Ganbarov, K., et al. (2020). Oral spirochetes: Pathogenic mechanisms in periodontal disease. Microbial Pathog. 144, 104193. doi: 10.1016/j.micpath.2020.104193
Zeng, H., Chan, Y., Gao, W., Leung, W. K., and Watt, R. M. (2021). Diversity of treponema denticola and other oral treponeme lineages in subjects with periodontitis and gingivitis. Microbiol. Spectr. 9, e0070121. doi: 10.1128/Spectrum.00701-21
Zhang, J., Luo, D., Lin, Z., Zhou, W., Rao, J., Li, Y., et al. (2020). Dysbiosis of gut microbiota in adult idiopathic membranous nephropathy with nephrotic syndrome. Microbial Pathog. 147, 104359. doi: 10.1016/j.micpath.2020.104359
Zhang, Z., Liu, L., Tang, H., Jiao, W., Zeng, S., Xu, Y., et al. (2018). Immunosuppressive effect of the gut microbiome altered by high-dose tacrolimus in mice. Am. J. Transplant. 18, 1646–1656. doi: 10.1111/ajt.14661
Zhou, C., Li, K., Zhao, L., Li, W., Guo, Z., Xu, J., et al. (2020). The relationship between urinary stones and gut microbiomeby 16S sequencing. BioMed Res. Int. 2020, 1582187. doi: 10.1155/2020/1582187
Zhou, R., Wen, W., Gong, X., Zhao, Y., and Zhang, W. (2022). Nephro-protective effect of Daphnetin in hyperoxaluria-induced rat renal injury via alterations of the gut microbiota. J. Food Biochem. 2022, e14377. doi: 10.1111/jfbc.14377
Zhou, S., Dong, J., Liu, Y., Yang, Q., Xu, N., Yang, Y., et al. (2021). Effects of acute deltamethrin exposure on kidney transcriptome and intestinal microbiota in goldfish (Carassius auratus). Ecotoxicol. Environ. Saf. 225, 112716. doi: 10.1016/j.ecoenv.2021.112716
Keywords: salivary microbiota, renal function, renal transplantation, perioperative period, 16S rRNA
Citation: Xiang X, Peng B, Liu K, Wang T, Ding P, Li H, Zhu Y and Ming Y (2023) Association between salivary microbiota and renal function in renal transplant patients during the perioperative period. Front. Microbiol. 14:1122101. doi: 10.3389/fmicb.2023.1122101
Received: 12 December 2022; Accepted: 07 March 2023;
Published: 29 March 2023.
Edited by:
Weiqi He, Soochow University, ChinaReviewed by:
Jozsef Soki, University of Szeged, HungaryZheng Sun, Harvard Medical School, United States
Copyright © 2023 Xiang, Peng, Liu, Wang, Ding, Li, Zhu and Ming. 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: Yingzi Ming, bWluZ3l6X2NoaW5hJiN4MDAwNDA7Y3N1LmVkdS5jbg==
†These authors have contributed equally to this work and share first authorship